<div>北京大学掲載のプログラムの練習問題。</div> <div>問題No1038を解く。</div> <div>動的計画法にBit演算まで動員して解いたのだけど速度が出なかった、コード実行速度下位ランクという結果に。</div> <div>2重再帰の深さ優先探索で良かったのだろうか?</div> <div>しかしそれだと縦に入れる横に入れるの分岐で組み合わせ数が凄いことになりそうなのだけど???</div> <div> </div> <div> </div> <div>#include<stdio.h></div> <div>#include<map></div> <div>#include<string.h></div> <div> </div> <div>//グローバル変数まとめ</div> <div>int map[153];//1bit一マスで盤面を表す</div> <div>const int okBlock=0;</div> <div>const int badBlock=1;</div> <div>int n,m,k;</div> <div>//まとめ終わり</div> <div> </div> <div>void lineCalc(int row,int col,int count,int kekka,std::map<int,int>& next){</div> <div>//一列だけ走査していく</div> <div>if(row==m){</div> <div>kekka=kekka>>m;//一行分削除する</div> <div>if(next.find(kekka)==next.end()){</div> <div>next[kekka]=count;//初めてのパタン</div> <div>}else if(next[kekka]<count){</div> <div>next[kekka]=count;//より大きいパタンが見つかった</div> <div>}</div> <div>return ;</div> <div>}</div> <div>//今のマスが不良品かチップで埋まっているなら</div> <div>if((kekka>>(row-1))&1){</div> <div>lineCalc(row+1,col,count,kekka,next);</div> <div>return;</div> <div>}</div> <div>//このマスに何も入れない</div> <div>lineCalc(row+1,col,count,kekka,next);</div> <div> </div> <div>//横を入れれるかチェックする</div> <div>int chip=3<<(row-1);</div> <div>int col0=kekka&chip;</div> <div>int col1=kekka&(chip<<m);</div> <div>int col2=kekka&(chip<<(m*2));</div> <div>if(col0==0 && col1==0 && col2==0){</div> <div> lineCalc(row+1,col,count+1,kekka|chip|(chip<<m)|(chip<<(m*2)),next);//横を入れた</div> <div>}</div> <div>if(row>=m-1)return;//縦を入れれないマスに到達。</div> <div>chip=7<<(row-1);</div> <div>col0=kekka&chip;</div> <div>col1=kekka&(chip<<m);</div> <div>//縦を入れれるかチェックする</div> <div>if(col0==0 && col1==0){</div> <div>lineCalc(row+1,col,count+1,kekka|chip|(chip<<m),next);//縦を入れた</div> <div>}</div> <div>}</div> <div> </div> <div> </div> <div>void setData(){</div> <div>//一行ずつ動的計画法で計算していく。</div> <div>memset(map,okBlock,sizeof(map));</div> <div>int x,y;</div> <div>scanf("%d %d %d",&n,&m,&k);</div> <div>while(k--){</div> <div>scanf("%d %d",&x,&y);</div> <div>map[x]|=badBlock<<(y-1);//1bit一マスで一列ずつ升目を格納する。</div> <div>}</div> <div>map[n+1]=~((int)0);//番兵</div> <div>std::map<int,int> memo,next;//</div> <div>std::map<int,int>::iterator it;</div> <div>memo[0]=0;</div> <div>for(int col=1;col<n;col++){</div> <div>next.clear();</div> <div>int badLines=map[col]|(map[col+1]<<m)|(map[col+2]<<(2*m));//邪魔なブロックを記録する</div> <div> </div> <div>for(it=memo.begin();it!=memo.end();it++){</div> <div> lineCalc(1,col,(*it).second,(*it).first|badLines,next);//今のラインで動的計画法を試す</div> <div>}</div> <div>memo.clear();</div> <div>memo.insert(next.begin(),next.end());</div> <div>}</div> <div>int ans=0;</div> <div>for(it=memo.begin();it!=memo.end();it++){</div> <div>if(ans==0||(*it).second>ans){</div> <div>ans=(*it).second;</div> <div>}</div> <div>}</div> <div>printf("%d\n",ans);</div> <div>}</div> <div> </div> <div>int main(){</div> <div>int d;</div> <div>scanf("%d",&d);</div> <div>while(d--){</div> <div>setData();</div> <div>}</div> <div>}</div>