登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 

标志
提示
(来自的问候整数序列在线百科全书!)
A057521号 n的强大(1)部分:如果n=Product_i(pi^ei),则a(n)=Product_{i:ei>1}(pi^ei);如果n=b*c^2*d^3,那么当b最小化时,a(n)=c^2*d^3。 70
1, 1, 1, 4, 1, 1, 1, 8, 9, 1, 1, 4, 1, 1, 1, 16, 1, 9, 1, 4, 1, 1, 1, 8, 25, 1, 27, 4, 1, 1, 1, 32, 1, 1, 1, 36, 1, 1, 1, 8, 1, 1, 1, 4, 9, 1, 1, 16, 49, 25, 1, 4, 1, 27, 1, 8, 1, 1, 1, 4, 1, 1, 9, 64, 1, 1, 1, 4, 1, 1, 1, 72, 1, 1, 25, 4, 1, 1, 1, 16, 81, 1, 1, 4, 1, 1, 1, 8, 1, 9, 1, 4, 1, 1, 1, 32, 1 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,4
链接
安蒂·卡图恩,n=1..16383时的n,a(n)表(T.D.Noe的前1000个术语)
维克托·乌夫纳罗夫斯基和博奥伦德,如何区分数字,J.整数序列。,第6卷(2003年),第03.3.4条。
配方奶粉
a(n)=n/A055231号(n) 。
对于e>1,与a(p)=1和a(p^e)=p^e相乘-弗拉德塔·约沃维奇2001年11月1日
发件人安蒂·卡图恩,2017年11月22日:(开始)
a(n)=A064549号(A003557号(n) )。
A003557号(a(n))=A003557号(n) 。
(结束)
a(n)=gcd(n,A003415号(n) ^k),对于所有k>=2。[这个公式是由Christian Krause的LODA miner以k=3的形式发现的。关于为什么成立,请参见Ufnarovski和Åhlander的论文,第4页的定理5。]-安蒂·卡图恩2021年3月9日
Dirichlet g.f.:zeta(s-1)*产品{p素数}-阿米拉姆·埃尔达尔2023年9月18日
例子
a(40)=8,因为40=2^3*5,所以强大的部分是2^3=8。
枫木
A057521号:=进程(n)
局部a、d、e、p;
a:=1;
对于ifactors(n)[2]do中的d
e:=d[1];
p:=d[2];
如果e>1,则
a:=a*p^e;
结束条件:;
结束do:
返回a;
结束进程:#R.J.马塔尔2016年6月9日
数学
rad[n_]:=倍@@First/@FactorInteger[n];a[n_]:=n/分母[n/rad[n]^2];表[a[n],{n,1,97}](*让-弗朗索瓦·奥尔科弗2013年6月20日*)
f[p_,e_]:=如果[e>1,p^e,1];a[1]=1;a[n_]:=倍@@f@@FactorInteger[n];数组[a,100](*阿米拉姆·埃尔达尔2020年9月21日*)
黄体脂酮素
(PARI)a(n)=我的(f=系数(n));prod(i=1,#f~,if(f[i,2]>1,f[i;1]^f[i、2],1))\\查尔斯·格里特豪斯四世2013年8月13日
(PARI)a(n)=我的(f=系数(n));对于(i=1,#f~,如果(f[i,2]==1,f[i、1]=1));因子回收(f)\\米歇尔·马库斯2021年1月29日
(Python)
来自症状输入因子int,prod
def a(n):如果n==1,则返回1 else prod(1 if e==1 else p**e for p,e in factorint(n).items()))
打印([a(n)代表范围(1,51)中的n])#因德拉尼尔·戈什2017年7月19日
(Python)
从数学导入prod
来自sympy导入因子
定义A057521号(n) :return n//prod(p代表p,e在因子(n).items()中,如果e==1)#柴华武2022年11月14日
交叉参考
关键字
非n,复数,容易的
作者
亨利·博托姆利2000年9月1日
状态
经核准的

查找|欢迎光临|维基|注册|音乐|地块2|演示|索引|浏览|更多|网络摄像头
贡献新序列。或评论|格式|样式表|变换|超级搜索|最近
OEIS社区|维护人OEIS基金会。

许可协议、使用条款、隐私政策。.

上次修改时间:美国东部夏令时2024年6月23日22:09。包含373661个序列。(在oeis4上运行。)