登录
无平方数的幂。
124

%I#38 2024年8月19日02:20:45

%S 1,2,3,4,5,6,7,8,9,10,11,13,14,15,16,17,19,21,2,23,25,26,27,29,30,31,

%电话:32、33、34、35、36、37、38、39、41、42、43、46、47、49、51、53、55、57、58、59、61、62、64,

%U 65、66、67、69、70、71、73、74、77、78、79、81、82、83、85、86、87、89、91、93、94、95、97

%N无平方数的幂。

%C a(n)=A072775(n)^ A072776(n);A059404的补充。

%C基本上与A062770相同_R.J.Mathar,2008年9月25日

%C数m,使得在正则素因式分解中,所有素数指数都相同:A124010(m,k)=A124010_Reinhard Zumkeller,2014年4月6日

%均匀分区的C-Heinz数。如果所有部分都以相同的多重性出现,则整数分区是一致的。整数分区的Heinz数(y_1,…,y_k)是素数(y_1)**素数(y_k).-_Gus Wiseman_,2018年4月16日

%H Reinhard Zumkeller,n的表,n=1..10000的a(n)</a>

%t选择[Range[100],Length[Union[FactorInteger[#][[All,2]]]==1&](*_Geoffrey Criter_,2015年3月30日*)

%o(哈斯克尔)

%o导入数据。映射(空,findMin,deleteMin,insert)

%o导入合格数据。地图。Lazy as Map(空)

%o a072774 n=a072774_列表!!(n-1)

%o(a072774_list、a072775_list和a072776_list)=解压缩3$

%o(1,1,1):f(尾部a005117_list)为空,其中

%o f vs'@(v:vs)m

%o|Map.null m||xx>v=(v,v,1):

%o f vs(插入(v^2)(v,2)m)

%o |否则=(xx,bx,ex):

%o f vs'(插入(bx*xx)(bx,ex+1)$删除最小值m)

%o其中(xx,(bx,ex))=findMin m

%o——Reinhard Zumkeller,2014年4月6日

%o(PARI)是(n)=功率(n,,&n);2015年10月16日,免发行(n)

%o(Python)

%o从数学导入isqrt

%o来自sympy import mobius,integer_nthroot

%o定义A072774(n):

%o定义g(x):返回int(范围(1,isqrt(x)+1)中k的总和(mobius(k)*(x//k**2))-1

%o定义f(x):对于范围(1,x.bit_length())中的k,返回n-2+x-sum(g(integer_ntroot(x,k)[0])

%o kmin,kmax=1,2

%o,而f(kmax)>=kmax:

%o kmax<<=1

%o为True时:

%o kmid=kmax+kmin>>1

%o如果f(kmid)<kmid:

%o kmax=公里ID

%o其他:

%o kmin=kmid

%o如果kmax-kmin<=1:

%o中断

%o 2024年8月19日归还kmax#_Chai Wah Wu_

%Y参见A072777(子序列)、A005117、A072778、A329332(表格安排)。

%Y A242414的子序列。

%Y参见A000009、A000837、A007916、A047966、A052409、A052410、A072774、A078374、A289023、A289509、A300486、A302491、A302796、A302979。

%K nonn公司

%O 1,2号机组

%A _Reinhard Zumkeller,2002年7月10日