%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日