（PARI）数字计数（n，基数=10）=本地（d）；如果（n==0，返回（1））；d=1+地板（对数（n）/对数（底座））；而（n>=基数^d，d++）；而（n<base^（d-1），d-）；d；
{
a（n）=
本地（r，num2，num5，d，M，pLeft，mainP，searchP，fixed，x，rNeeded，y，nextP）；
r=n；
d=数字计数（n）；
while（num2<d&&！（r%2），
num2++；
r=r/2
);
while（num5<d&&！（r%5），
数字5++；
r=r/5
);
M=10^d/2^num2/5^num5；
p左=n-num2-num5；
mainP=如果（num2==d，2，3）；
搜索P=最小值（4，pLeft）；
固定=2^num2*5^num5；
x=mainP^（pLeft-搜索P）；
r需要=升力（Mod（r，M）/Mod（x，M））；
while（bigomega（需要）！=搜索P，
r需要+=M
);
y=固定*x*r需要；
如果（mainP==2，
下一个P=3，
nextP=if（数字5==d，5，7）
);
while（searchP<pLeft&&fixed*x*nextP^（1+searchP）/mainP<y，
搜索P++；
x/=主P
);
r需要=升力（Mod（r，M）/Mod（x，M））；
while（bigomega（需要）！=搜索P，
r需要+=M
);
return（固定*x*rNeeded）；
}（结束）
