poj2983&&poj1364

poj 1364

题目大意：

题目大意：n个数的一个序列，m个约数，si, ni, oi, ki, 代表了序列中第si个数到第si+ni个数的和大于或小于ki，gt = 大于，lt = 小于

问是否存在相悖的约束

 将这些约束条件转化为差分约束，不妨设sum[x] = S[1]+S[2]+....S[x]，那么上面式子就可以转化为： 

     1. sum[si+ni] - T[si-1] > ki          2. sum[si+ni] - T[si-1] < ki  

     又差分约束系统的求解只能针对>=或<=，无法求解>的系统，还好ki是整数,可以很容易将<转化为<=  
     1. sum[si-1] - T[si+ni] <= -ki - 1    2. sum[si+ni] - T[si-1] <= ki - 1  

     剩下的就是普通的差分约束系统的求解了

#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>#include<queue>using namespace std;const int N=70010;int n,m;struct node{int x,y,z,next;}sa[N*4];int len=0,first[N];void ins(int x,int y,int z){len++;sa[len].x=x;sa[len].y=y;sa[len].z=z;sa[len].next=first[x];first[x]=len;}queue<int> q;int d[500010],ru[500010];bool u[500010];bool spfa(){    memset(d,0,sizeof(d));    memset(ru,0,sizeof(ru));    memset(u,true,sizeof(u));    for(int i=0;i<=2*m;i++) q.push(i);    while(!q.empty())    {        int x=q.front();        for(int i=first[x];i!=-1;i=sa[i].next)        {            int y=sa[i].y;            if(d[x]+sa[i].z>d[y])            {                d[y]=d[x]+sa[i].z;                ru[i]++;if(ru[i]>m) return false;                if(!u[y]) u[y]=true,q.push(y);            }        }        q.pop();        u[x]=false;    }    return true;}int main(){while(1){scanf("%d",&m);if(m==0) break;m++;scanf("%d",&n);int a,b,c;char s[7];memset(first,-1,sizeof(first));len=0;for(int i=1;i<=n;i++){scanf("%d%d%s%d",&a,&b,s,&c);a++;if(s[0]=='g'){ins(a-1,a+b,c+1);}else{ins(a+b,a-1,-c+1);}}if(!spfa()){printf("successful conspiracy\n");}else{printf("lamentable kingdom\n");}}}

poj2983

题意：给出一些不等式，求是否存在可行解。

从P  A  B  X得到的差分系统为

dist[A] - dist[B] >= X  &&  dist[A] - dist[B] <= X

等价于
dist[B] <= dist[A] - X  &&  dist[A] <= dist[B] + X

从V  A  B得到的差分系统为
dist[A] >= dist[B] +1

等价于
dist[B] <= dist[A] -1 

#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>#include<queue>using namespace std;const int N=2000005;int n,m;struct node{int x,y,z,next;}sa[N*4];int len=0,first[N];void ins(int x,int y,int z){len++;sa[len].x=x;sa[len].y=y;sa[len].z=z;sa[len].next=first[x];first[x]=len;}queue<int> q;int d[N],ru[N];bool u[N];bool spfa(){    memset(d,0,sizeof(d));    memset(ru,0,sizeof(ru));    memset(u,true,sizeof(u));    for(int i=0;i<=n;i++) q.push(i);    while(!q.empty())    {        int x=q.front();u[x]=false;        for(int i=first[x];i!=-1;i=sa[i].next)        {            int y=sa[i].y;u[x]=false;            if(d[x]+sa[i].z>d[y])            {                d[y]=d[x]+sa[i].z;                ru[i]++;if(ru[i]>n) return false;                if(!u[y]) u[y]=true,q.push(y);            }        }        q.pop();    }    return true;}int main(){while(scanf("%d %d",&n,&m)!=EOF){int a,b,c;char s[7];bool tf=true;memset(first,-1,sizeof(first));len=0;for(int i=1;i<=m;i++){scanf("%s",s);if(s[0]=='P'){scanf("%d%d%d",&a,&b,&c);ins(a,b,c);ins(b,a,-c);}else{if(a==b) tf=false;scanf("%d%d",&a,&b);ins(a,b,1);}}if(!spfa()||tf==false){printf("Unreliable\n");}else{printf("Reliable\n");}}}