[NOIP2011day2]观光公交 贪心

观光公交

题目描述

风景迷人的小城 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。共一行，包含一个整数，表示最小的总旅行时间。

样例输入

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

样例输出

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 分钟。
对于 100%的数据，1 <=n <=1,000，1 <=m <=10,000，0 <=k <=100,000，0 <=Di <=100，
0 <=Ti <=100,000

我的思路

• 先想了30分钟的贪心，发现了这样的式子

定义$arr[i]$为到达i点的时间，$last[i]$为最后一个到达i点人的时间
$Ans=\sum arr[Bi]-\sum Ti~~~~~~~~~ (1<=i<=m)$

$arr[i]=max(last[i-1]，arr[i-1])+Di-vi vi<=Di$

$Ans=\sum arr[i]*cnt[i]-\sum Tj (1<=i<=n,1<=j<=m)$

要求$\sum arr[i]*cnt[i]$最小

• 然后我们的决策就是在什么地方使用加速器了。。
• 由于贪心弄得有点烦了，先敲了一个$Dp$

struct P1{
    int dp[55][105][M];
    /*dp[i][j][p]表示到达i点，使用了j个加速器，到达i-1的时间为p的最小代价*/
    void work(){
        memset(dp,63,sizeof(dp));
        dp[1][0][0]=0;
        For(i,1,n+1)For(j,0,k)For(p,0,1000)For(c,0,D[i])if(dp[i][j][p]<1e9){
            int Arr=max(last[i-1],p)+D[i]-c;
            if(j+c<=k)Min(dp[i+1][j+c][Arr],dp[i][j][p]+Arr*cnt[i]);
        }
        int Res=1e9;
        For(j,0,k)For(p,0,1000)Min(Res,dp[n+1][j][p]);
        ptn(Res+Ans);
    }
}P1;

• 然后考虑怎么贪心。
• 发现每次我们只要知道贡献最大的位置放置一个加速器即可
• 我们需要完成的任务就是高效的计算出，一个$D[i]$减少后答案的变化
• 对于一个$D[i]$的减少会有两个贡献
  
  • $arr[i]$的减少。
  • 会有的$\sum cnt[j](j<limt),limt$表示最近的无法被影响的点
  • 然后就是高效的对每个点求一个$limt$。。。
• 然后认为很烦，想了很多乱七八糟的东西。。。™都堆到了这样还没有$Ac$

下面是比赛的时候想到的没有屁用的东西

Di最多可以减少min(last[a]-arr[a],Di)次

然后我们需要的就是

1.二分查询limt

2.区间最值的查询last[a]-arr[a]

3.区间更新arr[a]

写一个线段树，加一个$BIT$哈哈哈！
复杂度为$O(n*log(n)^2*k)$只有60分，不划算。。

• 然后没有搞清$arr[a]$，什么乱七八糟的。。。很烦。。。不想这样写。。
• 然后干脆打了一个$O(n^2*k)$的暴力，没有继续想了。。

思路更正

• ™的。。。没有发现吗，对$limt$的只要做到$O(1)$查询，$O(n)$维护就可以了！！
• 因为每次你需要枚举所有的点，复杂度也为$O(n)*O(1)+O(n)$！！！！
• 然后随便维护一下就可以了。。
• 我们用$nxt[i]$表示从$i$出发可以影响到的点。
  
  • $nxt[i]=nxt[i+1] ~~(arr[i]>last[i])$
  • $nxt[i]=i ~~~~~~~~~~~~~~(arr[i]<=last[i])$
  • 然后一个$D[i]$减少的贡献就是$\sum_{j=i}^{nxt[i]}cnt[j]$了，前缀和维护一下即可

  • struct P3{
        int nxt[M],arr[M],sum[M];
    
        ll calc(){
            ll res=0;
            arr[1]=0;
            For(x,2,n)arr[x]=max(arr[x-1],last[x-1])+D[x];
            For(x,1,n)res+=(ll)arr[x]*cnt[x];
            return res;
        }
    
        void init(){
            arr[1]=0;
            For(x,2,n)arr[x]=max(arr[x-1],last[x-1])+D[x];
            nxt[n+1]=n;
            Ror(i,2,n){
                if(arr[i]>last[i])nxt[i]=nxt[i+1];
                else nxt[i]=i;
            }
        }
    
        void work(){
            For(i,1,n)sum[i]=cnt[i]+sum[i-1];
            while(k--){
                init();
                int mx=0,id=0;
                For(i,2,n)if(D[i])if(mx<sum[nxt[i]]-sum[i-1])id=i,mx=sum[nxt[i]]-sum[i-1];
                D[id]--;
            }
            ptn(calc()+Ans);
        }
    }P3;