[CF662C]Binary Table

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

OneInDark

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

阅读量182

点赞数

分类专栏： C++ 数学

本文链接：https://blog.csdn.net/qq_42101694/article/details/108587729

版权

C++ 同时被 2 个专栏收录

547 篇文章 3 订阅

订阅专栏

数学

280 篇文章 1 订阅

订阅专栏

题目

传送门 to luogu

思路

$n$ 很小而 $m$ 很大，难道 $n$ 是可枚举的吗？

枚举之后每一列都重新看看不现实。因为我们是同时对 $m$ 列进行了操作。所以用桶排的思路，用 $C_x$ 表示列的状态是 $x$ 的数量。

然后用 $f (x)$ 表示列状态为 $x$ 时的最小 $1$ 数量。那么假如行的翻转情况为 $t$ ，则原来的 $x$ 就变成了 $x\oplus t$ ，所以答案是

$\sum_{x=0}^{2^n-1}C_{x}f(x\oplus t)$

注意到二者的下标的异或值固定，所以

$ans_t=\sum_{a\oplus b=t}C_af(b)$

这就是一个按位异或卷积。可以做到 $\mathcal O(nm+n2^n)$ 了。

代码

#include <cstdio>
#include <iostream>
#include <vector>
using namespace std;
typedef long long int_;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(int x){
	if(x > 9) writeint(x/10);
	putchar((x%10)^48);
}
inline int qkpow(int_ b,int q,int m){
	int ans = 1;
	for(; q; q>>=1,b=b*b%m)
		if(q&1) ans = ans*b%m;
	return ans;
}

const int Mod = 998244353;
const int inv2 = (Mod+1)>>1;
void FWT(int a[],int n,int opt){
	for(int l=1; l<n; l<<=1)
	for(int *p=a; p!=a+n; p+=(l<<1))
	for(int i=0; i<l; ++i){
		int x = p[i], y = p[i+l];
		p[i] = (x+y)%Mod;
		p[i+l] = (x-y+Mod)%Mod;
		if(opt == -1){
			p[i] = 1ll*p[i]*inv2%Mod;
			p[i+l] = 1ll*p[i+l]*inv2%Mod; 
		}
	}
}

const int MaxN = 20;
int cnt[1<<MaxN], f[1<<MaxN];
int a[100005]; // 每一列的状压
int main(){
	int n = readint(), m = readint();
	for(int i=1; i<=n; ++i){
		for(int j=1; j<=m; ++j)
			a[j] = (a[j]<<1)|(getchar()-'0');
		getchar(); // get '\n'
	}
	for(int i=1; i<=m; ++i)
		++ cnt[a[i]]; // 桶排数量
	for(int S=1; S<(1<<n); ++S)
		f[S] = f[S-(S&-S)]+1;
	for(int S=0; S<(1<<n); ++S)
		if(f[S] > n-f[S])
			f[S] = n-f[S];
	FWT(cnt,1<<n,1), FWT(f,1<<n,1);
	for(int i=0; i<(1<<n); ++i)
		f[i] = 1ll*f[i]*cnt[i]%Mod;
	FWT(f,1<<n,-1);
	int ans = n*m;
	for(int i=0; i<(1<<n); ++i){
		ans = min(ans,f[i]);
// printf("f[%d] = %d\n",i,f[i]);
	}
	printf("%d\n",ans);
	return 0;
}