A340585-OEIS公司

A340585型

非立方体完美力量。

4, 9, 16, 25, 32, 36, 49, 81, 100, 121, 128, 144, 169, 196, 225, 243, 256, 289, 324, 361, 400, 441, 484, 529, 576, 625, 676, 784, 841, 900, 961, 1024, 1089, 1156, 1225, 1296, 1369, 1444, 1521, 1600, 1681, 1764, 1849, 1936, 2025, 2048, 2116, 2187, 2209, 2304, 2401, 2500

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

抵消

1,1

这是对A239870型然而，该序列的真实定义似乎略有不同。

链接

雨果·普福尔特纳，n=1..10000时的n，a（n）表

配方奶粉

求和{n>=1}1/a（n）=1-zeta（3）+求和{k>=2}mu（k）*（1-zeta（k））=1-A002117号+A072102号= 0.6724074652... -阿米拉姆·埃尔达尔2021年1月12日

MAPLE公司

过滤器：=proc（n）局部g；

g： =igcd（op（ifactors（n）[2][..，2]））；

g>1和（g mod 3<>0）

结束进程：

选择（筛选器，[$1..10000]）#罗伯特·伊斯雷尔2021年1月12日

数学

选择[Range[2，2500]，（g=GCD@@FactorInteger[#][[；；，2]]）>1&&！可分[g，3]&](*阿米拉姆·埃尔达尔，2021年1月12日*）

黄体脂酮素

（PARI）对于（n=22500，如果（ispower（n）%3，print1（n，“，”））

（Python）

从数学导入isqrt

来自sympy import mobius，integer_nthroot

定义A340585型（n）：

定义f（x）：返回int（n+x-isqrt（x）+和（mobius（k）*（integer_nthroot（x，k）[0]-1），用于范围（5，x.bit_length（））中的k））

kmin，kmax=1,2

而f（kmax）>=kmax：

kmax≤1

为True时：

kmid=kmax+kmin>>1

如果f（kmid）<kmid：

kmax=kmid

其他：

kmin=kmid

如果kmax-kmin<=1：

打破

返回kmax#柴华武2024年8月14日

交叉参考

囊性纤维变性。A000578号,A001597号,A002117号,A007412号,A072102号,A076467号,A097054美元,A239728型.

上下文中的序列：A331219 A337534型 A126589号*A010409号 A010457号 A244833号

相邻序列：A340582型 A340583型 A340584型*A340586型 A340587 A340588型

关键词

非n

作者

雨果·普费尔特纳2021年1月12日

状态

经核准的