【CF662C】Binary Table(FWT)

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

原文地址:https://www.cnblogs.com/cjyyb/p/9065801.html

时间: 2024-11-02 23:43:50

【CF662C】Binary Table(FWT)的相关文章

【CF662C】Binary Table 按位处理

[CF662C]Binary Table 题意:给你一个$n\times m$的01网格,你可以进行任意次操作,每次操作是将一行或一列的数都取反,问你最多可以得到多少个1? $n\le 20,m\le 10^5$ 题解:我也不知道叫啥了,说状压也不对,说fwt也不太对,就叫按位处理得了. 显然有$O(2^nm)$暴力,先枚举每行是否取反,然后枚举每列,如果0多就取反,否则不取. 但我们发现我们完全可以将本质相同的列一起处理,什么叫本质相同的列呢?假如我们对每行是否取反的状态为S,则所有$xor

【HDU5909】Tree Cutting(FWT)

[HDU5909]Tree Cutting(FWT) 题面 vjudge 题目大意: 给你一棵\(n\)个节点的树,每个节点都有一个小于\(m\)的权值 定义一棵子树的权值为所有节点的异或和,问权值为\(0..m-1\)的所有子树的个数 题解 考虑\(dp\) 设\(f[i][j]\)表示以\(i\)为根节点的子树中,异或和为\(j\)的子树的个数 很显然,每次合并就是两个\(dp\)值做\(xor\)卷积 那么直接用\(FWT\)优化就行了 #include<iostream> #inclu

【CF662C】Binary Table

题目 好吧,我连板子都不会了 有一个非常显然的做法就是\(O(2^nm)\)做法就是枚举每一行的状态,之后我们贪心去看看每一列是否需要翻转就好啦 显然这个做法非常垃圾过不去 首先我们发现每一列都不超过\(20\),考虑把每一列都压成一个状态 我们考虑设一些奇怪的东西 设\(g_i\)表示行的翻转状态为\(i\)的最优解,\(f_i\)表示有多少列的状态为\(i\),\(dp_i\)表示\(i\)这个状态最少有多少个\(1\) 显然\(dp_i=min\{bit(i),n-bit(i)\}\) 我

【tornado】系列项目(二)基于领域驱动模型的区域后台管理+前端easyui实现

本项目是一个系列项目,最终的目的是开发出一个类似京东商城的网站.本文主要介绍后台管理中的区域管理,以及前端基于easyui插件的使用.本次增删改查因数据量少,因此采用模态对话框方式进行,关于数据量大采用跳转方式修改,详见博主后续博文. 后台界面展示: 地区管理包含省市县的管理.详见下文. 一.数据库设计 ? 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 class Province(Base):  

【Chrome】Chrome插件开发(一)插件的简单实现

不同浏览器插件开发比较 Chrome的插件开发起来最简单,总体上看没什么新的技术,开发语言就是javascript,web前端工程师能很快上手. Firefox的插件开发则复杂许多,涉及到环境的搭建和一些WEB以外的技术. IE的插件开发就更复杂了,需要熟悉C++和COM技术,当然还要装微软的Visual Studio. 这里有篇老外写的文章,对比Chrome.Opera和Firefox的插件开发的:http://blog.nparashuram.com/2011/10/writing-brow

【TYVJ】1307 联络员(最小生成树)

http://tyvj.cn/Problem_Show.aspx?id=1307 kruskal裸题.(水题红色警报) #include <cstdio> #include <cstring> #include <cmath> #include <string> #include <iostream> #include <algorithm> using namespace std; #define rep(i, n) for(int

【转】引导界面(五)实现应用程序只启动一次引导界面

这篇文章算是对整个引导界面开发专题的一个终结了吧,个人觉得大部分的引导界面基本上都是千篇一律的,只要熟练掌握了一个,基本上也就没什么好说的了,要是在今后的开发中遇到了更好玩,更有趣的引导界面,博主也会在这里及时的跟大家分享,今天的内容主要是教大家的应用程序只有在第一次启动的时候显示引导界面,以后在启动程序的时候就不再显示了. 其实要想实现这样的效果,只要使用SharedPreferences类,就会让程序变的非常简单,下面来详细介绍一下这个类的使用方法 一.SharedPreferences的详

【转】Robust regression(稳健回归)

Robust regression(稳健回归) 语法 b=robustfit(X,y) b=robustfit(X,y,wfun,tune) b=robustfit(X,y,wfun,tune,const) [b,stats]=robustfit(...) 描述 b=robustfit(X,y) 通过执行稳健回归来估计线性模型y=Xb,并返回一个由回归系数组成的向量b.X是一个n*p预测变量矩阵,y是一个n*1观测向量.计算使用的方法是加 上bisquare加权函数的迭代重加权最小二乘法.默认的

【转载】PHP的(&lt;&lt;&lt;EOT)在PHP中添加html

在W3school上学PHP,看到第一句就是"PHP 文件可包含文本.HTML 标签以及脚本"在后来的学习别人的代码,发现在需要HTML代码的PHP脚本中,多用这么几种方法第一种是在HTML中加PHP.大段大段的html代码中,在各个需要执行php的地方.这种方法在ASP的程序中比较常见.例子:     <head>    <meta http-equiv="Content-Type" content="text/html; charse