bzoj3436小K的农场
题意：
n个数，知道m条关系：a-b≥c、a-b≤c或a==b。问是否存在满足所有关系的情况。n≤10000，m≤10000。
题解：
差分约束。因为只要求是否满足，因此最短路最长路都可以。不过要注意如果是用spfa的bfs写法，每个点都必须作为源点判一次负环，因为图可能不连通。正因为如此，虽说加了SLF的bfs写法spfa能卡过，但比dfs写法慢不只10倍。
代码：
1 #include <cstdio> 2 #include <cstring> 3 #include <algorithm> 4 #include <queue> 5 #define ll long long 6 #define inc(i,j,k) for(int i=j;i<=k;i++) 7 #define maxn 10010 8 #define INF 0x3fffffff 9 using namespace std; 10 11 inline int read(){ 12 char ch=getchar(); int f=1,x=0; 13 while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();} 14 while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar(); 15 return f*x; 16 } 17 18 struct e{int t,w,n;}; e es[maxn*40]; int g[maxn],ess; 19 void pe(int f,int t,int w){es[++ess]=(e){t,w,g[f]}; g[f]=ess;} 20 int n,m,cnt[maxn],d[maxn]; bool inq[maxn]; deque<int>q; 21 ll spfa(int s){ 22 q.clear(); memset(inq,0,sizeof(inq)); memset(cnt,0,sizeof(cnt)); 23 q.push_back(s); inq[s]=1; d[s]=0; cnt[s]=1; 24 while(!q.empty()){ 25 int x=q.front(); q.pop_front(); inq[x]=0; 26 for(int i=g[x];i;i=es[i].n)if(d[es[i].t]>d[x]+es[i].w){ 27 d[es[i].t]=d[x]+es[i].w; 28 if(!inq[es[i].t]){ 29 if(!q.empty()&&d[es[i].t]<d[q.front()])q.push_front(es[i].t);else q.push_back(es[i].t); 30 inq[es[i].t]=1; cnt[es[i].t]++; if(cnt[es[i].t]>=n)return 0; 31 } 32 } 33 } 34 return 1; 35 } 36 int main(){ 37 n=read(); m=read(); 38 inc(i,1,m){ 39 int opt=read(),a,b,c; 40 if(opt==1)a=read(),b=read(),c=read(),pe(a,b,-c); 41 if(opt==2)a=read(),b=read(),c=read(),pe(b,a,c); 42 if(opt==3)a=read(),b=read(),pe(b,a,0),pe(a,b,0); 43 } 44 inc(i,1,n)if(!spfa(i)){puts("No"); return 0;} puts("Yes"); return 0; 45 }
20161018
