%I#41 2023年12月30日23:07:49

%S 1,2,3,4,6,13,23,7,43,73,9,5,12,8,10,14,16,15,18,21,20,59,22101,24，

%电话：27,19,25,71,30,26127,33,28,32,31,35,34,36,39109,38,40,11,89,42,44，

%U 45131,46,37,48262347,51,50,49,52151,54257,55,58,65

%对于N＜＝3，N a（N）＝N；当n>=4时，a（n）是该序列中以前没有出现过的最小正整数，与a（n-1）+a（n-2）+a（n-3）共用一个因子。

%C猜想：这个序列是自然数的排列。

%谢一凡，n的表，n的a（n）=1..2000</a>

%ta[1]=1；a[2]=2；a[3]=3；a[n]：=a[n]=模块[{s，i=1}，s=a[n-1]+a[n-2]+a[n-3]；而[MemberQ[a/@Range[1，n-1]，i]||GCD[s，i]==1，i++]；i] ；

%t表[a[n]，{n，1，65}]（*_Robert P.P.McGone_，2023年12月30日*）

%o（PARI）列表a（nn）={my（v=[1,2,3]）；对于（n=4，nn，my（t=1）；而（prod（X=1，n-1，v[X]-t）==0||gcd（v[n-3]+v[n-2]+v[n-1]，t）==1，t++）；v=concat（v，t））；v；}

%o（Python）

%o从数学导入gcd

%o a=[1,2,3]

%o t=设置（a）

%o定义next_element（）：

%o s=a[-1]+a[-2]+a[-3]

%o n=1

%o当n in t或gcd（s，n）==1时：

%o n+=1

%o返回n

%o定义a_seq（ul）：

%o对于范围（4，ul+1）内的_：

%o nn=下一个元素（）

%o a.附录（nn）

%o t.添加（nn）

%o返回a

%o打印（a_seq（65））#_Robert P.P.McKone_，2023年12月30日

%Y参考A064413，A337136。

%K nonn，简单

%O 1,2号机组

%阿一凡谢，2023年5月21日