Description
在一个r行c列的网格地图中有一些高度不同的石柱,一些石柱上站着一些蜥蜴,你的任务是让尽量多的蜥蜴逃
到边界外。 每行每列中相邻石柱的距离为1,蜥蜴的跳跃距离是d,即蜥蜴可以跳到平面距离不超过d的任何一个石
柱上。石柱都不稳定,每次当蜥蜴跳跃时,所离开的石柱高度减1(如果仍然落在地图内部,则到达的石柱高度不
变),如果该石柱原来高度为1,则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个
石柱上。
Input
输入第一行为三个整数r,c,d,即地图的规模与最大跳跃距离。以下r行为石竹的初始状态,0表示没有石柱
,1~3表示石柱的初始高度。以下r行为蜥蜴位置,“L”表示蜥蜴,“.”表示没有蜥蜴。
Output
输出仅一行,包含一个整数,即无法逃离的蜥蜴总数的最小值。
Sample Input
5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........
Sample Output
1
思路:
这个题是点上有权值,所以两个点之间的容量不好赋值,所以我们想到拆点,一个点分成两个点,然后这两个点之间的容量就是这个点的权值。
跑最大流就行了。
#include <bits/stdc++.h>
#define mem(x,v) memset(x,v,sizeof(x))
#define go(i,a,b) for (int i = a; i <= b; i++)
#define og(i,a,b) for (int i = a; i >= b; i--)
using namespace std;
typedef long long LL;
const double EPS = 1e-10;
const int INF = 0x3f3f3f3f;
const int N = 2e6+10;
struct Edge{
int v,flow,next;
}edge[N];
int Head[900],cnt;
int cur[900],dep[900],S,T;
int n,m,d,ans;
char s[30];
int mp[30][30],mark[30][30],tt;
void Add_edge(int u, int v, int flow){
edge[++cnt] = {v,flow,Head[u]};
Head[u] = cnt;
edge[++cnt] = {u,0,Head[v]};
Head[v] = cnt;
}
bool check(int i, int j, int ii, int jj){
if(ii < 1 || jj < 1) return 0;
if ((i-ii)*(i-ii)+(j-jj)*(j-jj) <= d*d && mp[i][j] && mp[ii][jj]) return 1;
return 0;
}
void init(){
cnt = -1; mem(Head,-1);
scanf("%d%d%d",&n,&m,&d);
S = 0; T = 801;
go(i,1,n){
go(j,1,m){
scanf("%1d",&mp[i][j]);
mark[i][j] = ++tt;
}
}
go(i,1,n){
scanf("%s",s);
go(j,1,m){
if (s[j-1] == 'L')Add_edge(S,mark[i][j],1),ans++;
}
}
go(i,1,d){
go(j,d+1,n-d){
Add_edge(mark[j][i]+400,T,INF);
Add_edge(mark[j][m-i+1]+400,T,INF);
}
}
go(i,1,d){
go(j,1,m){
Add_edge(mark[i][j]+400,T,INF);
Add_edge(mark[n-i+1][j]+400,T,INF);
}
}
go(i,1,n)
go(j,1,m)
go(ii,i-d,i+d)
go(jj,j-d,j+d){
if (check(i,j,ii,jj) && (i != ii || j != jj)) Add_edge(mark[i][j]+400,mark[ii][jj],INF);
}
go(i,1,n)
go(j,1,m) if (mp[i][j])
Add_edge(mark[i][j],mark[i][j]+400,mp[i][j]);
}
bool bfs(){
queue<int>Q; while(!Q.empty()) Q.pop();
mem(dep,0); dep[S] = 1;
Q.push(S);
while(!Q.empty()){
int u = Q.front(); Q.pop();
for (int i = Head[u]; i != -1; i = edge[i].next){
int v = edge[i].v,flow = edge[i].flow;
if (dep[v] == 0 && flow > 0){
dep[v] = dep[u] + 1;
if (v == T) return 1;
Q.push(v);
}
}
}
if (dep[T]) return 1;
return 0;
}
int dfs(int u, int flow){
if (u == T || flow == 0) return flow;
for (int &i = cur[u]; i != -1; i = edge[i].next){
int v = edge[i].v;
if ((dep[v] > dep[u]) && edge[i].flow>0){
int d = dfs(v,min(flow,edge[i].flow));
if (d){
edge[i].flow -= d;
edge[i^1].flow += d;
return d;
}
}
}
return 0;
}
int Dinic(){
int d;
while(bfs()){
go(i,0,T) cur[i] = Head[i];
while(d = dfs(S,INF)) ans-=d;
}
return ans;
}
int main(){
init();
Dinic();
cout<<ans<<endl;
return 0;
}
/*
5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........
*/