%I M0510#43 2024年8月2日04:02:32

%S 1,2,3,4,5,6,9,8,7,10,17,12,33,18,11,16,65,14129,20,19,34257,24,13，

%电话：66,15,36513,221025,32,35130,21,82049258,67,40497,388193,68，

%电话：2351416385,48,25,2613113232769,30,37,72259102665537,441310732050,39,64

%N Doudna序列A005940的逆序列。

%C a（2^k）=2^k.-罗贝特·威尔逊v_，2005年2月22日

%C固定点：A029747.-_Reinhard Zumkeller，2006年8月23日

%C问题：有没有一个简单的证据证明a（C）=C永远不会允许奇复合C作为解？另见A364551_Antti Karttunen_，2023年7月30日

%D J.H.Conway，个人沟通。

%D N.J.A.斯隆和Simon Plouffe，《整数序列百科全书》，学术出版社，1995年（包括该序列）。

%H R.J.Mathar，n的表，a（n）表示n=1，。。，5000个</a>

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

%F a（n）=h（g（n，1，1），0）/2+1，其中h（n，m）=如果n=0，则m else h（floor（n/2），2*m+n mod 2），g（n、i，x）=如果n=1，则x else（如果n mod prime（i）=0，那么g（n/prime（i_Reinhard Zumkeller_，2006年8月23日

%F a（n）=1+A156552（n）.-_Antti Karttunen，2014年6月26日

%p A005941:=程序（n）

%p局部k；

%从1do到k的p

%p如果A005940（k）=n，则#代码重用

%p返回k；

%p end if；

%p端do；

%结束程序：#_R.J.Mathar_，2010年3月6日

%t f[n_]：=块[{p=Partition[Split[Join[Integer Digits[n-1，2]，{2}]]，2]}，Times@@Flatten[Table[q=Take[p，-i]；素数[Count[Flatten[q]，0]+1]^q[[1，1]]，{i，Length[p]}]]；t=表格[f[n]，{n，10^5}]；扁平[表[位置[t，n，1，1]，{n，64}]]（*_Robert G.Wilson v_，2005年2月22日*）

%o（方案）（定义（A005941 n）（+1（A156552 n）））_Antti Karttunen，2014年6月26日

%o（Python）

%o来自sympy import primepi，factorint

%o def A005941（n）：返回和（（1<<素数pi（p）-1）<<i for i，p in enumerate（factorint（n，multiple=True））+1#_Chai Wah Wu_，2023年3月11日

%o（PARI）A005941（n）={my（f=因子（n），p，p2=1，res=0）；对于（i=1，#f~，p=1<<（素数（f[i，1]）-1）；res+=（p*p2*（2^（f[i,2]）-1）））；p2<<=f[i、2]）；（1+res）}；\\（继David A.Corneth的A156552计划之后）-安蒂·卡图内，2023年7月30日

%Y参考A103969。A005940的逆运算。比A156552多一个。

%Y参见A364559[=a（n）-n]，A364557（莫比乌斯变换），A3645508。

%Y参见A029747[其中a（n）=n的已知位置]，A364560[其中a。

%K nonn公司

%O 1,2号机组

%A _N.J.A.斯隆_

%E更多条款摘自2005年2月22日_Robert G.Wilson v_

%E a（61）由R.J.Mathar_插入，2010年3月6日