OFFSET
0,4
COMMENTS
Self-convolution is A105614.
EXAMPLE
Given this sequence:
A = [1, 1, 1, 2,3, 1, 6,11,12, 1, 24,48,71,78, 1, 155, ...],
the self-convolution yields:
A^2 = [1, 2,3, 6,11,12, 24,48,71,78, 155,272,420,652,903, ...],
giving this sequence with 1's omitted at positions {0,2,5,9,14,...}.
PROG
(PARI) a(n)=local(V, t); if(n==0 || n==1, 1, t=ceil((sqrt(8*n+9)-1)/2); if(ceil((sqrt(8*n+17)-1)/2)-t==1, 1, V=vector(n-t+2, k, a(k-1)); polcoeff(Ser(V)^2+O(x^(n-t+2)), n-t+1)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Apr 16 2005
STATUS
approved