题面
题意:
有一个 n n nx n n n的区域。黑色为墙,白色为空白,你现在要在空白区域安装大炮。大炮的可以摧毁同行和同列的所有物品,但是大炮无法摧毁墙。求为了避免大炮之间两两攻击,最多放几门大炮。
IDEA:
我们构造二分图,左面的 n n n个点为行,右面的 n n n个点为列。如果对于第 i i i行有 n u m num num个不连续的空白区域,说明第 i i i行最多可以和 n u m num num列进行匹配。如果第 j j j列有 n u m num num个不连续的空白区域,那么说明第 j j j列最多可以和 n u m num num行进行匹配。建好图用 D i n i c Dinic Dinic跑二分图最大匹配即可。
code
#include<bits/stdc++.h>
using namespace std;
const int INF = 0x3f3f3f3f;
typedef long long ll;
//Dinic最大流, 节点编号从0开始
struct MaxFlow{
const static ll MAX_V = 1005;
ll V;
//终点、容量、反向边
struct edge{
ll to, cap, rev;
};
vector<edge> G[MAX_V];
ll level[MAX_V];//顶点到源点的距离标号
ll iter[MAX_V];// 当前弧,在其之前的边已经没有用了
void add_edge(ll from, ll to, ll cap){
G[from].push_back((edge){to, cap, (ll)G[to].size()});
G[to].push_back((edge){from, 0, (ll)G[from].size()-1});
}
// 通过BFS计算从源点出发的距离标号
void bfs(ll s){
fill(level, level + V, -1);
queue<ll> que;
level[s] = 0;
que.push(s);
while (!que.empty()){
ll v = que.front();
que.pop();
for (ll i=0; i< G[v].size(); i++){
edge& e = G[v][i];
if (e.cap > 0 && level[e.to] < 0){
level[e.to] = level[v] + 1;
que.push(e.to);
}
}
}
}
//通过DFS寻找增广路
ll dfs(ll v, ll t, ll f){
if (v == t)
return f;
for (ll &i = iter[v]; i < G[v].size(); i++){
edge& e = G[v][i];
if (e.cap > 0 && level[v] < level[e.to]){
ll 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;
}
//求解从s到t的最大流
ll max_flow(ll s, ll t){
ll flow = 0;
for (;;){
bfs(s);
if(level[t] < 0)
return flow;
fill(iter, iter + V, 0);
ll f;
while ((f = dfs(s, t, INF)) > 0){
flow += f;
}
}
}
void init(ll n = 0){
for (ll i = 0; i < V; i++){
G[i].clear();
}
V = n;
}
}mf;
char maze[10][10];
int main(){
int n;
while(scanf("%d",&n)&&n){
mf.init();
mf.V=2*n+2;
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) cin>>maze[i][j];
}
for(int i=1;i<=n;i++){
int num=0;
for(int j=1;j<=n;){
if(maze[i][j]=='.'){
num++;
while(maze[i][j]=='.') j++;
}
else{
while(maze[i][j]=='X') j++;
}
}
mf.add_edge(0,i,num);
mf.add_edge(i,0,0);
}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
if(maze[i][j]=='.'){
mf.add_edge(i,j+n,1);
mf.add_edge(j+n,i,0);
}
}
}
for(int j=1;j<=n;j++){
int kk=0;
for(int i=1;i<=n;){
if(maze[i][j]=='.'){
kk++;
while(maze[i][j]=='.') i++;
}
else{
while(maze[i][j]=='X') i++;
}
}
mf.add_edge(j+n,2*n+1,kk);
}
int ans=mf.max_flow(0,2*n+1);
printf("%d\n",ans);
}
return 0;
}