尤拉計劃712:指數差

這個p至少要是D的一個引數的質因數才有意義。 比如14，24只要算2、3、7。
D(14,24)
=|v₂(14)-v₂(24)| +|v₃(14)-v₃(24)| +|v₇(14)-v₇(24)|
=|1-3|+|0-1|+|1-0|
=4

測試程式碼

def v(p,n):
 n1=n
 t=0
 while n1%p==0: n1/=p;t+=1
 return t 

primes=[2,3,5,7,11,13,17]
def prime(n):
 primes=[]
 h  = [True] * int(n+1)
 h[:2] = [False, False]
 for i in range(2, int(n ** (1.0*0.5)) + 1):
    if h[i]:
        h[i*i::i] = [False] * len(h[i*i::i])
 for i, e in enumerate(h):
    if e:
        primes.append(i)
 return primes

primes=prime(1000)

def D(n,m):
 s=0
 for p in primes:
  if p>max(n,m):break
  s+=abs(v(p,n)-v(p,m))
 return s 

def S(N):
 s=0
 for n in range(1,N+1):
  for m in range(n+1,N+1):
   s+=D(n,m)
 return s*2