f: =proc(n)使用priqueue;局部t,x,p,i;
初始化(pq);
插入([-3^n,3$n],pq);
做
t: =提取物(pq);
x: =-t[1];
如果numtheory:-bigomega(x-1)=n,则返回x-1
elif numtheory:-bigomega(x+1)=n,然后返回x
fi;
p: =下一素数(t[-1]);
当t[i]=t[-1]do时,i从n+1到2乘-1
插入([t[1]*(p/t[-1])^(n+2-i),op(t[2..i-1]),p$(n+2-i)],pq)
od;
日
结束进程: