【NOIP2013模拟11.7A组】图书馆题解

Description

圣玛格丽特大图书馆是一座由石材砌成的角柱型高塔，是欧洲屈指可数的巨大书库。图书馆整面墙壁都是巨大的书架，书架与书架之间就像巨大的迷宫一般，以细窄的木制楼梯连结。大图书馆的最高处是一个绿意盎然的植物园，维多利加正在那无聊地看着书。今天，一如往常地，久城要爬上这迷宫般的楼梯给维多利加送讲义。

图书馆墙壁上有N个平台，编号为1到N，入口为1号，植物园为N号。有M个连接两个不同平台的楼梯，爬每个楼梯需要消耗一定的体力值。楼梯一定是由低处通往高处的，为了省时间，久城只能选择上楼梯而不能下楼梯，也就是说，楼梯之间不会形成环路。而且，出于人性化考虑，不管久城选择哪条路线上楼，他爬的楼梯数量一定小于20。

为了使体力消耗尽量平稳，久城需要选择一条“每个楼梯消耗体力值的方差最小”的路径上楼。请帮助久城计算出这个最小方差。

Solution

首先，让我们来推导方差的公式。
大家如果不会看我的另外一篇博客：方差公式化简
化简完之后，我们发现我们只需要记录路程总和还有经过的边数，所以我们可以设一个dp。
设$f[i][j][k]$表示当前走过了$i$条边，现在在$j$这个平台，经过的路程为$k$的最小平方和。
那么很容易可以得到状态转移方程：
$f[i][j][k]=min(f[i][j][k],f[i-1][v][k-map[v,j]]+map[v,j{]}^2)$
我们只需要枚举v进行转移，时间复杂度是$O(200000m)$

Code

#include<cstdio>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;
int n,m,cnt;
int head[51],f[20][51][1001];
double ans=1000000.0;
struct node {
	int to,len,next;
}edge[301];
inline int read() {
	int s=0;
	char ch=getchar();
	while(ch<48||ch>57)ch=getchar();
	while(ch>=48&&ch<=57)
		s=(s<<1)+(s<<3)+(ch^48),ch=getchar();
	return s;
}
inline void add(int x,int y,int z) {
	edge[cnt]=(node){y,z,head[x]};
	head[x]=cnt++;
}
inline int min(int x,int y) {return x<y?x:y;}
int main() {
	freopen("library.in","r",stdin);
	freopen("library.out","w",stdout);
	n=read();m=read();
	memset(head,-1,sizeof(head));
	for(int i=1;i<=m;++i) {
		int u=read(),v=read(),w=read();
		add(v,u,w);
	}
	for(int i=0;i<20;++i)
		for(int j=1;j<=n;++j)
			for(int k=0;k<=1000;++k)
				f[i][j][k]=1e5;
	f[0][1][0]=0;
	for(int i=1;i<20;++i)
		for(int j=2;j<=n;++j)
			for(int k=0;k<=1000;++k) {
				for(int t=head[j];t!=-1;t=edge[t].next) {
					int to=edge[t].to,len=edge[t].len;
					if(len>k)continue;
					if(f[i-1][to][k-len]<1e5)
						f[i][j][k]=min(f[i][j][k],f[i-1][to][k-len]+len*len);
				}
			}
	for(int i=1;i<20;++i)
		for(int k=0;k<=1000;++k) {
			if(f[i][n][k]==1e5)continue;
			double t=((1.0*f[i][n][k]-(1.0*(k*k))/(1.0*i))/(1.0*i));
			if(abs(t)<ans)ans=abs(t);
		}
	printf("%0.4lf",ans);
	fclose(stdin);
	fclose(stdout);
	return 0;
}