zoukankan      html  css  js  c++  java
  • 差分约束

    众所周知,差分约束有m个不等式,要找到一组解(本题中是非负解)满足所有约束,因为不等式都是差分形式,又要满足所有约束,所以叫差分约束。

    我们采取这样的方式建边:对于ai - aj <= b,从j向i连一条边权为b的边

    原因:ai - aj <= b十分形似 dis[i] - dis[j] <= b,即dis[i] <= dis[j] + b,在一个有最短路的图中,对于每一个点均满足这个等式,所以我们想把ai转化为i的最短路,最后求出的最短路就是答案

    关于无解:当图中存在负环即无解

    一方面,如果有负环,原图中没有最短路,因此无解。从另一个角度考虑,如果有负环,则意味着永远不可能满足不等式(dis[i]和dis[j]都可以无穷小),所以无解(不懂的话可以画一个有两个点的图试一试

    细节:防止图不联通的情况,我们建一个超级源点S,由S向其他点连一条边权为0的边,图即联通

    spfa判负环就不说了

    #include<cstdio>
    #include<algorithm>
    #include<iostream>
    #include<queue>
    #include<cstring>
    using namespace std;
    inline int read()
    {
        int ans = 0,op = 1;
        char ch = getchar();
        while(ch < '0' || ch > '9')
        {
            if(ch == '-') op = -1;
            ch = getchar();
        }
        while(ch >= '0' && ch <= '9')
        {
            (ans *= 10) += ch - '0';
            ch = getchar();
        }
        return ans * op;
    }
    const int maxn = 2005;
    struct egde
    {
        int to,next,cost;
    }e[maxn * 10];
    int fir[maxn],alloc;
    inline void adde(int u,int v,int w)
    {
        e[++alloc].next = fir[u];
        fir[u] = alloc;
        e[alloc].to = v;
        e[alloc].cost = w;
    }
    int n,m;
    int dis[maxn],popst[maxn],minm;
    bool instack[maxn];
    void spfa(int s)
    {
        queue<int> q;
        memset(dis,0x3f,sizeof(dis));
        dis[s] = 0;
        q.push(s);
        instack[s] = 1;
        while(q.size())
        {
            int u = q.front();
            q.pop();
            popst[u]++;
            if(popst[u] > n - 1) { printf("NO SOLUTION"); return;}
            instack[u] = 0;
            for(int i = fir[u];i;i = e[i].next)
            {
                int v = e[i].to,w = e[i].cost;
                if(dis[v] > dis[u] + w)
                {
                    dis[v] = dis[u] + w;
                    if(!instack[v]) q.push(v),instack[v] = 1;
                }
            }
        }
        for(int i = 1;i <= n;i++) minm = min(minm,dis[i]);
        for(int i = 1;i <= n;i++) printf("%d
    ",dis[i] - minm);
    }
    int main()
    {
        n = read(),m = read();
        for(int i = 1;i <= m;i++)
        {
            int u = read(),v = read(),w = read();
            adde(v,u,w);    
        }
        for(int i = 1;i <= n;i++) adde(n + 1,i,0);
        spfa(n + 1);
    }
    View Code
  • 相关阅读:
    flask的简单使用
    Android Studio 使用Lambda表达式
    安卓SDK更新 国内镜像站点记录
    REACT-NATIVE入门(环境搭建及HELLOWORLD)
    深入理解Android 读感(一)
    JNI初入手
    (转)c++ 中的using namespace std是什么意思,什么时候用
    (转)const char to LPCTSTR不能转化问题
    Git常用操作之add操作
    Git日志操作
  • 原文地址:https://www.cnblogs.com/LM-LBG/p/10645034.html
Copyright © 2011-2022 走看看