%I #9 Aug 12 2015 02:00:17
%S 1,-1,6,-12,-60,720,-2520,-20160,362880,-1814400,-19958400,479001600,
%T -3113510400,-43589145600,1307674368000,-10461394944000,
%U -177843714048000,6402373705728000,-60822550204416000,-1216451004088320000,51090942171709440000,-562000363888803840000
%N Real term generated from a complex product operation.
%C The signed sequence is a(n)=n!*T_n(-1/2) for n>1 where T are the Chebyshev polynomials. Therefore a(n)=n!=A000142(n) if 3 divides n, else a(n)=-n!/2=-A001710(n) (n>1). - _R. J. Mathar_, Aug 25 2006
%F Let a(1) = (1 Angle 2Pi/3); for n>1, a(n) = Re:[Product(1,n):(n Angle 2Pi/3)].
%e a(3) = 6 = Re:[(1 Angle 2Pi/3)*(2 Angle 2Pi/3)*(3 Angle 2Pi/3).
%p with(orthopoly) ; A121735 := proc(n) if n= 1 then RETURN(1) ; else RETURN( n!*T(n,-1/2)) ; fi ; end: for n from 1 to 25 do print(A121735(n)) ; od ; # _R. J. Mathar_, Aug 25 2006
%o (PARI) A121735(n)={ local(T) ; if(n==1, return(1), x=-1/2 ; T=poltchebi(n) ; return(n!*eval(T)) ; ) ; } { for(n=1,25, print1(A121735(n),",") ; ) ; } \\ _R. J. Mathar_, Nov 07 2006
%K sign
%O 1,3
%A _Gary W. Adamson_, Aug 18 2006
%E More terms from _R. J. Mathar_, Nov 07 2006
|