%I#6 2012年3月30日18:57:32

%S 1,3,2,6,4,5,10,7,9,8,16,12,15,13,11,25,19,24,21,18,14,39,30,37,33,28，

%电话：22,17,60,46,57,51,43,34,27,20,91,70,87,78,66,52,32,3138106132，

%电话：118100,79,64,48,36,26208160199178151120,97,73,55

%N（2,5,8,11,14,17，…）的反对偶扩散。

%C第1行：A152009。

%C有关分散体的背景讨论，请参见A191426。

%C。。。

%C序列（3n，n>0）、（3n+1，n>O）、（3d+2，n>=0）中的每一个序列都会产生色散。每个补码（以第一项>1开始）也会产生一个离散。这里列出了六个序列和分散体：

%C。。。

%C A191449=A008583的分散度（0 mod 3）

%C A191451=A016777的分散度（1 mod 3）

%C A191450=A016789的分散度（2 mod 3）

%C A191656=A001651的离散度（1或2 mod 3）

%C A083044=A007494的离散度（0或2 mod 3）

%C A191655=A032766的离散度（0或1 mod 3）

%C。。。

%C除了最多2个初始项（因此第1列总是以1开头）：

%C A191449具有第1列A001651，其他所有列A008583

%C A191451具有第一列A007494，其他所有列为A016777

%C A191450具有第一列A032766，所有其他列为A016789

%C A191656第1列为A008583，其他所有为A001651

%C A083044具有第一列A016777，所有其他列为A083044

%C A191655具有第一列A016789，所有其他列为A032766

%C。。。

%C“（a或b mod m）”类型的序列有一个公式（如下面的Mathematica程序所示）：

%C如果f（n）=（n mod 2），则（a、b、a、b，a、b…）由下式给出

%C a*f（n+1）+b*f（n），因此“（a或b mod m）”由

%C a*f（n+1）+b*f（n）+m*楼面（（n-1）/2）），对于n>=1。

%e西北角：

%e 1…3…6…10…16

%e 2…4…7…12…19

%e 5…9…15…24…37

%e 8…13…21…33…51

%e 11..18..28…43…66

%t（*程序生成递增序列f[n]*的色散阵列t）

%t r=40；r1=12；c=40；c1=12；

%t a=3；b=4；m[n_]：=如果[Mod[n，2]==0，1，0]；

%tf[n_]：=a*m[n+1]+b*m[n]+3*层[（n-1）/2]

%t表[f[n]，{n，1，30}]（*A032766:（3+5k，4+5k，k>=0）*）

%t mex[list_]：=NestWhile[#1+1&，1，并集[list][[#1]]<=#1&，1、长度[Union[list]]]

%t行={NestList[f，1，c]}；

%t Do[rows=Append[rows，NestList[f，mex[Flatten[rows]]，r]]，{r}]；

%t[i_，j_]：=行[[i，j]]；

%t表格形式[表格[t[i，j]，{i，1,10}，{j，1,10{]]（*A191655数组*）

%t压扁[表[t[k，n-k+1]，{n，1，c1}，{k，1，n}]]（*A191655序列*）

%Y参考A016789、A032766、A191426。

%K nonn，表格

%O 1,2号机组

%A_Clark Kimberling_，2011年6月10日