- f[n_]:=Module[{P=x,Q=1,s},
- Do[s=Exp[I i Pi/(n+1)];
- {P,Q}={(P/. x->x-s)(Q/. x->x+s),(P/. x->x+s)(Q/. x->x-s)},{i,1,n}];
- (P+Q)/2]
- m[n_]:=Module[{poly=x,s},Do[s=Exp[I i Pi/n];
- poly=(poly/. x->x-s)(poly/. x->x+s),{i,1,n-1}];
- poly]
- Table[FullSimplify[Coefficient[f[n]-m[n],x^(2^(n-1)-2 n)]],{n,5,7}]
Copy the Code {-10240,1712128,-345899008}
有何规律? |