%I#23 2020年3月26日07:11:17
%S 1,1,8,1,16,29,1,32,95,80,1,64313426195,1128103323201745444,1,
%电话:25634112706163476838969,1512112656966215425911257225897,
%电话:2056,110243720538196414544951859660752245952924279
%N方阵T(N,k),N>=1,k>=2,由反对角线读取,其中T(N,k)是N×k网格图中自回避行走的次数,自回避行走从图左侧的N个顶点中的任何一个开始,到右侧的N个顶点中的任何一个结束。
%e方阵T(n,k)开始:
%e 1,1,1,1。。。
%e 8、16、32、64、128。。。
%e 29、95、313、1033、3411。。。
%e 80、426、2320、12706、69662。。。
%e 195、1745、16347、154259、1454495。。。
%电子邮箱4446838112572185966030549774。。。
%o(Python)
%o#使用石墨
%o从graphillion导入GraphSet
%o导入graphillion.tutorial作为tl
%o定义A(开始,目标,n,k):
%o宇宙=tl.grid(n-1,k-1)
%o GraphSet.set_universe(universe)
%o paths=GraphSet.paths(开始,目标)
%o返回路径.len()
%o定义A333509(n,k):
%o如果n=1:返回1
%o s=0
%o对于范围(1,n+1)中的i:
%对于范围内的j(k*n-n+1,k*n+1):
%o s+=A(i,j,k,n)
%o返回s
%o打印([A333509(j+1,i-j+2)表示范围(9)中的i,表示范围(i+1)中的j])
%Y列k=2-3给出:A333510、A333511。
%Y行n=1-3给出:A000012、A000079(n+1)、2*A082574(n+1”)+1。
%Y T(n,n)表示A121785(n-1)。
%K nonn,表格
%氧1,3
%A _Seiichi Manyama,2020年3月25日
|