%I#18 2024年2月26日01:22:47
%S 1,2,5,4,3,10,7,8,25,6,13,20,11,14,15,16,17,50,19,12,35,26,23,40,9,22,
%电话:125,28,29,30,31,32,65,34,21100,37,38,55,24,41,70,43,52,75,46,47,80,
%U 49、18、85、44、53250、39、56、95、58、61、60、59、62175、64、33
%交换N的素因式分解中形式为(4*k+3,4*(k+1)+1)的孪生素对。
%H Amiram Eldar,n的表,n=1..10000的a(n)</a>
%H<a href=“/index/Per#IntegerPermutation”>为自然数置换序列的索引条目</a>。
%F与a(p)=相乘(如果p mod 4=3且p+2是素数,则p+2其他(如果p mod 4=1且p-2是素数则p-2其他p),则p为素数。
%F a(a(n))=n,自然数的自反转排列。
%F和{k=1..n}a(k)~c*n^2/2,其中c=Product_{(p<q)交换对}((p^2-p)*(q^2-q)/_Amiram Eldar,2024年2月26日
%e a(42)=a(2*3*7)=a。
%t a[n_]:=乘积[{p,e}=pe;哪个[
%t Mod[p,4]==3&&PrimeQ[p+2],p+2,
%t Mod[p,4]==1&&PrimeQ[p-2],p-2,
%t True,p]^e,{pe,FactorInteger[n]}];
%t数组[a,100](*_Jean-François Alcover_,2021年11月21日*)
%o(PARI)a(n)={my(f=因子(n));prod(i=1,#f~,p=f[i,1]);if(p==2,p,if(p%4==3&isprime(p+2),p+2,if
%Y参见A061898、A064505、A071698、A0711699、A072026、A072027、A072028。
%K nonn,多个
%O 1,2号机组
%A _Reinhard Zumkeller,2002年6月7日
|