首页 > 代码库 > 【网络流】【BZOJ1006】【SCOI2007】蜥蜴
【网络流】【BZOJ1006】【SCOI2007】蜥蜴
学弟@lher在周末训练赛中出的题目的原题(这个人拿省选题来当作提高组模拟,太丧了。。。)
链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1066
题意简析:看题目:)
解题思路:题目显然是最大流。
首先拆点将点权变为边权,然后按照题意对于所有有跳板的点向可以跳到的点连一条权值为inf的边,对于能够跳出地图边界的点,将它与汇点连一条权值为inf的边。对于有蜥蜴的点,从源点向这个点连一条权值为1的边,然后跑一次最大流,答案就是蜥蜴数-最大流。
AC代码:
#include<stdio.h> #include<iostream> #include<cmath> #include<cstring> #include<string> #include<algorithm> #define ll long long #define inf 0x7fffffff using namespace std; struct zxy{ int next,to,v; }edge[1000005]; int n,m,d,head[10005],iter[10005],lev[10005],que[100005],cnt=1,mz=0; bool b[55][55]; inline int getno(int x,int y){ return (x-1)*m+y;} inline int getdis(int x,int y,int a,int b){ return ceil(sqrt((x-a)*(x-a)+(y-b)*(y-b)));} inline void ins(int x,int y,int l){edge[++cnt].next=head[x],head[x]=cnt,edge[cnt].to=y,edge[cnt].v=l;} inline int in(){ int x=0,f=1; char ch=getchar(); while(ch<‘0‘||ch>‘9‘) {if(ch==‘-‘) f=-1; ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘) x=x*10+ch-‘0‘,ch=getchar(); return x*f; } inline bool bfs(int s,int e){ memset(lev,-1,sizeof(lev)); int t=1,h=0; que[1]=s; lev[s]=0; do{ h++; int k=head[que[h]]; while(k){ if (lev[edge[k].to]==-1&&edge[k].v){ lev[edge[k].to]=lev[que[h]]+1; que[++t]=edge[k].to; } k=edge[k].next; } }while(h<t); return lev[e]!=-1; } inline int dfs(int u,int v,int f){ if (u==v) return f; int used=0,k=head[u]; while(k){ if (edge[k].v&&lev[edge[k].to]==lev[u]+1){ int w=dfs(edge[k].to,v,min(edge[k].v,f-used)); used+=w; edge[k].v-=w; edge[k^1].v+=w; if (used==f) return used; } k=edge[k].next; } return used; } int dinic(int s,int t){ int flow=0; while(bfs(s,t)) flow+=dfs(s,t,inf); return flow; } void init(){ n=in(),m=in(),d=in(); for (int i=1; i<=n; ++i){ for (int j=1; j<=m; ++j){ int x=getchar()-‘0‘,no=getno(i,j); if (x) ins(no,no+n*m,x),ins(no+n*m,no,0),b[i][j]=1; } getchar(); } for (register int i=1; i<=n; ++i){ for(register int j=1; j<=m; ++j) if (getchar()==‘L‘) ++mz,ins(0,getno(i,j),1),ins(getno(i,j),0,0); getchar(); } for (register int i=1; i<=n; ++i) for (register int j=1; j<=m; ++j) if(b[i][j]){ if (i-d<1||i+d>n||j-d<1||j+d>m) ins(getno(i,j)+n*m,n*m*2+1,inf),ins(n*m*2+1,getno(i,j)+n*m,0); for (int k=(i-d<1?1:i-d); k<=(i+d>n?n:i+d); ++k) for (int t=(j-d<1?1:j-d); t<=(j+d>m?m:j+d); ++t) if (b[k][t]&&getdis(i,j,k,t)<=d&&(i!=k||j!=t)) ins(getno(i,j)+n*m,getno(k,t),inf),ins(getno(k,t),getno(i,j)+n*m,0); } } int main(){ init(); printf("%d",mz-dinic(0,2*n*m+1)); return 0; }
【网络流】【BZOJ1006】【SCOI2007】蜥蜴
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。