OFFSET
3,1
COMMENTS
a(8) = 70244258770074672.
PROG
(Python)
# Using graphillion
from graphillion import GraphSet
def make_CnXPk(n, k):
grids = []
for i in range(1, k + 1):
for j in range(1, n):
grids.append((i + (j - 1) * k, i + j * k))
grids.append((i + (n - 1) * k, i))
for i in range(1, k * n, k):
for j in range(1, k):
grids.append((i + j - 1, i + j))
return grids
def A(start, goal, n, k):
universe = make_CnXPk(n, k)
GraphSet.set_universe(universe)
paths = GraphSet.paths(start, goal)
return paths.len()
def A338963(n):
m = n * n
s = 0
for i in range(1, m):
for j in range(i + 1, m + 1):
s += A(i, j, n, n)
return s
print([A338963(n) for n in range(3, 7)])
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Seiichi Manyama, Dec 18 2020
STATUS
approved