%I#10 2022年12月1日11:04:01
%S 1,4,9,2,25,36,49,8,3100121,18169196225,32289,12361,50441484,
%电话:529,72,5676,27,98841900961,1610891156125,613691441521,
%电话:200168117641849242,7521162209288,7,20260133828091083025392
%N a(N)是序列中尚未出现的最小正整数,因此N*a(N”)是一个立方体。
%这是正整数的自反转置换。
%H Amiram Eldar,n的表,n的a(n)=1..10000</a>
%F与a(p^{3n})=p^{2n},a(p*3n+1})=p*3n+2),以及a(p*3n+2))=p*3n+1)的乘法。
%F和{k=1..n}a(k)~c*n^3,其中c=(zeta(6)/(3*zeta(3)))*Product_{p素数}(1-1/p^2+1/p^3)=A347328*A330596/3=0.211705….-_Amiram Eldar,2022年12月1日
%t f[p_,e_]:=如果[(r=模式[e,3])==0,p^e,p^(e-(-1)^r)];a[n_]:=倍@@(f@@@FactorInteger[n]);阵列[a,100](*_Amiram Eldar_,2019年8月29日*)
%o(PARI)a(n)={my(f=因子(n),r);prod(i=1,#f~,r=f[i,2]%3;f[i、1]^if(r==0,f[i;2],f[i,2]-(-1)^r));}\\阿米拉姆·埃尔达尔,2022年12月1日
%Y参考A064429(指数函数)
%Y参见A330596、A347328。
%K mult,非n
%O 1,2号机组
%A _弗兰克林·T·亚当斯-沃特斯,2009年4月7日
|