似乎没什么值得一提的，用筛法求出范围内所有数的质因数，若
$$n=p_1^{a_1}\ldots p_k^{a_k},$$
是n的质因数分解，则显然其约数总数为$(1+a_1)\ldots(1+a_k)$。直接计算和比较即得结果$986262$。
target=10**7
subTarget=3200
d=[[] for i in range(target)]
n=2
while n < subTarget:
for i in range(n+n,target,n):
d[i].append(n)
n+=1
while n < subTarget and d[n]:
n+=1
r=1000
s=1
for n in range(8985,target):
if d[n]:
m=n
t=1
for p in d[n]:
i=1
while m%p==0:
m/=p
i+=1
t*=i
if m>1:
t*=2
if s==t:
r+=1
s=t
else:
s=2
print r
|
