【POJ1201】Intervals——差分约束
我的第一道差分约束系统......刚开始有点摸不着头脑,搞清楚后发现完全就是套路...
题目大意:在线段[ai,bi]上至少选ci个点,使被选出的点的个数最少而且满足所有的限制条件,输出这个最小值。
分析:差分约束系统我就不在这里介绍啦~\(≧▽≦)/~(这种东西网上一搜一大堆)
我们设s[i]表示区间[0,i]的整数个数,这道题的差分约束系统为:
s[b]-s[a-1]>=ci;
s[b]-a[b-1]>=0;
s[b-1]-s[b]>=-1。
然后我们再以ai-1或bi的值为结点,ci为权值连有向边,用SPFA求出从min跑到max的最长路就是答案啦。
至于为什么是最长路,不理解请自行跳转到此链接
ps.注意有向边的方向!!!
代码:
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
#include<cstdio> #include<algorithm> #include<cstring> const int maxn=5e5+5,inf=0x3f3f3f3f; using namespace std; struct point{ int next,w,to; }e[maxn*3]; int s,t,tot=0,first[maxn],minn=inf,maxx=0,q[maxn]; int dis[maxn]; bool vis[maxn]; int read() { int ans=0,f=1;char c=getchar(); while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();} while(c>='0'&&c<='9'){ans=ans*10+c-48;c=getchar();} return ans*f; } void add(int u,int v,int wi) { tot++;e[tot].next=first[u];first[u]=tot;e[tot].to=v;e[tot].w=wi; } void spfa() { int head=0,tail=1; q[head]=minn;dis[minn]=0;vis[minn]=1; while(head!=tail){ int x=q[head];head++;if(head>=50000)head=0; for(int i=first[x];i;i=e[i].next){ int to=e[i].to; if(dis[to]<dis[x]+e[i].w){ dis[to]=dis[x]+e[i].w; if(!vis[to]){ vis[to]=1; q[tail]=to; tail++; if(tail>=50000)tail=0; } } } vis[x]=0; } } int main() { int n=read(),a,b,c; for(int i=1;i<=n;i++){ a=read();b=read();c=read(); add(a-1,b,c); minn=min(minn,a-1); maxx=max(maxx,b); } for(int i=minn;i<=maxx;i++){ add(i-1,i,0); add(i,i-1,-1); } memset(dis,-127,sizeof(dis)); memset(vis,0,sizeof(vis)); spfa(); printf("%d",dis[maxx]); return 0; }