差分约束系统

差分约束系统

 

一：基础知识

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

        求解差分约束系统，可以转化成图论的 单源最短路径（或最长路径）问题。

        观察xj-xi<=bk，会发现它类似最短路中的三角不等式d[v]<=d[u]+w[u,v]，即d[v]-d[u]<=w[u,v]。因此，以每个变量xi为结点，对于约束条件xj-xi<=bk，连接一条边(i,j)，边权为bk。我们再增加一个源点s,s与所有定点相连，边权均为0。对这个图，以s为源点运行Bellman-ford算法（或SPFA算法），最终{d[ i]}即为一组可行解。

                          

         例如，考虑这样一个问题，寻找一个5维 向量x=(xi)以满足：

         这一问题等价于找出未知量xi，i=1,2,…,5，满足下列8个差分约束条件： 

  

x1-x2≤0

                       x1-x5≤-1

                       x2-x5≤1

                       x3-x1≤5

                       x4-x1≤4

                       x4-x3≤-1

                       x5-x3≤-3

                       x5-x4≤-3

        该问题的一个解为x=(-5,-3,0,-1,-4)，另一个解y=(0,2,5,4,1)，这2个解是有联系的：y中的每个元素比x中相应的元素大5。

        引理：设x=(x1,x2,…,xn)是差分约束系统Ax≤b的一个解，d为任意常数。则x+d=(x1+d,x2+d,…,xn+d)也是该系统Ax≤b的一个解。

        bellman-ford算法 伪代码：

[cpp]  view plain copy

1. for each v V do d[v] <-- 无限大; d[s] <-- 0  
2. Relaxation  
3. for i =1,...,|V|-1 do  
4. for each edge (u,v) 属于 E do  
5. d[v] <-- min{d[v], d[u]+w(u,v)}  
6. Negative cycle checking  
7. for each v 属于V do if d[v]> d[u] + w(u,v) then no solution  
8. <span style="font-family:Arial;BACKGROUND-COLOR: #ffffff"></span>  

 

 

      在实际的应用中，一般使用 SPFA(Shortest Path Fast Algorithm)算法来实现。

      差分约束系统中源点到每个点的距离确定

      关于Dist[]的初始化化

            1.如果将源点到各点的距离初始化为0，最终求出的最短路满足 它们之间相互最接近了

            2.如果将源点到各点的距离初始化为INF(无穷大)，其中之1为0，最终求出的最短路满足 它们与该点之间相互差值最大。

            3.差分约束系统的确立要根据自己确定的约束条件，从约束点走向被约束点

     连边一般有两种方法，第一种是连边后求最长路的方法，第二种是连边后求最短路的方法。

     例：d[x]-d[y]>=Z

           如果想连边后求最长路 那么将不等式变形为这种形式 d[x]>=d[y]+z y---x连一条权值为z的边

           求最短路则变形成d[y]<=d[x]-z x---y连一条权值为-z的边。

           如果是别的不等式，也可以根据情况变形。但是要保证的是 两个变量（x,y)的系数一定要是正的。而 常量则不一定。

      参考资料：百度百科

      笔者经验：差分约束系统就是给定一组不等式组，然后对于每个不等式<x1-x2≤k>中x2到x1连一条有向边，权值为k，然后建立一个超级汇点<定为0>，那么0号节点到其他点的距离为0，注意负权边的spfa，若某一点入队超过N<即所有点数目>，则无解

 

二：例题<POJ1364>

King

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 9040		Accepted: 3395

Description

对题目中给定的si,ni,ki,和一个给定的序列S[1....N]，如果格式为(si,ni,gt,ki),意思就是新增一约束条件S[si]+S[si+1]+...S[si+ni]>ki，如果格式为(si,ni,lt,ki),意思就是新增一约束条件S[si]+S[si+1]+...S[si+ni]<ki，判断所给的约束条件有符合要求的序列S[1....N]，有解就输出lamentable kingdom，无解就输出successful conspiracy。

Sample Input

4 2
1 2 gt 0
2 2 lt 2
1 2
1 0 gt 0
1 0 lt 0
0

Sample Output

lamentable kingdom
successful conspiracy

【分析】：

    设a[i]表示∑a[k]，k∈[1,i]且k∈Z*，那么对于约束条件S[si]+S[si+1]+...S[si+ni]>ki可转变为a[si+ni]-a[si-1]>ki(因为要包含[si,si+ni]，因此根据a数组的定义，下限为a[si-1])，而差分约束系统只能解形如Ax≤b的不等式组，因此转化给出的约束条件S[si]+S[si+1]+...S[si+ni]>ki→a[si+ni]-a[si-1]>ki→a[si-1]-a[si+ni]<=-ki-1，另一个同理：a[si+ni]-a[si-1]<=ki-1，然后跑一边差分约束即可<这里用的是spfa+邻接表存储>

【基础理论】：差分约束系统，用最短路算法实现

【代码】：

[cpp]  view plain copy

1. /* 
2.    Problem: 1364  User: csyzcyj  
3.    Memory: 192K  Time: 0MS  
4.    Language: C++  Result: Accepted  
5. */  
6. #include<stdio.h>  
7. #include<string.h>  
8. #include<stdlib.h>  
9. #include<algorithm>  
10. #include<iostream>  
11. #include<queue>  
12. using namespace std;  
13. #define MAX 201  
14. #define MAXP 1001  
15. #define IMAX 21474836  
16. struct EDGE{int to,next,v;};  
17. EDGE a[MAXP];  
18. int N,M,tot=0,last[MAX],dist[MAX],tot1[MAX];       
19. bool vis[MAX];  
20. void add(int x,int y,int value)  
21. {  
22.       a[++tot].to=y;  
23.       a[tot].next=last[x];  
24.       a[tot].v=value;  
25.       last[x]=tot;  
26. }  
27. bool check()  
28. {  
29.       queue<int> Q;  
30.       memset(dist,0,sizeof(dist));  
31.       for(int i=0;i<=N;i++)  
32.             Q.push(i);  
33.       while(!Q.empty())  
34.       {  
35.             int now=Q.front();  
36.             Q.pop();  
37.             vis[now]=false;  
38.             for(int i=last[now];i!=-1;i=a[i].next)  
39.             {  
40.                   int now_to=a[i].to;  
41.                   if(dist[now]+a[i].v<dist[now_to])  
42.                   {  
43.                         dist[now_to]=dist[now]+a[i].v;  
44.                         if(!vis[now_to])  
45.                         {  
46.                               vis[now_to]=true;  
47.                               tot1[now_to]++;  
48.                               if(tot1[now_to]>N)   return false;     
49.                               Q.push(now_to);  
50.                         }  
51.                   }   
52.             }  
53.       }  
54.       return true;  
55. }  
56. int main()  
57. {  
58.       //freopen("input.in","r",stdin);  
59.       //freopen("output.out","w",stdout);   
60.       while(scanf("%d%d",&N,&M)!=EOF && N!=0 && M!=0)  
61.       {  
62.               memset(last,-1,sizeof(last));        
63.               memset(vis,true,sizeof(vis));  
64.               memset(tot1,0,sizeof(tot1));  
65.               memset(a,0,sizeof(a));  
66.               tot=0;  
67.               for(int i=1;i<=M;i++)  
68.               {  
69.                     int A,B,C;  
70.                     char ch1,ch2;  
71.                     scanf("%d%d",&A,&B);  
72.                     scanf("%c%c%c",&ch1,&ch1,&ch2);  
73.                     scanf("%d",&C);  
74.                     if(ch1=='g')  
75.                           add(A+B,A-1,-C-1);  
76.                     if(ch1=='l')  
77.                           add(A-1,A+B,C-1);  
78.               }  
79.               if(check())   printf("lamentable kingdom\n");  
80.               else    printf("successful conspiracy\n");  
81.       }  
82.       //system("pause");  
83.       return 0;  
84. }  

 

转载注明出处：http://blog.csdn.net/u011400953