• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

CF662C Binary Table 【状压 + FWT】

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

CF662C

题解

行比较少,容易想到将每一列的状态压缩
在行操作固定的情况下,容易发现每一列的操作就是翻转\(0\)\(1\),要取最小方案,方案唯一
所以我们只需求出每一种操作的答案
如果操作的行的集合为\(S\),那么对于状态为\(e\)的列,将会变成\(e \; xor \; S\),同时产生\(e \; xor \; S\)的答案
如果\(s\)的答案记为\(b[s]\),状态为\(s\)的列数量为\(a[s]\)
那么对于操作\(S\),最后的答案为

\[\sum\limits_{i \; xor \; j = S}a[i] \centerdot b[j] \]

\(b[s]\)\(a[s]\)数组都可以预处理出来

\(N = 2^{20}\)
复杂度为\(O(NlogN)\)

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#include<map>
#define LL long long int
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define cls(s,v) memset(s,v,sizeof(s))
#define mp(a,b) make_pair<int,int>(a,b)
#define cp pair<int,int>
using namespace std;
const int maxn = 2100005,maxm = 100005,INF = 0x3f3f3f3f;
inline int read(){
	int out = 0,flag = 1; char c = getchar();
	while (c < 48 || c > 57){if (c == '-') flag = 0; c = getchar();}
	while (c >= 48 && c <= 57){out = (out << 1) + (out << 3) + c - 48; c = getchar();}
	return flag ? out : -out;
}
char S[22][100005];
int n,m;
LL A[maxn],B[maxn];
void fwt(LL* a,int n,int f){
	for (int i = 1; i < n; i <<= 1)
		for (int j = 0; j < n; j += (i << 1))
			for (int k = 0; k < i; k++){
				LL x = a[j + k],y = a[j + k + i];
				a[j + k] = x + y,a[j + k + i] = x - y;
				if (f == -1) a[j + k] /= 2,a[j + k + i] /= 2;
			}
}
int main(){
	n = read(); m = read();
	REP(i,n) scanf("%s",S[i] + 1);
	REP(j,m){
		int s = 0;
		REP(i,n) s = s << 1 | (S[i][j] - '0');
		A[s]++;
	}
	int maxv = (1 << n) - 1;
	for (int s = 0; s <= maxv; s++){
		int cnt = 0;
		for (int i = s; i; i >>= 1) cnt += (i & 1);
		B[s] = min(cnt,n - cnt);
	}
	int deg = 1;
	while (deg <= maxv) deg <<= 1;
	fwt(A,deg,1); fwt(B,deg,1);
	for (int i = 0; i < deg; i++) A[i] = A[i] * B[i];
	fwt(A,deg,-1);
	LL ans = INF;
	for (int i = 0; i <= maxv; i++) ans = min(ans,A[i]);
	printf("%lld\n",ans);
	return 0;
}


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
C语言位域发布时间:2022-07-14
下一篇:
C#堆与栈、值类型与引用类型、可空类型发布时间:2022-07-14
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap