首页 > 代码库 > AC日记——[SCOI2007]蜥蜴 bzoj 1066
AC日记——[SCOI2007]蜥蜴 bzoj 1066
1066
思路：
网络流最大流；
拆点，每个点拆成两个，流量为这个点的高度；
注意，文中说的距离是曼哈顿距离（劳资以为开根号wa了不知道多少次）；
每两个距离不大于d的点连边，流量inf；
如果距离能够延伸到边界外，就将这个点连向t；
最后输出，蜥蜴个数减去最大流；
来，上代码：
#include <cmath>#include <cstdio>#include <cstring>#include <iostream>#include <algorithm> using namespace std; #define maxn 905#define INF 0x7ffffff struct NodeType { int x,y,hi;};struct NodeType ai[maxn]; int n,m,cnt=1,tot,map[maxn][maxn],s,t,que[maxn*maxn*4],tott,d;int head[maxn],deep[maxn],F[maxn*maxn*2],E[maxn*maxn*2],V[maxn*maxn*2];inline void edge_add(int u,int v,int f){ E[++cnt]=head[u],V[cnt]=v,F[cnt]=f,head[u]=cnt; E[++cnt]=head[v],V[cnt]=u,F[cnt]=0,head[v]=cnt;} bool bfs(){ for(int i=s;i<=t;i++) deep[i]=-1; int h=0,tail=1;que[h]=s,deep[s]=0; while(h<tail) { int now=que[h++]; for(int i=head[now];i;i=E[i]) { if(deep[V[i]]<0&&F[i]>0) { deep[V[i]]=deep[now]+1; if(V[i]==t) return true; que[tail++]=V[i]; } } } return false;} int flowing(int now,int flow){ if(now==t||flow<=0) return flow; int oldflow=0; for(int i=head[now];i;i=E[i]) { if(deep[V[i]]==deep[now]+1&&F[i]>0) { int pos=flowing(V[i],min(flow,F[i])); F[i]-=pos,F[i^1]+=pos; flow-=pos,oldflow+=pos; if(flow==0) return oldflow; } } if(oldflow==0) deep[now]=-1; return oldflow;} int main(){ scanf("%d%d%d",&n,&m,&d); char ch[101]; for(int i=1;i<=n;i++) { scanf("%s",ch+1); for(int j=1;j<=m;j++) { if(ch[j]!=‘0‘) ai[++tot].x=i,ai[tot].y=j,ai[tot].hi=ch[j]-‘0‘,map[i][j]=tot; } } t=tot*2+1; for(int i=1;i<=tot;i++) { edge_add(i,i+tot,ai[i].hi); for(int j=1;j<=tot;j++) { if(j==i) continue; if(abs(ai[i].x-ai[j].x)+abs(ai[i].y-ai[j].y)<=d) edge_add(i+tot,j,INF); } if(ai[i].x<=d||ai[i].y<=d||n-ai[i].x+1<=d||m-ai[i].y+1<=d) edge_add(i+tot,t,INF); } for(int i=1;i<=n;i++) { scanf("%s",ch+1); for(int j=1;j<=m;j++) if(ch[j]==‘L‘) edge_add(s,map[i][j],1),tott++; } int ans=0; while(bfs()) ans+=flowing(s,INF); cout<<tott-ans; return 0;}
AC日记——[SCOI2007]蜥蜴 bzoj 1066
声明：以上内容来自用户投稿及互联网公开渠道收集整理发布，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任，若内容有误或涉及侵权可进行投诉： 投诉/举报 工作人员会在5个工作日内联系你，一经查实，本站将立刻删除涉嫌侵权内容。