|
%I #10 Jul 30 2015 22:41:38
%S 1,1,1,2,1,12,1,16,6,88,1,960,1,656,264,896,1,48960,1,53504,1968,
%T 36544,1,2795520,44,272768,2448,2980864,1,1547335680,1,2781184,109632,
%U 15196672,14432,8635760640,1,113429504,818304,8677064704,1,4808968273920,1,9252356096,415337472,6319476736,1,26795484119040,328,3584860454912
%N Sylvester dividends for A002605.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/SylvesterCyclotomicNumber.html">Sylvester Cyclotomic Number</a>.
%F a(n) = A002605(n)/A105607(n); a(n) = Product_{k=1..n-1, if(gcd(n, k)>1, (1+sqrt(3))-(1-sqrt(3))*exp(2*Pi*I*k/n), 1)}, I=sqrt(-1),
%t f[n_] := FullSimplify[ Expand[ Times @@ ((1+Sqrt[3])-(1-Sqrt[3])*Exp[2Pi*I*Select[Range[n], GCD[ #, n] == 1 &]/n])]]; Table[ f[n], {n, 2, 32] (* _Robert G. Wilson v_, Aug 02 2005 *)
%Y Cf. A105607, A105609.
%K easy,nonn
%O 1,4
%A _Paul Barry_, Apr 15 2005
%E More terms from _David Wasserman_, May 06 2008
|
