OFFSET
1,1
FORMULA
a(n)= (sqrt(11)-1)^degree(cyclotomic(n,x),x)*cyclotomic(n,(6+sqrt(11)/5) L(n)=10*F(n-1)+F(n+1) F(2n)=Product(d|2n) a(d), F(2n+1)=Product(d|2n+1) a(2d). L(2n+1)=Product(d|2n+1, a(d)), for k>0: L(2^k*(2n+1))=Product(d|2n+1, a(2^(k+1)*d)). for odd prime p, a(p)=L(p)/2, a(2p)=f(p) a(1)=2, a(2)=1; a(2^(k+1))=L(2^k);
EXAMPLE
F(12)=a(1)*a(2)*a(3)*a(4)*a(6)*a(12)=2*1*34*24*14*276=6306048
F(9)=a(2)*a(6)*a(18)= 1*14*5624=78736
L(12)=a(8)*a(24)=376*111376=41877376
L(21)=a(1)*a(3)*a(7)*a(21)=2*34*13784*23235136=21778571794432
MAPLE
with(numtheory): a[1]:=2:a[2]:=1:for n from 3 to 60 do a[n]:=round(evalf((sqrt(11)-1)^degree(cyclotomic(n, x), x) *cyclotomic(n, (6+sqrt(11))/5), 30)) od: seq(a[n], n=1..60);
CROSSREFS
KEYWORD
nonn
AUTHOR
Miklos Kristof, Apr 03 2007
STATUS
approved