%I#9 2017年7月29日19:25:23
%S 1,4,9,25,49,8,62115,77,35169,1621,10,12187,21209,27,91,65,
%电话:361,20289,32,55,18,14,33161,3913329119,95,85247,34125,22,45,
%U 44,69,26,24203,81217323259377155287,51115143,38145,36205,46,57,52
%N表示1<=k<=j的三角形T(j,k)按以下方式逐行填充:T(1,1)=1;T(j,k)是尚未使用的最小非素数,它与它的左上、左上、上和右上邻居互素,前提是存在这样的邻居(即属于三角形)。序列包含按行排列的三角形。
%C这会产生一个三角形,其中每个元素与其所有(最多八个)相邻元素互素。
%C非质数的重排。
%e三角形的前几行是
%e 1;
%e 4、9;
%e 25、49、8;
%e 6、121、15、77;
%e第35、169、16、221、10页;
%e第12、187、21、209、27、91页;
%e T(4,2)=121是与邻域25、49、8、6、15、35、169和16的互素。
%o(PARI){未使用的(m,T,j,k)=b=1;对于(jj=1,j-1,对于(kk=1,jj,if(T[jj,kk]==m,b=0));对于(kk=1,k-1,if
%o(PARI){下一个数(T,j,k)=T=1;m=4;而(T>0,如果(!i素数(m)&&未使用的(m,T,j,k)&&if(k>1,gcd(m,T[j-1,k-1])==1,1)&&如果(k<j,gcd j,k-1])==1,1),T=0,m++);m}
%o(PARI){n=50;T=矩阵(n,n
%Y参考A082201、A082202、A082203。
%K nonn,表
%O 1,2号机组
%A _Amarnath Murthy,2003年4月7日
%E由Klaus Brockhaus编辑、更正和扩展,2003年5月14日
|