|
%I#102 2021年1月9日11:06:21
%S 0,1,2,3,4,7,6,14,5,15,27,12,24,10,19,30,8,31,43,28,39,13,35,45,48,62，
%电话：20,57,37,63,60,79,9126,91,11,86,29,56,23,54,75,26,51,70,46,27,22,89，
%U 21,93,83,40,61114,78,38,18,71,87,77,42124127,16254187,92151,90,44,58117单位
%N a（N）是尚未使用的最小数字k，使得k的二进制表示中的1位数字等于k*N的二进制表示的1位数量。
%C我会把这个序列称为“邪恶的野兽”，因为它显示了许多模式，但对于每个模式，似乎都有一个n值，在这个值处规则突然改变或出现一些意想不到的异常。
%C如果n是2的幂，则任何数字都满足条件，因为1位的数量不会因乘以2的幂而改变。正因为如此，每个数字最终都有机会出现在这个序列中；这证明了这个序列是非负整数的置换。
%C此序列可能有助于找到b和C的小对，以便A000120（b）=A000120。
%C A340100中描述了所有固定点n=a（n），它们是A077436的子集。
%C在n=0..100000范围内，最大值a（n）为131072=a（32769），但在n=30000..40000范围内的最小值为a（32778）=137。
%C如果A000120（b）=A000120；这导致这个序列中的一些模式可能在n的有限范围内有效。我们能找到一个对n的大范围有效的模式吗？
%如果a（n）是2的幂，那么n也是2的幂。但如果n是2的幂，a（n）并不总是2的幂。
%C在A000120（C）=A000120（C*b）形式的方程中，对于所有A000120（C）=2，我们发现b的所有解为b=0，b=2^d或b=（2^d）*（1+2^（（（C-1）/2）+e*（C-1）））/C，如果C是奇数。对于偶数c，用2的最大可能幂除以c。c=3的示例是b=A263132。
%C a（n）>=A292849（n）。这个下限是这个序列中一些峰值的原因。
%H Thomas Scheurele，n的表，n=0..10000的a（n）</a>
%H Thomas Scheurele，<a href=“/A340069/A340069.svg”>此序列显示了一个极端混沌图。
%H<a href=“/index/Per#IntegerPermutation”>自然数排列序列的索引项</a>
%如果n<5，F a（n）=n。
%如果n<5，则F a（2^（2*n））=2^（1+n）。
%如果n<5，F a（2^（2*n+1））=2^（1+n）+1。
%如果n>0且<4，则F a（3*2^n）=3*2^（n+1）。
%o（MATLAB）
%o函数a=A340069（max_n）
%o a（1）=1；
%o n=2；
%o t=1；
%o，而n<=max_n
%o%搜索尚未在中使用的下一个数字
%o当~为空时（查找（a==t，1））
%o t=t+1；
%o端
%o比特1=长度（查找（比特（t，1:32）==1））；
%o比特2=长度（查找（比特（t*n，1:32）==1））；
%o if（位1==位2）
%o%我们找到了候选人
%o a（n）=t；
%o t=1；
%o n=n+1；
%o其他
%o%number t尚不适合
%o t=t+1；
%o端
%o端
%o结束
%o（PARI）列表a（nn）={my（va=向量（nn，k，-1））；对于（n=0，nn-1，my（k=0）；而（！（（hammingweight（k*n）==hammingweight（k）））&&！（#select（x->（x==k），va）），k++）；va[n+1]=k；）；va；}\\_Michel Marcus_，2020年12月30日
%o（Python）
%o定义binwt（n）：返回bin（n）。计数（'1'）
%o定义缺陷（n）：
%另外，设置=[]，设置（）
%对于范围（n+1）中的k，为o：
%o ak=0
%o当为True时：
%o当ak处于aset:ak+=1时
%o如果binwt（ak）==binwt
%o ak+=1
%o附加（ak）
%o基准添加（ak）
%o返回alst
%o打印（aupto（72））#_Michael S.Branicky_，2021年1月2日
%Y参见A000120、A077436、A340100、A263132（数字A000120（3）=A000120。
%K nonn，基础
%0、3
%A Thomas Scheuere，2020年12月28日
| 