https://codeforces.com/problemset/problem/1391/D
每个偶数长度的正方形都有奇数个1,那么存在长度为4的正方形肯定不行
那么把<=3的那一维状压,然后扫描另外一维进行状态转移就行了,根据状态算出这几个变成这个状态s的代价,枚举上一个位置的几个的状态t,如果满足2*2的正方形奇数个1,那么就可以从上一个位置的t状态转移到当前的s状态
#include<bits/stdc++.h>
#define pb push_back
using namespace std;
typedef long long ll;
const int maxl=1e6+10;
int n,m,cas,k,cnt,tot,ans;
vector<int> a[maxl];
int dp[maxl][9];
char ch[maxl];
bool in[maxl];
inline void prework()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
a[i].resize(m+2);
scanf("%s",ch+1);
for(int j=1;j<=m;j++)
a[i][j]=ch[j]-'0';
}
}
inline bool jug(int s,int t,int l)
{
int tmp1=0,tmp2=0;
if((s>>0)&1) tmp1++;
if((s>>1)&1) tmp1++;
if((t>>0)&1) tmp1++;
if((t>>1)&1) tmp1++;
if(tmp1%2==0) return false;
if(l>2)
{
if((s>>1)&1) tmp2++;
if((s>>2)&1) tmp2++;
if((t>>1)&1) tmp2++;
if((t>>2)&1) tmp2++;
if(tmp2%2==0) return false;
}
return true;
}
inline void mainwork()
{
int tmp=0;
if(n>=4 && m>=4)
{
ans=-1;
return;
}
if(n==1 || m==1)
{
ans=0;
return;
}
ans=n*m;
memset(dp,0x3f,sizeof(dp));
if(n<4)
{
for(int s=0;s<(1<<n);s++)
{
tmp=0;
for(int i=1;i<=n;i++)
if(a[i][1]!=((s>>(i-1))&1))
tmp++;
dp[1][s]=tmp;
}
for(int j=2;j<=m;j++)
for(int s=0;s<(1<<n);s++)
{
tmp=0;
for(int i=1;i<=n;i++)
if(a[i][j]!=((s>>(i-1))&1))
tmp++;
for(int t=0;t<(1<<n);t++)
if(jug(t,s,n))
dp[j][s]=min(dp[j-1][t]+tmp,dp[j][s]);
if(j==m)
ans=min(dp[j][s],ans);
}
}
else
{
for(int s=0;s<(1<<m);s++)
{
tmp=0;
for(int j=1;j<=m;j++)
if(a[1][j]!=((s>>(j-1))&1))
tmp++;
dp[1][s]=tmp;
}
for(int i=2;i<=n;i++)
for(int s=0;s<(1<<m);s++)
{
tmp=0;
for(int j=1;j<=m;j++)
if(a[i][j]!=((s>>(j-1))&1))
tmp++;
for(int t=0;t<(1<<m);t++)
if(jug(t,s,m))
dp[i][s]=min(dp[i-1][t]+tmp,dp[i][s]);
if(i==n)
ans=min(dp[i][s],ans);
}
}
}
inline void print()
{
printf("%d\n",ans);
}
int main()
{
int t=1;
//scanf("%d",&t);
for(cas=1;cas<=t;cas++)
{
prework();
mainwork();
print();
}
return 0;
}