%I #4 Mar 30 2012 18:36:50
%S 1,2,3,3,5,3,7,2,9,10,11,9,13,7,15,4,17,18,19,15,21,11,23,6,25,26,27,
%T 21,29,15,31,8,33,34,35,27,37,19,39,10,41,42,43,33,45,23,47,12,49,50,
%U 51,39,53,27,55,14,57,58,59,45,61,31,63,16,65,66,67,51,69,35,71,18,73,74,75
%N Fourth column of A109626.
%F G.f.: (1 + 2x + 3x^2 + 3x^3 + 5x^4 + 3x^5 + 7x^6 + 2x^7 + 7x^8 + 6x^9 +5x^10 + 3x^11 + 3x^12 + x^13 + x^14)/(1-x^8)^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, 128}]]; g[n_, m_] := f[n][[m]]; Table[g[n, 4 + 1], {n, 75}]
%Y Cf. A109626.
%K nonn
%O 1,2
%A _Paul D. Hanna_ and _Robert G. Wilson v_, Aug 01 2005
|