|
|
A109409
|
|
Coefficients of polynomials triangular sequence produced by removing primes from the odd numbers in A028338.
|
|
0
|
|
|
1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 9, 10, 1, 0, 0, 0, 0, 9, 10, 1, 0, 0, 0, 0, 0, 9, 10, 1, 0, 0, 0, 0, 0, 135, 159, 25, 1, 0, 0, 0, 0, 0, 0, 135, 159, 25, 1, 0, 0, 0, 0, 0, 0, 0, 135, 159, 25, 1, 0, 0, 0, 0, 0, 0, 0, 2835, 3474, 684, 46, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,19
|
|
COMMENTS
|
The row sums also appear to be new: b = Flatten[Join[{{1}}, Table[Apply[Plus, Abs[CoefficientList[Product[x + g[n], {n, 0, m}], x]]], {m, 0, 10}]]] {1, 2, 2, 2, 2, 20, 20, 20, 320, 320, 320, 7040} Since the row sum of A028338 is the double factorial A000165: this result seems to be a factorization of the double factorial numbers by relatively sparse nonprime odd numbers. It might be better to reverse the order of the coefficients to get the higher powers first.
|
|
LINKS
|
Table of n, a(n) for n=1..78.
|
|
FORMULA
|
p[n]=Product[If[PrimeQ[2*n+1]==false,x+(2*n+1),x] a(n) =CoefficientList[p[n],x]
|
|
EXAMPLE
|
{1},
{1, 1},
{0, 1, 1},
{0, 0, 1, 1},
{0, 0, 0, 1, 1},
{0, 0, 0, 9, 10, 1},
{0, 0, 0, 0, 9, 10, 1},
{0, 0, 0, 0, 0, 9, 10, 1}
|
|
MATHEMATICA
|
a = Join[{{1}}, Table[CoefficientList[Product[x + g[n], {n, 0, m}], x], {m, 0, 10}]]; Flatten[a]
|
|
CROSSREFS
|
Cf. A028338, A000165, A039757.
Sequence in context: A220450 A318147 A324663 * A262551 A160563 A158286
Adjacent sequences: A109406 A109407 A109408 * A109410 A109411 A109412
|
|
KEYWORD
|
nonn,uned
|
|
AUTHOR
|
Roger L. Bagula, May 19 2007
|
|
STATUS
|
approved
|
|
|
|