|
Search: id:A147851
|
|
|
| A147851 |
|
Coefficient expansion of : 1/(1 - x^3 - x^4 - x^5 + x^8)^2= 1/(-1+ 2 x^3 + 2 x^4 + 2 x^5 - x^6 - 2 x^7 - 5x^8 - 2 x^9 - x^10 + 2 x^11 + 2 x^12 + 2 x^13 - x^16) |
|
+0 2
|
|
| 1, 0, 0, 2, 2, 2, 3, 6, 7, 10, 15, 18, 27, 38, 50, 66, 92, 126, 165, 224, 300, 400, 536, 714, 948, 1258, 1676, 2218, 2932, 3882, 5128, 6768, 8924, 11760, 15479, 20366, 26780, 35174, 46182, 60602, 79473, 104158, 136445, 178654, 233797, 305834, 399881
(list; graph; listen)
|
|
|
OFFSET
|
0,4
|
|
|
COMMENT
|
Base Polynomial found using census program: ( also in the Salem list as #23): Clear[h, k, a, b, c, d, e, f, g, h, j, k, l, m, p, x]; h[x_] = a + b*x + c*x^2 + d*x^3 + e*x^4 + f*x^5 + g*x^6 + j*x^7 + x^8; n0 = 8; Union[Flatten[Table[If [Abs[x /. NSolve[h[x] == 0, x][[n0]]] - Max[Abs[Table[ x /. NSolve[h[x] == 0, x][[n]], {n,1, n0}]]] == 0 && Im[x /. NSolve[h[x] == 0, x][[n0]]] ==0 && Sort[Abs[Table[x /.NSolve[h[x] == 0, x][[n]], {n, 1, n0}]]][[n0 - 1]] <= 1 &&Sort[Abs[Table[x /. NSolve[h[x] == 0,x][[n]], {n, 1, n0}]]][[1]] > 0, {Sort[Abs[Table[x /.NSolve[h[x] == 0, x][[n]], {n, 1, n0}]]], h[x]}, {}], {a, -1,1}, {b, -1, 1}, {c, -1, 1}, {d, -1, 1}, {e, -1, 1}, {f, -1,1}, {g, -1, 1}, {j, -1, 1}], 7]]
|
|
REFERENCES
|
http://www.cecm.sfu.ca/~mjm/Lehmer/lists/SalemList.html
|
|
FORMULA
|
Coefficient expansion of : 1/(1 - x^3 - x^4 - x^5 + x^8)^2= 1/(-1+ 2 x^3 + 2 x^4 + 2 x^5 - x^6 - 2 x^7 - 5x^8 - 2 x^9 - x^10 + 2 x^11 + 2 x^12 + 2 x^13 - x^16)
|
|
MATHEMATICA
|
f[x_] = 1 - x^3 - x^4 - x^5 + x^8; g[x] = ExpandAll[f[x]*x^8*f[1/x]]; a = Table[SeriesCoefficient[Series[1/g[x], {x, 0, 50}], n], {n, 0, 50}]
|
|
CROSSREFS
|
Sequence in context: A038715 A057040 A096235 this_sequence A143596 A091712 A125721
Adjacent sequences: A147848 A147849 A147850 this_sequence A147852 A147853 A147854
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Nov 15 2008
|
|
|
Search completed in 0.002 seconds
|