分割表を代数統計する

• 2x2表をSingular でやってみよう
• x+y=p1,z+w=p2,x+z=q1,y+w=q2とする(p1,p2,q1,q2が周辺度数)

ring p = 0,(p1,p2,q1,q2),lp;
ring x = 0,(x,y,z,w),lp;
setring p;
map f = x,p1*q1,p1*q2,p2*q1,p2*q2;
ideal i0 = 0;
setring x;
preimage(p,f,i0);

• これは前記事の連鎖平衡とまったく同じ
• 2x3表をやってみよう

ring p =0,(p1,p2,q1,q2,q3),lp;
ring x =0,(x1,x2,x3,y1,y2,y3),lp;
setring p;
map f = x,p1*q1,p1*q2,p1*q3,p2*q1,p2*q2,p2*q3;
ideal i0=0;
setring x;
preimage(p,f,i0);

• が返る