A1251. 设计铁路(钱桥) DP+斜率优化

最新推荐文章于 2019-06-03 13:52:02 发布

ws_fqk

最新推荐文章于 2019-06-03 13:52:02 发布

阅读量667

点赞数

CC 4.0 BY-SA版权

分类专栏： My Code

本文链接：https://blog.csdn.net/Phenix_2015/article/details/51227814

My Code 专栏收录该内容

349 篇文章

订阅专栏

明明是一道裸题，我却做了这么久。感觉自己斜率优化真是弱啊，省选就没做出来，现在做还是这么费力。
我们用 $d_i,a_i$ 分别表示第 $i$ 个车站到终点的距离和人数。
首先按照 $d_i$ 降序排序，然后我们用 $f_i$ 表示在第i个村庄出修建一个车站，且 $1~i$ 的所有人都开车到车站，此时的总花费。
然后我们维护两个前缀和 $s_i,p_i$ 。

s i = \sum j = i + 1 n d j * a j

$s_i=\sum_{j=i+1}^n d_j*a_j$

p i = \sum j = i + 1 n a j

$p_i=\sum_{j=i+1}^n a_j$
那么状态转移方程如下：

f i = m i n {f j + m + s j - s i - (p j - p i) \times d i}

$f_i=min\{f_j+m+s_j-s_i-(p_j-p_i)\times d_i\}$
把与

j $j$ 无关的量提到外面：

f i = m i n {f j + s j - p j \times d i} + m + p i \times d i - s i

$f_i=min\{f_j+s_j-p_j\times d_i\}+m+p_i\times d_i-s_i$
那么考虑若存在

j<k $j<k$ ，且

k $k$ 比

j $j$ 更优，那么有

f k + s k - p k \times d i < f j + s j - p j \times d i

$f_k+s_k-p_k\times d_i<f_j+s_j-p_j\times d_i$
化简可得

( f k + s k ) - ( f j + s j ) p k - p j > d i

$\frac{(f_k+s_k)-(f_j+s_j)}{p_k-p_j}>d_i$
然后就可以斜率优化了。

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;

const int N=40005;
int n,m;
struct node
{
    int d,r;
}p[N];
long long s[N],g[N],f[N];
int q[N];

inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}

inline bool operator<(node a,node b)
{
    return a.d>b.d;
}

inline long long Y(int j,int k)
{
    return f[j]+s[j]-f[k]-s[k];
}

inline long long X(int j,int k)
{
    return g[j]-g[k];
}

int main()
{
    n=read(); m=read();
    for (int i=1;i<=n;i++)
        p[i].d=read(),p[i].r=read();
    sort(p+1,p+n+1);
    for (int i=n-1;~i;i--)
        s[i]=s[i+1]+p[i+1].d*p[i+1].r,g[i]=g[i+1]+p[i+1].r;
    int t=0,w=0;
    for (int i=1;i<=n+1;i++)
    {
        while (t<w&&Y(q[t+1],q[t])<X(q[t+1],q[t])*p[i].d) t++;
        f[i]=f[q[t]]+s[q[t]]-s[i]+(g[i]-g[q[t]])*p[i].d+m;
        while (t<w&&Y(i,q[w])*X(q[w],q[w-1])>Y(q[w],q[w-1])*X(i,q[w])) w--;
        q[++w]=i;
    }
    cout << f[n+1]-m << endl;
    return 0;
}