I think there is something wrong in the formulation of this question.
With Impulse defined as
#vec J = int_(t = a)^b vec F(t) \ dt #
#= int_(t = a)^b vec dot p (t) \ dt = vec p(b) - vec p(a)#
then the Impulse on the object at t= 1 is
#vec J = int_(t = 1)^1 vec F(t) \ dt = vec p(1) - vec p(1) = 0#
It may be that you want the total impulse applied for #t in [0,1]# which is
#vec J = int_(t = 0)^1 vec F(t) \ dt = vec p(1) - vec p(0) qquad star#
To evaluate #star# that we note that if the rate of change of kinetic energy #T# is constant, ie:
#(dT)/(dt) = const #
then
#T= alpha t + beta#
#T(0) = 8 implies beta = 8#
#T(4) = 136 = alpha(4) + 8 implies alpha = 32#
#T= 32 t + 8#
Now #T = abs(vec p)^2/(2m)#.
#implies (vec p * vec p)= 4(32 t + 8)#
#vec p = 2sqrt( (32 t + 8)) hat p#
and
#vec p(1) - vec p (0)#
# = (2sqrt( (32 + 8)) - 2sqrt( 8) )hat p#
# = 4(sqrt( 10) - sqrt( 2) )hat p# N s