codeforses 662C. Binary Table

题意：

给出一个01矩阵，可以将若干列或行取反，是最后1的个数尽量小。
n<=20,m<=100000

题解：

假如取反那些行确定的话，列的操作是唯一的。
所以列的状态确定了，他对答案的贡献也确定了。
考虑暴力，显然就是枚举选那些行。
用fwt优化即可。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
LL bin[22];
void fwt(LL *f,LL n,LL op)
{
    for(LL i=1;i<n;i<<=1)
        for(LL j=0;j<n;j+=i<<1)
            for(LL k=0;k<i;k++)
            {
                LL x=f[j+k],y=f[i+j+k];
                f[j+k]=x+y;f[i+j+k]=x-y;
                if(op==-1) f[j+k]>>=1,f[i+j+k]>>=1;
            }
}
LL n,m,a[200010];
LL A[1100000],B[1100000];
char s[200010];
LL calc(LL x)
{
    LL k[2];k[0]=k[1]=0;
    for(LL i=0;i<n;i++) k[x&bin[i]?1:0]++;
    return k[0]<k[1]?k[0]:k[1];
}
int main()
{
    bin[0]=1;for(LL i=1;i<=20;i++) bin[i]=bin[i-1]<<1;
    scanf("%lld %lld",&n,&m);
    for(LL i=1;i<=n;i++)
    {
        scanf("%s",s+1);
        for(LL j=1;j<=m;j++) a[j]+=s[j]=='0'?0:bin[i-1];
    }
    for(LL i=0;i<bin[n];i++) A[i]=calc(i);
    for(LL i=1;i<=m;i++) B[a[i]]++;
    fwt(A,bin[n],1);fwt(B,bin[n],1);
    for(LL i=0;i<bin[n];i++) A[i]*=B[i];
    fwt(A,bin[n],-1);
    LL ans=1LL<<60;
    for(LL i=0;i<bin[n];i++) ans=min(ans,A[i]);
    printf("%lld\n",ans);
}