

A174065


Convolved with its aerated variant = A000041.


9



1, 1, 1, 2, 3, 4, 5, 7, 9, 12, 15, 19, 25, 31, 38, 48, 60, 73, 89, 109, 133, 161, 193, 232, 279, 333, 395, 470, 558, 658, 775, 912, 1071, 1254, 1464, 1708, 1991, 2313, 2681, 3107, 3595, 4149, 4782, 5506, 6331, 7268, 8330, 9538, 10912, 12462, 14213, 16199
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

0,4


COMMENTS

A000041 = (1, 1, 2, 3, 5, 7, ...) = (1, 1, 1, 2, 3, 4, 5, 7, ...) * (1, 0, 1, 0, 1, 0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 7, 0, 9, 0, ...).
The sequence diverges from A100853 after 16 terms; and is a conjectured Euler transform of A035263: (1, 0, 1, 1, 1, 0, 1, 0, 1, ...).


LINKS

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000 (terms 0..5000 from Alois P. Heinz)


FORMULA

Aerate and convolve sequences are generated by triangles (in this case A174066) in which ongoing terms are placed in the left column and at the top as a heading. Columns >1 are shifted down k times (k=2) in this case corresponding to (k1) interpolated zeros. Next term in left column = nth term in the "target sequence" S(n) (in this case A000041) minus (sum of terms in nth row for columns >1). Place the latter term in the heading filling in missing terms.
G.f.: Product_{i>=1, j>=0} (1 + x^(i*4^j)).  Ilya Gutkovskiy, Sep 23 2019
a(n) ~ exp(2*Pi*sqrt(n)/3) / (2^(11/8) * 3^(3/4) * n^(7/8)).  Vaclav Kotesovec, Sep 24 2019


EXAMPLE

Heading at top, with triangle A174066 underneath (the generator for A174065):
1, 1, 1, 2, 3, 4,.... = heading
1;................... = 1
1;................... = 1
1, 1;................ = 2
2, 1;................ = 3
3, 1, 1;............. = 5
4, 2, 1;............. = 7
5, 3, 1, 2;.......... = 11
7, 4, 2, 2;.......... = 15
9, 5, 3, 2, 3;....... = 22
...
... where terms in the left column are the result of the two rules: multiply heading * left column, and row sums = partition numbers.
Thus leftmost term in column 8 must be 7 = 15  (4 + 2 + 2). Then the 7 is placed in its spot in the left column and as the next heading term.


MAPLE

p:= combinat[numbpart]:
a:= proc(n) option remember; `if`(n=0, 1, p(n)add(a(j)*
`if`(irem(nj, 2, 'r')=1, 0, a(r)), j=0..n1))
end:
seq(a(n), n=0..61); # Alois P. Heinz, Jul 27 2019


MATHEMATICA

nmax = 60; CoefficientList[Series[Product[QPochhammer[1, x^(4^j)]/2, {j, 0, Log[nmax]/Log[4]}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Sep 24 2019 *)


CROSSREFS

Cf. A000041, A174066, A174067.
Sequence in context: A280909 A003413 A100853 * A121659 A096814 A039861
Adjacent sequences: A174062 A174063 A174064 * A174066 A174067 A174068


KEYWORD

nonn


AUTHOR

Gary W. Adamson, Mar 06 2010


EXTENSIONS

More terms from R. J. Mathar, Mar 18 2010
Offset corrected by Alois P. Heinz, Jul 27 2019


STATUS

approved



