HDU(3592)---POJ (3169) Layout 差分约束系统

题目链接（POJ）

HDU

题意简述：

当排队等候喂食时，奶牛喜欢和它们的朋友站得靠近些。FJ有N（2<=N<=1000）头奶牛，编号从1到N，沿一条直线站着等候喂食。奶牛排在队伍中的顺序和它们的编号是相同的。因为奶牛相当苗条，所以可能有两头或者更多奶牛站在同一位置上。即使说，如果我们想象奶牛是站在一条数轴上的话，允许有两头或更多奶牛拥有相同的横坐标。

一些奶牛相互间存有好感，它们希望两者之间的距离不超过一个给定的数L。另一方面，一些奶牛相互间非常反感，它们希望两者间的距离不小于一个给定的数D。给出ML条关于两头奶牛间有好感的描述，再给出MD条关于两头奶牛间存有反感的描述。（1<=ML,MD<=10000，1<=L,D<=1000000）

你的工作是：如果不存在满足要求的方案，输出-1；如果1号奶牛和N号

奶牛间的距离可以任意大，输出-2；否则，计算出在满足所有要求的情况下，1号奶牛和N号奶牛间可能的最大距离。

分析：

如果当前问题比较复杂，我们应该学会“退一步”思考，由简单到复杂。

求最大值不知从何下手，我们先从容易的开始分析。

我们先研究，如果不要求输出1和N的最大距离，而只需一个可行的距离，应该如何操作。

我们用D[i]表示I号奶牛和1号奶牛间的距离。

因为在队伍中的顺序必须和编号相同，所以对于任意I号奶牛，1<=I<N，在距离上应该满足：

D[I+1] - D[I] >= 0

对于每个好感的描述(i，j，k)，假设i<=j，体现到距离上的要求就是：

D[j] - D[I] <= k

对于每个反感的描述(i，j，k)，假设i<=j，体现到距离上的要求就是：

D[j] - D[I] >= k

这时的模型有一个名称，叫作：差分约束系统。

为了方便起见，我们将每种不等式写成我们约定的形式：

D[I] <= D[I+1]

D[j] <= D[I] + k

D[I] <= D[j] - k

看到这些不等式，你想到了什么？

没错，在求顶点间地最短路问题中，我们有这样的不等式：

若顶点u到顶点v有边e=uv，且边权为w(e)，设d(u),d(v)为源点到顶点u和顶点v的最短路长，

则 d(v) <= d(u) + w(e)

这个不等式和前面的条件形式十分相似，这就启发我们用构图用最短路做。（想想为什么求最长距离，正好是最短路？）

具体步骤是：

作有向图G=(V,E)，V={ v1，v2，v3，…，vn}，E={e1，e2，e3，…}，对于相邻两点i和(i+1)，对应的顶点vi+1向vi引一条边，费用为0；对于每组好感描述(ai，bi，di)，我们假设有ai<bi，否则ai和bi交换，则顶点vai向vbi引一条边，费用为di；对于每组反感描述(ai，bi，di)，我们假设有ai<bi，否则ai和bi交换，则顶点vbi向vai引一条边，费用为-di。

于是问题变为在G中求v1到其它所有顶点的最短路。我们证明若G中无负权回路，则问题有解，即存在满足条件的数列，若G中有负权回路，则问题无解，即不存在满足条件的数列。

定理5 问题是否有解等价于图G是否没有负权回路。

证明：若G中无负权回路，我们可以求出v1其他顶点u的最短路长，设为d(u)。由于是最短路，因此对于任意边eE，e=uv，有d(u)+w(e)>=d(v)，从而所有的约束条件都被满足，问题一定有解。若G中有负权回路，说明在任何时刻，G中至少有一个点v的最短路长可以更新，因此必须存在一条边e=uv，使得d(u)+w(e)<d(v)。所以无论何时，都会有某个约束条件不被满足，问题无解。（证毕）

检测负权回路，可以用Bellman-Ford算法。

回到原问题。

先说说考试时我的做法。

将所有点的最短路估计值设为一个充分大的值，v1的最短路估计值设为0。然后运行一次Bellman-Ford。

如果图G中有负权回路，那么输出-1；

否则，如果标号为N的顶点的最短路估计仍为一个充分大的值，那么它和标号为1的顶点间的距离可以任意大，这时输出-2；

如果以上两种情况都不满足，那么输出标号为N的点的最短路径估计值。

什么是充分大呢？

只需要比原图中最大的边权×顶点数还大就行了。

因为只要无圈，每个顶点最多只会被经过一次，所以肯定比这个值小，所以在该图中，我们可以把这个值看作是无穷大。

该方法可以通过竞赛的所有测试数据。

考试的时候，我没有去刻意证明这个方法的正确性，只是感觉它应该可行。

现在让我们从理论上证明它。

定理6 若运行Bellman-Ford后，标号为N的顶点的最短路估计值仍为充分大，那么N和1的距离可以任意大。

证明：从刚才的操作可以看出，到了这一步，已经把含有负权回路的情况排除掉了。在图G中，该充分大的值比可能得到的最大距离大，因此，它和任意大的值对于G的效果都是一样的（同样大于合法的最大距离）。由于充分大的值在G中满足约束，所以任意大的值亦满足约束，从而距离可以任意大。（证毕）

剩下还有一个问题。

定理7 若运行Bellman-Ford后，标号为N的顶点的最短路估计值比充分大小，那么它是N和1可能的最大距离。

证明：设D[i]是顶点i和1的最短路径估计值，d[i]是顶点i和1可能的最大距离。

我们首先证明，d[n]<=D[n]，运用反证法。

假如d[n]>D[n]，那么在Bellman-Ford运行之前，将赋予每个顶点i的充分大的值换成对应的d[i]。由于d本身满足所有约束条件，所以运行后，得出D'=d。由于充分大的值比所有d[i]都大，而求最短路运用的是逐步松弛操作，我们设立一个更大的初值不可能导致我们的终值反而更小。所以对于任意i，必定有D[i]>=D'[i]，即有D[n]>=d[n]，这与我们的假设矛盾。

然后我们证明，d[n]>=D[n]。

根据d[i]的定义，它是i和1的可能最大距离。由于D[i]是满足题目的所有约束的，所以D[i]是顶点i和1可能的距离。如果D[i]>d[i]，那么与d[i]的定义矛盾。

综合上述，有D[i]=d[i]。从而D[n]是n和1可能的最大距离。（证毕）

运行Bellman-Ford最坏情况下的复杂度是O((ML+MD)*N)=O(2*10 7 )，可以在规定时间内出解。另外，实际运行的速度相当理想，大部分数据运行基本上不需要时间。

#include<iostream>
#include<queue>
#include<cstring>
using namespace std;
const int maxn=2000;
const int maxm=40000;
const int inf=0xffffff;
int p[maxn],eid,n,ml,md;
struct edge{
    int v,next,w;
}e[maxm];
void init(){
    memset(p,-1,sizeof(p));
    eid=0;
}
void insert(int a,int b,int w)
{
    e[eid].v=b;
    e[eid].next=p[a];
    e[eid].w=w;
    p[a]=eid++;
}
bool intq[maxn];
int dist[maxn],cnt[maxn];
bool spfa(int v){
    for(int i=1;i<=n;i++)dist[i]=inf;
    memset(intq,0,sizeof(intq));
    memset(cnt,0,sizeof(cnt));
    dist[v]=0;
    intq[v]=true;
    queue<int>q;
    q.push(v);
    while(!q.empty()){
        int j=q.front();
        if(cnt[j]>n)return false;
        q.pop();
        intq[j]=false;
        for(int i=p[j];i!=-1;i=e[i].next){
            int x=e[i].v;
            if(dist[x]>e[i].w+dist[j]){
                dist[x]=e[i].w+dist[j];
                if(!intq[x]){
                    q.push(x);
                    intq[x]=true;
                    cnt[x]++;
                }
            }
        }
    }
    return true;
}
int main(){
    int t;
    cin>>t;
    while(t--){
    init();
    cin>>n>>ml>>md;
    for(int i=0;i<ml;i++){
        int a,b,c;
        cin>>a>>b>>c;
        insert(a,b,c);
    }
    for(int i=0;i<md;i++){
        int a,b,c;
        cin>>a>>b>>c;
        insert(b,a,-c);
    }
    for(int i=1;i<=n;i++)insert(i,i-1,0);
    for(int i=1;i<=n;i++)insert(0,i,0);//插入一个超级源！！！。
    if(!spfa(1))cout<<"-1\n";
    else if(dist[n]==inf)cout<<"-2\n";
    else cout<<dist[n]<<endl;
    }
}