%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 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(rNeeded)!=搜索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的更多条款
|