結果
問題 | No.421 しろくろチョコレート |
ユーザー |
![]() |
提出日時 | 2015-02-12 05:47:00 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 4 ms / 2,000 ms |
コード長 | 2,789 bytes |
コンパイル時間 | 1,982 ms |
コンパイル使用メモリ | 87,356 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-09-23 07:03:54 |
合計ジャッジ時間 | 3,235 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 65 |
ソースコード
#include <iostream>#include <vector>#include <cstdio>#include <sstream>#include <map>#include <string>#include <algorithm>#include <queue>#include <cmath>using namespace std;typedef struct{int to;int cap;int rev;}edge;// node v : distance from s => level[v]void bfs(vector<vector<edge> > &G, vector<int> &level, int s){fill(level.begin(), level.end(), -1);queue<int> q;q.push(s);level[s] = 0;while(!q.empty()){int e=q.front(); q.pop();for(int i=0; i<G[e].size(); i++){if(G[e][i].cap > 0 && level[G[e][i].to] < 0){level[G[e][i].to] = level[e] + 1;q.push(G[e][i].to);}}}}int dfs(vector<vector<edge> > &G, vector<int> &level, vector<bool> &used, vector<int> &iter, int s, int f, int t){if(s==t) return f;else{//iter[e] : done searching from v[0] to v[ iter[e]-1 ]for(int &i=iter[s]; i<G[s].size(); i++){//distance from s to v[e][i].to must be longer than dist from s to vif(G[s][i].cap > 0 && level[s] < level[ G[s][i].to ]){int d = dfs(G, level, used, iter, G[s][i].to, min(f, G[s][i].cap), t);if(d>0){G[s][i].cap -= d;G[ G[s][i].to ][ G[s][i].rev ].cap += d;return d;}}}return 0;}}int dinic_maxflow(vector<vector<edge> > &G, int s, int t){const int INF = 100000000;int flow=0;while(true){vector<int> level(G.size(), -1);vector<int> iter(G.size(), 0);vector<bool> used(G.size(), false);bfs(G, level, s);if(level[t] < 0) return flow; //unable to achieve to twhile(true){int f = dfs(G, level, used, iter, s, INF, t);if(f==0) break;else flow += f;}}}void add_edge(vector<vector<edge> > &G, int from, int to, int cap){G[from].push_back((edge){to, cap, (int)G[to].size()});G[to].push_back((edge){from, 0, (int)G[from].size() - 1});}int main(){int n,m;cin >> n >> m;vector<string> v(n);for(int i=0; i<n; i++){cin >> v[i];}// s -> w -> b -> tvector< vector<edge> > G(n*m + 2);int s = n*m;int t = n*m+1;int cnt_w = 0;int cnt_b = 0;for(int i=0; i<n; i++){for(int j=0; j<m; j++){if(v[i][j] == 'w'){cnt_w++;add_edge(G, s, i*m + j, 1);if(j+1<m && v[i][j+1] == 'b'){add_edge(G, i*m + j, i*m + j+1, 1);}if(i+1<n && v[i+1][j] == 'b'){add_edge(G, i*m + j, (i+1)*m + j, 1);}}else if(v[i][j] == 'b'){cnt_b++;add_edge(G, i*m + j, t, 1);if(j+1<m && v[i][j+1] == 'w'){add_edge(G, i*m + j+1, i*m + j, 1);}if(i+1<n && v[i+1][j] == 'w'){add_edge(G, (i+1)*m + j, i*m + j, 1);}}}}int bw = dinic_maxflow(G, s,t);int ans = 0;ans += bw*100;cnt_w -= bw;cnt_b -= bw;ans += min(cnt_w, cnt_b)*10;ans += max(cnt_w, cnt_b) - min(cnt_w, cnt_b);cout << ans << endl;return 0;}