POJ 1201&&HDU 1384 Intervals(差分制約システム)
タイトルアドレス:POJ 1201 HDU 1384
タイトルの意味によって、不等式を以下のようにリストすることができます.
Sj-Si>=c;
Si-S(i-1)>=0;
S(i-1)-Si>=-1;
次いで,この不等式を最短spfaで解決した.maxをソースポイントとし、0を終点とします.最終的な-d[0]が答えです.
コードは次のとおりです.
タイトルの意味によって、不等式を以下のようにリストすることができます.
Sj-Si>=c;
Si-S(i-1)>=0;
S(i-1)-Si>=-1;
次いで,この不等式を最短spfaで解決した.maxをソースポイントとし、0を終点とします.最終的な-d[0]が答えです.
コードは次のとおりです.
#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <ctype.h>
#include <queue>
#include <map>
#include <set>
#include <algorithm>
using namespace std;
const int INF=0x3f3f3f3f;
int d[60000], vis[60000], head[60000], cnt, q[5000000];
struct node
{
int u, v, w, next;
} edge[1000000];
void add(int u, int v, int w)
{
edge[cnt].v=v;
edge[cnt].w=w;
edge[cnt].next=head[u];
head[u]=cnt++;
}
void spfa(int s)
{
memset(d,INF,sizeof(d));
d[s]=0;
memset(vis,0,sizeof(vis));
int ss=0, ee=0;
q[ss++]=s;
while(ss>ee)
{
int u=q[ee++];
vis[u]=0;
for(int i=head[u]; i!=-1; i=edge[i].next)
{
int v=edge[i].v;
if(d[v]>d[u]+edge[i].w)
{
d[v]=d[u]+edge[i].w;
if(!vis[v])
{
vis[v]=1;
q[ss++]=v;
}
}
}
}
}
int main()
{
int a, b, c, n, i, max1, min1;
while(scanf("%d",&n)!=EOF)
{
memset(head,-1,sizeof(head));
max1=-1;
min1=INF;
cnt=0;
while(n--)
{
scanf("%d%d%d",&a,&b,&c);
if(max1<b)
max1=b;
if(min1>a)
min1=a;
add(b,a-1,-c);
}
for(i=1; i<=max1; i++)
{
add(i-1,i,1);
add(i,i-1,0);
}
spfa(max1);
printf("%d
",-d[0]);
}
return 0;
}