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

Description

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

Input

输入的第一行有两个整数，n和maxc。
之后的n行，每行两个整数，分别表示0到n-1号车站的c和v.

Output

输出一行n个整数，其中第i个整数代表从0号车站到i号车站的最小花费，若不能从0号车站到达i号车站，则在i号车站的位置输出-1。

Sample Input

输入1：
1 1
1 1
输入2：
9 5
2 5
5 2
5 14
1 18
4 13
3 17
1 16
1 7
5 4

Sample Output

输出1：
1
输出2：
-1 5 -1 10 -1 15 19 20 33

Data Constraint

对于30%的数据满足，1<=n<=5000
对于60%的数据满足，1<=n<=10^5
对于另20%的数据满足，maxc=1
对于100%的数据满足，1<=n<=10^6,1<=maxc<=10,1<=ci<=maxc,1<=vi<=1000
数据存在梯度。

Hint

样例1说明：从0号车站坐1站地，到达1号车站，花费为1，可以发现这是从0号车站到1号车站的最小花费。

题解

用$f_i$ 表示到达点i的最小值。
我们可以发现，每一辆公交车的收费都是一个一次函数，斜率是$v_i$，截距与$f_i$ 有关。
对于所有的公交车还可以按照出发点i%$c_i$ 和间隔$c_i$ 进行分类，因为他们经过的路线都是相同的。这样就可以将公交车分为55类。

对于一个点i的转移，它可以选择1~maxc这些间隔的公交车。
枚举一个间隔j，可以到达i的是i%j，j这一类公交车。

考虑如何决策最优，对于每一类的公交车，维护一个栈。
栈里面存着这一类型的公交车，
而它们的票价都是一条直线，就维护凸包。
首先，斜率一定是单调递增的，因为对于一条新插入的直线i，$f_i$ 一定是小于等于当前栈顶转移到i的费用，如果栈顶的斜率更大，显然没有新加入的直线i优。
第二，两条相邻的直线交点也要单调的。因为两条直线的交点，就是它们的谁更优分界点。
第三，当处理的点已经超过了这两条直线的交点，分界点之后，前面的那一条直线就不是最优的，就应该把它弹出栈。

code

#include<queue>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include <cstring>
#include <string.h>
#include <cmath>
#include <math.h>
#define ll long long
#define N 1000333
#define db double
#define P putchar
#define G getchar
#define mo 998244353
using namespace std;
char ch;
void read(int &n)
{
    n=0;
    ch=G();
    while((ch<'0' || ch>'9') && ch!='-')ch=G();
    ll w=1;
    if(ch=='-')w=-1,ch=G();
    while('0'<=ch && ch<='9')n=(n<<3)+(n<<1)+ch-'0',ch=G();
    n*=w;
}

int max(int a,int b){return a>b?a:b;}
int min(int a,int b){return a<b?a:b;}
ll abs(ll x){return x<0?-x:x;}
ll sqr(ll x){return x*x;}
void write(ll x){if(x>9) write(x/10);P(x%10+'0');}
void writeln(ll x){write(x);P('\n');}

int n,m,c[N],v[N],f[N],x,y;
int d[10][11][N],top[10][11];
db g[N];

bool pd(int x1,int x2,int x3)
{
    db t1=(g[x2]-g[x1])/(v[x1]-v[x2]);
    db t2=(g[x3]-g[x1])/(v[x1]-v[x3]);
    return t1>=t2;
}

int S(int x,int y)
{
    return f[y]+(x-y)/c[y]*v[y];
}

void ins(int x,int y,int z)
{
    while((top[x][y] && v[d[x][y][top[x][y]]]>=v[z])||
    (top[x][y]>1 && pd(z,d[x][y][top[x][y]],d[x][y][top[x][y]-1])))top[x][y]--;
    d[x][y][++top[x][y]]=z;
}

void del(int x,int y,int z)
{
    while(top[x][y]>1 && S(z,d[x][y][top[x][y]])>=S(z,d[x][y][top[x][y]-1]))top[x][y]--;
}

int main()
{
    freopen("bus.in","r",stdin);
    freopen("bus.out","w",stdout);
    memset(f,127,sizeof(f));f[0]=0;
    read(n);read(m);
    for(int i=0;i<n;i++)
        read(c[i]),read(v[i]);
    ins(0,c[0],0);

    for(int i=1;i<=n;i++)
    {
        if(i==78)
            x=1;
        for(int j=1;j<=m;j++)
        {
            x=i%j;y=j;
            del(x,y,i);
            if(top[x][y])f[i]=min(f[i],S(i,d[x][y][top[x][y]]));
        }
        g[i]=f[i]-1.00*i/c[i]*v[i]; 
        if(i!=n && f[i]<=1000000000)del(i%c[i],c[i],i),ins(i%c[i],c[i],i);
        if(f[i]>1000000000)P('-'),P('1');else write(f[i]);
        P(' ');
    }
}