（Python）
从sympy导入除数
将numpy导入为np
对于范围（3，10**4，2）中的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：
………….d2=[2*x，对于d中的x，如果n>2*x和n%（2*x）]+\
………….[x代表除数（2*n-1）中的x，如果n>x>=2且n%x]+\
………….[x代表除数（2*n+1）中的x，如果n>x>=2且n%x]
…………..s，dmax=总和（d2），最大值（d2
………….如果不是s%2和2*dmax<=s：
…………..d2.删除（dmax）
s2，ld=int（s/2-dmax），len（d2）
z=np.零（（ld+1，s2+1），数据类型=int）
对于范围（1，ld+1）中的i：
y=最小值（d2[i-1]，s2+1）
…………..z[i，范围（y）]=z[i-1，范围（y）]
…………..z[i，范围（y，s2+1）]=np.最大值（z[i-1，范围（y,s2+1）]，z[i-l，范围（0，s2+1-y）]+y）
如果z[i，s2]==s2：
............................A243104型.append（n）
…………..断裂
…………..断裂
