%I#39 2022年9月8日08:45:07
%S 1,9,25,27,49,81121125169225243289343361441529625729841,
%电话9611089122513311369152116811849202521872197220924012601,
%电话2809302531253249337534813769422544894761491350 4153295625
%N次奇完全幂(1加上数字m^k,m奇数,k>=2)。
%H Reinhard Zumkeller,n的表,n=1..10000的a(n)</a>
%F总和{n>=1}1/a(n)=1+总和{k>=2}mu(k)*(1-泽塔(k)x(2^k-1)/2^k)=1.2890375574…-阿米拉姆·埃尔达尔,2020年12月19日
%pq:=n->n=1或n::奇数和igcd(seq(i[2],i=ifactors(n)[2]))>1:
%p选择(q,[1.6000]美元)[];#_阿洛伊斯·海因茨,2022年5月4日
%t Take[Union[Flatten[Table[a^b,{a,1,99,2},{b,2,15}]],50](*_Alonso del Arte_2011年11月22日*)
%o(哈斯克尔)
%o a075109 n=a075109_列表!!(n-1)
%o a075109_list=奇数过滤器a001597_list--_Reinhard Zumkeller_,2012年10月4日
%o(Magma)[1]cat[n:n in[3..6000 by 2]| IsPower(n)];//_文森佐·利班迪(Vincenzo Librandi),2014年3月31日
%o(PARI)isok(m)=(m==1)||((m%2)&&ispower(m));\\_米歇尔·马库斯,2022年5月4日
%A001597和A005408的Y交点。
%Y参考A008683、A075090。
%K容易,不是
%O 1,2号机组
%A _Zak Seidov,2002年10月11日
%E定义由N.J.A.Sloane澄清,2021年12月25日
|