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

题意不说了。

有点鬼畜，比赛只会50分，花了30分钟思考一无所获，连按模分组都没想到迟早退役。
然后花了十分钟打50分扔了这题。巨亏。
其实比赛还剩一个小时的时候发现是斜率优化，但是没时间打所以就放了。。
事实上这并不是简单的斜率优化，，因为如果直接做不能保证他的斜率是单调递增的。
首先先把点按模分组，只有%c相同的点才能被到达。
然后就可以直接斜率优化了，每次先枚举c，对于栈中的点，我们把栈顶和前一个取出来，如果栈顶到i的答案不如前一个到i的答案优，那肯定要把栈顶踢出。。据说只做这个优化能拿80。。
然后就可以直接单调递增维护斜率了。
注意实现的话最好用链表会比较方便。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<cstring>
#include<queue>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e6+5;
typedef double db;
int c[N],v[N],n,maxc,cnt,q[N];
typedef long long ll;
int f[N],sta[N],h[15][15],t[15][15],last[N],next[N];
int top,inf=1e9;
db K[N];
inline int get(int j,int i)
{
    return (i-j)/c[j]*v[j]+f[j];
}
double getk(int j,int k)
{
    return (double)1.0*(k/c[k]*v[k]-j/c[j]*v[j]+f[j]-f[k])/(v[k]-v[j]);
}
int main()
{
    freopen("bus.in","r",stdin);
    freopen("bus.out","w",stdout);
    scanf("%d%d",&n,&maxc);
    fo(i,1,n)
    {
        scanf("%d%d",&c[i],&v[i]);
        f[i+1]=inf;
    }
    h[c[1]][1%c[1]]=t[c[1]][1%c[1]]=1;
    fo(i,2,n+1)
    {
        fo(j,1,maxc)
        {
            int x=i%j;
            if (!h[j][x])continue;
            while (h[j][x]!=t[j][x]&&t[j][x]&&get(t[j][x],i)>get(last[t[j][x]],i))
            t[j][x]=last[t[j][x]],next[t[j][x]]=0;
            f[i]=min(f[i],get(t[j][x],i));
        }
        if(i==n+1||f[i]==inf)continue;
        int x=c[i],y=i%x;
        while (t[x][y]&&v[t[x][y]]>=v[i])
            t[x][y]=last[t[x][y]],next[t[x][y]]=0;
            if (!t[x][y]) h[x][y]=0;
        while (h[x][y]!=t[x][y]&&t[x][y]&&getk(last[t[x][y]],t[x][y])<getk(t[x][y],i))
        t[x][y]=last[t[x][y]],next[t[x][y]]=0;
        if (t[x][y])next[t[x][y]]=i,last[i]=t[x][y];
        if (!h[x][y])h[x][y]=i;
        t[x][y]=i;
    }
    fo(i,2,n+1)
    if (f[i]<inf)printf("%d ",f[i]);
    else printf("%d ",-1);
    return 0;
}