結果
| 問題 | No.421 しろくろチョコレート |
| コンテスト | |
| ユーザー |
 tsutaj
|
| 提出日時 | 2019-02-16 14:13:27 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.89.0) |
| 結果 |
AC
|
| 実行時間 | 19 ms / 2,000 ms |
| コード長 | 7,899 bytes |
| 記録 | |
| コンパイル時間 | 1,043 ms |
| コンパイル使用メモリ | 112,768 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-09-23 07:22:53 |
| 合計ジャッジ時間 | 2,711 ms |
|
ジャッジサーバーID （参考情報） |
judge5 / judge3 |
（要ログイン）
| ファイルパターン | 結果 |
|---|---|
| other | AC * 65 |
ソースコード
// #define _GLIBCXX_DEBUG // for STL debug (optional)
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <string>
#include <cstring>
#include <deque>
#include <list>
#include <queue>
#include <stack>
#include <vector>
#include <utility>
#include <algorithm>
#include <map>
#include <set>
#include <complex>
#include <cmath>
#include <limits>
#include <cfloat>
#include <climits>
#include <ctime>
#include <cassert>
#include <numeric>
#include <fstream>
#include <functional>
#include <bitset>
using namespace std;
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define int long long int
template<typename T> void chmax(T &a, T b) {a = max(a, b);}
template<typename T> void chmin(T &a, T b) {a = min(a, b);}
template<typename T> void chadd(T &a, T b) {a = a + b;}
typedef pair<int, int> pii;
typedef long long ll;
int dx[] = {0, 0, 1, -1};
int dy[] = {1, -1, 0, 0};
const ll INF = 1001001001001001LL;
const ll MOD = 1000000007LL;
// Ford-Fulkerson 法による 最大流 O( F |E| )
// Bellman-Ford 法による 最小費用流 O( F |V| |E| )
// [条件に注意] Dijkstra 法による 最小費用流 O( F |E| log |V| )
// Verified: AOJ GRL_6_A (Maximum Flow)
// Verified: AOJ GRL_6_B (Minimum Cost Flow)
// 行き先と容量と逆辺のインデックスを記録する構造体
// 通常のグラフの辺の構造体と異なるため注意
template <typename Type>
struct Edge {
int to; Type cap, cost; int rev;
Edge(int t, Type c, int r) : to(t), cap(c), rev(r) {}
Edge(int t, Type ca, Type co, int r) : to(t), cap(ca), cost(co), rev(r) {}
};
template <typename Type>
using Graph = vector< vector< Edge<Type> > >;
template <typename Type>
struct Flow {
vector< vector< Edge<Type> > > G;
const Type MAXC = 1 << 30;
int n;
vector<bool> used;
vector<int> prevv, preve, dist;
Flow(int _n) : G(_n), n(_n), used(_n, false),
prevv(_n), preve(_n), dist(_n, MAXC) {}
// G[e.to][e.rev] で逆辺を操作できる
void add_edge(int from, int to, Type cap) {
G[from].push_back(Edge<Type>(to, cap, G[to].size()));
G[to].push_back(Edge<Type>(from, 0, G[from].size() - 1));
}
void add_edge(int from, int to, Type cap, Type cost) {
G[from].push_back(Edge<Type>(to, cap, cost, G[to].size()));
G[to].push_back(Edge<Type>(from, 0, -cost, G[from].size() - 1));
}
Type dfs(int v, int t, Type f) {
if(v == t) return f;
used[v] = true;
for(int i=0; i < G[v].size(); i++) {
Edge<Type> &e = G[v][i];
if(!used[e.to] && e.cap > 0) {
Type d = dfs(e.to, t, min(f, e.cap));
if(d > 0) {
e.cap -= d;
G[e.to][e.rev].cap += d;
return d;
}
}
}
return 0;
}
Type max_flow(int s, int t) {
Type flow = 0;
while(1) {
fill(used.begin(), used.end(), false);
Type f = dfs(s, t, INF);
if(f == 0) return flow;
flow += f;
}
}
// ベルマンフォード法で最小費用流を解く
// 負辺や負閉路があるときは基本的にこちらのほうが確実
Type mincost_flow(int s, int t, Type f) {
Type res = 0;
Type ma = MAXC;
while(f > 0) {
fill(dist.begin(), dist.end(), ma);
dist[s] = 0;
bool update = true;
while(update) {
update = false;
for(int v = 0; v < n; v++) {
if(dist[v] == ma) continue;
for(int i=0; i<G[v].size(); i++) {
Edge<Type> &e = G[v][i];
if(e.cap>0 && dist[e.to] > dist[v] + e.cost) {
dist[e.to] = dist[v] + e.cost;
prevv[e.to] = v;
preve[e.to] = i;
update = true;
}
}
}
}
if(dist[t] == ma) return -1;
Type d = f;
for(int v = t; v != s; v = prevv[v]) {
d = min(d, G[prevv[v]][preve[v]].cap);
}
f -= d;
res += d * dist[t];
for(int v = t; v != s; v = prevv[v]) {
Edge<Type> &e = G[prevv[v]][preve[v]];
e.cap -= d;
G[v][e.rev].cap += d;
}
}
return res;
}
// ポテンシャルの導入により、ダイクストラ法で最小費用流を解く
// [仮定している条件]
// 1. グラフに負の閉路が存在しない (流量の 0 初期化のため)
// もし存在するならベルマンフォードで負の閉路を見つけ
// そこに流せるだけ流してスタート
// 2. グラフに負の辺が存在しない (pot_0 の計算可能性)
// もし存在する場合は最初のみベルマンフォードを使う必要あり
Type fast_mincost_flow(int s, int t, Type f) {
Type res = 0, ma = MAXC;
vector<Type> pot(n);
while(f > 0) {
using PT = pair<Type, int>;
priority_queue< PT, vector<PT>, greater<PT> > que;
fill(dist.begin(), dist.end(), ma);
dist[s] = 0;
que.push(make_pair(0, s));
while(!que.empty()) {
PT cur = que.top(); que.pop();
int v = cur.second;
if(dist[v] < cur.first) continue;
for(size_t i=0; i<G[v].size(); i++) {
Edge<Type>& e = G[v][i];
if(e.cap > 0 and dist[e.to] > dist[v] + e.cost + pot[v] - pot[e.to]) {
dist[e.to] = dist[v] + e.cost + pot[v] - pot[e.to];
prevv[e.to] = v;
preve[e.to] = i;
que.push(make_pair(dist[e.to], e.to));
}
}
}
if(dist[t] == ma) {
return -1;
}
for(int v=0; v<n; v++) pot[v] += dist[v];
Type d = f;
for(int v=t; v!=s; v=prevv[v]) {
d = min(d, G[prevv[v]][preve[v]].cap);
}
f -= d;
res += d * pot[t];
for(int v=t; v!=s; v=prevv[v]) {
Edge<Type>& e = G[prevv[v]][preve[v]];
e.cap -= d;
G[v][e.rev].cap += d;
}
}
return res;
}
};
char board[55][55];
signed main() {
int N, M; cin >> N >> M;
for(int i=0; i<N; i++) {
for(int j=0; j<M; j++) {
cin >> board[i][j];
}
}
Flow<int> fl(N*M + 2);
int source = N*M, sink = source + 1;
int cntW = 0, cntB = 0;
for(int i=0; i<N; i++) {
for(int j=0; j<M; j++) {
if(board[i][j] == '.') continue;
int u = i * M + j;
if((i + j) % 2 == 1) {
// 'b'
fl.add_edge(u, sink, 1);
cntB++;
}
else {
// 'w'
fl.add_edge(source, u, 1);
cntW++;
for(int k=0; k<4; k++) {
int ni = i + dx[k];
int nj = j + dy[k];
if(ni < 0 or ni >= N or nj < 0 or nj >= M) continue;
if(board[ni][nj] == '.') continue;
int v = ni * M + nj;
fl.add_edge(u, v, 1);
}
}
}
}
int cntPair = fl.max_flow(source, sink);
cntB -= cntPair, cntW -= cntPair;
int cntApartPair = min(cntB, cntW);
cntB -= cntApartPair, cntW -= cntApartPair;
int ans = cntPair * 100 + cntApartPair * 10 + cntB + cntW;
cout << ans << endl;
return 0;
}