差分约束:如果一个系统由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]}即为一组可行解。
解题思路:
1 2 gt 0
a1+a2+a3>0
2 2 lt 2
a2+a3+a4<2
a1+a2+a3>0
2 2 lt 2
a2+a3+a4<2
首先设Si=a1+a2+a3+...+ai
那么根据样例可以得出
S3-S0>0---->S0-S3<=-1
S4-S1<2---->S4-S1<=1
那么根据差分约束建图,加入这些有向边
gt: <a+b,a-1>=-ki-1
lt: <a-1,a+b>=ki-1
再根据bellman_ford判断是否有无负环即可
若出现负环了则这个序列不满足所有的不等式
源码:
#include <iostream>
#include <memory.h>
using namespace std;
typedef struct
{
int s,t,w;
}Edge;
Edge edge[110];
int n,m;
bool bellman_ford()
{
int dist[110];
memset(dist,0,sizeof(dist));
for(int i=1;i<=n;i++)
for(int j=0;j<m;j++)
if(dist[edge[j].s]+edge[j].w<dist[edge[j].t])
dist[edge[j].t] = dist[edge[j].s]+edge[j].w;
for(int j=0;j<m;j++)
if(dist[edge[j].s]+edge[j].w<dist[edge[j].t])
return 0;
return 1;
}
int main()
{
int a,b,c;
char s[3];
while(cin>>n)
{
if(!n)
break;
cin>>m;
for(int i=0;i<m;i++)
{
cin>>a>>b>>s>>c;
if(s[0]=='g')
{
edge[i].s=b+a;
edge[i].t=a-1;
edge[i].w=-c-1;
}
else
{
edge[i].s=a-1;
edge[i].t=b+a;
edge[i].w=c-1;
}
}
if(!bellman_ford())
cout<<"successful conspiracy"<<endl;
else
cout<<"lamentable kingdom"<<endl;
}
return 0;
}