【GDOI2014模拟】服务器

前言

直到比赛最后几分钟，才发现60%数据居然是一个水dp，结果没打完。

题目

我们需要将一个文件复制到n个服务器上，这些服务器的编号为S1, S2, …, Sn。
首先，我们可以选择一些服务器，直接把文件复制到它们中；将文件复制到服务器Si上，需要花费ci > 0的置放费用。对于没有直接被复制文件的服务器Si来说，它依次向后检查Si+1, Si+2, …直到找到一台服务器Sj：Sj中的文件是通过直接复制得到的，于是Si从Sj处间接复制得到该文件，这种复制方式的读取费用是j – i（注意j>i）。另外，Sn中的文件必须是通过直接复制得到的，因为它不可能间接的通过别的服务器进行复制。我们设计一种复制方案，即对每一台服务器确定它是通过直接还是间接的方式进行复制（Sn只能通过直接方式），最终使每一台服务器都得到文件，且总花费最小。

分析

60%的数据

水dp，设\(f[i]\)表示直接复制文件给\(i\)，把文件复制到i~n中的所有服务器最小的花费。
显然，转移为\(f[i]=c[i]+\min(f[j]+(1+j-(i+1))*(j-(i+1))/2)\).

100%的数据

现在就要想办法给dp加斜率优化了。
有两个位置\(j和k(j!=k 并且 i<j、k)\)，要使选\(j\)比\(k\)优，
那么$$\dfrac{f[j]+(1+j-(i+1))(j-(i+1))}{2}<\dfrac{f[k]+(1+k-(k+1))(j-(k+1))}{2}$$
移项得$$\dfrac{f[j]−f[k]+(j^{2}−k+k−j)/2}{j-k}<i$$
接着就可以打斜率优化dp了。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const long long maxlongint=214748364700;
using namespace std;
long long f[200000],n,m,c[200000],d[200000];
long long sum(long x,long y)
{
    return y*(y-x)-(y-x)*(x+y-1)/2;
}
double slope(long long x,long long y)
{
	return (f[x]-f[y]+(x*x-y*y+y-x)/2)*1.0/(x-y);
}
int main()
{
	scanf("%lld",&n);
	for(long long i=1;i<=n;i++)
	{
		scanf("%lld",&c[i]);
	}
	f[n]=c[n];
	long long l,r;
	d[l=r=1]=n;
	for(long long i=n-1;i>=1;i--)
	{
		while(l<r && slope(d[l],d[l+1])>=i) l++;
		f[i]=c[i]+f[d[l]]+sum(i+1,d[l]);
		while(l<r && slope(d[r-1],d[r])<slope(d[r],i)) r--;
		d[++r]=i;
	}
	long long ans=maxlongint;
	for(long long i=1;i<=n;i++)
		ans=min(ans,f[i]+sum(1,i));
	printf("%lld",ans);
}