首页 > 代码库 > BZOJ 1266 上学路线route(最小割)
BZOJ 1266 上学路线route(最小割)
题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=1266
题意:给出一个无向图,每条边有长度和代价。求出1到n的最短路。之后删掉一些边使得1到n的最短路变大?在此情况下使得删掉边的代价之和最小。
思路:首先求出每个点到1和n的最短路。之后可以确定每条边是否为关键边(就是最短路上的边)。将关键边建立网络流图,求最小割即可。
struct node{ int v,cap,next;}; node edges[N];int head[N],e; void add(int u,int v,int cap){ edges[e].v=v; edges[e].cap=cap; edges[e].next=head[u]; head[u]=e++;} void Add(int u,int v,int cap){ add(u,v,cap); add(v,u,0);} int pre[N],cur[N],num[N],h[N]; int Maxflow(int s,int t,int n){ int i; for(i=0;i<=n;i++) cur[i]=head[i],num[i]=h[i]=0; int u=s,Min,k,v; int ans=0; while(h[u]<n) { if(u==t) { Min=INF; for(i=s;i!=t;i=edges[cur[i]].v) { k=cur[i]; if(edges[k].cap<Min) Min=edges[k].cap,v=i; } ans+=Min; u=v; for(i=s;i!=t;i=edges[cur[i]].v) { k=cur[i]; edges[k].cap-=Min; edges[k^1].cap+=Min; } } for(i=cur[u];i!=-1;i=edges[i].next) { if(edges[i].cap>0&&h[u]==h[edges[i].v]+1) break; } if(i!=-1) { cur[u]=i; pre[edges[i].v]=u; u=edges[i].v; } else { if(--num[h[u]]==0) break; k=n; cur[u]=head[u]; for(i=head[u];i!=-1;i=edges[i].next) { if(edges[i].cap>0&&h[edges[i].v]<k) { k=h[edges[i].v]; } } num[k+1]++; h[u]=k+1; if(u!=s) u=pre[u]; } } return ans;}struct Node{ int u,v,w,c,next;};Node edge[300000];int head1[505],e1;void add1(int u,int v,int w,int c){ edge[e1].u=u; edge[e1].v=v; edge[e1].w=w; edge[e1].c=c; edge[e1].next=head1[u]; head1[u]=e1++;}int n,m,dis1[505],dis2[505];int inq[505];void BFS(int dis[],int s){ int i; FOR1(i,n) dis[i]=INF,inq[i]=0; queue<int> Q; Q.push(s); dis[s]=0; int u,v; while(!Q.empty()) { u=Q.front(); Q.pop(); inq[u]=0; for(i=head1[u];i!=-1;i=edge[i].next) { v=edge[i].v; if(dis[v]>dis[u]+edge[i].w) { dis[v]=dis[u]+edge[i].w; if(!inq[v]) inq[v]=1,Q.push(v); } } }}void build(){ clr(head,-1); int i,u,v,c; FOR0(i,e1) { u=edge[i].u; v=edge[i].v; c=edge[i].c; if(dis1[u]+edge[i].w+dis2[v]==dis1[n]) { if(u!=n&&v!=1) Add(u,v,c); } }}int main(){ RD(n,m); clr(head1,-1); int i,u,v,w,c; FOR1(i,m) { RD(u,v); RD(w,c); add1(u,v,w,c); add1(v,u,w,c); } BFS(dis1,1); BFS(dis2,n); build(); PR(dis1[n]); PR(Maxflow(1,n,n));}
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。