%I#10 2014年8月27日14:54:09

%S 1,2,4,3,6,7,8,5,9,13,10,17,20,19,12,11,46,23166,41,15,29858,59,14，

%电话：71,16,676186,3758645,31,18199,22,83705348983,3217910428487,47，

%电话：184718194109、2166593840230006277、27、43、65353

%N自然数的置换：a（N）=A227413（A005941（N））。

%H<a href=“/index/Per#IntegerPermutation”>自然数排列序列的索引项</a>

%F a（n）=A227413（A005941（n））=A227413（1+A156552（n））。

%F其他身份：

%F对于所有n>=1，A010051（a（n））=1-A000035（n）。[此置换以某种顺序将偶数映射到素数，将奇数映射到非素数，因为置换A227413具有相同的属性，而A005941保留奇偶性]。

%o（方案）（定义（A246368 n）（A227413（A005941 n））

%o（PARI）

%o默认（质数限制，（2^31）+（2^30））；

%o A002808（n）={my（k=-1）；while（-n+n+=-k+k=素数（n），）；n}；\\此函数来自_M.F.Hasler_

%o A064989（n）={my（f）；f=因子（n）；如果（n>1&&f[1,1]==2），f[1,2]=0）；对于（i=1，#f~，f[i，1]=前一素数（f[i、1]-1））；因子回退（f）}；

%o A156552（n）=如果（1==n，0，如果（！（n%2），1+（2*A156552（n/2）），2*A156552（A064989（n）））；

%o A005941（n）=A156552（n）+1；

%o A227413（n）=如果（1==n，1，如果（！（n%2），素数（A227412（n/2）），A002808（A2274013（（n-1）/2）））；

%o A246368（n）=A227413（A005941（n））；

%o表示（n=1,52，写入（“b246368.txt”，n，“”，A246368（n））；

%Y反转：A246367。

%Y类似或相关排列：A005941、A156552、A227413、A246364、A246365。

%Y参考A000035、A010051。

%K nonn公司

%O 1,2号机组

%2014年8月26日，安蒂·卡图内