NOIP 2011 观光公交

题目描述

风景迷人的小城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该如何安排使用加速器，才能使所有乘客的旅行时间总和最小？

输入输出格式

输入格式：

输入文件名为bus.in。

第1 行是3 个整数n, m, k ，每两个整数之间用一个空格隔开。分别表示景点数、乘客数和氮气加速器个数。

第2 行是n-1 个整数，每两个整数之间用一个空格隔开，第i 个数表示从第i 个景点开往第i+1 个景点所需要的时间，即 Di 。

第3 行至m+2 行每行3 个整数 Ti, Ai, Bi，每两个整数之间用一个空格隔开。第 i+2 行表示第i 位乘客来到出发景点的时刻，出发的景点编号和到达的景点编号。

输出格式：

输出文件名为bus.out。共一行，包含一个整数，表示最小的总旅行时间。

输入输出样例

输入样例#1：

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

输出样例#1：

10

说明

【输入输出样例说明】

对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% 的数据，k=1 ；

对于40% 的数据，2 ≤ n ≤ 50，1 ≤ m ≤ 1,000，0 ≤ k ≤ 20，0 ≤ Di ≤ 10，0 ≤ T i ≤ 500；

对于60% 的数据，1 ≤ n ≤ 100，1 ≤ m ≤ 1,000，0 ≤ k ≤ 100 ，0 ≤ Di ≤ 100，0 ≤ T i ≤ 10,000 ；

对于100%的数据，1 ≤ n ≤ 1,000，1 ≤ m ≤ 10,000 ，0 ≤ k ≤ 100,000，0 ≤ Di ≤ 100 ，0 ≤ T i ≤ 100,000。

令car[i]为车到i的时间,st[i]为i位置最大的t，sum[i]为下车人数的前缀和

显然car[i]=max(car[i-1],st[i-1])+d[i-1];

ans=∑car[b[i]]-t[i]

t是固定的，只要使∑car[b[i]]最小就行

但是不能简单的取下车人数最大的路加速

因为可能会有一条路径，它的值减小可能会影响多条路径

令f[i]表示i能影响到的路径的终点

f[i]=f[i+1]      st[i]<car[i]

    =i+1          st[i]>=car[i]

然后枚举路径，取出路径i，使sum[f[i]]-sum[i]最大

给它一个加速器，然后d值变化

又使car变化，所以要重算一边

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct Peo
{
  int t,l,r;
}a[100001];
long long ans;
int d[10001],n,m,k,car[10001],st[10001],b[10001],f[10001];
int main()
{int i,j;
  cin>>n>>m>>k;
  for (i=1;i<=n-1;i++)
    {
      scanf("%d",&d[i]);
    }
  for (i=1;i<=m;i++)
    {
      scanf("%d%d%d",&a[i].t,&a[i].l,&a[i].r);
      st[a[i].l]=max(st[a[i].l],a[i].t);
      b[a[i].r]++;
    }
  for (i=2;i<=n;i++)
    b[i]+=b[i-1];
  for (i=1;i<=n-1;i++)
    car[i+1]=max(car[i],st[i])+d[i];
  for (i=1;i<=m;i++)
    ans+=car[a[i].r]-a[i].t;
  while (k--)
    {
      f[n-1]=n;
      for (i=n-2;i>=1;i--)
    {
      if (car[i+1]<=st[i+1])
        {
          f[i]=i+1;
        }
      else f[i]=f[i+1];
    }
      int maxx=0;j=0;
      for (i=n-1;i>=1;i--)
    if (b[f[i]]-b[i]>maxx&&d[i])
      {
        maxx=b[f[i]]-b[i];
        j=i;
      }
      if (maxx==0) break;
      ans-=maxx;
      d[j]--;
       for (i=1;i<=n-1;i++)
        car[i+1]=max(car[i],st[i])+d[i];
    }
  cout<<ans;
}