|
程序来自大卫·瓦瑟曼2008年9月19日：（开始）
（PARI）位数（n，基数=10）=局部（d）；如果（n==0，返回（1））；基底（d+1）/log）；而（n>=基^d，d++）；而（n<base^（d-1），d--）；d；
{
a（n）=
本地（r、num2、num5、d、M、pLeft、mainP、searchP、固定、x、rNeeded、y、nextP）；
r=n；
d=数字计数（n）；
当（num2<d&！（r%2）时，
num2++；
r=r/2
);
当（num5<d&！（r%5）时，
num5++；
r=r/5
);
M=10^d/2^num2/5^num5；
pLeft=n-num2-num5；
mainP=if（num2==d，2，3）；
searchP=最小值（4，pLeft）；
固定=2^num2*5^num5；
x=mainP^（pLeft-searchP）；
rNeeded=升力（Mod（r，M）/Mod（x，M））；
while（bigomega（needed）！=搜索P，
rNeeded+=米
);
y=固定*x*rNeeded；
如果（mainP==2，
nextP=3，
nextP=if（num5==d，5，7）
);
而（searchP<pLeft&&fixed*x*nextP^（1+searchP）/mainP<y，
searchP++；
x/=主管道
);
rNeeded=升力（Mod（r，M）/Mod（x，M））；
while（bigomega（needed）！=搜索P，
rNeeded+=米
);
返回（固定*x*rNeeded）；
}（结束）
| 