第一个(n)={my(res=向量(n,i,oo),todo=n-1);res[1]=1;对于(i=2,oo,c=最大功率(i);如果(c<=n&&res[c]==oo,待办事项--;res[c]=i; 如果(todo==0,返回(res))))}最大功率(n)={my(f=因子(n),res=oo);对于(i=1,#f~,res=最小值(res,val(n,f[i,1])\f[i,2])); 物件}addhelp(val,“n中素数p的指数!”)val(n,p)=我的(r=0);而(n,r+=n \=p);第页