%I#19 2023年12月29日10:58:54

%S 3,6,9,18,12,16,8,24,32,34,30,48,15,14,2,17,11,20,4,26,7,21,22,5,27，

%电话：10,13,25,40,19,28,33,37,23,42,38,44,35,31,41,29,46,49,39,50,36,43,47，

%U 45,51,54,55,53,52,56,57,62,61,60,64,68,67,58,63,70,69,71,65,77,66,72

%N a（1）=3。a（n）是序列中没有出现的最小正整数，因此（和{k=1到n}a（k））将乘积{j=1到n}a（j）除以。

%这是正整数的置换吗？根据汉斯·哈弗曼（Hans Havermann）的说法，该序列与序列A127562匹配，从每个序列的第104项开始。与序列A127562不同，此序列的前n项之和将前n项的乘积除以每个正整数n。

%p N:=1000:#得到a（1）到a（m-1），其中a（m）是第一项>N

%p a[1]：=3：

%p R：={$1..N}减去{3}：

%p p：=3：

%p S：=3：

%p成功：=真：

%从2开始p代表n，而成功和R代表{}

%p成功：=假；

%r do中r的p

%p如果类型为（（p*r）/（S+r），整数），则

%pa[n]：=r；

%p nmax:=n；

%p R：=R减去{R}；

%p成功：=真；

%p p：=p*r；

%p S：=S+r；

%p断点

%功率因数

%操作说明：

%日期：

%p序列（a[i]，i=1..nmax）；#_罗伯特·伊斯雷尔，2014年12月13日

%t f[l_List]：=块[{k=1，s=Plus@@l，p=Times@@l}，而[MemberQ[l，k]||Mod[k*p，k+s]>0，k++]；附加[l，k]]；巢穴[f，{3}，75]（*雷·钱德勒，2007年1月22日*）

%o（PARI）v=[3]；打印1（3，“，”）；n=1；而（n<100，p=prod（i=1，v，v[i]）；如果（p*n\（vecsum（v）+n）==p*n/（vecsum+n）&&！vecsearch（vecsort（v），n），v=concat（v，n）；打印1（n，“，”）；n=0）；n++）\\_Derek Orr_，2014年12月13日

%Y参见A127562、A127645和A127646。

%K nonn公司

%O 1,1号机组

%A _罗伊查询，2007年1月22日

%E由雷·钱德勒扩展，2007年1月22日