A111617 - OEIS

%I #4 Mar 30 2012 18:36:50

%S 6,10,18,20,22,30,34,36,38,44,50,52,56,58,62,72,78,84,86,90,92,96,98,

%T 100,102,104,110,116,134

%N Where A111616(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]]; Select[ Range[141], # / g[ #, # + 3] == 2 &]

%Y Cf. A109626, A111615.

%K nonn

%O 1,1

%A _Paul D. Hanna_ and _Robert G. Wilson v_, Jul 30 2005