|
MATHEMATICA
|
a[n_] := a[n] = Block[{s = Sum[ a[i]*x^i, {i, 0, n - 1}]}, If[ IntegerQ@ Last@ CoefficientList[ Series[ Sqrt[s + x^n], {x, 0, n}], x], 1, 2]]; t = Union@ Table[ If[ a[n] == 1, n, 0], {n, 0, 5506}]; f[n_] := Block[{k = n, m = 2}, If[n == 1, 0, While[m < Length@t && t[[m - 1]] + 2 == t[[m]] || t[[m]] + 2k != t[[m + k]] || t[[m + k]] + 2 == t[[m + k + 1]], m++ ]; If[m < Length@t - k, t[[m]], 0]]]; Table[ f[n], {n, 0, 9}]
|