[BZOJ2330]SCOI2011糖果|差分约束

差分约束题。。学了一下差分约束，我觉得还是挺简单的，考虑f[u]-f[v]<=c，发现和最短路的松弛操作神似，最短路跑完之后对于一条边（v,u），显然有d[u]<=d[v]+c，不就是上面差分约束的式子吗。。那就转化成最短(长)路做咯，以最短路为例了，对于每个f[u]-f[v]<=c，连一条v到u权值为c的边，然后跑最短路，如果有负环那就说明无解，没有的话跑完之后各点的d值就是解。。

对于本题，d[a]==d[b]，有 d[a]-d[b]>=0, d[b]-d[a]>=0

d[a]<d[b]，有 d[b]-d[a]>=1

d[a]>=d[b]，有 d[a]-d[b]>=0

d[a]>d[b]，有 d[a]-d[b]>=1

d[a]<=d[b]，有 d[b]-d[a]>=0

这里要求一组和最小的解，而每个点要>=1，所以用最长路做，每个点初始赋为1，然后一开始把它们都扔进队列，跑最长路即可。。

对于差分约束来说，无论跑最长路还是最短路都可以出解，但是我们求出的都只是一组相对大小，要利用这个相对大小来求一组符合题目的解。。像这题，每个点的值要求大于0，所以我们赋初值1，要保证在运行算法的过程中每个点的值始终大于0，那么就是跑最长路了。。

[cpp]  view plain copy print ?

1. #include<iostream>  
2. #include<cstdio>  
3. #include<memory.h>  
4. #define N 100005  
5. using namespace std;  
6. struct edge{  
7.     int e,q,next;  
8. }ed[N*2];  
9. int n,k,i,s,t,opt,ne=0,a[N],dis[N],que[N*2],u[N],inq[N];  
10. long long ans=0ll;  
11. void add(int s,int e,int q)  
12. {  
13.     ed[++ne].e=e;ed[ne].q=q;  
14.     ed[ne].next=a[s];a[s]=ne;  
15. }  
16. bool spfa()  
17. {  
18.     int head=1,tail=n,get,hh=1,tt=n,j,i,to;  
19.     for (i=1;i<=n;i++) inq[i]=dis[i]=u[i]=1,que[i]=i;  
20.     while (hh<=tt)  
21.     {  
22.         get=que[head++];hh++;  
23.         if (head>200000) head=1;  
24.         for (j=a[get];j;j=ed[j].next)  
25.             if (dis[get]+ed[j].q>dis[to=ed[j].e])  
26.             {  
27.                 dis[to]=dis[get]+ed[j].q;  
28.                 if (!inq[to])  
29.                 {  
30.                     if (++u[to]>=n) return false;  
31.                     tail++;tt++;  
32.                     if (tail>200000) tail=1;  
33.                     que[tail]=to;  
34.                     inq[to]=1;  
35.                 }  
36.             }  
37.         inq[get]=0;  
38.     }  
39.     return true;  
40. }  
41. int main()  
42. {  
43.     scanf("%d%d",&n,&k);  
44.     for (i=1;i<=n;i++) a[i]=0;  
45.     for (i=1;i<=k;i++)  
46.     {  
47.         scanf("%d%d%d",&opt,&s,&t);  
48.         switch (opt)  
49.         {  
50.             case 1:add(s,t,0);add(t,s,0);break;  
51.             case 2:if (s==t) {puts("-1");return 0;}  
52.                 add(s,t,1);break;  
53.             case 3:add(t,s,0);break;  
54.             case 4:if (s==t) {puts("-1");return 0;}  
55.                 add(t,s,1);break;  
56.             case 5:add(s,t,0);break;  
57.         }  
58.     }  
59.     if (spfa())  
60.     {  
61.         for (i=1;i<=n;i++) ans+=(int)dis[i];  
62.         printf("%lld",ans);  
63.     }  
64.     else printf("-1");  
65. }