N: =10^9:#获取所有术语<N
n0:=mul(i),i=1..8):
素数:=选择(isprime,[1..floor(N/n0)]):
nTimes:=nops(素数):
因为我从1到9岁
将j从1改为nTimes do
M[i,j]:=转换(素数[1..min(j,i)],`*`);
日期:
A: ={}:
对于i9,从9到nTimes do
m9:=素数[i9];
对于select中的i8(t->M[7,t-1]*素数[t]*m9<=N,[$8..i9-1])do
m8:=m9*素数[i8];
对于选择中的i7(t->M[6,t-1]*素数[t]*m8<=N,[$7..i8-1])do
m7:=m8*素数[i7];
对于选择中的i6(t->M[5,t-1]*素数[t]*m7<=N,[$6..i7-1])do
m6:=m7*素数[i6];
对于选择中的i5(t->M[4,t-1]*素数[t]*m6<=N,[$5.i6-1])do
m5:=m6*素数[i5];
对于选择中的i4(t->M[3,t-1]*素数[t]*m5<=N,[$4..i5-1])do
m4:=m5*素数[i4];
对于选择中的i3(t->M[2,t-1]*素数[t]*m4<=N,[$3.i4-1])do
m3:=m4*底漆[i3];
对于选择中的i2(t->M[1,t-1]*素数[t]*m3<=N,[2..i3-1])do
m2:=m3*底漆[i2];
对于select中的i1(t->素数[t]*m2<=N,[$1..i2-1])do
A: =并集{m2*素数[i1]};
ododododod ododod: