Codeforces 2011-2012 Waterloo Local Contest, 24 September, 2011 E.Harmonious Matrices

链接：http://codeforces.com/gym/100008

题意：给定矩阵的长宽n,m。要求在每个格子填0/1并且每个格子自己+周围4个格子为偶数，要求1尽可能多。

分析：很明显我们只需要知道第一行的填写情况就能将整个矩阵填出来，但是暴力2^m填写是不可能的，我们设这m个位置为x1,x2...xm。xi=0/1。那么我们可以得到第n+1行每一个位置的情况（一个关于x的集合）且n+1行应该每一个位置的值都是0。这样我们就能得到m个异或方程，然后异或高斯消元就能求解啦，因为要求1尽可能的多，所以我们在自由元的位置全部放1。

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<vector>
#include<string>
#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=50;
const int mod=100000000;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000007;
const int INF=1000000010;
const ll MAX=1000000000000;
const double pi=acos(-1.0);
typedef double db;
typedef unsigned long long ull;
int a[N][N];
ll f[N][N];
void gauss(int n,int m) {
    int i=1,j,k,h,w=0;//i-1:矩阵的秩
    for (j=1;j<m;j++,w=0) {
        for (k=i;k<=n;k++)
        if (a[k][j]) w=k;
        if (w) {
            for (k=j;k<=m;k++) swap(a[i][k],a[w][k]);
            for (k=i+1;k<=n;k++)
            if (a[k][j]) {
                for (h=j;h<=m;h++) a[k][h]^=a[i][h];
            }
            i++;
        }
        if (i>n) break ;
    }
    for (j=1;j<m;j++) f[1][j]=1;
    for (j=i-1;j;j--) {
        for (k=1;k<m;k++)
        if (a[j][k]) break ;
        f[1][k]=f[j][m];
        for (h=k+1;h<m;h++)
        if (f[1][h]&&a[j][h]) f[1][k]^=1;
    }
}
int main()
{
    int i,j,n,m,t;
    scanf("%d", &t);
    while (t--) {
        scanf("%d%d", &n, &m);
        memset(a,0,sizeof(a));
        memset(f,0,sizeof(f));
        for (i=1;i<=m;i++) f[1][i]=1ll<<i;
        for (i=2;i<=n+1;i++)
            for (j=1;j<=m;j++) f[i][j]=f[i-1][j-1]^f[i-1][j]^f[i-1][j+1]^f[i-2][j];
        for (i=1;i<=m;i++)
            for (j=1;j<=m;j++)
            if (f[n+1][i]&(1ll<<j)) a[i][j]=1;
        gauss(m,m+1);
        for (i=2;i<=n;i++)
            for (j=1;j<=m;j++) f[i][j]=f[i-1][j-1]^f[i-1][j]^f[i-1][j+1]^f[i-2][j];
        for (i=1;i<=n;i++) {
            for (j=1;j<m;j++) printf("%d ", f[i][j]);
            printf("%d\n", f[i][m]);
        }
    }
    return 0;
}