去看闪闪发光的砖头,rosettacode.org上的强力手镯test for primes页面defExpand_x_1(p ) :
forIinrange(p ) :
return ex[:-1]
efaks_test(p ) :
if p 2:返回假
returnnotany (多% pformultinex [ 0:-1 ]
打印(' # p : (x-1 ) ^p for small p () ) ) ) ) ) )。
forpinrange(12 ) :
print('%3I:%s'%(p,“. joinzqdyi%s' % (e ) e,) (x ^ % I“% n”if n else“”) ) )。
for n,einenumerate(Expand_x_1) p ) )
print ((n # smallprimesusingtheakstest ) )。
print([pforpinrange(101 ) if aks_test(p ) p ) ]
0: 1
1: -1 1x^1
2: 1 -2x^1 1x^2
3: -1 3x^1 -3x^2 1x^3
4: 1 -4x^1 6x^2 -4x^3 1x^4
3:-15x ^1- 10x ^ 210 x ^3- 5x ^ 41x ^ 5
6:1-6x ^ 115 x ^2- 20x ^ 315 x ^4- 6x ^ 51x ^ 6
7:-17x ^1- 21x ^ 235 x ^3- 35x ^ 421 x ^5- 7x ^ 61x ^ 7
8:1-8x ^ 128 x ^2- 56x ^ 370 x ^4- 56x ^ 528 x ^6- 8x ^ 71x ^ 8
93360-19x ^1- 36x ^ 284 x ^3- 126 x ^ 4126 x ^5- 84x ^ 636 x ^7- 9x ^ 81x ^ 9
10:1-10x ^ 145 x ^2- 120 x ^ 3210 x ^4- 252 x ^ 5210 x ^6- 120 x ^ 745 x ^8- 10x ^ 91x ^ 10
1:-111 x ^1- 55x ^ 2165 x ^3- 330 x ^ 4462 x ^5- 462 x ^ 6330 x ^7- 165 x ^ 855 x ^9- 11x ^ 101 x ^ 11
#小primes using the aks test
[2、3、5、7、11、13、17、19、29、31、37、43、47、53、59、61、67、71、73、79、83、89、97]