%I#38 2022年10月2日10:33:30

%S 1,2,3,4,5,6,9,8,7,10,15,12,25,18,27,16,11,14,21,20,35,30,45,24,49,50，

%电话：75,36125,54,81,32,13,22,33,28,55,42,63,40,77,70105,60175,90135，

%U 48,99,98147100245150225,72343250375108625162243,64,17单位

%N如果N是2的幂，a（N）=N。否则，如果2^j是2的最大幂，但不超过N，如果k=N-2^j，那么a（N）是尚未发生的最小m*a（k），其中m是奇数。

%这个定义意味着，如果n不是2的幂，那么也不是a（n）。

%C与Doudna序列（A005940）类似，只是这里用于计算a（n）的a（k）的倍数是最小奇数（而不是最小奇数素数），因此a（n）是一个新术语。术语与A005940中的术语相同，直到a（49）=99（而不是121），然后再加上奇数非素数倍数，与A00594产生更多差异；其次是a（71）=117（而不是99）。

%C正整数的一种置换，其中素数按自然顺序出现，而奇数则不按自然顺序（9先于7，25先于21，以此类推）。

%H Michael De Vlieger，n的表，n的a（n）=1..16384</a>

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

%Fa（2^n+1）是尚未出现的最小奇数。

%e n=49=2 ^ 5+17，a（17）=11，因此a（49）是之前没有出现的最小m*a（17，其中m是奇数。到目前为止，我们已经看到了3*11、5*11、7*11，但没有看到9*11。因此，a（49）=9*11=99（与A005940（71）=99相比）。

%t nn=65；m=1；c[_]=错误；做[Set[{m，k}，{1，n-2^Floor[Log2[n]]}]；如果[k==0，集[{a[n]，c[n]}，{n，True}]，而集[t，ma[k]]；或[m==1，c[t]]，m+=2]；设置[｛a[n]，c[t]｝，｛t，True｝]]，｛n，nn｝]；数组[a，nn]（*_迈克尔·德弗里格，2022年9月21日*）

%o（Python）

%o从sympy导入nextprime

%o来自症状合成导入数字

%o从itertools导入计数，islice

%o定义b（n）：返回n-2**（len（bin（n）[2:]）-1）

%o定义代理（）：

%o aset，alst=set（），[无]

%计数（1）中n的o：

%o k=b（n）

%o如果k==0:an=n

%o其他：

%o ak，p=阿尔斯特[k]，3

%o，而p*ak在aset:p+=2

%o an=p*ak

%o产量；增加（a）；附加（an）

%o打印（列表（islice（agen），65））#_Michael S.Branicky_，2022年9月21日

%o（PARI）f（n）=n-2^（logint（n，2））；\\A053645号

%o列表（nn）={my（va=向量（nn，sa=集合（va））；对于（n=1，nn，my（x=f（n）））；如果（x==0，va[n]=n，my

%Y参考A005940、A356867、A356886。

%Y参考A053644、A053645。

%K nonn公司

%O 1,2号机组

%2022年9月21日，一位狂热的詹姆斯·西卡莫雷（A _David James Sycamore）