Problem:
おせんべい

問題

IOI製菓では，創業以来の伝統の製法で煎餅（せんべい）を焼いている．この伝統の製法は，炭火で一定時間表側を焼き，表側が焼けると裏返して，炭火で一定時間裏側を焼くというものである．この伝統を守りつつ，煎餅を機械で焼いている．この機械は縦 R (1 ≤ R ≤ 10) 行, 横 C (1 ≤ C ≤ 10000) 列の長方形状に煎餅を並べて焼く．通常は自動運転で，表側が焼けたら一斉に煎餅を裏返し裏側を焼く．

ある日，煎餅を焼いていると，煎餅を裏返す直前に地震が起こり何枚かの煎餅が裏返ってしまった．幸いなことに炭火の状態は適切なままであったが，これ以上表側を焼くと創業以来の伝統で定められている焼き時間を超えてしまい，煎餅の表側が焼けすぎて商品として出荷できなくなる．そこで，急いで機械をマニュアル操作に変更し，まだ裏返っていない煎餅だけを裏返そうとした．この機械は，横の行を何行か同時に裏返したり縦の列を何列か同時に裏返したりすることはできるが，残念なことに，煎餅を１枚ごと裏返すことはできない．

裏返すのに時間がかかると，地震で裏返らなかった煎餅の表側が焼けすぎて商品として出荷できなくなるので，横の何行かを同時に１回裏返し，引き続き，縦の何列かを同時に１回裏返して，表側を焼きすぎずに両面を焼くことのできる煎餅，つまり，「出荷できる煎餅」の枚数をなるべく多くすることにした．横の行を１行も裏返さない，あるいは，縦の列を１列も裏返さない場合も考えることにする．出荷できる煎餅の枚数の最大値を出力するプログラムを書きなさい．

地震の直後に，煎餅が次の図のような状態になったとする．黒い丸が表側が焼ける状態を，白い丸が裏側が焼ける状態を表している．

1行目を裏返すと次の図のような状態になる．

さらに， 1列目と5列目を裏返すと次の図のような状態になる．この状態では，出荷できる煎餅は9枚である．

ヒント

R の上限 10 は C の上限 10000 に比べて小さいことに注意せよ．

入力

入力は複数のデータセットからなる．各データセットは以下の形式で与えられる．

入力の1行目には2つの整数 R, C (1 ≤ R ≤ 10, 1 ≤ C ≤ 10 000) が空白を区切りとして書かれている．続く R 行は地震直後の煎餅の状態を表す． (i+1) 行目 (1 ≤ i ≤ R) には， C 個の整数 ai,1, ai,2, ……, ai,C が空白を区切りとして書かれており， ai,j は i 行 j 列 の煎餅の状態を表している. ai,j が 1 なら表側が焼けることを， 0 なら裏側が焼けることを表す．

C, R がともに 0 のとき入力の終了を示す.　データセットの数は 5 を超えない．

出力

データセットごとに，出荷できる煎餅の最大枚数を1行に出力する．

入出力例

入力例

2 5
0 1 0 1 0
1 0 0 0 1
3 6
1 0 0 0 1 0
1 1 1 0 1 0
1 0 1 1 0 1
0 0

出力例

9
15

上記問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C++ (judge verdict: Time Limit Exceeded)

#include<bits/stdc++.h>
#define MAX_C 10000
using namespace std;
int r,c;
int senbei[1<<10][MAX_C+1];

int bitCnt(int n){
  int cnt=0;
  while(n>0){
    if(n%2==1)cnt++;
    n/=2;
  }
  if(r-cnt<cnt){
    return r-cnt;
  }else{
    return cnt;
  }
}

int rec(int s,int n){
  if(senbei[s][n]!=-1){
    return senbei[s][n];
  }
  int tmp=s;
  for(int i=0;tmp>0;i++){
    if(tmp%2==1){
      for(int j=0;j<c;j++){
	int bitpos=1<<i;
	senbei[s][j]=rec(s^bitpos,j)^bitpos;
      }
    }
    tmp/=2;
  }
  return senbei[s][n];
}

int minimam(){
  int ans=1<<30;
  int tmp=0;
  for(int i=0;i<1<<r;i++){
    tmp=0;
    cout << "reset" << endl;
    for(int j=0;j<c;j++){
      tmp+=bitCnt(senbei[i][j]);
      cout<<"tmp is "<<tmp<<endl;
    }
    if(ans>tmp)ans=tmp;
  }
  return ans;
}

int main(){
  while(cin>>r>>c){
    if(c==0)break;
    int tmp[c];
    memset(tmp,0,sizeof(tmp));
    memset(senbei,-1,sizeof(senbei));
    for(int i=0;i<r;i++){
      for(int j=0;j<c;j++){
	tmp[j]=tmp[j]<<1;
	int bit;
	cin>>bit;
	if(bit==1){
	  tmp[j]=tmp[j]|1;
	}
      }
    }
    for(int i=0;i<c;i++){
      senbei[0][i]=tmp[i];
      //cout << tmp[i] << endl;
    }
    rec((1<<r)-1,0);
    for(int i=0;i<1<<r;i++){
      for(int j=0;j<c;j++){
	cout<<senbei[i][j];
      }
      cout<<endl;
    }
    cout<<r*c-minimam()<<endl;
  }
  return 0;
}