差分约束系统很棒的应用,首先运用前缀和的思想,令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; }