|
|
A060719
|
|
a(0) = 1; a(n+1) = a(n) + Sum_{i=0..n} binomial(n,i)*(a(i)+1).
|
|
6
|
|
|
1, 3, 9, 29, 103, 405, 1753, 8279, 42293, 231949, 1357139, 8427193, 55288873, 381798643, 2765917089, 20960284293, 165729739607, 1364153612317, 11665484410113, 103448316470743, 949739632313501, 9013431476894645, 88304011710168691
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
a(3) = 29 = (30 - 1) = A186021(4) - 1
|
|
MAPLE
|
A060719 := proc(n) option remember; local i; if n=0 then 1 else A060719(n-1)+add(binomial(n-1, i)*(A060719(i)+1), i=0..n-1); fi; end;
|
|
MATHEMATICA
|
|
|
PROG
|
(PARI) vector(26, n, my(m=n-1); 2*sum(k=0, m+1, stirling(m+1, k, 2)) -1 ) \\ G. C. Greubel, Feb 12 2020
(Magma) [2*Bell(n+1) -1: n in [0..25]]; // G. C. Greubel, Feb 12 2020
(Sage) [2*bell_number(n+1)-1 for n in (0..25)] # G. C. Greubel, Feb 12 2020
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|