标签：void space div || ble pre 取反 line href
你有一个\(n*m\)的\(01\)矩阵。你可以把任意一行或者一列的\(01\)取反。求矩阵中最少的\(1\)的数量。
\(n\le20,m\le10^5\)
很自然地有一个\(O(2^nm)\)的暴力：枚举横行的取反情况，然后纵列就取\(01\)数量较少的一者。
我们记状态\(x\)在原矩阵中的出现次数为\(a[x]\)，状态中\(01\)较少一者的数量为\(b[x]\)。
会发现当最终的取反状态为\(i\)时，会有\(ans_i=\sum_{j\otimes k=i}a_j*b_k\)。
就是一个异或卷积。FWT即可。
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int gi(){
int x=0,w=1;char ch=getchar();
while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if (ch=='-') w=0,ch=getchar();
while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return w?x:-x;
}
#define ll long long
const int N = 1<<20;
const int M = 1e5+5;
int n,m,len;ll a[N],b[N];char s[20][M];
void fwt(ll *P,int opt){
for (int i=1;i<len;i<<=1)
for (int p=i<<1,j=0;j<len;j+=p)
for (int k=0;k<i;++k){
ll x=P[j+k],y=P[j+k+i];
P[j+k]=(x+y)/opt;P[j+k+i]=(x-y)/opt;
}
}
int main(){
n=gi();m=gi();len=1<<n;
for (int i=0;i<n;++i) scanf("%s",s[i]+1);
for (int i=1;i<=m;++i){
int x=0;
for (int j=0;j<n;++j) x=(x<<1)+s[j][i]-'0';
++a[x];
}
for (int i=1;i<len;++i) b[i]=b[i>>1]+(i&1);
for (int i=1;i<len;++i) b[i]=min(b[i],n-b[i]);
fwt(a,1);fwt(b,1);
for (int i=0;i<len;++i) a[i]*=b[i];
fwt(a,2);
ll ans=a[0];
for (int i=1;i<len;++i) ans=min(ans,a[i]);
printf("%I64d\n",ans);return 0;
}标签：void space div || ble pre 取反 line href
原文地址：https://www.cnblogs.com/zhoushuyu/p/9068280.html
