51nod 1781 Pinball

来源：互联网发布：abbyy12序列号mac 编辑：程序博客网时间：2024/06/10 03:39

昨天看到时，这题还是三级题的。然后各种想不出。大佬说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);elsechange(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);elsereturn 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");elseprintf("%lld\n",ans);}}

0 0