昨天看到时,这题还是三级题的。然后各种想不出。大佬说DP+线段树+离散,想了想可以做。然后今天改成了五级题了。。。
L[i]表示球从1位置落下,最后从第i个漏斗出去的情况下,的最小费用。
R[i]表示球从n位置落下,最后从第i个漏斗出去的情况下,的最小费用。
那么答案就是min(L[i]+R[i]-d[i])。
动态转移方程:L[i]=min{L[j] | a[i]≤c[j]≤b[i]} + d[i]
由于转移时是求区间最小值,所以可以用线段树优化。
n十分大,所以要离散。
然后被一点坑到了,一直RE。。。。共M个漏斗,离散后得到3M个点,所以线段树要至少开9M大小。。。。
#include<bits/stdc++.h>
using namespace std;
const int MAXN=100100;
const long long inf=1ll<<60;
struct node
{
int l,r;
long long nsum;
}segtree[MAXN<<4];
void build(int i,int l,int r)
{
segtree[i].l=l;
segtree[i].r=r;
if(l==r)
{
segtree[i].nsum=inf;
return;
}
int mid=(l+r)>>1;
build(i<<1,l,mid);
build(i<<1|1,mid+1,r);
segtree[i].nsum=min(segtree[i<<1].nsum,segtree[i<<1|1].nsum);
}
void change(int i,int a,long long c)
{
if(segtree[i].l==a&&segtree[i].r==a)
{
segtree[i].nsum=min(segtree[i].nsum,c);
return;
}
int mid=(segtree[i].l+segtree[i].r)>>1;
if(a<=mid)
change(i<<1,a,c);
else
change(i<<1|1,a,c);
segtree[i].nsum=min(segtree[i<<1].nsum,segtree[i<<1|1].nsum);
}
long long query(int i,int a,int b)
{
if(segtree[i].l==a&&segtree[i].r==b)
return segtree[i].nsum;
int mid=(segtree[i].l+segtree[i].r)>>1;
if(b<=mid)
return query(i<<1,a,b);
else if(a>mid)
return query(i<<1|1,a,b);
else
return min(query(i<<1,a,mid),query(i<<1|1,mid+1,b));
}
long long a[MAXN],b[MAXN],c[MAXN],d[MAXN],sca[MAXN<<2],l[MAXN],r[MAXN];
void read(long long &a){
char ch;while(!((ch=getchar())>='0')&&(ch<='9'));
a=ch-'0';while(((ch=getchar())>='0')&&(ch<='9'))a*=10,a+=ch-'0';
}
int main()
{
long long n,m,i,scanum,ltmp,rtmp,ans;
while(~scanf("%lld%lld",&m,&n))
{
scanum=1;
for(i=1;i<=m;i++)
{
read(a[i]); read(b[i]); read(c[i]); read(d[i]);
sca[scanum++]=a[i];
sca[scanum++]=b[i];
sca[scanum++]=c[i];
}
sca[scanum++]=1;
sca[scanum++]=n;
sort(sca+1,sca+scanum);
scanum=unique(sca+1,sca+scanum)-sca;
for(i=1;i<=m;i++)
{
a[i]=lower_bound(sca+1,sca+scanum,a[i])-sca;
b[i]=lower_bound(sca+1,sca+scanum,b[i])-sca;
c[i]=lower_bound(sca+1,sca+scanum,c[i])-sca;
}
n=scanum-1;
build(1,1,n);
change(1,1,0);
for(i=1;i<=m;i++)
{
l[i]=query(1,a[i],b[i])+d[i];
change(1,c[i],l[i]);
}
build(1,1,n);
change(1,n,0);
for(i=1;i<=m;i++)
{
r[i]=query(1,a[i],b[i])+d[i];
change(1,c[i],r[i]);
}
ans=inf;
for(i=1;i<=m;i++)
ans=min(ans,l[i]+r[i]-d[i]);
if(ans==inf)
printf("-1\n");
else
printf("%lld\n",ans);
}
}