不妨规定先翻转某些行,再翻转某些列。
设(F[x]=min(mathbb{pop}(x),mathbb{pop}((2^n-1)veebar x))),(S_i)为第(i)列的表格状态。
枚举每一行的翻转情况(p),此时的答案为(sum_i F[S_iveebar p])。
整体答案
[ans
=min_{p=0}^{2^n-1} sum_{i=0}^{m-1} F[S_iveebar p]
=min_{p=0}^{2^n-1} sum_{q=0}^{2^n-1}F[q]sum_{i=0}^{m-1} [S_iveebar p=q]
]
记(T[x]=sum_{i=0}^{m-1}[S_i=x]),则
[ans
=min_{p=0}^{2^n-1} sum_{q=0}^{2^n-1}F[q]T[pveebar q]
=min_{p=0}^{2^n-1} sum_{x=0}^{2^n-1}sum_{y=0}^{2^n-1}[xveebar y=p]F[x]T[y]
]
因此直接对(F,T)做XOR卷积,拿出结果的最小值即可。
我受够了你谷公式渲染, 求您们修修MathJax吧
#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int L=1<<20;
const int N=1e5;
int n,m,s[N];
ll F[L],T[L];
void FWT_XOR(ll a[],int len) {
for(int m=1; m<len; m<<=1)
for(int i=0,s=m<<1; i<len; i+=s)
for(int j=0; j<m; ++j) {
ll x=a[i+j], y=a[i+m+j];
a[i+j]=x+y, a[i+m+j]=x-y;
}
}
void IFWT_XOR(ll a[],int len) {
for(int m=1; m<len; m<<=1)
for(int i=0,s=m<<1; i<len; i+=s)
for(int j=0; j<m; ++j) {
ll x=a[i+j], y=a[i+m+j];
a[i+j]=(x+y)/2;
a[i+m+j]=(x-y)/2;
}
}
int main() {
scanf("%d%d",&n,&m);
int len=1<<n,col;
for(int i=0; i<n; ++i)
for(int j=0; j<m; ++j)
scanf("%1d",&col),s[j]|=col<<i;
for(int i=0; i<len; ++i)
F[i]=min(__builtin_popcount(i),__builtin_popcount((len-1)^i));
for(int i=0; i<m; ++i) T[s[i]]++;
FWT_XOR(F,len);
FWT_XOR(T,len);
for(int i=0; i<len; ++i) F[i]*=T[i];
IFWT_XOR(F,len);
printf("%lld
",*min_element(F,F+len));
return 0;
}