poj3169 差分约束

【问题描述】

　　如果一个系统由n个变量和m个约束条件组成，其中每个约束条件形如xj-xi≤bk(i,j∈[1,n],k∈[1,m]),则称其为差分约束系统。亦即，差分约束系统是求解关于一组变量的特殊不等式组的方法。

　　下面是关于差分约束系统一个实际问题的描述：

　　FJ有 n 头奶牛（编号从 1 到 n），沿一条直线站着等候喂食，由于奶牛们身材都比较苗条，所以可能有多头奶牛站在同一个位置的情况。设第i头奶牛站队位置为d[i]，在安排每头奶牛站队位置d[i]时需要满足如下的条件：br>
　　1、奶牛应按照编号顺序站队，即d[i]-d[i+1]<=0。
　　2、m对友好的牛希望彼此站的近些。即对于友好关系(i,j,D)有d[j]-d[i]≤D（1 ≤ i < j ≤ n）。
　　3、p对敌对的牛希望彼此站的远些。即对于敌对关系(i,j,L)有d[j]-d[i]≥L（1 ≤ i < j ≤ n）。

　　现在请你帮助FJ判断是否存在满足这些条件的一种站队方案，如果存在请计算1号奶牛与n号奶牛之间的最大距离。如果不存在满足要求的方案，输出-1；如果1号奶牛和N号奶牛间的距离可以任意大，输出-2；

【输入格式】

　　第 1 行输入 n,m,p，接下来的 m 行，每行三个整数 i,j,D，表示牛 i 和牛 j 的距离最多为 D，再接下来的 p 行，每行三个整数 i,j,L，表示牛 i 和牛 j 的距离最少为 L。

【输出格式】

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

【输入样例】

4 2 1
1 3 10
2 4 20
2 3 3

【输出样例】

27

【样例解释】

四只奶牛的位置分别为：0,7,10,27

【数据范围】

2≤n≤1000 , 1≤m,p≤10000 , 1≤D,L≤1000000

实在想不到方法暴力求解可得一定分数。
也可以用动态规划：d[i]表示1到i的最长距离
但最好还是用图论求最优路径。
有序：d[i]<=d[i+1]
朋友：d[j]<=d[i]+D
敌人：d[i]<=d[j]-L
根据不等式连接各点
首先有n个点，于是把每条边向前一个点连一条长为0的边
互为朋友的奶牛从编号小的往大的连一条长为D的边
互为敌人的奶牛从编号大的往小的连一条长为-L的边
（把不等式转换为图的边）
找图上1到n的最短路径（因为每次查找敌人的时候都让他尽量与前面的远，所以1到n的最短路径就是最远距离）

#include<cstdio>
#include<map>
#include<vector>
#define maxn 10005
#define oo 100000000
using namespace std;
int n,m,p,x,y,z,front,rear;
int q[maxn*maxn],vis[maxn]={0},d[maxn];
bool inq[maxn]={false}; 
vector<int>g[maxn],w[maxn];
void in()
{
    scanf("%d%d%d",&n,&m,&p);

    for(int i=1;i<=n;i++)
    {
        g[i].push_back(i-1);
        w[i].push_back(0);
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        if(x>y) swap(x,y);
        g[x].push_back(y);
        w[x].push_back(z);
    }
    for(int i=1;i<=p;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        if(x>y) swap(x,y);
        g[y].push_back(x);
        w[y].push_back(-z);
    }

}

int spfa()
{
    for(int i=1;i<=n;i++) d[i]=oo;

    front=rear=1;
    q[rear++]=1;
    vis[1]++;
    inq[1]=1;
    d[1]=0;
    while(front!=rear)
    {
        int i=q[front++];
        inq[i]=0;
        for(int k=0;k<g[i].size();k++)
        {
            int j=g[i][k],l=w[i][k];
            if(d[i]+l>=d[j]) continue;
            d[j]=d[i]+l;
            if(inq[j]) continue;
            q[rear++]=j;
            inq[j]=1;
            vis[j]++;
            if(vis[j]>n) return -1;
        }
    }
    if(d[n]==oo) return -2;
    return d[n];
}

int main()
{
    //freopen("in.txt","r",stdin);
    in();
    if(g[1].size()==0)
    {
        printf("-2");
        return 0;
    }
    printf("%d",spfa());
    return 0;
}