【CF662C】Binary Table（FWT）

题面

洛谷
CF
翻译：
有一个$n*m$的表格($n<=20,m<=10^5$)，
每个表格里面有一个$0/1$，
每次可以将一行或者一列的$01$全部翻转
回答表格中最少有多少个$1$

题解

发现$n$很小，$m$很大
状压是跑不掉了
如果我们确定翻转哪些行，那么答案唯一确定（贪心的选每一列中$0/1$的较小值）
相同的列显然可以合并，
把每一列按照$01$状压，记$a[i]$为状态为$i$的列的个数
记$f[i]$表示翻转状态为$i$的那些行的结果
假设翻转的行为$S$，这一列的状态为$i$，显然最终就变成了$i\oplus S$
而对于每一列的任意一种状态$i$
答案显然是$min(Numberof(0),Numberof(1))$
预处理出来
这样子的话，我们就知道了$f[i]=\sum_{j\oplus k=i}a[j]*b[k]$
显然这是一个$xor$卷积，用$FWT$优化即可。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
inline int read()
{
    RG int x=0,t=1;RG char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}
int N,n,m;
ll a[1<<20],b[1<<20];
int g[21][100001];
void FWT(ll *P,int opt)
{
    for(int i=2;i<=N;i<<=1)
        for(int p=i>>1,j=0;j<N;j+=i)
            for(int k=j;k<j+p;++k)
            {
                ll x=P[k],y=P[k+p];
                P[k]=x+y;P[k+p]=x-y;
                if(opt==-1)P[k]/=2,P[k+p]/=2;
            }
}
int main()
{
    n=read();m=read();N=1<<n;
    char ch[100001];
    for(int i=1;i<=n;++i)
    {
        scanf("%s",ch+1);
        for(int j=1;j<=m;++j)g[i][j]=ch[j]-48;
    }
    for(int i=1;i<=m;++i)
    {
        int x=0;
        for(int j=1;j<=n;++j)x=(x<<1)|g[j][i];
        a[x]++;
    }
    for(int i=0;i<N;++i)b[i]=b[i>>1]+(i&1);
    for(int i=0;i<N;++i)b[i]=min(b[i],n-b[i]);
    FWT(a,1);FWT(b,1);
    for(int i=0;i<N;++i)a[i]*=b[i];
    FWT(a,-1);
    ll ans=n*m;
    for(int i=0;i<N;++i)ans=min(ans,a[i]);
    printf("%I64d\n",ans);
    return 0;
}