题目:https://vijos.org/p/1901

题解:这题比较有意思。

经过一番思考之后我想出了下面的算法:

我们反着来推,按i从大到小

f[i]表示从>=m到 i 需要多长时间,则如果v[j]=i,则我们可以用f[i]+t[j] 去更新 r[j]-inf 的所有f值

由于f[i]是一个单调递减的函数,则我们只需要用它去更新r[j]-i-1即可

然后问题就变成了区间修改,单点查询,线段树即可

代码:

 #include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 100000000000000ll
#define maxn 250000
#define maxm 500+100
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define mod 1000000007
using namespace std;
inline int read()
{
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}
return x*f;
}
int n,m,a[*maxn],b[*maxn],c[*maxn];
ll ti[maxn];
struct seg{int k,l,r;ll mi,tag;}t[*maxn];
inline void build(int k,int l,int r)
{
t[k].l=l;t[k].r=r;int mid=(l+r)>>;
t[k].mi=inf;t[k].tag=inf;
if(l==r)return;
build(k<<,l,mid);build(k<<|,mid+,r);
}
inline void update(int k,ll z)
{
t[k].tag=min(t[k].tag,z);
t[k].mi=min(t[k].mi,z);
}
inline void pushdown(int k)
{
if(t[k].tag==inf)return;
update(k<<,t[k].tag);
update(k<<|,t[k].tag);
t[k].tag=inf;
}
inline void pushup(int k)
{
t[k].mi=min(t[k<<].mi,t[k<<|].mi);
}
inline void change(int k,int x,int y,ll z)
{
int l=t[k].l,r=t[k].r,mid=(l+r)>>;
if(l==x&&r==y){update(k,z);return;}
pushdown(k);
if(y<=mid)change(k<<,x,y,z);
else if(x>mid)change(k<<|,x,y,z);
else change(k<<,x,mid,z),change(k<<|,mid+,y,z);
pushup(k);
}
inline ll query(int k,int x)
{
int l=t[k].l,r=t[k].r,mid=(l+r)>>;
if(l==r)return t[k].mi;
pushdown(k);
if(x<=mid)return query(k<<,x);else return query(k<<|,x);
}
inline bool cmp1(int x,int y){return a[x]<a[y];}
inline bool cmp2(int x,int y){return b[x]>b[y];}
int main()
{
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
int cs=read(),mm=;
while(cs--)
{
n=read();m=read();
for1(i,n)a[i]=read(),a[n+i]=read(),ti[i]=read();a[*n+]=m;a[*n+]=;
for1(i,*n+)c[i]=i;
sort(c+,c+*n+,cmp1);
int tot=;
for1(i,*n+)
{
if(i==||a[c[i]]!=a[c[i-]])tot++;
b[c[i]]=tot;
}
for1(i,n)c[i]=i;
sort(c+,c+n+,cmp2);
build(,,tot);
change(,b[*n+],tot,);
for1(i,n)
{
ll x=query(,b[c[i]]);if(x==inf)continue;
if(b[c[i]]>b[n+c[i]])change(,b[n+c[i]],b[c[i]]-,x+ti[c[i]]);
}
ll ans=query(,b[*n+]);
printf("Case #%d: ",++mm);
if(ans==inf)printf("-1\n");else printf("%lld\n",ans);
}
return ;
}
05-11 14:02