CF662C Binary Table

一、题目

点此看题

二、解法

观察到$n$很小，我们肯定要从这里入手，设$i$列的状态压缩为$s_i$，我们先枚举行的翻转选择$x$，那么列翻转后的状态是$y_i=s_i\oplus x$，设$f_i$为二进制$i$最少有多少个$1$（可翻转），那么答案就是：
$$\sum _{i=1}^m{f}_{s_i\oplus x}$$我们枚举这个$s_i$和$y_i$，那么$x$就不需要了，设$q_i$为有多少列的初始状态是$i$：
$$\sum _{s=0}^{2^n-1}\sum _{y=0}^{2^n-1}[s\oplus y=x]f_y\times q_s$$改一种写法，就会清晰很多：
$$ans(x)=\sum _{s\oplus y=x}{f}_y\times q_s$$直接搞$xor$卷积就可以了，时间复杂度$O(n{2}^n)$

#include <cstdio>
#include <iostream>
using namespace std;
const int M = 110005;
#define int long long
int read()
{
    int num=0,flag=1;
    char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,lim,a[25][M],s[M],f[M*10],g[M*10];
void fwt(int *a,int n,int op)
{
    for(int i=1; i<n; i<<=1)
        for(int p=i<<1,j=0; j<n; j+=p)
            for(int k=0; k<i; k++)
            {
                int x=a[j+k],y=a[i+j+k];
                a[j+k]=x+y;
                a[i+j+k]=x-y;
                if(op==-1)
                {
                    a[j+k]=a[j+k]/2;
                    a[i+j+k]=a[i+j+k]/2;
                }
            }
}
signed main()
{
    n=read();
    m=read();
    lim=1<<n;
    for(int i=1; i<=n; i++)
        for(int j=1; j<=m; j++)
            scanf("%1lld",&a[i][j]);
    for(int i=1; i<=m; i++)
        for(int j=1; j<=n; j++)
            s[i]|=(a[j][i]<<j-1);
    for(int i=1; i<=m; i++)
        f[s[i]]++;
    for(int i=0; i<lim; i++)
        g[i]=g[i>>1]+(i&1);
    for(int i=0; i<lim; i++)
        g[i]=min(g[i],n-g[i]);
    fwt(f,lim,1);
    fwt(g,lim,1);
    for(int i=0; i<lim; i++) f[i]=f[i]*g[i];
    fwt(f,lim,-1);
    int ans=1e9;
    for(int i=0; i<lim; i++) ans=min(ans,f[i]);
    printf("%lld\n",ans);
}