D. 505(状压dp)

https://codeforces.com/problemset/problem/1391/D

---

思路：我记得这个题我似乎去年碰到了。当时看了一眼就扔了。今年给补上了。

这个题有个结论。就是n>=4&&m>=4的时候是-1的。所以不用存图。

能找到的最小的子矩阵应该是2 * 2 ，如果一个2 * 2矩阵内含有奇数个1，那么四个2 * 2矩阵拼接形成一个4 * 4矩阵中肯定就含有偶数个1，所以如果n和m同时大于等于4，那么一定无法构造出符合题意的矩阵.

那么手动模拟一下可以发现，固定第一列的形式

1的个数:                    0个 1个          2个

第二列必须                1个  0个/2个  1个

所以可以暴力枚举改成第一列的情况来算最终答案。

我是把n=2和n=3一起做了。当成状压去做感觉方便一点。

我们设 dpi,j​ 为第 i行在 j 状态下时需要修改的最小次数。

统计修改次数只要算当前的数列和目标序列第一行+第二行+（第三行）不同的数。

判断前一列和当前列可否满足可以麻烦点讨论。但是位运算^更加方便。只要正方形区间^答案为1说明就是奇数个，不然就不行。一下子大大简化了代码量。

#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<cstdio>
#include<algorithm>
#define debug(a) cout<<#a<<"="<<a<<endl;
#define getlow(x,k)  ( (x>>k)&1 )
using namespace std;
const int maxn=1e6+100;
typedef int LL;
inline LL read(){LL x=0,f=1;char ch=getchar();	while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;}
LL a[10][maxn];
LL dp[10][maxn];
bool ok2(LL x,LL y){
     if(getlow(x,0)^getlow(y,0)^getlow(x,1)^getlow(y,1)==0) return 0;
     else return 1;
}
bool ok3(LL x,LL y){
     if(getlow(x,0)^getlow(y,0)^getlow(x,1)^getlow(y,1)==0) return 0;
     if(getlow(x,1)^getlow(y,1)^getlow(x,2)^getlow(y,2)==0) return 0;
     else return 1;
}
LL cal2(LL id,LL y){
     LL sum=0;
     if(a[1][id]!=getlow(y,0)) sum++;
     if(a[2][id]!=getlow(y,1)) sum++;
     return sum;
}
LL cal3(LL id,LL y){
     LL sum=0;
     if(a[1][id]!=getlow(y,0)) sum++;
     if(a[2][id]!=getlow(y,1)) sum++;
     if(a[3][id]!=getlow(y,2)) sum++;
     return sum;
}
int main(void){
   cin.tie(0);std::ios::sync_with_stdio(false);
   LL n,m;cin>>n>>m;
   if(n>=4&&m>=4){
      cout<<"-1"<<"\n";return 0;
   }
   for(LL i=1;i<=n;i++){
       for(LL j=1;j<=m;j++) {
            char c;cin>>c;a[i][j]=c-'0';
       }
   }
   if(n==1){
      cout<<"0"<<"\n";
   }
   else if(n==2){
      memset(dp,0x3f,sizeof(dp));
      for(LL j=0;j<4;j++)  dp[j][1]=cal2(1,j);
      for(LL i=2;i<=m;i++){
          for(LL j=0;j<4;j++){
              for(LL k=0;k<4;k++){
                  if(!ok2(j,k)) continue;
                  dp[k][i]=min(dp[k][i],dp[j][i-1]+cal2(i,k));
              }
          }
      }
      LL ans=0x3f3f3f3f;
      for(LL j=0;j<4;j++){
          ans=min(ans,dp[j][m]);
      }cout<<ans<<"\n";
   }
   else if(n==3){
      memset(dp,0x3f,sizeof(dp));
      for(LL j=0;j<8;j++) dp[j][1]=cal3(1,j);
      for(LL i=2;i<=m;i++){
          for(LL j=0;j<8;j++){
              for(LL k=0;k<8;k++){
                  if(!ok3(j,k)) continue;
                  dp[k][i]=min(dp[k][i],dp[j][i-1]+cal3(i,k));
              }
          }
      }
      LL ans=0x3f3f3f3f;
      for(LL j=0;j<8;j++){
          ans=min(ans,dp[j][m]);
      }cout<<ans<<"\n";

   }
   return 0;
}