（纪中）3528. 图书馆【DP】【数学】【方差】

---

解题思路

首先我们知道方差的公式是：
$K=({\sum }_{i=1}^m(x_i-p{)}^2)/m$

把 $(x_i-p{)}^2$ 拆开，可以得到
$K=({\sum }_{i=1}^m(x_i^2-x_i\ast p-x_i\ast p+p^2))/m$也就是

$K=({\sum }_{i=1}^m(x_i^2-2\ast x_i\ast p+p^2))/m$
再把 ∑拆开，就可以得到

$K={\sum }_{i=1}^m{x}_i^2/m-2\ast ({\sum }_{i=1}^m{x}_i)\ast p/m+p^2$

令 $A={\sum }_{i=1}^m{x}_i^2,B={\sum }_{i=1}^m{x}_i$
又知道 $p=\frac{{\sum }_{i=1}^m{x}_i}{m}$ ,即 $p=\frac{B}{m}$

再把这个式子带入，得

$K=A/m-2\frac{B^2}{m^2}+(\frac{B}{m}{)}^2$

再化简，得

$K=A/m-2\frac{B^2}{m^2}+\frac{B^2}{m^2}$

即$K=A/m-\frac{B^2}{m^2}$

---

死了死了，打到手费，来吧，回到此题：
设$f[i][j][k]$表示当前走过了 i 条边，现在在 j 这个平台，经过的路程为 k的上文中A，
那么很容易可以得到状态转移方程：
$f[i][j][k]=min(f[i][j][k],f[i-1][y][k-w[y,j]]+w[y,j{]}^2)$
我们需要枚举y进行转移。

---

代码

#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
#include<algorithm>
#include<iomanip>
#include<cmath>
using namespace std;

const int maxn=2147483600;
int n,m,x,y,z,k,head[610],f[25][310][1020];
double ans;

struct c{
	int x,next,w;
}a[610];

void add(int x,int y,int w){
	a[++k].x=y;
	a[k].next=head[x];
	a[k].w=w;
	head[x]=k;
}

int main(){
	freopen("library.in","r",stdin);
	freopen("library.out","w",stdout);
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		add(y,x,z);//因为后面要枚举一个所有能走到下文中j的点，所以把每个点与能走向他的点连边
	}
	for(int i=0;i<20;i++)
		for(int j=1;j<=n;j++)
			for(int k=0;k<=1010;k++)
				f[i][j][k]=maxn;
	f[0][1][0]=0;
	for(int i=1;i<20;i++)
	{
		for(int j=2;j<=n;j++)
		{
			for(int k=0;k<=1010;k++)
			{
				for(int t=head[j];t;t=a[t].next)
				{
					int y=a[t].x;
					if(k>=a[t].w&&f[i-1][y][k-a[t].w]<maxn) 
						f[i][j][k]=min(f[i][j][k],f[i-1][y][k-a[t].w]+a[t].w*a[t].w);
				}
			}
		}
	}
	ans=maxn;
	for(int i=1;i<20;i++)
	{
		for(int k=0;k<=1010;k++)
			if(f[i][n][k]!=maxn)
				ans=min(ans,abs(f[i][n][k]*1.0/(i*1.0)-((k*k)*1.0)/(i*1.0)/(i*1.0)));//套个方差公式
	}
	printf("%.4lf",ans);
}