POJ 1201 Intervals（差分约束）

我做的第一道差分约束题，把求满足不等式的解的个数最值问题转化成求最短路真是神转化，看半天才看明白这到底是什么意思。这题中给的不等式是1 <= ci <= bi - ai+1，问要满足这些条件至少需要多少个数。通过不等式我们能得到一个方程，令s(i)表示从最左边的点到第i个点满足约束条件的情况下至少需要多少个点，因为要保证所有子区间也要满足条件，所以应取最大值，这样就转化成了一个SPFA求最长路的问题。SPFA最长路的松弛条件为dis[j]>=dis[i]+v，那么我们就把题目中给出的不等式条件也转化成三角不等式的形式，即为bi+1>=ai+ci，那么类比过来就是dis[bi+1]>=dis[ai]+map[ai][bi+1]，那么我们就只需要在ai和bi+1这两个顶点之间连一条长度为ci的边，再加上题目中的隐含条件1>=s[i+1]-s[i]>=0，这样建图之后，然后把dis数组初始化为负无穷，最左边的点距离初始化为0，跑一次SPFA，然后输出到达最右边的点的最长路径即为要求解。

#include <iostream>
#include <string.h>
#include <cstdio>
#include <vector>
#include <queue>
using namespace std;
#define MAXN 50000
struct Node
{
    int v,len;
};
vector<Node> g[MAXN+5];
int dis[MAXN+5];
bool vis[MAXN+5];
const int INF=10000000;
void spfa(int s,int t)     //SPFA最长路
{
    memset(vis,0,sizeof(vis));
    for(int i=s; i<=t; i++)
        dis[i]=-INF;
    dis[s]=0;
    queue<int> q;
    q.push(s);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        vis[u]=false;
        for(int i=0; i<g[u].size(); i++)
        {
            int j=g[u][i].v;
            if(dis[j]<dis[u]+g[u][i].len)
            {
                dis[j]=dis[u]+g[u][i].len;
                if(!vis[j])
                {
                    vis[j]=true;
                    q.push(j);
                }
            }
        }
    }
}
int main()
{
    //freopen("input.txt","r",stdin);
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<=MAXN;i++)
            g[i].clear();
        int start=100000,end=-1;
        for(int i=0; i<n; i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            start=min(start,x);        //维护最左边的点的坐标
            end=max(end,y+1);     //维护做右边的点的坐标
            Node temp;
            scanf("%d",&temp.len);
            temp.v=y+1;
            g[x].push_back(temp);           //从ai向bi+1连一条长度为ci的边
        }
        for(int i=start; i<=end; i++)
        {
            Node temp;
            temp.v=i;
            temp.len=-1;
            g[i+1].push_back(temp);     //从i+1向i连一条长度为-1的边（1>=s[i+1]-s[i] => s[i]>=s[i+1]-1）
            temp.v=i+1;
            temp.len=0;
            g[i].push_back(temp);         //从i向i+1连一条长度为0的边（s[i+1]-s[i]>=0  => s[i+1]>=s[i]+0）
        }
        spfa(start,end);
        printf("%d\n",dis[end]);
    }
}