CF 662C Binary Table

Binary Table

题目描述
照例贴传送门：http://codeforces.com/problemset/problem/662/C

题解

暴力做法
显然每一行每一列只会被修改一次或者不被修改
2^n枚举每一行是否被修改
然后O(m)扫一遍每一列，每一列初始状态为a[i]，在经过状态为S的修改之后的状态为T=S^a[i]
ans=min(Σmin(f(T),n-f(T)))，f(T)表示T中1的个数
所以对于T相同的状态对答案的贡献是一样的
T=S^a[i] a[i]^T=S
枚举S，ans=Σsum[a[i]]*min(n-f[T],f[T])，sum[a[i]]表示a[i]这种状态的个数，其中a[i]^T=S
这东西可以用fwt搞
fwt就是fft把i+j换成i^j或i|j或i&j
并不会写fwt，去抄了个板子（都变成acm选手了抄抄板子好啦）
和fft有点像但还是有很大差别的

代码

#include<bits/stdc++.h>
#define ll long long
#define N 22
#define M 100005
#define P 2100000
using namespace std;
int n,m,S,s[N][M];ll a[P],b[P],ans;

void FWT(ll *a,int n)
{
  for(int d=1;d<n;d<<=1)
    for(int m=d<<1,i=0;i<n;i+=m)
      for(int j=0;j<d;j++)
      {
        ll x=a[i+j],y=a[i+j+d];
        a[i+j]=x+y;a[i+j+d]=x-y;
        //xor:a[i+j]=x+y,a[i+j+d]=x-y;  
        //and:a[i+j]=x+y;  
        //or:a[i+j+d]=x+y;  
      }
}

void UFWT(ll *a,int n)
{
  for(int d=1;d<n;d<<=1)
    for(int m=d<<1,i=0;i<n;i+=m)
      for(int j=0;j<d;j++)
      {
        ll x=a[i+j],y=a[i+j+d];
        a[i+j]=x+y>>1,a[i+j+d]=x-y>>1;
        //xor:a[i+j]=(x+y)/2,a[i+j+d]=(x-y)/2;  
        //and:a[i+j]=x-y;  
        //or:a[i+j+d]=y-x; 
      }
}

int main()
{
  char c;
  scanf("%d%d",&n,&m);S=(1<<n)-1;ans=n*m;
  for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
      scanf(" %c",&c),s[i][j]=c-48;
  for(int i=1;i<=m;i++)
  {
    int num=0;
    for(int j=1;j<=n;j++)
      num=(num<<1)+s[j][i];
    a[num]++;
  }
  for(int i=0;i<=S;i++)
  {
    int cnt=0;
    for(int x=i;x;x-=x&-x)cnt++;
    b[i]=min(cnt,n-cnt);
  }
  FWT(a,S);FWT(b,S);
  for(int i=0;i<=S;i++)a[i]*=b[i];
  UFWT(a,S);
  for(int i=0;i<=S;i++)ans=min(ans,a[i]);
  printf("%I64d\n",ans);
  return 0;
}