OFFSET
0,3
COMMENTS
a(n) = Phi_21(n) where Phi_k(x) is the k-th cyclotomic polynomial.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Cyclotomic Polynomial
Index entries for linear recurrences with constant coefficients, signature (13, -78, 286, -715, 1287, -1716, 1716, -1287, 715, -286, 78, -13, 1).
FORMULA
G.f.: (1 - 12*x + 2424*x^2 + 337214*x^3 + 8182695*x^4 + 58741344*x^5 + 156377856*x^6 + 168607380*x^7 + 73943271*x^8 + 12191420*x^9 + 612600*x^10 + 5406*x^11 + x^12)/(1-x)^13.
Sum_{n>=0} 1/a(n) = 2.00042670913...
MATHEMATICA
Table[Cyclotomic[21, n], {n, 0, 17}]
CoefficientList[Series[(1 -12x +2424x^2 +337214x^3 +8182695x^4 +58741344 x^5 +156377856x^6 +168607380x^7 +73943271x^8 +12191420x^9 + 612600 x^10 +5406x^11 +x^12)/(1-x)^13, {x, 0, 33}], x] (* Vincenzo Librandi, Feb 28 2016 *)
PROG
(Magma) [n^12-n^11+n^9-n^8+n^6-n^4+n^3-n+1: n in [0..20]]; // Vincenzo Librandi, Feb 28 2016
(Python)
A269483_list, m = [], [479001600, -2674425600, 6386688000, -8501915520, 6889478400, -3482100720, 1080164160, -194177280, 17948256, -666714, 5418, 0, 1]
for _ in range(10**2):
A269483_list.append(m[-1])
for i in range(12):
m[i+1] += m[i] # Chai Wah Wu, Feb 28 2016
(PARI) a(n) = polcyclo(21, n); \\ Michel Marcus, Feb 29 2016
(Sage) [n^12-n^11+n^9-n^8+n^6-n^4+n^3-n+1 for n in (0..20)] # G. C. Greubel, Apr 24 2019
(GAP) List([0..20], n-> n^12-n^11+n^9-n^8+n^6-n^4+n^3-n+1) # G. C. Greubel, Apr 24 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Ilya Gutkovskiy, Feb 27 2016
STATUS
approved