|
|
A144048
|
|
Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is Euler transform of (j->j^k).
|
|
19
|
|
|
1, 1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 5, 6, 5, 1, 1, 9, 14, 13, 7, 1, 1, 17, 36, 40, 24, 11, 1, 1, 33, 98, 136, 101, 48, 15, 1, 1, 65, 276, 490, 477, 266, 86, 22, 1, 1, 129, 794, 1828, 2411, 1703, 649, 160, 30, 1, 1, 257, 2316, 6970, 12729, 11940, 5746, 1593, 282, 42, 1, 1, 513
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,6
|
|
COMMENTS
|
In general, column k > 0 is asymptotic to (Gamma(k+2)*Zeta(k+2))^((1-2*Zeta(-k)) /(2*k+4)) * exp((k+2)/(k+1) * (Gamma(k+2)*Zeta(k+2))^(1/(k+2)) * n^((k+1)/(k+2)) + Zeta'(-k)) / (sqrt(2*Pi*(k+2)) * n^((k+3-2*Zeta(-k))/(2*k+4))). - Vaclav Kotesovec, Mar 01 2015
|
|
LINKS
|
|
|
FORMULA
|
G.f. of column k: Product_{j>=1} 1/(1-x^j)^(j^k).
|
|
EXAMPLE
|
Square array begins:
1, 1, 1, 1, 1, 1, ...
1, 1, 1, 1, 1, 1, ...
2, 3, 5, 9, 17, 33, ...
3, 6, 14, 36, 98, 276, ...
5, 13, 40, 136, 490, 1828, ...
7, 24, 101, 477, 2411, 12729, ...
|
|
MAPLE
|
with(numtheory): etr:= proc(p) local b; b:= proc(n) option remember; `if`(n=0, 1, add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n) end end: A:= (n, k)-> etr(j->j^k)(n); seq(seq(A(n, d-n), n=0..d), d=0..13);
|
|
MATHEMATICA
|
etr[p_] := Module[{ b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[j]}]*b[n - j], {j, 1, n}]/n]; b]; A[n_, k_] := etr[Function[j, j^k]][n]; Table[Table[A[n, d - n], {n, 0, d}], {d, 0, 13}] // Flatten (* Jean-François Alcover, Dec 27 2013, translated from Maple *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|