CF662C Binary Table

最新推荐文章于 2022-01-04 21:00:13 发布

C202044zxy

最新推荐文章于 2022-01-04 21:00:13 发布

阅读量210

点赞数

分类专栏： FWT

本文链接：https://blog.csdn.net/C202044zxy/article/details/106738575

版权

FWT 专栏收录该内容

12 篇文章 0 订阅

订阅专栏

该博客介绍了CF662C题目解法，通过观察n值较小的特点，采用状态压缩和位运算的方法求解。通过枚举列的翻转状态和计算二进制中1的个数，利用XOR卷积达到O(n^2^n)的时间复杂度。

摘要由CSDN通过智能技术生成

一、题目

点此看题

二、解法

观察到 $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$ 直接搞 $x o r$ 卷积就可以了，时间复杂度 $O(n2^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);
}