> with (numtheory); -------------------------------------------------------------------------------- > fft_param:=proc(n) > pms:=max(op(invphi(n))); > d:=gcd(pms,2*3*5); > for i from pms while i<2*pms by d do > if phi(i)<=n then > pms:=i; > # print(phi(pms),pms,ifactor(pms)); > fi; > od; > RETURN(n,phi(pms),pms,ifactor(pms),n/2*pms); > end; -------------------------------------------------------------------------------- > for k from 8 to 20 do fft_param(2^k); od; 2 256, 240, 1050, (2) (3) (5) (7), 134400 512, 480, 2310, (2) (3) (5) (7) (11), 591360 2 1024, 960, 4620, (2) (3) (5) (7) (11), 2365440 3 2048, 1920, 9240, (2) (3) (5) (7) (11), 9461760 2 4096, 4032, 19110, (2) (3) (5) (7) (13), 39137280 8192, 7680, 39270, (2) (3) (5) (7) (11) (17), 160849920 16384, 16128, 79170, (2) (3) (5) (7) (13) (29), 648560640 2 32768, 31680, 159390, (2) (3) (5) (7) (11) (23), 2611445760 2 65536, 63360, 330330, (2) (3) (5) (7) (11) (13), 10824253440 131072, 126720, 690690, (2) (3) (5) (7) (11) (13) (23), 45265059840 2 262144, 253440, 1381380, (2) (3) (5) (7) (11) (13) (23), 181060239360 2 524288, 518400, 2852850, (2) (3) (5) (7) (11) (13) (19), 747857510400 2 2 1048576, 1036800, 5705700, (2) (3) (5) (7) (11) (13) (19), 2991430041600 -------------------------------------------------------------------------------- >