nextp:=proc(P)局部k,m,newP,PL,PG,iv,i;
m: =nops(P);
从m-1到-1 do的k
如果P[k]>P[k+1],则
PL,PG:=选择删除(`<`,P[k+1..m],P[k]);
iv:=最大[指数](PL);
return[op(P[1..k-1]),PL[iv],op(sort([op(subsop(iv=P[k],PL)),op(PG)],`>`))]
fi(菲涅耳)
操作系统
结束进程:
五: =数组(0..21):V[0]:=1:
P: =[seq(i$2,i=9..0,-1)]:计数:=1:
当计数<3 do时
x: =加(P[i]*10^(19-i),i=1..19):
w: =理论数量:-bigomega(x);
如果w<=2且V[w]=0,则V[w]:=x;计数:=计数+1;fi;
P: =下一个P(P);
日期:
P: =[seq(i$2,i=9..0,-1)]:
当计数<22时
x: =加(P[i]*10^(20-i),i=1..20):
w: =理论数量:-bigomega(x);
如果w<=21且V[w]=0,则V[w]:=x;计数:=计数+1;fi;
P: =nextp(P);
日期:
转换(V,列表);
|