【BZOJ2300】【SCOI2011】糖果

最新推荐文章于 2018-06-16 15:42:01 发布

转载最新推荐文章于 2018-06-16 15:42:01 发布 · 45 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/JSL2018/p/6539363.html

本文介绍了一道关于幼儿园小朋友分配糖果的问题，通过差分约束系统与SPFA算法解决该问题。讨论了如何建立差分约束模型，并利用SPFA算法求解最短路径以确保所有小朋友都能得到糖果并满足特定条件。

差点就忘了还有差分约束这个东西……看见了就要学习一个

原题：

幼儿园里有N个小朋友，lxhgww老师现在想要给这些小朋友们分配糖果，要求每个小朋友都要分到糖果。但是小朋友们也有嫉妒心，总是会提出一些要求，比如小明不希望小红分到的糖果比他的多，于是在分配糖果的时候，lxhgww需要满足小朋友们的K个要求。幼儿园的糖果总是有限的，lxhgww想知道他至少需要准备多少个糖果，才能使得每个小朋友都能够分到糖果，并且满足小朋友们所有的要求。

N<=100000，K<=100000，1<=A, B<=N

差分约束模板题，学习一个！

以操作3为例，如果a不少于b，那么a-b>=0，a>=b+0

所以连b到a权值为0的边，根据三角形不等式，要跑最长路来使条件满足

具体原理我也想不太清楚，反正记住搞出三角形不等式后跑符号相反的spfa就对了 _(:3 」∠)_

因为写出了像"q[tl]=false"酱sb的东西+没对拍所以没1A……

注意longlong

代码：

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 #include<cmath>
 6 using namespace std;
 7 int rd(){int z=0,mk=1;  char ch=getchar();
 8     while(ch<'0'||ch>'9'){if(ch=='-')mk=-1;  ch=getchar();}
 9     while(ch>='0'&&ch<='9'){z=(z<<3)+(z<<1)+ch-'0';  ch=getchar();}
10     return z*mk;
11 }
12 struct ddd{int nxt,y,v;}e[210000];  int lk[110000],ltp=0;
13 inline void ist(int x,int y,int z){  e[++ltp].nxt=lk[x],lk[x]=ltp,e[ltp].y=y,e[ltp].v=z;}
14 int n,m;
15 int dstc[110000],cnt[110000];
16 int q[110000],hd=0,tl=0,tp=100001;  bool vstd[110000];
17 bool spfa(){
18     memset(vstd,0,sizeof(vstd));
19     //memset(dstc,0,sizeof(dstc));
20     for(int i=1;i<=n;++i)  q[++hd]=i,dstc[i]=1;
21     while(tl!=hd){
22         tl=(tl==tp ? 0 : tl+1);
23         for(int i=lk[q[tl]];i;i=e[i].nxt)if(dstc[q[tl]]+e[i].v>dstc[e[i].y]){
24             dstc[e[i].y]=dstc[q[tl]]+e[i].v;
25             if(++cnt[e[i].y]==n)  return false;
26             if(!vstd[e[i].y])  q[hd=(hd==tp ? 0 : hd+1)]=e[i].y,vstd[e[i].y]=true;
27         }
28         vstd[q[tl]]=false;
29     }
30     return true;
31 }
32 int main(){//freopen("ddd.in","r",stdin);
33     cin>>n>>m;
34     int mk,l,r;
35     while(m--){
36         mk=rd(),l=rd(),r=rd();
37         if(mk==1)  ist(l,r,0),ist(r,l,0);
38         else if(mk==2){
39             if(l==r){  cout<<-1<<endl;  return 0;}
40             ist(l,r,1);
41         }
42         else if(mk==3)  ist(r,l,0);
43         else if(mk==4){
44             if(l==r){  cout<<-1<<endl;  return 0;}
45             ist(r,l,1);
46         }
47         else if(mk==5)  ist(l,r,0);
48     }
49     if(!spfa()){  cout<<-1<<endl;  return 0;}
50     long long ans=0;
51     for(int i=1;i<=n;++i)  ans+=dstc[i];
52     cout<<ans<<endl;
53     return 0;
54 }