【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

Solution

看这题的转移比较像斜率
但是直接斜率又会出错
那么分组
按照他跳的步数，能跳到的分在一组
那么有100组
做100个单调栈

Code

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define N 1000100
#define INF 2000000000
#define db double
using namespace std;
int n,mc,a[11][11],d[11][11][2010],k[N],b[N],tot=0;
int xj(int i,int j,int y,int z)
{
    return (b[i]-b[j])*1.0/(k[j]-k[i])<z/y;
}
int pd(int k1,int b1,int i,int j)
{
    return (b1-b[i])*(k[j]-k[i])>(b[i]-b[j])*(k[i]-k1);
}
int yi(int x,int i)
{
    return x*k[i]+b[i];
}
void write(int x)
{
    if(x==INF||x==0) printf("-1 ");else printf("%d ",x);
}
int main()
{
    freopen("bus.in","r",stdin);
    freopen("bus.out","w",stdout);
    scanf("%d%d",&n,&mc);
    fo(i,0,n-1)
    {
        int c,v,ans=INF;scanf("%d%d",&c,&v);
        if(i)
        fo(j,1,mc)
        {
            int q=i%j,w=j;
            while(a[q][w]>1&&xj(d[q][w][a[q][w]],d[q][w][a[q][w]-1],j,i)) a[q][w]--;
            if(a[q][w]) ans=min(ans,yi(i/j,d[q][w][a[q][w]]));
        }
        int q=i%c,w=c;
        if(i) write(ans);
        if(i&&ans==INF||ans==0) continue;
        if(!i) ans=0;
        int k1=v,b1=ans-v*(i/c);
        while(a[q][w]>0&&(k1<k[d[q][w][a[q][w]]]||k1==k[d[q][w][a[q][w]]]&&b1<=b[d[q][w][a[q][w]]])) a[q][w]--;
        while(a[q][w]>1&&pd(k1,b1,d[q][w][a[q][w]],d[q][w][a[q][w]-1])) a[q][w]--;
        k[++tot]=k1;b[tot]=b1;
        d[q][w][++a[q][w]]=tot;
    }
    int ans=INF;
    fo(j,1,mc)
    {
        int q=n%j,w=j;
        while(a[q][w]>1&&xj(d[q][w][a[q][w]],d[q][w][a[q][w]-1],j,n)) a[q][w]--;
        if(a[q][w]) ans=min(ans,yi(n/j,d[q][w][a[q][w]]));
    }
    write(ans);
    printf("\n");
}