|
EXAMPLE
|
1+x = 1 + 1*x/(1+x) + 1*x^2/(1+x)^4 + 3*x^3/(1+x)^9 + 18*x^4/(1+x)^16 + 172*x^5/(1+x)^25 + 2313*x^6/(1+x)^36 +...
Also forms the final terms in rows of the triangle where row n+1 equals the partial sums of row n with the final term repeated 2n+1 times, starting with a '1' in row 0, as illustrated by:
1;
1, 1, 1;
1, 2, 3, 3, 3, 3, 3;
1, 3, 6, 9, 12, 15, 18, 18, 18, 18, 18, 18, 18;
1, 4, 10, 19, 31, 46, 64, 82, 100, 118, 136, 154, 172, 172, 172, 172, 172, 172, 172, 172, 172;
1, 5, 15, 34, 65, 111, 175, 257, 357, 475, 611, 765, 937, 1109, 1281, 1453, 1625, 1797, 1969, 2141, 2313, 2313, 2313, 2313, 2313, 2313, 2313, 2313, 2313; ...
|
|
PROG
|
(PARI) {a(n)=local(F=1/(1+x+x*O(x^n))); polcoeff(1+x-sum(k=0, n-1, a(k)*x^k*F^(k^2)), n)}
(PARI) {A=[1, 1]; for(i=1, 40, A=concat(A, -Vec(sum(n=0, #A-1, A[n+1]*x^n/(1+x+x*O(x^#A))^(n^2)))[#A+1])); for(n=0, #A-1, print1(A[n+1], ", "))}
|