题意:告诉你一些关于两地距离的条件,然后判断满足所有条件时的情况是否真实
差分约束系统简析:
如果一个系统由n个变量和m个约束条件组成,其中每个约束条件形如xj-xi<=bk(i,j∈[1,n],k∈[1,m]),则称其为差分约束系统(system of difference constraints)。亦即,差分约束系统是求解关于一组变量的特殊不等式组的方法。
可以看到,xj-xi<=bk可以转化成 xj<=xi+bk 这是最短路中的三角不等式,即为求得最短路的约束条件,如果满足所有约束条件,则存在最短路,不等式有解,如果不满足约束条件,则存在负圈,不等式也无解。
那么将差分约束系统中的不等式转化成图论中的边
例如xj-xi<=bk
则转化成图中的约束条件 xj<=xi+bk 那么即为 xi->xj的一条边,当xj>xi+bk时 xj=xi+bk 进行松弛
注意求解时引入一个超级源点(可能转化成的图会不连通)跑一遍SPFA或Bellman-Ford。
#include<stdio.h>
#include<string.h>
#define INF 1000000000;
typedef struct{
int v,w;
int Val;
}Edge;
Edge E[100000*3];
int Head[1001];
int Next[100000*3],Enum;
int Queue[100000*3],Front,Rear;
int Dist[1001],InQ[1001],QCount[1001];
int N,M;
void AddEdge(int v,int w,int val){
E[Enum].v=v;
E[Enum].w=w;
E[Enum].Val=val;
Next[Enum]=Head[v];
Head[v]=Enum;
Enum++;
}
int SPFA(){
int V,e,i;
memset(InQ,0,sizeof(InQ));
memset(QCount,0,sizeof(QCount));
Front=Rear=0;
for(i=1;i<=N;i++)
Dist[i]=INF;
Dist[0]=0;
Queue[Rear++]=0;
while(Front!=Rear){
V=Queue[Front++];
Front=Front%300000;
InQ[V]=0;
for(e=Head[V];e!=-1;e=Next[e]){
if(Dist[E[e].w]>E[e].Val+Dist[V]){
Dist[E[e].w]=E[e].Val+Dist[V];
if(!InQ[E[e].w]){
InQ[E[e].w]=1;
Queue[Rear++]=E[e].w;
Rear=Rear%300000;
if(++QCount[E[e].w]>N)
return 0;
}
}
}
}
return 1;
}
int main(){
int i,j;
char op[3];
int a,b,c;
while(~scanf("%d %d",&N,&M)){
memset(Head,-1,sizeof(Head));
memset(Next,-1,sizeof(Next));
Enum=0;
for(i=0;i<M;i++){
scanf("%s",op);
if(op[0]=='P'){
scanf("%d %d %d",&a,&b,&c);
AddEdge(a,b,-c);
AddEdge(b,a,c);
}
else{
scanf("%d %d",&a,&b);
AddEdge(a,b,-1);
}
}
for(i=1;i<=N;i++){
AddEdge(0,i,0); //超级源点
}
if(SPFA())
printf("Reliable\n");
else
printf("Unreliable\n");
}
return 0;
}