poj 1275(约束差分求解：汇总)

有是一道差分约束题，虽然啊以前看过，但不是很了解，现在明白了。

转：对于差分约束系统的一个资料参考：

1) 建图后求最短路（对应 <= 的差分约束）

从求最短路后分析开始。求最短路后一定有：d(u) + w(u,v) >= d(v) 转换为：d(v) – d(u) <= w(u,v)。然后这个式子与我们差分约束条件得到的不等式（ X(a) - X(b) <= c ）相似。

所以得到建图的规则为：一条b到a的边权值为c。

建图后求最短路即为一组解。

2) 建图后求最长路（对应 >= 的差分约束）

同上。求最长路后一定得到：d(u) + w(u,v) <= d(v) 转换为：d(v)-d(u)>=w(u,v)。然后对比差分约束条件的不等式（ X(a) - X(b) >= c ）。

所以得到建图的规则为：一条b到a的边权值为c。

建图后求最长路即为一组解。

如果用SPFA做的话，必须加一个超级源点，令它到每个点距离<=0就好，然后该怎么做就怎么做，存在负环表示无解，否则存在正解。

如果不加源点的话可以用Bellman_ford()做，否则有些不在统一图上的点是搜不到的

此题正解：http://blog.csdn.net/morgan_xww/article/details/6370743

Source Code

Problem: 1275  User: 1013101127 
Memory: 392K  Time: 0MS 
Language: G++  Result: Accepted 

Source Code 
#include <stdio.h>   
#include <stdlib.h>   
#include <string.h>   
  
#define MAXV 30   
#define MAXE 1000   
#define INF  100000000   
  
struct  
{  
    int  s, e, v;  
} edge[MAXE];  
  
int neg;  
int dis[MAXV];  
int r[30], t[30], sum;  
  
void add(int s, int e, int v)  
{  
    edge[neg].s = s;  
    edge[neg].e = e;  
    edge[neg].v = v;  
    neg++;  
}  
  
void change(int s)  
{  
    int i, j;  
    neg = 48;  
    for (j=1; j<=24; j++)  
    {  
        i = (j+7)%24 + 1;  
        if (i > j)  
            add(i, j, -r[i]);  
        else  
            add(i, j, s-r[i]);  
    }  
    add(24, 0, -s);  
}  
  
int relax(int s, int e, int v)  
{  
    if (dis[s]+v < dis[e])  
    {  
        dis[e] = dis[s]+v;  
        return 1;  
    }  
    return 0;  
}  
  
int BellmanFord()  
{  
    int i, j;  
    for (i=0; i<=24; i++)  
        dis[i] = INF;  
    dis[0] = 0;  
    for (i=0; i<=24; i++)  
        for (j=0; j<neg; j++)  
            relax(edge[j].s, edge[j].e, edge[j].v);  
    for (j=0; j<neg; j++)  
        if (relax(edge[j].s, edge[j].e, edge[j].v))  
            return 0;  
    return 1;  
}  
  
void Bsearch(int low, int high)  
{  
    int mid = (low+high)/2;  
  
    if (low > high) return ;  
    change(mid);  
    if (BellmanFord())  
    {  
        sum = mid;  
        Bsearch(low, mid-1);  
    }  
    else  
        Bsearch(mid+1, high);  
}  
  
int main()  
{  
    int  nc, n, k, i;  
  
    scanf("%d", &nc);  
    while (nc--)  
    {  
        neg = 0;  
        memset(t, 0, sizeof(t));  
        for (i=1; i<=24; i++)  
            scanf("%d", &r[i]);  
        scanf("%d", &n);  
        for (i=0; i<n; i++)  
        {  
            scanf("%d", &k);  
            t[k+1]++;  
        }  
        for (i=1; i<=24; i++)  
        {  
            add(i-1, i, t[i]);  
            add(i, i-1, 0);  
        }  
        sum = -1;  
        Bsearch(0, n);  
        if (sum == -1)  
            printf("No Solution\n");  
        else  
            printf("%d\n", sum);  
    }  
}