A113877-OEIS公司

A113877号

半素数到半素数幂。

256, 1296, 4096, 6561, 10000, 38416, 46656, 50625, 194481, 234256, 262144, 390625, 456976, 531441, 1000000, 1048576, 1185921, 1336336, 1500625, 2085136, 2313441, 4477456, 5764801, 6765201, 7529536, 9150625, 10077696, 10556001, 11316496, 11390625, 14776336

(列表;图表;参考;听;历史;文本;内部格式)

抵消

1,1

这是的半素数模拟A053810号.

链接

查尔斯·格里塔斯四世，n=1..10000时的n，a（n）表

配方奶粉

{a（n）}={a^b其中a和b是A001358号}.

{a（n）}={（p*q）^（r*s）=（p^（r*s））*（q^r*s ^（p*q）=（p^（p*q））*（q^（p*q）。

a（n）~（n log n/log log n）^4。 -查尔斯·格里特豪斯四世2013年6月5日

例子

a（1）=256=4^4=半素数（1）^半素数。

a（2）=1296=6^4=半素数（2）^半素数。

a（3）=4096=4^6=半素（1）^半素（2）。

a（4）=6561=9^4=半素数（3）^半素数。

a（5）=10000=10^4。

数学

lim=10^8；s=选择[Range[lim^（1/4）]，Total[Transpose[FactorInteger[#]][2]]==2&]；t={}；j=1；而[b=s[[j]]；i=1；而[a=s[i]]；e=a ^b；如果[e<=lim，追加到[t，e]]；e<lim&&i<长度[s]，i++]；i>1，j++]；t=联管节[t](*T.D.诺伊2013年6月5日*）

黄体脂酮素

（PARI）是（n）=我的（b，e=功率（n，&b），o）；如果（e==0，返回（0））；o=大ω（e）；（o==2&&bigomega（b）==2）||（e%2==0&&o==3&&i素数（b））\\查尔斯·格里特豪斯四世2013年6月5日

（PARI）列表（lim）=我的（v=列表（））；对于（e=4，log（lim\=1+.5）\log（4），if（bigomega（e）！=2，下一个）；对于（b=4，（lim+.5）^（1/e），如果（bigomega（b）==2，listput（v，b^e）））；向量排序（Vec（v））\\查尔斯·格里特豪斯四世2013年6月5日

（Python）

从数学导入isqrt

从sympy导入primepi，primerange，integer_nthroot，factorint

定义A113877号（n）：

定义A072000型（n）：return int（-（（t:=primepi（s:=isqrt（n）））*（t-1）>>1）+总和

定义f（x）：返回int（n+x-sum(A072000型（integer_nthroot（x，p）[0]），对于范围（4，x.bit_length（））中的p，如果sum（factorint（p）.values（））==2））

定义平分（f，kmin=0，kmax=1）：

而f（kmax）>kmax:kmax<<=1

当kmax-kmin>1时：

kmid=kmax+kmin>>1

如果f（kmid）<=kmid：

kmax=kmid