zoukankan      html  css  js  c++  java
  • ZOJ2760_How Many Shortest Path

    给一个图,求从某个点到另一个点的最短路有多少条?所有的路都不共边。

    首先从终点开始Spfa标记最短距离,然后建图。

    建图的时候,如果满足两点之间的最短路只差为两点之间的边长,那么在网络流的模型中连接一条边。

    最终也只需要跑最大流即可。

    注意此题没有要求不能经过同一个点,所有不需要拆点,由于我们在网络流的模型中间加边的时候边容量为1,也就保证了每条边只遍历一边了。

    注意,有可能两个不同点之间的距离也为0,真是深坑啊,无法直视。

    召唤代码君:

    #include <iostream>
    #include <cstdio>
    #include <cstring>
    #define maxn 111
    #define maxm 44444
    typedef long long ll;
    using namespace std;
    
    ll inf=~0U>>2;
    ll to[maxm],next[maxm],c[maxm],first[maxn],edge;
    ll d[maxn],tag[maxn],TAG=520;
    ll Q[maxm],bot,top;
    ll ans,n,s,t;
    ll dis[maxn][maxn],f[maxn];
    bool can[maxn];
    
    void _input()
    {
        for (ll i=1; i<=n; i++)
        {
            f[i]=inf,first[i]=-1;
            for (ll j=1; j<=n; j++) scanf("%lld",&dis[i][j]);
        }
        scanf("%lld%lld",&s,&t);
        s++,t++;
        Q[bot=top=1]=t,f[t]=0;
        while (bot<=top)
        {
            ll cur=Q[bot++];
            for (ll i=1; i<=n; i++)
                if (dis[i][cur]>=0 && f[cur]+dis[i][cur]<f[i])
                    f[i]=f[cur]+dis[i][cur],Q[++top]=i;
        }
    }
    
    void addedge(ll U,ll V)
    {
        edge++;
        to[edge]=V,c[edge]=1,next[edge]=first[U],first[U]=edge;
        edge++;
        to[edge]=U,c[edge]=0,next[edge]=first[V],first[V]=edge;
    }
    
    void build_graph()
    {
        edge=-1;
        for (int i=1; i<=n; i++)
            for (int j=1; j<=n; j++)
                if (i!=j && dis[i][j]>=0 && f[i]==f[j]+dis[i][j])
                    addedge(i,j);
    }
    
    bool bfs()
    {
        Q[bot=top=1]=t,d[t]=0,tag[t]=++TAG,can[t]=false;
        while (bot<=top)
        {
            ll cur=Q[bot++];
            for (ll i=first[cur]; i!=-1; i=next[i])
                if (c[i^1]>0 && tag[to[i]]!=TAG)
                {
                    tag[to[i]]=TAG,can[to[i]]=false;
                    d[to[i]]=d[cur]+1,Q[++top]=to[i];
                    if (to[i]==s) return true;
                }
        }
        return false;
    }
    
    ll dfs(ll cur,ll num)
    {
        if (cur==t) return num;
        ll tmp=num,k;
        for (ll i=first[cur]; i!=-1; i=next[i])
            if (c[i]>0 && tag[to[i]]==TAG && d[to[i]]==d[cur]-1 && !can[to[i]])
            {
                k=dfs(to[i],min(num,c[i]));
                if (k) num-=k,c[i]-=k,c[i^1]+=k;
                if (num==0) break;
            }
        if (num) can[cur]=true;
        return tmp-num;
    }
    
    int main()
    {
        inf*=inf;
        while (scanf("%lld",&n)!=EOF)
        {
            _input();
            if (s==t)
            {
                puts("inf");
                continue;
            }
            build_graph();
            for (ans=0; bfs(); ) ans+=dfs(s,inf);
            printf("%lld
    ",ans);
        }
        return 0;
    }
    如有转载,请注明出处(http://www.cnblogs.com/lochan)
  • 相关阅读:
    存储过程访问外部服务器
    MVC4 code first 增加属性,对应自动修改列的方法笔记
    摘抄
    一条命令使win7可以直接运行.net3.5程序
    工作心得
    删除一个不存在的东西可以抛出异常吗
    洛谷P3379 【模板】最近公共祖先(LCA)
    洛谷P1967 货车运输
    洛谷P1653 猴子
    洛谷P2278 [HNOI2003]操作系统
  • 原文地址:https://www.cnblogs.com/lochan/p/3854865.html
Copyright © 2011-2022 走看看