Codeforces #662C Binary Table

时间:2023-03-09 08:49:56
Codeforces #662C Binary Table

听说这是一道$ Tourist$现场没出的题

Codeforces #662C

题意:

给定$n*m的 01$矩阵,可以任意反转一行/列($0$变$1$,$1$变$0$),求最少$ 1$的数量

$ n<=20 \ m<=100000$


$ Solution$

考虑暴力

枚举每一行反转/不反转

预处理$ g(s)$表示某状态为$ s$的列的最少$ 1$的数量

显然$ g(s)=min(popcount(s),n-popcount(s))$

枚举每行是否反转之后直接$ O(m)$计算即可

时间复杂度$ O(2^n m)$,无法通过这题

容易发现瓶颈在于暴力枚举行状态之后无法快速计算答案

我们令$ f(s)$表示列状态为$ s$的列的出现次数,$ F(s)$表示行反转状态为$ s$的时候的答案

转移有$ F(s)=\sum\limits_{i=0}^{2^n-1}f(i)g(i \ xor \ s)$

由于$ i \ xor \ i \  xor \  s = s$

所以可以化简为$ F(s)=\sum\limits_{i \ xor \ j =s}f(i)g(j)$

是一个$ FWT$卷积的形式

直接$ FWT$优化

时间复杂度:$ O(nm+2^n n)$

注意$ FWT$过程中可能要开$ long \ long$


$ my \ code:$

#include<ctime>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#define rt register int
#define ll long long
using namespace std;
inline ll read(){
ll x = ; char zf = ; char ch = getchar();
while (ch != '-' && !isdigit(ch)) ch = getchar();
if (ch == '-') zf = -, ch = getchar();
while (isdigit(ch)) x = x * + ch - '', ch = getchar(); return x * zf;
}
void write(ll y){if(y<)putchar('-'),y=-y;if(y>)write(y/);putchar(y%+);}
void writeln(const ll y){write(y);putchar('\n');}
int i,j,k,m,n,x,y,z,cnt,invn;
void fwt(int n,ll *a,int fla){
for(rt i=;i<n;i<<=)
for(rt j=;j<n;j+=i<<)
for(rt k=;k<i;k++){
ll x=a[j+k],y=a[i+j+k];
a[j+k]=x+y;a[i+j+k]=x-y;
}
if(fla==-)for(rt i=;i<n;i++)a[i]/=n;
}
char c[][];
int s[];ll f[],g[];
#define cnt(x) __builtin_popcount(x)
int main(){
n=read();m=read();
for(rt i=;i<=n;i++)scanf("%s",c[i]+);
for(rt i=;i<=n;i++)
for(rt j=;j<=m;j++)s[j]=s[j]<<|(c[i][j]=='');
for(rt i=;i<=m;i++)g[s[i]]++;
for(rt i=;i<(<<n);i++)f[i]=min(cnt(i),n-cnt(i));
fwt(<<n,f,);fwt(<<n,g,);
for(rt i=;i<<<n;i++)f[i]=f[i]*g[i];
fwt(<<n,f,-);
cout<<*min_element(f,f+(<<n));
return ;
}