%I#12 2024年7月23日18:53:29

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

%电话：4050761432683621580281610010754117249670981823202182619，

%电话：14417920886805499321767092568221474909801623251658243922632451125

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

%H Robert Israel，n的表，n的a（n）=1..2078</a>

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

%pf:=proc（n）使用priqueue；

%p局部pq，d，td，t，v，p，i，w2，w5，more2，more5；

%p初始化（pq）；

%p d：=1+ilog10（n）；td:=10^d；

%pw2：=padic:-ordp（n，2）；

%p更多2:=（w2>=d）；

%pw2：=最小值（w2，d）；

%pw5:=padic:-ordp（n，5）；

%p更多5:=（w5>=d）；

%pw5：=最小值（w5，d）；

%p如果更多2，则插入（[-5^w5*2^（n-w5），2，n-w5-w2]，pq）

%p其他插入（[-2^w2*5^w5*3^（n-w5-w2），3，n-w5-w2]，pq）

%p fi；

%p do公司

%p t：=提取物（pq）；

%p v：=-t[1]；

%p如果v mod td＝n，则返回v fi；

%p:=下一素数（t[2]）；

%p如果p=5且不大于5，则p:=7fi；

%p代表i从1到t[3]do

%p插入（[t[1]*（p/t[2]）^i，p，i]，pq）；

%p od；

%p od；

%p端程序：

%p映射（f，[1..30]美元）；#_Robert Israel，2024年7月23日

%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（bigmomega（r需要）！=搜索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的更多条款