hdu 5294 最短路+最大流

来源:互联网 发布:iface102下载数据 编辑:程序博客网 时间:2024/06/05 15:17

点击打开链接

题意:题目要求1走到n的路程为最短路的路径,但是路径可能有多个,都要用,如果不是最短路径上的边就不走它了,第一个结果是问最少删除几条边1走不到n,第二个是最多删除几条边后,1仍能够走到n

思路:我们要先求出最短路,并将路径上的边挑出来,然后建图,此时的图流量和时间都是1,因为最后的时候我们只要的是边的个数,这样求最小割=最大流就是第一个答案,然后求最短路ans[n],用总变数m减去就是第二个答案

#include <queue>#include <vector>#include <stdio.h>#include <string.h>#include <stdlib.h>#include <iostream>#include <algorithm>#include <functional>using namespace std;typedef long long ll;const int inf=0x3f3f3f3f;const int maxn=2010;struct edge{    int to,cap,rev;    edge(){}    edge(int a,int b,int c){to=a;cap=b;rev=c;}};vector<edge> G[maxn];int level[maxn],iter[maxn];void add_edge(int from,int to,int cap){    G[from].push_back(edge(to,cap,G[to].size()));    G[to].push_back(edge(from,0,G[from].size()-1));}void bfs(int s){    memset(level,-1,sizeof(level));    queue<int>que;    level[s]=0;que.push(s);    while(!que.empty()){        int v=que.front();que.pop();        for(unsigned int i=0;i<G[v].size();i++){            edge &e=G[v][i];            if(e.cap>0&&level[e.to]<0){                level[e.to]=level[v]+1;                que.push(e.to);            }        }    }}int dfs(int v,int t,int f){    if(v==t) return f;    for(int &i=iter[v];i<G[v].size();i++){        edge &e=G[v][i];        if(e.cap>0&&level[v]<level[e.to]){            int d=dfs(e.to,t,min(f,e.cap));            if(d>0){                e.cap-=d;                G[e.to][e.rev].cap+=d;                return d;            }        }    }    return 0;}int max_flow(int s,int t){    int flow=0;    while(1){        bfs(s);        if(level[t]<0) return flow;        memset(iter,0,sizeof(iter));        int f;        while((f=dfs(s,t,inf))>0) flow+=f;    }}struct edge1{    int to,cost;    edge1(int a,int b){to=a;cost=b;}};typedef pair<int,int>P;vector<edge1>G1[maxn];int dis[maxn];void dijkstra(int s){    priority_queue<P,vector<P>,greater<P> >que;    fill(dis,dis+maxn,inf);    dis[s]=0;que.push(P(0,s));    while(!que.empty()){        P p=que.top();que.pop();        int v=p.second;        if(dis[v]<p.first) continue;        for(unsigned int i=0;i<G1[v].size();i++){            edge1 e=G1[v][i];            if(dis[e.to]>dis[v]+e.cost){                dis[e.to]=dis[v]+e.cost;                que.push(P(dis[e.to],e.to));            }        }    }}int A[60010],B[60010],C[60010];int main(){    int n,m;    while(scanf("%d%d",&n,&m)!=-1){        for(int i=0;i<maxn;i++) G[i].clear();        for(int i=0;i<maxn;i++) G1[i].clear();        add_edge(0,1,inf);        add_edge(n,n+1,inf);        for(int i=0;i<m;i++){            scanf("%d%d%d",&A[i],&B[i],&C[i]);            G1[A[i]].push_back(edge1(B[i],C[i]));            G1[B[i]].push_back(edge1(A[i],C[i]));        }        dijkstra(1);        for(int i=0;i<maxn;i++) G1[i].clear();        for(int i=0;i<m;i++){            if(dis[A[i]]>dis[B[i]]) swap(A[i],B[i]);            if(dis[B[i]]-dis[A[i]]==C[i]){                G1[A[i]].push_back(edge1(B[i],1));                G1[B[i]].push_back(edge1(A[i],1));                add_edge(A[i],B[i],1);            }        }        int ans=max_flow(0,n+1);        dijkstra(1);        int ans1=m-dis[n];        printf("%d %d\n",ans,ans1);    }    return 0;}

0 0
原创粉丝点击