洛谷2472 蜥蜴(網路流)
阿新 • • 發佈:2018-10-31
傳送門
【題目分析】
令人智熄的字串讀入操作。。。。。在BZOJ上過了然後氵谷全T?emmm。。。。。。
網路流的題難點就在於建圖,這道題還是比較明顯,首先每個點還是要拆點限制流量,上限設為石柱高度表示最多可以跳過這麼多蜥蜴。
建立一個起點s和終點t,每個有蜥蜴的石柱對應的入點就和起點s連線,s與入點間的弧容量為1(含義是該石柱只能有一個蜥蜴,如果容量大於1,就有可能在跑最大流時出現一個石柱有多個蜥蜴的情況),每個與邊界距離足夠近的石柱對應的出點和終點t連線,t與出點間的弧容量為INF,這樣t與出點間的弧容量不會影響到結果,此題的圖就建立完成了,再跑一遍最大流,最大流等於能逃脫的蜥蜴數量,輸出時用總蜥蜴數減去能逃脫的蜥蜴數即可。
【程式碼~】
#include<bits/stdc++.h> using namespace std; const int MAXN=1e4+10; const int MAXM=5e5+10; const int INF=0x3f3f3f3f; int n,m,cnt,s,t,mind,ans; int head[MAXN],cur[MAXN],depth[MAXN]; int nxt[MAXM],to[MAXM],w[MAXM]; int mx[MAXN],my[MAXN],tot; char x; int Read() { int i=0,f=1; char c; for(c=getchar();(c>'9'||c<'0')&&c!='-';c=getchar()); if(c=='-') f=-1,c=getchar(); for(;c>='0'&&c<='9';c=getchar()) i=(i<<3)+(i<<1)+c-'0'; return i*f; } void Add(int x,int y,int z) { nxt[cnt]=head[x]; head[x]=cnt; to[cnt]=y; w[cnt]=z; cnt++; } void add(int x,int y,int z) { Add(x,y,z); Add(y,x,0); } bool bfs() { queue<int> q; memset(depth,0,sizeof(depth)); depth[s]=1; q.push(s); while(!q.empty()) { int u=q.front(); q.pop(); for(int i=head[u];i!=-1;i=nxt[i]) { int v=to[i]; if(!depth[v]&&w[i]) { depth[v]=depth[u]+1; q.push(v); } } } return depth[t]!=0; } int dfs(int u,int dist) { if(u==t) return dist; for(int &i=cur[u];i!=-1;i=nxt[i]) { int v=to[i]; if(depth[v]==depth[u]+1&&w[i]) { int di=dfs(v,min(w[i],dist)); if(di>0) { w[i]-=di; w[i^1]+=di; return di; } } } return 0; } int dinic() { int ans=0; while(bfs()) { for(int i=s;i<=t;++i) cur[i]=head[i]; while(int d=dfs(s,INF)) ans+=d; } return ans; } int dist(int x1,int y1,int x2,int y2) { return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2); } int zh(int x,int y) { return (x-1)*m+y; } int main() { memset(head,-1,sizeof(head)); scanf("%d%d%d\n",&n,&m,&mind); s=0,t=n*m*2+1; for(int i=1;i<=n;++i) { for(int j=1;j<=m;++j) { cin>>x; if(x!='0') { add(zh(i,j),zh(i,j)+n*m,x-'0'); if(i<=mind||j<=mind||i+mind>n||j+mind>m) add(zh(i,j)+n*m,t,INF); tot++; mx[tot]=i,my[tot]=j; } } } for(int i=1;i<=n;++i) { for(int j=1;j<=m;++j) { cin>>x; if(x=='L') ans++,add(s,zh(i,j),1); } } for(int i=1;i<=tot;++i) for(int j=1;j<=tot;++j) if(i!=j) { if(dist(mx[i],my[i],mx[j],my[j])<=mind*mind) add(zh(mx[i],my[i])+n*m,zh(mx[j],my[j]),INF); } cout<<ans-dinic(); return 0; }