(PARI)
r=平方(2);
z=100;
s(n)=如果(n<1,1,s(n-1)+1+楼层(n*r));
p(n)=n+1+总和(k=0,n,floor((n-k)/r));
u=v=矢量(z+1);
对于(n=1101,(v[n]=s(n-1)));
对于(n=1101,(u[n]=p(n-1)));
w(i,j)=u[i]+v[j]+(i-1)*(j-1)-1;
表(nn)={表示(n=1,nn,表示(k=1,n,打印1(w(k,n-k+1),“,”););打印();};
(Python)
从sympy导入sqrt
导入数学
定义s(n):如果n<1其他s(n-1)+1,则返回1+
int(math.floor(n*sqrt(2))
定义p(n):返回n+1+sum([int(math.floor((n-k)/sqrt(2))for k in range(0,n+1)])
v=[范围(0,101)中n的s(n)]
u=[p(n),对于范围(0,101)中的n]
定义w(i,j):返回u[i-1]+v[j-1]+(i-1)*(j-1)-1
对于范围(1,11)中的n:
….打印[w(k,n-k+1)for k in range(1,n+1)]#因德拉尼尔·戈什2017年3月21日
|