%I#14 2019年12月7日12:33:53
%S 1,2,3,4,5,6,7,8,9,11,10,12,13,15,17,14,16,18,20,22,25,19,21,23,26,28,
%电话:31,34,24,27,29,32,35,38,41,44,30,33,36,39,42,46,49,52,56,37,40,43,47,
%U 50,54,58,61,65,69,45,48,51,55,59,63,67,71,75,79,84,53
%N sqrt签名序列的间隔(1/2)。
%C每一行分散所有其他行,每一列分散所有其他列。数组是补码的离散值(第1列=A022776)。
%C R(n,m)=当所有数字k*R+h联合排序时,n*R+m的位置,其中R=sqrt(2),k>=1,h>=0_克拉克·金伯利(Clark Kimberling),2017年10月6日
%H Clark Kimberling,反对角线n=1..60,扁平</a>
%H Clark Kimberling和John E.Brown,<a href=“http://www.cs.uwaterloo.ca/journals/JIS/VOL7/Kimberling/kimber67.html“>部分互补和转座分散</a>,《整数序列杂志》,2004年第7卷。
%F R(i,j)=R(i、0)+R(0、j)+i*j-1,对于i>=1,j>=1。
%e R西北角:
%电子1 2 4 7 10 14 19 24 30
%电子邮箱3 5 8 12 16 21 27 33 40
%电子6 9 13 18 23 29 36 43 51
%电子11 15 20 26 32 39 47 44 64
%电子17 22 28 35 42 50 59 68 78
%电子邮箱:25 31 38 46 54 63 73 83 94
%t r=平方英尺[1/2];z=100;
%t s[0]=1;s[n]:=s[n]=s[n-1]+1+楼层[n*r];
%t u=表[n+1+总和[下限[(n-k)/r],{k,0,n}],{n,0,z}](*A022775,A283962的第1列*)
%t v=表[s[n],{n,0,z}](*A022776,A283962的第1行*)
%tw[i,j]:=u[i]]+v[[j]]+(i-1)*(j-1)-1;
%t网格[表[w[i,j],{i,1,10},{j,1,10}]](*A283962,数组*)
%t压扁[表[w[k,n-k+1],{n,1,20},{k,1,n}]](*A283962,序列*)
%o(PARI)
%o r=平方(1/2);
%o z=100;
%o s(n)=如果(n<1,1,s(n-1)+1+楼层(n*r));
%o p(n)=n+1+总和(k=0,n,楼层((n-k)/r));
%o u=v=矢量(z+1);
%o表示(n=1101,(v[n]=s(n-1));
%o表示(n=1101,(u[n]=p(n-1));
%o w(i,j)=u[i]+v[j]+(i-1)*(j-1)-1;
%o表(nn)={对于(n=1,nn,对于(k=1,n,打印1(w(k,n-k+1),“,”););打印();};
%o表(10),2017年3月21日
%o(Python)
%o来自sympy导入sqrt
%o导入数学
%o定义s(n):如果n<1其他s(n-1)+1,则返回1+
%o int(数学楼层(n*sqrt(1/2))
%o定义p(n):对于k in,返回n+1+总和([int(math.floor((n-k)/sqrt(1/2)))
%o范围(0,n+1)])
%o v=[范围(0,101)内n的s(n)]
%o u=[p(n)表示范围(0,101)内的n]
%o定义w(i,j):返回u[i-1]+v[j-1]+(i-1)*(j-1)-1
%o表示范围(1,11)中的n:
%o。。。。print[w(k,n-k+1)for k in range(1,n+1)]#_Indranil Ghosh,2017年3月21日
%Y参见A010503、A022775、A022766、A283939。
%K non,tabl,简单
%O 1,2号机组
%A _百灵鸟金伯利,2017年3月19日
|