OFFSET
0,4
LINKS
Paul D. Hanna, Table of n, a(n), n = 0..500.
EXAMPLE
G.f.: A(x) = 1 + x + x^2 + 2*x^3 + 3*x^4 + 4*x^5 + 5*x^6 + 6*x^7 +...
The square of the g.f. A(x) begins:
A(x)^2 = 1 + 2*x + 3*x^2 + 6*x^3 + 11*x^4 + 18*x^5 + 28*x^6 + 42*x^7 + 61*x^8 + 88*x^9 + 124*x^10 + 168*x^11 + 227*x^12 + 302*x^13 + 394*x^14 +...
The integer part of the square-root of the coefficients of the square of the g.f. A(x) equals the respective coefficients of A(x), as illustrated by:
[sqrt(1)]=1; [sqrt(2)]=1; [sqrt(3)]=1; [sqrt(6)]=2; [sqrt(11)]=3;
[sqrt(18)]=4; [sqrt(28)]=5; [sqrt(42)]=6; [sqrt(61)]=7; [sqrt(88)]=9;
[sqrt(124)]=11; [sqrt(168)]=12; [sqrt(227)]=15; [sqrt(302)]=17; ...
PROG
(PARI) {a(n)=local(A=[1, 1, 1, t], T); for(i=1, n-2, for(N=A[ #A-1], 2*A[ #A-1], T=subst(Vec(Ser(A)^2)[ #A], t, N); if(sqrtint(T)==N, A[ #A]=N; A=concat(A, t); break))); A[n+1]}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Sep 22 2010
STATUS
approved