[POJ-1201]Intervals

最新推荐文章于 2022-01-26 15:56:35 发布

转载最新推荐文章于 2022-01-26 15:56:35 发布 · 128 阅读

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/szmssf/p/11085434.html

差分约束系统很棒的应用，首先运用前缀和的思想，令s[k]为0~k中被选出来数的个数，则s[b[i]]>=s[a[i]-1]+c,这样会有一个问题，那就是下标可能会出现负数，因此我们可以改变前缀和表示的方法，设s[k]为0~k-1中被选出来的数，则约束条件可以转化为：s[b[i]+1]-s[a[i]]>=c;简单处理后可得到s[b[i]+1]>=s[a[i]+c;但是只有这一个约束条件远远不够，这只是数据硬性规定的，而还有一些条件是隐含的，有下面两种隐含条件也需要加入约束系统：1）s[k]是前缀和的形式，所以s[k]不可能比s[k-1]小，可得到s[k]>=s[k-1]+0;（2）s[k]最多比s[k-1]多选k这一个数，所以s[k]<=s[k-1]+1，简单处理后得到s[k-1]>=s[k]-1,加入以上两种边后SPFA跑最长路求解即可

下面给出参考代码：

#include<iostream>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;
int t,n,a,b,c,v[550005],w[550005],head[550005],nxt[550005],cnt,dist[500005],maxn;
bool vis[500005];
void add(int a,int b,int c)
{
	v[++cnt]=b;
	w[cnt]=c;
	nxt[cnt]=head[a];
	head[a]=cnt;
}
void spfa(int node)
{
    queue<int>q;
    q.push(node);
    vis[node]=1;
    while(!q.empty())
    {
        int c=q.front();
        q.pop();
        vis[c]=0;
        for(int i=head[c];i;i=nxt[i])
        {
            int y=v[i];
            if(dist[y]<dist[c]+w[i])
            {
                dist[y]=dist[c]+w[i];
                if(!vis[y])
                {
                    q.push(y);
                    vis[y]=1;
                }
            }
        }
    }
}
int main()
{
	scanf("%d\n",&t); 
	while(t--)
	{
		scanf("%d\n",&n);
		memset(v,0,sizeof(v));
		memset(w,0,sizeof(w));
		memset(head,0,sizeof(head));
		memset(nxt,0,sizeof(nxt));
		cnt=0;maxn=-21374404;
		for(int i=1;i<=n;i++)
		{
			scanf("%d%d%d",&a,&b,&c);
			//dist[b+1]>=c+dist[a]
			add(a,b+1,c);
			maxn=max(maxn,b+1);
		}
		for(int i=1;i<=maxn;i++)
		{
			//dist[i]>=0+dist[i-1]
			add(i-1,i,0);
			//dist[i-1]>=dist[i]-1
			add(i,i-1,-1);
		}
		for(int i=1;i<=maxn;i++)
		{
			dist[i]=-21374404;
			add(0,i,0);
		}
		spfa(0);
		printf("%d",dist[maxn]);
		if(t>0)printf("\n");
	}
	return 0;
}

转载于:https://www.cnblogs.com/szmssf/p/11085434.html