登录
A277802型
最小k>0,这样k*A004709号(n) 是一个立方体。
2
1, 4, 9, 2, 25, 36, 49, 3, 100, 121, 18, 169, 196, 225, 289, 12, 361, 50, 441, 484, 529, 5, 676, 98, 841, 900, 961, 1089, 1156, 1225, 6, 1369, 1444, 1521, 1681, 1764, 1849, 242, 75, 2116, 2209, 7, 20, 2601, 338, 2809, 3025, 3249, 3364, 3481, 450, 3721, 3844
抵消
1,2
评论
这是立方数的排列(A004709号).
a(n)<=A004709号(n) ^2,带等式iffA004709号(n) 是自由的-罗伯特·伊斯雷尔2016年11月9日
链接
配方奶粉
a(n)=A048798号(A004709号(n) )。
和{k=1..n}a(k)~c*zeta(3)^3*n^3/3,其中c=Product_{p素数}(1-1/p^2+1/p^5-1/p^6)=0.36052971192705404983-阿米拉姆·埃尔达尔2024年2月20日
例子
a(8)=3,因为3*A004709号(8) = 3 * 9 = 3^3.
a(16)=12,因为A004709号(16) = 18 = 2^1 * 3^2. 使k*2^1*3^2成为立方体的最小k是2^(3-(1mod3))*3^(3-(2mod3),)=12-大卫·A·科内斯2016年11月1日
MAPLE公司
f: =proc(n)局部f,E;
F: =系数(n)[2];
E: =F[..,2];
如果max(E)>=3,则返回NULL fi;
mul(F[i,1]^(3-E[i]),i=1..nos(F));
结束过程:
地图(f,[1..1000]美元)#罗伯特·伊斯雷尔2016年11月9日
数学
表[k=1;当[!IntegerQ[(k#)^(1/3)],k++]&@#[n]];k、 {n,53}]&@Select[Range[10^4],FreeQ[FactorInteger@#,{_,k_/;k>2}]&](*迈克尔·德弗利格2016年11月1日之后简·曼加尔丹A004709号*)
f[p_,e_]:=如果[e>2,0,p^(模态[-e,3])];s[n_]:=倍@@f@@FactorInteger[n];选择[Array[s,100],#>0&](*阿米拉姆·埃尔达尔2024年2月20日*)
黄体脂酮素
(PARI)\\大约n个术语的列表(可能性稍大一些)。
lista(n)={n=ceil(1.21*n);my(l=List([1]),f);forprime(p=2,n,for(i=1,#l,if(l[i]*p<=n,listput(l,l[i]*p);if(l[i]*p^2<=n,listput(l,l[i]*p^2)));listsort(l);for(i=2,#l,f=因子(l[i]);f[,2]=向量(#f[,2],i,3-(f[i,2]%3))~;l[i]=因子返回(f));l}\\大卫·A·科内斯2016年11月1日
(Python)
从数学导入prod
来自sympy import mobius,factorint,integer_nthroot
定义A277802型(n) :
定义f(x):对于范围(1,integer_ntroot(x,3)[0]+1)中的k,返回n+x-sum(mobius(k)*(x//k**3))
m、 k=n,f(n)
而m!=克:
m、 k=k,f(k)
factorint(m).items()中p,e的return prod(p**(-e%3))#柴华武2024年8月5日
关键词
非n,
作者
彼得·卡吉2016年10月31日
状态
经核准的