|
%I #4 Mar 30 2012 18:36:50
%S 1,3,5,7,11,15,17,19,21,23,27,31,33,35,39,43,53,57,59,61,63,67,79,95,
%T 97,101,115,127,129,131,135,137,139,143
%N Where A111621(n)=2.
%t f[n_] := f[n] = Block[{a}, a[0] = 1; a[l_] := a[l] = Block[{k = 1, s = Sum[ a[i]*x^i, {i, 0, l - 1}]}, While[ IntegerQ[ Last[ CoefficientList[ Series[(s + k*x^l)^(1/n), {x, 0, l}], x]]] != True, k++ ]; k]; Table[a[j], {j, 0, 144}]]; g[n_, m_] := f[n][[m]]; t = Table[ n/g[n, n - 1], {n, 2, 144}]; Select[ Range[143], t[[ # ]] == 2 &]
%Y Cf. A109626, A111618.
%K nonn
%O 1,2
%A _Paul D. Hanna_ and _Robert G. Wilson v_, Aug 01 2005
|
