%I#21 2020年12月22日03:47:12
%S 1,2,4,8,16,6,9,12,18,5,10,7,24,15,11,14,20,28,13,25,17,19,21,32,22,
%电话:26,40,23,27,30,35,29,42,31,36,34,33,39,38,44,46,48,54,52,50,45,56,49,
%U 51,57,55,58,63,60,62,66,37,69,72,75,64,78,81,84,90,87,80,68,76,88,96,92104,65100,70,74,85108单位
%N从第一个数字开始,每次读取一个数字d;用序列中尚未出现的数据的最小倍数扩展S。不读取零。
%这是大于0的整数的置换(因为素数是1的倍数)。
%H Carole Dubois,<a href=“/A39853/b339853.txt”>n,a(n)表,n=1..5000</a>
%H<a href=“/index/Per#IntegerPermutation”>自然数排列序列的索引项</a>
%e第一个d是1;由于1已经在序列中,我们用2对其进行了扩展;
%e下一个d现在是2;由于序列中已经有2,我们将其扩展为4(序列中没有2的最小倍数);
%e下一个d是4;由于序列中已经有4,我们将其扩展为8(最小倍数);
%e下一个d是8;由于序列中已经有8,我们将其扩展为16(最小倍数);
%e下一个d是1;我们用3扩展了序列,因为3是序列中尚未出现的1的最小倍数;
%e下一个d是6;由于6还没有出现,我们用6扩展了序列;
%e下一个d是3;我们用9扩展了序列,因为9是3的最小倍数;
%e下一个d是6;我们用12扩展了序列,因为12是6的最小倍数;等。
%由于10的零不会被读取,我们将用7的最小倍数来扩展这个序列,这个倍数是14。
%o(Python)
%o定义缺陷(n):
%o alst,astr,used,字符串=[1],“1”,{1},0
%o对于范围(1,n)中的k:
%o而astr[strind]==“0”:strind+=1
%o ak=数字=int(astr[字符串])
%o使用ak时:ak+=数字
%附加附加(ak);astr+=str(ak);used.add(ak);字符串+=1
%o return alst#对a(n)使用alst[n-1]
%o打印(aupto(71))#_Michael S.Branicky_,2020年12月19日
%Y参考A316749。
%K基,nonn
%O 1,2号机组
%A _Eric Angelini和Carole Dubois,2020年12月19日
|