观光公交 [NOIP 2011] [思维推导]

Description

风景迷人的小城Y 市，拥有n 个美丽的景点。由于慕名而来的游客越来越多，Y 市特意安排了一辆观光公交车，为游客提供更便捷的交通服务。观光公交车在第0 分钟出现在1号景点，随后依次前往2、3、4……n 号景点。从第i 号景点开到第i+1 号景点需要Di 分钟。
任意时刻，公交车只能往前开，或在景点处等待。
设共有m 个游客，每位游客需要乘车1 次从一个景点到达另一个景点，第i 位游客在Ti 分钟来到景点Ai，希望乘车前往景点Bi（Ai＜bi）。为了使所有乘客都能顺利到达目的地，公交车在每站都必须等待需要从该景点出发的所有乘客都上车后才能出发开往下一景点。
假设乘客上下车不需要时间。
一个乘客的旅行时间，等于他到达目的地的时刻减去他来到出发地的时刻。因为只有一辆观光车，有时候还要停下来等其他乘客，乘客们纷纷抱怨旅行时间太长了。于是聪明的司机ZZ
给公交车安装了k 个氮气加速器，每使用一个加速器，可以使其中一个Di 减1。对于同一个Di 可以重复使用加速器，但是必须保证使用后Di
大于等于0。
那么ZZ 该如何安排使用加速器，才能使所有乘客的旅行时间总和最小？

Input

第1 行是3 个整数n, m, k，每两个整数之间用一个空格隔开。分别表示景点数、乘客数和氮气加速器个数。
第2 行是n-1 个整数，每两个整数之间用一个空格隔开，第i 个数表示从第i 个景点开往第i+1 个景点所需要的时间，即Di。
第3 行至m+2 行每行3 个整数Ti, Ai, Bi，每两个整数之间用一个空格隔开。第i+2 行表示第i 位乘客来到出发景点的时刻，出发的景点编号和到达的景点编号。

Output

共一行，包含一个整数，表示最小的总旅行时间。

Sample Input

3 3 2
1 4
0 1 3
1 1 2
5 2 3

Sample Output

Hint

对D2 使用2 个加速器，从2 号景点到3 号景点时间变为2 分钟。
公交车在第1 分钟从1 号景点出发，第2 分钟到达2 号景点，第5 分钟从2 号景点出发，第7 分钟到达3 号景点。
第1 个旅客旅行时间 7-0 = 7 分钟。
第2 个旅客旅行时间 2-1 = 1 分钟。
第3 个旅客旅行时间 7-5 = 2 分钟。
总时间 7+1+2 = 10 分钟。
数据范围：
对于10%的数据，k=0；
对于20%的数据，0≤k≤1；
对于40%的数据，2 ≤ n ≤ 50，1 ≤ m≤ 1,000，0 ≤ k ≤ 20，0 ≤ Di ≤ 10，0 ≤ Ti ≤ 500；
对于60%的数据，1 ≤ n ≤ 100，1 ≤ m≤ 1,000，0 ≤ k ≤ 100，0 ≤ Di ≤ 100，0 ≤ Ti ≤ 10,000；
对于100%的数据，1 ≤ n ≤ 1,000，1 ≤ m ≤ 10,000，0 ≤ k ≤ 100,000，0 ≤ Di ≤ 100，0 ≤ Ti ≤ 100,000。

思路

第一眼看到题，我认为，我是没有思路的这应该是一道思维推导题

我们应该从何下手？

首先列出几个显而易见的观点

　　1.一个点的出发时间 leave_i = max{ arrive_i , max{ get_on_person_i } };

　　2.一个人的在车上时间=get_off_i - arrive_station_i

我们抓住观点2，显然到达时间是固定的，我们要做的是尽量把上车时间提前

如何将氮气加速和提前时间搭建起关系？

对比每一个点到达时间和最后一个人上车时间，定义晚点为到达时间晚于最后一个人上车时间，于在i晚点后会对后面一段区间产生延误，直到某一个点j的最后一个人上车时间晚于列车到达时间，则不会继续对后面产生影响。

对于i设它对后面的波及范围为f(i)

则它的递推式为　　f ( i )= arrive(i+1)<=last_get_on(i+1) ? i+1 : f ( i+1)

所以f函数我们要从后往前推

那么f（i）表示了影响区间，那么用前缀和求出区间中的个数，最后选择最大的贪心，循环K次后即可AC

代码

 #include<set>

 #include<map>

 #include<stack>

 #include<queue>

 #include<cstdio>

 #include<cstring>

 #include<iostream>

 #include<algorithm>

 #define RG register int

 #define rep(i,a,b)    for(RG i=a;i<=b;i++)

 #define per(i,a,b)    for(RG i=a;i>=b;i--)

 #define ll long long

 #define inf (1<<30)

 #define maxn 10005

 using namespace std;

 int n,m,k,ans;

 int val[maxn],st[maxn],ed[maxn],up[maxn],lat[maxn],downs[maxn];

 int arrive[maxn],range[maxn];

 inline int read()

 {

     int x=,f=;char c=getchar();

     while(c<''||c>''){if(c=='-')f=-;c=getchar();}

     while(c>=''&&c<=''){x=x*+c-'';c=getchar();}

     return x*f;

 }

 int main()

 {

     n=read(),m=read(),k=read();

     rep(i,,n-) val[i]=read();

     rep(i,,m)    up[i]=read(),st[i]=read(),ed[i]=read(),lat[st[i]]=max(lat[st[i]],up[i]),downs[ed[i]]++;

     rep(i,,n)    downs[i]+=downs[i-];

     arrive[]=lat[];

     rep(i,,n)    arrive[i]=max(arrive[i-],lat[i-])+val[i-];

     rep(i,,m)    ans+=arrive[ed[i]]-up[i];

     while(k--)

     {

         range[n]=range[n-]=n;

         per(i,n-,)    range[i]=arrive[i+]<=lat[i+]?i+:range[i+];

         int mx=,mxid;

         rep(i,,n-)

             if(downs[range[i]]-downs[i]>mx&&val[i])    mx=downs[range[i]]-downs[i],mxid=i;

         ans-=mx,val[mxid]--;

         rep(i,,n) arrive[i]=max(arrive[i-],lat[i-])+val[i-];

     }

     printf("%d",ans);

     return ;

 }