【gmoj】 【NOIP2013模拟11.7A组】 【DP】【数学】图书馆

【gmoj】 【NOIP2013模拟11.7A组】 【DP】【数学】图书馆

题目

---

解题思路

雾~ 求分差是什么鬼

由初二知识可得
方差公式长这个亚子

将方差先乘上n 化开 变成这样

因为平均数等于ai的和/n 那么 ai的和=n*平均数
最后再合并一下就变成了这个亚子
s²=∑(a[i])²-(x)²

设f[i][j][k]为走了i步到达第j个平台用了k的力气
里面存的是（a[i]）²

---

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
struct lzf{
	int to,next,q;
}f[320];
int n,m,x,y,z,t,head[52],d[21][52][1001];
double ans=1000000;
void add(int x,int y,int q)
{
	 f[++t].q=q;
	 f[t].to=y;
	 f[t].next=head[x];
	 head[x]=t;
}
int main()
{
	freopen("library.in","r",stdin);
	freopen("library.out","w",stdout);
	memset(d,0x7f,sizeof(d));
	scanf("%d%d",&n,&m);
	for (int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		add(y,x,z);
	}
	d[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 z=head[j];z;z=f[z].next)
	                 if (f[z].q<=k&&d[i-1][f[z].to][k-f[z].q]!=d[0][0][0])
	                 {
	                 	d[i][j][k]=min(d[i][j][k],d[i-1][f[z].to][k-f[z].q]+f[z].q*f[z].q);
					 }
	 for (int i=1;i<20;i++)
	     for (int j=0;j<=1000;j++)
	     if (d[i][n][j]!=d[0][0][0])
	     {
	     	 double da=(1.0*d[i][n][j]/i)-(1.0*j*j/i/i);
	     	 if (abs(da)<ans) ans=abs(da);
		 }
	  printf("%.4lf",ans);
	  fclose(stdin);
	  fclose(stdout);
	  return 0;
}