hdu 6582path （2019 Multi-University Training Contest 1 1005）

题意：jerry要去见女朋友，但他的好朋友tom觉得jerry因此忽视了他，现在tom想要通过破坏某些路让jerry不能在最短时间见到jerry的女朋友，破坏每条路的花费就是这条路的长度(单向边)，求tom的最小花费

思路：最短路加最小割，我们先求出起点到所有点的最短距离和终点到所有点的最短距离，然后通过遍历所有的边找出最短路径的所有边，( disu[ edge[ i ].u ]+edge[ i ].w+disv[ edge[ i ].v ]==disu[ v ] )，disu[ i ](起点到i的最短距离) disv[ i ](终点到i的最短距离)

有这些边之后建一个新图跑一边最小割就好了

AC代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
#include<cmath>
#include<stack>
#include<map>
#include<set>
using namespace std;
#define LL long long
const int MOD=1e9+7;
const int inf=0x3f3f3f3f;
const LL inff=0x3f3f3f3f3f3f3f3f;
const LL N=205;
const LL M=5005;
const double pi=acos(-1);
#define MEF(x) memset(x,-1,sizeof(x))
#define ME0(x) memset(x,0,sizeof(x))
#define MEI(x) memset(x,inf,sizeof(x))
struct Edge
{
    int u,v,next;
    LL w;
}edge1[20005],edge2[20005],edge3[20005];
struct Node
{
    int i;
    LL dis;
    bool operator < (const Node &a)const
    {
        return dis>a.dis;
    }
};
int cnt,cnt2,cnt3,first[10005],first2[10005],vis[10005];
LL disu[10005],disv[10005],dis[10005],Max;
void init()
{
    cnt=cnt3=0;
    MEF(first);
    MEF(first2);
}
inline void add(int u,int v,int w)
{
    edge1[cnt].u=u;
    edge1[cnt].v=v;
    edge1[cnt].w=w;
    edge1[cnt].next=first[u];
    first[u]=cnt++;
}
inline void add2(int u,int v,int w)
{
    edge2[cnt2].u=u;
    edge2[cnt2].v=v;
    edge2[cnt2].w=w;
    edge2[cnt2].next=first[u];
    first[u]=cnt2++;
}
inline void add3(int u,int v,int w)
{
    edge3[cnt3].u=u;
    edge3[cnt3].v=v;
    edge3[cnt3].w=w;
    edge3[cnt3].next=first2[u];
    first2[u]=cnt3++;
}
void dijstra(int s,int t,LL di[],Edge edge[],int first1[])
{
    priority_queue<Node> pq;
    ME0(vis);
    di[s]=0;
    Node l;
    l.i=s;
    l.dis=0;
    pq.push(l);
    while(!pq.empty())
    {
        l=pq.top();
        pq.pop();
        int x=l.i;
        if(vis[x])
        {
            continue;
        }
        vis[x]=1;
        for(int i=first1[x];i!=-1;i=edge[i].next)
        {
            int lv=edge[i].v;
            if(di[lv]>di[x]+edge[i].w)
            {
                di[lv]=di[x]+edge[i].w;
                l.i=lv;
                l.dis=di[lv];
                pq.push(l);
            }
        }
    }
}
int n,m;
inline int bfs(int s,int t)//建层次图
{
    queue<int> q;
    MEF(dis);
    q.push(s);
    dis[s]=0;
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        if(x==t)
        {
            return 1;
        }
        for(int i=first[x];i!=-1;i=edge2[i].next)
        {
            int v=edge2[i].v;
            if(edge2[i].w&&dis[v]==-1)
            {
                dis[v]=dis[x]+1;
                q.push(v);
            }
        }
    }
    if(dis[t]==-1)
    {
        return 0;
    }
    else
    {
        return 1;
    }
}
inline LL Dfs(int s,LL flow)//flow 到大某一节点的最大流
{
    if(s==n)
    {
        return flow;
    }
    LL ans=0;
    for(int i=first[s];i!=-1;i=edge2[i].next)
    {
        int v=edge2[i].v;
        if(edge2[i].w&&dis[v]==dis[s]+1)
        {
            int f=Dfs(v,min(flow-ans,edge2[i].w));//(flow-ans) 某一节点的最大流-流过该节点的流量
            edge2[i].w-=f;
            edge2[i^1].w+=f;
            ans+=f;
            if(ans==flow)//某一节点流量达到最大流
            {
                return ans;
            }
        }
    }
    return ans;//某一节点的最终流量
}
LL dinc(int s,int t)
{
    LL flow=0;
    while(bfs(s,t))
    {
        flow+=Dfs(s,inf);
    }
    return flow;
}
int main()
{
    int t;
    scanf("%d",&t);
    for(int t1=1;t1<=t;t1++)
    {
        scanf("%d%d",&n,&m);
        init();
        for(int m1=1;m1<=m;m1++)
        {
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            if(a==b)
            {
                continue;
            }
            add(a,b,c);
            add3(b,a,c);
        }
        MEI(disu);
        MEI(disv);
        dijstra(1,n,disu,edge1,first);
        //cout<<'&'<<endl;
        dijstra(n,1,disv,edge3,first2);
        Max=disu[n];
        if(Max==inf)
        {
            printf("0\n");
        }
        else
        {
            cnt2=0;
            MEF(first);
            for(int i=0;i<cnt;i++)
            {
                int lu=edge1[i].u,lv=edge1[i].v;
                if(disu[lu]+disv[lv]+edge1[i].w==Max)
                {
                    add2(lu,lv,edge1[i].w);
                    add2(lv,lu,0);
                }
            }
            LL ANS=dinc(1,n);
            printf("%lld\n",ANS);
        }
    }
    return 0;
}