51nod 1781 Pinball

昨天看到时，这题还是三级题的。然后各种想不出。大佬说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);
	}
}