%I#2 2012年3月30日17:31:20
%编号:8,454123918389913950854046429416573622929292956573,
%电话:45851698382473238120286488029224195550280826966354
%第10^N-第3个几乎素数。
%t ThreeAlmostPrimePi[n_]:=总和[PrimePi[n/(素数@i*素数@j)]-j+1,{i,PrimePi[n^(1/3)]},{j,i,PricePi@Sqrt[n/素数@i]}]; ThreeAlmostPrime[n_]:=块[{e=楼层[Log[2,n]],a,b},a=2^e;做[b=2^p;而[ThreeAlmostPrimePi[a]<n,a=a+b];a=a-b/2,{p,e,0,-1}];a+b/2];执行[打印@ThreeAlmostPrime[10^n],{n,0,13}]
%t三素数[n_]:=块[{e=楼层[Log[2,n]+2],a,b},a=2^e;Do[b=2^p;While[ThreePrimePi[a]<n,a=a+b];a=a-b/2,{p,e,0,-1}];a+b/2];表[ThreePrime[n],{n,0,13}]
%Y参见A109251、A006988、A114125、A120045、A120046。
%K更多,nonn
%0、1
%A _Robert G.Wilson v_,2006年2月15日
|