【JZOJ 5415】 公交运输

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

Analysis

考虑dp
$f[i]=min(f[j]+(i-j)*v[j]/c[j]),(j\mod c[j]=i\mod c[j])$
为了把后面的转移限制去掉，我们将决策i扔进i%ci这个队列，那么只有10个队列
考虑到maxc只有10，我们可以枚举c[j]来转移
那么转移式形如$f[i]=min(f[j]+(i-j)*v[j])$
这是一次函数的形式，可以用线段树维护直线带log碾过去，但是有O(nmxc)做法

凸壳

没错，只是这次维护的凸壳长得很奇怪
将所有的决策的作用画成直线，可以发现我们只需要维护一个上凸壳，具体而言就是维护一个斜率递增，且 直线之间的交点不断左移的单调栈
综合上面，我们可以开10*10个单调栈，里面存ci，i%ci都相同的决策
注意弹栈的顺序一定是先保证斜率递增，再保证交点左移

Code

#include<vector>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,b,a) for(int i=b;i>=a;i--)
#define efo(i,v,u) for(int i=last[v],u=to[i];i;i=next[i],u=to[i])
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define mset(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
typedef double db;
char ch;
void read(int &n)
{
    n=0;int p=1;
    for(ch=getchar();ch<'0' || ch>'9';ch=getchar())
        if(ch=='-') p=-1;
    for(;'0'<=ch && ch<='9';ch=getchar()) n=n*10+ch-'0';
    n*=p;
}
const int N=1e6+5,INF=2139062143;
int n,mxc,c[N],w[N],f[N],top[11][11],pre[N];
db F(int j,int i)
{
    if(j==n+1) j=0;
    return f[j]+1.0*(i-j)*w[j]/c[j];
}
db crossx(int j,int k)
{
    if(j==n+1) j=0;if(k==n+1) k=0;
    db dj=1.0*w[j]/c[j],dk=1.0*w[k]/c[k];
    return 1.0*(dj*j-dk*k+f[k]-f[j])/(dj-dk);
}
int main()
{
    freopen("bus.in","r",stdin);
    freopen("bus.out","w",stdout);
    read(n),read(mxc);
    fo(i,0,n-1) read(c[i]),read(w[i]);
    mset(f,127);f[0]=0;
    top[0][c[0]]=n+1,w[n+1]=w[0],c[n+1]=c[0];
    fo(i,1,n)
    {
        fo(j,1,10) if(i>=j)
        {
            int ic=i%j;
            while(pre[top[ic][j]] && F(pre[top[ic][j]],i)<F(top[ic][j],i)) top[ic][j]=pre[top[ic][j]];
            if(top[ic][j]) f[i]=min(f[i],F(top[ic][j],i));
        }
        if(i==n || f[i]==INF) continue;
        int ic=i%c[i],tp=top[ic][c[i]];
        while(tp && w[tp]>=w[i]) tp=pre[tp];
        while(pre[tp] && crossx(pre[tp],tp)<=crossx(tp,i)) tp=pre[tp];
        pre[i]=tp,top[ic][c[i]]=i;
    }
    fo(i,1,n)
        if(f[i]==INF) printf("-1 ");
        else printf("%d ",f[i]);
    return 0;
}