尤拉計劃717:取模公式之和

lt發表於2020-05-24
公式

enter image description here

測試程式:

import time
P=[3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97]
def f(p):
 return (2**(2**p)//p)%(2**p)

def g(p):
 return (2**(2**p)//p)%(2**p)%p

def G(N):
 s=0
 for p in P:
  if p>=N:break
  s+=g(p) 
 return s

t=time.time();print(G(10));print(time.time()-t)

相關文章