sage - 用矩阵组合集合：sagemath
问题描述
我希望每个人都做得很好。我尝试做一些组合，但我不知道如何做。如果有人能指出我的任何方向，那就太好了。这是问题所在：
M = FiniteEnumeratedSet({1, 2})
P1 = cartesian_product([M]*2)
P1 = P.list()
P2 = cartesian_product([M]*4)
P2 = [p for p in P2 if p[0] <= p[1]]
C = P1 + P2 #len(C) total columns
R = P1*3 #len(R) total rows
我想做的是构造一个矩阵，其中列数表示集合 C（16 列），行数表示 R（12 行）。每个入口的规则将是（一些功能entries）：
if (i,j) in R and [(i,j) in C:
return -1
elif if (i,j) in R and [(i,u,j,v) in C or (u,i,v,j) in C]:
return 1
else:
return 0
不幸的是，我不知道如何定义有效的功能entries。谢谢你。