BZOJ1066: [SCOI2007]蜥蜴

xiaoxiao2021-02-28  29

题目描述:

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 HINT 100%的数据满足:1<=r, c<=20, 1<=d<=4

题解:

此题是网络流的又一个模型,用DINIC求解。 因为每一个柱子被跳了一次都会变短一个长度,这就可以看做是流量。但是这样只是点权,而网络流需要的是边权,怎么办呢? 我们可以把一个点拆成两个点,可以理解为一个石柱的顶部和底部。链接它们的这条路的权值就是柱子高度。 然后就是常规套路了。建立超级源和超级汇,连接超级源和所有有蜥蜴的柱子,边权为1;链接所有能一步跳出地图的柱子和超级汇,边权为INF;连接所有可以一步到达的柱子(两个点的顶部和底部分别相连),边权为INF。可以理解为只有蜥蜴走的路线才有真实流量,限制了最大流。 ps:DINIC复杂度玄学。。。 代码如下:

#include<cstdio> #include<string> #include<cstring> using namespace std; const int maxn=1005,maxm=160005,inf=1<<30; int n,m,ans,d,S,T,tot,lnk[maxn],son[maxm],nxt[maxm],w[maxm],dep[maxm],cur[maxm],que[maxm]; inline int read(){ int x=0; char ch=getchar(); while (ch<'0'||ch>'9') ch=getchar(); while (ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar(); return x; } int getid(int x,int y,int z) {return (x-1)*m+y+z*n*m;} void add(int x,int y,int z) { son[++tot]=y,w[tot]=z,nxt[tot]=lnk[x],lnk[x]=tot; son[++tot]=x,w[tot]=0,nxt[tot]=lnk[y],lnk[y]=tot; } int bfs(){ memset(dep,0,sizeof(dep)); int head=0,tail=1; que[tail]=S,dep[S]=1; while (head!=tail) { head++; for (int j=lnk[que[head]];j!=-1;j=nxt[j]) if (dep[son[j]]==0&&w[j]>0) dep[son[j]]=dep[que[head]]+1,que[++tail]=son[j]; } return dep[T]; } int dfs(int x,int sum){ if (x==T) return sum; for (int& j=cur[x];j!=-1;j=nxt[j]) if (w[j]>0&&dep[son[j]]==dep[x]+1) { int sum1=dfs(son[j],min(sum,w[j])); if (sum1>0) {w[j]-=sum1,w[j^1]+=sum1; return sum1;} } return 0; } int dinic(){ while (bfs()) { for (int i=1;i<=T;i++) cur[i]=lnk[i]; ans-=dfs(S,inf); } return ans; } int main(){ n=read(),m=read(),d=read(); tot=-1; memset(lnk,255,sizeof(lnk)); memset(nxt,255,sizeof(nxt)); S=2*n*m+1,T=S+1; for (int i=1;i<=n;i++) for (int j=1;j<=m;j++){ char ch=getchar(); while (ch<'0'||ch>'3') ch=getchar(); if (ch>'0') add(getid(i,j,0),getid(i,j,1),ch-48); if (min(min(i,j),min(n-i+1,m-j+1))<=d) add(getid(i,j,1),T,inf); } for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) { char ch=getchar(); while (ch!='L'&&ch!='.') ch=getchar(); if (ch=='L') {ans++; add(S,getid(i,j,0),1);} } for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) for (int i1=1;i1<=n;i1++) for (int j1=1;j1<=m;j1++) if ((i!=i1||j!=j1)&&(i-i1)*(i-i1)+(j-j1)*(j-j1)<=d*d) add(getid(i,j,1),getid(i1,j1,0),inf); printf("%d\n",dinic()); return 0; }
转载请注明原文地址: https://www.6miu.com/read-2626132.html

最新回复(0)