題意:
給出n個點的m條約束信息。每條信息表述爲(P a b c)表示a在b北方距離c的位置,或者(V a b) 表示a在b北方1單位距離或者更遠的位置。問是否可能存在符合以上m個要求的點。
解題思路:
差分約束一般給我的印象都是一個範圍性的式子才能寫出不等式,但是這道題,固定的距離可以用兩個只有不等式號相反的式子表示,這樣子可以將值限定在“=”的那個值上。
把dis[i]設爲其到始點的距離。第二個條件很簡單dis[a]-dis[b]>=1 也就是dis[b]<=dis[a]-1。對於第一個,帶等於號的條件dis[a]-dis[b]==c,我們可以轉化爲dis[a]-dis[b]>=c和dis[a]-dis[b]<=c
兩個不等式,然後再轉化爲最短路三角不等式。由於可能出現多個不連通圖的情況,所以要設一個虛擬始點,與所有點相連,長度設爲0。然後用spfa來判定是否有可行解。
代碼:
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<queue>
#include<vector>
#include<set>
#include<stack>
#include<map>
#include<ctime>
#include<bitset>
#define LL long long
#define db double
#define EPS 1e-8
#define inf 1e9
using namespace std;
const int nMax=200050;
const int mMax=1000050;
struct{
int u,v,nxt;
int w;
}e[mMax];
int n,cnt,head[nMax];
int dis[nMax];
int que[nMax],m,sum[nMax];
bool vis[nMax];
void add(int u,int v,int w){
e[cnt].w=w;
e[cnt].u=u; e[cnt].v=v;
e[cnt].nxt=head[u];
head[u]=cnt; cnt++;
}
bool spfa(int s){
int i,hhead=0,tail=1;
for (i=0;i<=n;i++){
dis[i]=inf;
vis[i]=0;
}
dis[s]=0; vis[s]=1; que[0]=s;
while (hhead!=tail){
int u=que[hhead];
vis[u]=0;
for (int p=head[u];p!=0;p=e[p].nxt){
int v=e[p].v;
if (dis[v]>dis[u]+e[p].w){
dis[v]=dis[u]+e[p].w;
if (!vis[v]){
vis[v]=1;
que[tail++]=v;
if (tail==nMax) tail=0;
if (++sum[v]>n) return 0;
}
}
}
hhead++;
if (hhead==nMax) hhead=0;
}
return 1;
}
int main(){
int m,a,b,c,s;
char str[20];
while (cin>>n>>m){
s=0; cnt=1;
memset(sum,0,sizeof(sum));
memset(head,0,sizeof(head));
memset(vis,0,sizeof(vis));
while (m--){
scanf("%s",str);
if (str[0]=='P'){
scanf("%d%d%d",&a,&b,&c);
add(b,a,c);
add(a,b,-c);
}
else {
scanf("%d%d",&a,&b);
add(a,b,-1);
}
}
for (int i=1;i<=n;i++){
add(s,i,0);
}
if (spfa(s))
printf("Reliable\n");
else printf("Unreliable\n");
}
return 0;
}