JZOJ 3823. 【NOIP2014模拟9.9】遇见

题目

Description
Zyh独自一人在街上漫步。Zyh相信不久后应该就可以和她一起漫步，可是去哪里寻找那个她呢？Zyh相信每个人都有一个爱情的号码牌，这个号码牌是一个n*n的矩阵。
每个人都要在矩阵中选择若干个元素，使得每行每列都有奇数个数被选中，且选中的数字的乘积是完全平方数。每当选出了这若干个元素，他/她就能找到那个她/他。
Zyh想知道对于一个号码牌有多少种选择的方法，使得zyh能够不再孤独。由于这个数字很大，只要输出对1,000,000,007取模后的余数即可。
Sample Input
输入1：
2
1 1
1 2
输入2：
2
620 620
620 620
Sample Output
输出1：
1
输出2：
2
Data Constraint
对于100%的数据 n<=30 Aij<=1000000000

题解

我们将这道题分为2部分。
第一部分：预处理。
首先，我们将每一个格子上的数分解质因数。然后分别求出第p个质因数分别各个格子里的指数。
设$a[p][i][j]$表示在第$(i,j)$格中第p个质数的指数(0/1)，因为完全平方数只需要满足各质因数指数为偶数。
设$x[i][j]$表示选或者不选第$(i,j)$格的数。
因为要使得每行每列都有奇数个数被选中，所以我们可以得到2n个方程：

$$x[i][1]异或x[i][2]异或...异或x[i][n]=1(1≤i≤n)$$
$$x[1][j]异或x[2][j]异或...异或x[n][j]=1(1≤j≤n)$$
还有P（P为质数的个数）个方程：
$$a[p][i][j]*x[i][j]异或....=0(1≤i≤n,1≤j≤n)$$
第二部分：高斯消元
设正在消第i个元，现在我们找到那些方程能够通过异或将第i个元消掉，那么就将两个方程的对应的系数一一异或。在找到第一个合法方程时，将它移到最上面，这样方便操作。
然后我们求自由元。（自由元：值不是固定的未知数）
产生自由元的情况有2种：
第一种，没有其他方程能够与之一起消掉。（显然这个元能够取任何值）
第二种，方程解完了，但是，剩下了一些元。那么这些元统统是自由元。（这些元的取值不影响方程）
答案为2的自由元个数次方。（自由元可选可不选，但是其他的必选或必不选）。