%I#39 2022年2月10日22:01:42

%S 0,3,342992625233782100351904324174272581607885361493776443，

%电话1395999034213112601717812370880486531171715902308080，

%电话：111329817298881106105729282726910139482913717352971123037685177087932300026230174178896660584964121902286389956293761485464

%N a（N）=半素数<10^N。

%C除第一个非零项外，序列与A036352相同_雨果·普福尔特纳，2003年7月22日

%H Eric Weisstein的数学世界，<a href=“http://mathworld.wolfram.com/Semiprime.html“>半素数</a>

%H<a href=“/index/Pri#primepop”>为与不同范围内的素数数量相关的序列索引条目</a>

%F（1/2）*（pi（10^（n/2））+求和{i=1..pi（10 ^n）}pi（（10 ^n-1）/P-i））=求和{i=1..pi（sqrt（10^n））}π（（10^n-1_Robert G.Wilson v_，2005年5月16日

%在10以下有三个半素数：4（2*2），6（2*3）和9（3*3）。

%t f[n_]：=和[PrimePi[（10^n-1）/Prime[i]]，{i，PrimePi[Sqrt[10^n]]}]-二项式[PrimePi[Sqrt[10 ^n]]，2]；Do[打印[f[n]]，{n，0，14}]（*_Robert G.Wilson v_，2005年5月16日*）

%t半素数Pi[n_]：=和[PrimePi[n/Prime@i]-i+1，{i，PrimePi@Sqrt@n}]；数组[SemiPrimePi[10^#-1]&，14，0]（*_Robert G.Wilson v_，2015年1月21日*）

%o（PARI）a（n）=我的；对于素数（p=2，sqrt（10^n），s+=primepi（（10^n-1）\p））；s-二项式（素数（sqrt（10^n）），2）\\-Charles R Greathouse IV_，2012年4月23日

%o（Perl）使用数学：：Prime:：Util qw/：all/；使用整数；sub-countsp{my（$k，$sum，$pc）=（$_[0]-1,0,1）；prime_precalc（60_000_000）；对于primes{$sum+=prime_count（$k/$_）+1-$pc++；}int（sqrt（$k））$总和；}foreach my$n（0..16）{say“$n:”，countsp（10**$n）；}#_Dana Jacobsen_，2014年5月11日

%Y参见A001358、A064911、A072000、A036352（从a（2）开始相同）、A220262、A292785。

%K非n

%0、2

%A _帕特里克·德·吉斯特，2001年12月10日

%E更多术语摘自Hugo Pfoertner，2003年7月22日

%E a（14）摘自_Robert G.Wilson v_，2005年5月16日

%E a（15）-a（16）摘自Donovan Johnson，2010年3月18日

%E a（17）-a（18）摘自_Dana Jacobsen，2014年5月11日

%E a（19）-a（21），摘自HHerri Lifchitz，2015年7月4日