%I#27 2023年8月2日10:02:03

%S 1,2,3,4,5,6,7,8,9,10,12,14,16,18,20,24,28,32,36,40,48,56,64,72,80,96，

%电话：11212814416019222425628832038444485125766407688961024，

%电话：11521280153617922048230425603072358440964608512061447168

%N位置N，使A005941（N）=A005940（N）。

%n>=5的C序列似乎是用A000079（二次幂）的五次剖分倍数进行的四次剖分：对于m>0，a（5m）=5,10,20,40…=5*2^（m-1）；a（5m+1）=6,12,24,48，…=6*2^（m-1）；同样，a（5m+2）=7*2^（m-1）；a（5m+3）=8*2^（m-1）；a（5m+4）=9*2^（m-1）_拉尔夫·斯蒂芬，2010年11月13日

%C来自_Antti Karttune_，2023年8月1日：（开始）

%C数k，其中A005940（A005940-（k））=k，或相等，其中A00.5941。

%C如果n是一个项，那么也存在2*n，反之亦然。

%C问题：1、3、5、7和9是这个序列中唯一的奇数项吗？

%C（结束）

%H R.J.Mathar，n表，n=1..60的a（n）。

%F经验g.F.：x*（1+x+x^2+x^3+x^4）^2/（1-2*x^5）。-_科林·巴克（Colin Barker），2016年11月18日

%e 56的顺序为A005940（56）=A005941（56）=72。

%e 7的顺序是因为A005940（7）=9，A005940-（9）=7，因此A005941（7）=9，以及A005941--_Antti Karttunen，2023年8月1日

%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^4}]；u=压扁[表[位置[t，n，1，1]，{n，10^4}]]；做[如果[u[[n]]=={}，u[[n]]={0}]，{n，10^4}]；压扁[位置[Take[t，10^4]-压扁[u]，0]]

%o（Python）

%o来自math导入prod

%o从itertools导入累加、计数、岛屿

%o来自sympy import prime，primepi，factorint

%o从集合导入计数器

%o定义A103969_gen（startvalue=1）：#术语生成器>=startvalue

%o返回过滤器（lambda n:sum（（1<<primepi（p）-1）<<i代表i，p代表枚举（factorint（n，multiple=True）））+1==prod（prime（len（a）+1）**b代表a，b代表计数器（累加（bin（n-1）[2:].split（'1'）[：0:-1]））.items（）），count（max（startvalue，1））））

%o A103969_list=列表（岛屿（A103969-gen（），20））#_Chai Wah Wu_，2023年3月11日

%Y参见A005940、A005941、A029747（子序列）。

%K nonn公司

%O 1,2号机组

%A _Robert G.Wilson v _，2005年2月22日

%2010年3月6日，R.J.Mathar修正了定义并更新了示例