(Python)
从sympy导入除数
将numpy导入为np
对于范围(1,10**6)中的n:
….d=除数(n)
….s=总和(d)
….如果不是s%2和2*n<=s:
……..删除(n)
……..s2,ld=int(s/2-n),len(d)
…….z=np.零((ld+1,s2+1),数据类型=int)
……..对于范围(1,ld+1)内的i:
…………y=分钟(d[i-1],s2+1)
…………z[i,范围(y)]=z[i-1,范围(y)]
…………z[i,范围(y,s2+1)]=np.最大值
…………如果z[i,s2]==s2:
…………..断裂
|