CodeForces - 662C Binary Table （FWT）

题目链接：http://codeforces.com/problemset/problem/662/C

题目大意：给出一个n*m的01矩阵，你可以对矩阵一行或者一列进行变换，每次变换会使0变成1，1变成0。你可以进行无限次操作，现在问你怎么样操作才能使得矩阵内1的个数最少。

题目思路：由于本题的n比较小（<=20），所以我们可以用二进制状压来表示每一列的状态。就是表示第j列的状态。

我们先考虑变换列对于答案的影响，我们令 表示某一列的状态为 j 时，它选择变换或者不变换最少的1的个数。那么

。其中 表示 j 的二进制表示中1的个数。

接着我们考虑变换行对答案的影响。我们可以用一个数来表示行的变换状态，如果的二进制第 i 位为1的话，就说明变换第 i 行。那么现在的变换状态对第 j 列的影响就是 令 变为 了 。(表示异或运算)

那么每一个变换状态的最终情况就为。

由于所有  相同的列对于答案的贡献是相同的，那么这个式子就可以变化为

（变化后的式子表示，原状态为，经过变换后变为了的状态）。

接着我们再令，式子就变为了

。这个就可以借助FWT来求解了，最后再利用O(2^n)的枚举，找出最优的状态即可。

具体实现看代码：

#include <bits/stdc++.h>
#define fi first
#define se second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define lowbit(x) x&-x
#define pb push_back
#define MP make_pair
#define clr(a) memset(a,0,sizeof(a))
#define _INF(a) memset(a,0x3f,sizeof(a))
#define FIN freopen("in.txt","r",stdin)
#define fuck(x) cout<<"["<<x<<"]"<<endl
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int>pii;
typedef vector<int> VI;
const int MX = (1<<20) + 7;

void FWT_xor(ll *a,int N,int opt){
	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){
				ll X=a[j+k],Y=a[i+j+k];
				a[j+k]=(X+Y);a[i+j+k]=(X-Y);
				if(opt==-1){
					a[j+k]=1ll*a[j+k]/2;
					a[i+j+k]=1ll*a[i+j+k]/2;
				}
			}
		}
	}
}

int n,m;
char str[MX];
ll a[MX],num[MX],f[MX];

int main(){
	//FIN;
	ios::sync_with_stdio(false);
	cin>>n>>m;
	for(int i = 0;i < n;i++){
		cin>>str;
		for(int j = 0;j < m;j++){
			int x = str[j]-'0';
			a[j] = (a[j]<<1) + x;
		}
	}
	int N = (1<<n);
	for(int i = 0;i < m;i++) num[a[i]]++;
	for(int i = 1;i < N;i++){
		int cnt = __builtin_popcount(i);
		f[i] = min(cnt,n - cnt);
	}
	FWT_xor(num,N,1);FWT_xor(f,N,1);
	for(int i = 0;i < N;i++) f[i] = f[i] * num[i];
	FWT_xor(f,N,-1);
	ll ans = f[0];
	for(int i = 1;i < N;i++) ans = min(ans,f[i]);
	cout<<ans<<endl;
	return 0;
}