POJ2983 差分约束系统

题意：告诉你一些关于两地距离的条件，然后判断满足所有条件时的情况是否真实

 

差分约束系统简析：

如果一个系统由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;
}