%I#5 2012年3月31日10:26:37

%编号：11,22,63,24405,96218714081246596561042712113891237614213，

%电话：4050761432683621580281610010754117249670981823202182619，

%电话：14417920886805499321767092568221474909801623251658243922632451125

%N以N的数字结尾的最小数字，正好有N个素因子（以重数计算）。

%e a（5）=405，因为405=3*3*3x3*5（5个因子），以5结尾，是最小的数字

%o David Wasserman的计划，2008年9月19日：（开始）

%o（PARI）数字计数（n，基数=10）=本地（d）；如果（n==0，返回（1））；d=1+地板（对数（n）/对数（底座））；而（n>=基数^d，d++）；而（n<基数^（d-1），d--）；d；

%o个{

%o a（n）=

%o本地（r，num2，num5，d，M，pLeft，mainP，searchP，fixed，x，rNeeded，y，nextP）；

%o r=n；

%o d=数字计数（n）；

%o while（num2<d&&！（r%2），

%o数量++；

%o r=r/2

%o）；

%o while（num5<d&&！（r%5），

%o数字5++；

%o r=r/5

%o）；

%o M=10^d/2^num2/5^num5；

%o pLeft=n-num2-num5；

%o mainP=如果（num2==d，2，3）；

%o搜索P=最小值（4，pLeft）；

%o固定=2^num2*5^num5；

%o x=mainP^（pLeft-搜索P）；

%o r所需=升力（Mod（r，M）/Mod（x，M））；

%o while（bigomega（需要）！=搜索P，

%o r需要+=M

%o）；

%o y=固定*x*r需要；

%o如果（mainP==2，

%o下一个P=3，

%o nextP=if（数字5==d，5，7）

%o）；

%o while（searchP<pLeft&&fixed*x*nextP^（1+searchP）/mainP<y，

%o搜索P++；

%o x/=主P

%o）；

%o r所需=升力（Mod（r，M）/Mod（x，M））；

%o while（bigomega（需要）！=搜索P，

%o r需要+=M

%o）；

%o返回（固定*x*rNeeded）；

%o}（结束）

%Y参考A109665[摘自David Wasserman，2008年9月30日]

%K基，nonn

%O 1，1

%A _弗里德曼，2005年8月7日

%E 2008年9月18日来自_David Wasserman的更多条款