%I#36 2015年1月15日10:18:49

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

%电话：55,26,33,28,45,32,51,38,17,19,30119,36,65,34,57,40,63,44,69,50,23，

%U 42,85,46,75,52,81,56,87,62,29,31,48

%如果N<=3，则N a（N）=N，否则，前面未出现的最小数与a（N-2）*a（N-3）至少有一个公因数，但与a（N-1）没有公因数。

%C从a（11）开始，序列与A098550不同。

%序列是自然数的排列。该证明与A098550的证明相似（有微小变化）_Vladimir Shevelev，2015年1月14日

%H Robert Israel，<a href=“/A247942/b247942.txt”>n，a（n）表，n=1.-10000</a>

%H David L.Applegate、Hans Havermann、Bob Selcoe、Vladimir Shevelev、N.J.A.Sloane和Reinhard Zumkeller，<A href=“http://arxiv.org/abs/1501.01669“>《黄石置换》，arXiv预印本arXiv:1501.016692015。

%p代表n从1到3 do

%p a[n]：=n：

%p b[n]：=1：

%日期：

%p表示n从4到1000 do

%pq:=a[n-2]*a[n-3]；

%从4do到k的p

%p如果未赋值（b[k]）且igcd（k，q）>1且igcd（k，a[n-1]）=1，则

%p a[n]：=k；

%p b[k]：=1；

%p断点

%功率因数

%日期：

%日期：

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

%t a[n/；n<=3]：=n；a[n_]：=a[n]=对于[aa=表[a[j]，{j，1，n-1}]；k=4，真，k++，如果[FreeQ[aa，k]&&！互质Q[k，a[n-2]*a[n-3]]&&互质Q[k，a[n-1]]，返回[k]]；表[a[n]，{n，1，60}]（*Jean-François Alcover_，2015年1月12日*）

%Y参见A098550、A249167、A251604、A247225。

%K nonn公司

%O 1,2号机组

%2015年1月11日，A_Vladimir Shevelev

%E更多来自Jean-François Alcover的条款，2015年1月12日