[JZOJ5415]【NOIP2017提高A组集训10.22】公交运输

Description

城市中有一条长度为n的道路，每隔1的长度有一个公交车站，编号从0到n，学校在0号车站的位置。其中每个公交车站(除了n号车站)有两个属性ci和vi，代表从这个公交车站出发的公交车的性质。ci代表这个从i出发的公交车，相邻两个停靠站之间的距离。vi表示每坐1站的花费。
注意，一辆公交车出发后会向n号车站的方向行进。同时，一名乘客只能从起点站上车，但可以从任意停靠站下车。校庆志愿者小Z为了帮助校友查询有关城市交通费用的问题，想知道从0号车站(也就是学校)出发，到达每个公交车站的最小花费，于是他找到了你。
对于100%的数据满足，1<=n<=10^6,1<=maxc<=10,1<=ci<=maxc,1<=vi<=1000

Solution

可以按c和%c来分组，分成100组，找的时候枚举它被转移过来的c，在对应的组中找

考虑优化

发现转移式是一个二次函数的形式。

那么在某个组插入相当于插入一条直线，而我们需要维护凸包

线段树维护。

然而还有更优的做法。

可以朝着斜率优化的方向想

但是尴尬的问题在于插入的$v_i$并不单调，导致没有办法斜率优化

但实际上，假如在同一组内有$i,j$,$i<j,v_i>=v_j$

那i是没有什么用的
因为在j以后i和j能转移到的集合相同，且i一定能转移到j
后面就可以直接用j转移了

所以在组内只需要存下v单调递增的就可以了，类似于一个单调栈

在这个单调栈的基础上再进行斜率优化，最优答案就是栈顶，发现它不优了就弹栈

Code

#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define N 1000005
using namespace std;
int n,m,f[N],t[N][2],lt[10*N],ls[11][11];
int get(int x,int w)
{
    return f[x]+(w-x)/t[x][0]*t[x][1];
}
bool pd(int x,int y,int w)
{
    return get(x,w)<get(y,w);
}
double xl(int x,int y)
{
    double qx=t[x][1]*1.0/t[x][0],qy=t[y][1]*1.0/t[y][0];
    return (f[x]-qx*x-f[y]+qy*y)/(qy-qx);
}
int main()
{
    cin>>n>>m;
    fo(i,0,n-1) scanf("%d%d",&t[i][0],&t[i][1]);
    memset(f,107,sizeof(f));
    f[0]=0;
    memset(ls,255,sizeof(ls));
    memset(lt,255,sizeof(lt));
    fo(i,0,n)
    {
        if(i!=0)
        {
            fo(j,1,m)
            {
                int q=i%j;
                if(ls[j][q]==-1) continue;
                while(lt[ls[j][q]]!=-1&&pd(lt[ls[j][q]],ls[j][q],i)) ls[j][q]=lt[ls[j][q]];
                if(lt[ls[j][q]]==-1) ls[j][q]=ls[j][q];
                f[i]=min(f[i],get(ls[j][q],i));
            }
            printf("%d ",(f[i]==1802201963)?-1:f[i]);
        }
        if(i!=n)
        {
            if(f[i]==1802201963) continue;
            int q=i%t[i][0],j=t[i][0],k=ls[j][q];
            while(k!=-1&&t[k][1]>=t[i][1]) k=lt[k];
            while(k!=-1&&lt[k]!=-1&&xl(k,i)>=xl(lt[k],k)) k=lt[k];
            lt[i]=k;
            ls[j][q]=i;
        }
    }
}