2020年3月23日月曜日

200323

Python


A003763(1)

Graphillion を使って計算してみた。

from graphillion import GraphSet
import graphillion.tutorial as tl

def A(n):
    universe = tl.grid(n - 1, n - 1)
    GraphSet.set_universe(universe)
    start, goal = 1, 2
    paths = GraphSet.paths(start, goal)
    for i in range(3, n * n + 1):
        paths = paths.including(i)
    return paths.len()

def A003763(n):
    return A(2 * n)

print([A003763(n) for n in range(1, 5)])

出力結果
[1, 6, 1072, 4638576]

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。