|
（PARI）我的（N=66，x='x+O（'x^N））；向量（（x*（1+x））^2/（1-2*x-4*x^2-2*x^3））
（Python）
#使用石墨
从graphillion导入GraphSet
定义make_T_nk（n，k）：
网格=[]
对于范围（1，k+1）中的i：
对于范围（1，n）中的j：
网格.附加（（i+（j-1）*k，i+j*k））
如果i<k：
网格.附加（（i+（j-1）*k，i+j*k+1））
对于范围（1，k*n，k）中的i：
对于范围（1，k）中的j：
网格.附加（（i+j-1，i+j））
回流格栅
universe=make_T_nk（n，k）
GraphSet.set_universe（宇宙）
cycles=GraphSet.cycles（is_hamilton=True）
返回周期.len（）
| 