seq[kmax_] := Module[{s0 = 1, s1 = 2, s2 = 4, s3, s = {}}, Do[s3 = (3*(k+1)*s2 + (k-4)*s1 - 3*(k-1)*s0)/(k+2); If[EvenQ[k] && !PrimeQ[k+1] && Divisible[s3 - 2*(-1)^(k/2), k+1], AppendTo[s, k+1]]; s0 = s1; s1 = s2; s2 = s3, {k, 3, kmax}]; s]; seq[500]