3076 翻翻乐

Task
翻转01矩阵，使对任意一个点，周围点的异或和为0.每个点对应一个代价，有的点不能翻转。求最小代价。
数据范围：矩阵的大小为n*n,n<=16,如果不可能完成输出-1.

Solution
终态枚举
对于翻转问题，同一个格子翻转2次会恢复原状。翻转格子集合是相同时，次序是不重要的。
不妨先假定第一行的01状态，那么第二行的状态是是由上一行唯一决定的，因此整个图的状态和代价都可以计算。
不要忘记输出-1。（10分）

const int M=18;
int n,ans=-1;
int A[M][M],val[M][M],a[M][M];
bool cnt[M][M];
int rx[]={0,0,-1},ry[]={1,-1,0};
inline void input(){
    int i,j,k;
    rd(n);
    rep(i,1,n)
        rep(j,1,n)rd(A[i][j]);
    rep(i,1,n)
        rep(j,1,n)rd(val[i][j]);
}
inline bool add(int x,int y){return x>=1&&x<=n&&y>=1&&y<=n&&a[x][y]==1;}
inline int solve(){
    int i,j,k,sum=0;
    memset(cnt,0,sizeof(cnt));
    rep(i,1,n){
        if(A[1][i]!=a[1][i])sum+=val[1][i];
        rep(j,0,2){
            cnt[1][i]^=add(1+rx[j],i+ry[j]);
        }
    }
    rep(i,2,n){
        rep(j,1,n){
            a[i][j]=cnt[i-1][j];
            if(A[i][j]!=a[i][j]){
                if(val[i][j]==-1)return 1e9;
                sum+=val[i][j];
            }
        }
        rep(j,1,n)
            rep(k,0,2)
                cnt[i][j]^=add(i+rx[k],j+ry[k]);
    }
    return sum;
}
inline void dfs(int id,int c){
    if(val[1][id]==-1&&A[1][id]!=c)return;//不能改变 
    if(id==n+1){
        MIN(ans,solve());
        return;
    }
    a[1][id]=c;
    dfs(id+1,0);
    dfs(id+1,1);
}
int main(){
    input();
    dfs(0,0);
    if(ans==1e9)ans=-1;
    printf("%d\n",ans);
    return 0;
}