Path 【多校1 HDU 6582】【最短路+最小割】

题目链接

题目大意

有n个城市，m条路，有个人每天都沿着最短路径从城市1到城市n，现在让你堵上一些路使得这个人走的路比之前长，堵上一个路的代价是这个路的长度，输出最小代价。

解题思路

先计算出来最短路，找出所有在最短路上的路，然后再这些路上跑最小割。
找bug找到吐血（QAQ）

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<queue>
using namespace std;

const int maxn=1e5+5;
struct node
{
    int v,nex;
    long long w;
};
node e[maxn*8];
int a[maxn],b[maxn],f[maxn],tot,book[maxn],n,m,le[maxn];
long long c[maxn],dis[maxn];
const long long  inf=0x3f3f3f3f3f3f3f3f;
void init()
{
    memset(f,-1,sizeof(f));
    tot=0;
}
void add(int u,int v,long long w)
{
    e[tot].v=v;
    e[tot].w=w;
    e[tot].nex=f[u];
    f[u]=tot++;
}
void add1(int u,int v,long long w)
{
    e[tot].v=v;
    e[tot].w=w;
    e[tot].nex=f[u];
    f[u]=tot++;

    e[tot].v=u;
    e[tot].w=0;
    e[tot].nex=f[v];
    f[v]=tot++;
}
void spfa()
{
    for(int i=1;i<=n;i++)
        dis[i]=inf;
    memset(book,0,sizeof(book));
    dis[1]=0;
    book[1]=1;
    queue<int>q;
    q.push(1);
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        book[x]=0;
        for(int i=f[x];i!=-1;i=e[i].nex)
        {
            //printf("*");
            int v=e[i].v;
            if(dis[v]>dis[x]+e[i].w)
            {
                dis[v]=dis[x]+e[i].w;
                if(!book[v])
                {
                    book[v]=1;
                    q.push(v);
                }
            }
        }
    }
}

int bfs()
{
    queue<int>q;
    memset(le,0,sizeof(le));
    q.push(1);
    le[1]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        if(u==n)
            return 1;
        for(int i=f[u]; i!=-1; i=e[i].nex)
        {
            int v=e[i].v;
            if(!le[v]&&e[i].w>0)
            {
                le[v]=le[u]+1;
                q.push(v);
            }
        }
    }
    return 0;
}
long long dfs(int u,int l,long long maxf)
{
    if(u==n||maxf==0)
        return maxf;
    int ret=0,tmp;
    int v;
    for(int i=f[u]; i!=-1; i=e[i].nex)
    {
        v=e[i].v;
        if(le[v]==le[u]+1&&e[i].w>0&&(tmp=dfs(v,n,min(e[i].w,maxf-ret)))>0)
        {
            ret+=tmp;
            e[i].w-=tmp;
            e[i^1].w+=tmp;
            if(ret==maxf)
                return ret;
        }
    }
    le[u]=0;
    return ret;
}
long long dinic()
{
    int ans=0;
    while(bfs())
    {
        ans+=dfs(1,n,inf);
    }
    return ans;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        init();
        scanf("%d %d",&n,&m);
        for(int i=1;i<=m;i++)
        {
            scanf("%d %d %lld",&a[i],&b[i],&c[i]);
            add(a[i],b[i],c[i]);
        }

        spfa();
        init();
        for(int i=1;i<=m;i++)
        {
            int u=a[i],v=b[i];
            long long w=c[i];
            if(dis[v]==dis[u]+w)
            {
                //printf("%d %d\n",u,v);
                add1(u,v,w);
            }

        }
        printf("%lld\n",dinic());
    }
    return 0;
}