[CF662C]Binary Table

luogu

题意

你有一个\(n*m\)的\(01\)矩阵。你可以把任意一行或者一列的\(01\)取反。求矩阵中最少的\(1\)的数量。
\(n\le20,m\le10^5\)

sol

很自然地有一个\(O(2^nm)\)的暴力：枚举横行的取反情况，然后纵列就取\(01\)数量较少的一者。
我们记状态\(x\)在原矩阵中的出现次数为\(a[x]\)，状态中\(01\)较少一者的数量为\(b[x]\)。
会发现当最终的取反状态为\(i\)时，会有\(ans_i=\sum_{j\otimes k=i}a_j*b_k\)。
就是一个异或卷积。FWT即可。

code

#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;
}