HDU5492——合肥1009

题目意思是找一条方差最小的路径
记得以前做过最小方差生成树，这题要解最小方差路径也是差不多的。方差与期望相关，而期望和选择的路径相关所以没法直接用期望来求，我用的方法和最小方差生成树的方法一样，就是枚举所有可能的期望值来计算最小的方差。
可以证明，当枚举的期望使得方差最小时，这条路径一定是方差最小的路径，然而我并不会证明，我只是记住了结论。
S=∑N1(ai−x¯)2<∑N1(ai−x)2(x≠x¯)

所以做这个题目的时候先走两遍DP把最大路径和最小路径求出来，然后从最小路径到最大路径枚举平均值，每次的间隔取（max+min）/(m+n-1),每次枚举按照枚举的EX来生成一个新的方差图，找出所有方差图中最小的，就是最小方差路径。

————————————————————————————
今天还看到一种用E(x^2)-(E(x))^2来DP的方法http://blog.csdn.net/Baileys0530/article/details/48768123
感觉都能做

#include <iostream>
#include <cstring>
using namespace std;
double f[35][35],a[35][35],b[35][35],ans,Max,Min,Exp,ave;
int i,j,n,m;
int main()
{
    int c=0,T;
    cin>>T;
    while (T--)
{
    c++;
    cin>>n>>m;
    Max=-1;
    Min=99999999;
    memset(f,0,sizeof(f));
    memset(a,0,sizeof(f));


    for (i=1; i<=n; i++)
        for (j=1; j<=m; j++)
            {
                cin>>a[i][j];
                f[i][j]=a[i][j];
            }
     for (i=1; i<=n; i++)
        for (j=1; j<=m; j++)
            if (f[i-1][j]>f[i][j-1])
                f[i][j]=f[i][j]+f[i-1][j];
            else f[i][j]=f[i][j]+f[i][j-1];

    Max=(f[n][m])/(m+n-1);

    for (i=1; i<=n; i++)
        for (j=1; j<=m; j++)
            f[i][j]=a[i][j];
   for (i=1;i<=n;i++) f[i][0]=99999999;
   for (j=1;j<=m;j++) f[0][j]=99999999;
    for (i=1; i<=n; i++)
        for (j=1; j<=m; j++)
            {
                if (f[i-1][j]<f[i][j-1])
                f[i][j]=f[i][j]+f[i-1][j];
            else f[i][j]=f[i][j]+f[i][j-1];
            f[1][1]=a[1][1];
            }
    Min=f[n][m]/(m+n-1);
    // cout<<Max<<" ..."<<Min<<endl;
    Exp=1/(float)(m+n-1);
   // cout<<exp<<endl;
    ans=99999999;
    ave=Min;
    while (ave<=Max+0.00000001)
    //for (ave=Min; ave<=Max; ave=ave+exp)
    {
        double temp;
        for (i=1; i<=n; i++)
            for (j=1; j<=m; j++)
                f[i][j]=(a[i][j]-ave)*(a[i][j]-ave);
        temp=f[1][1];
        for (i=1; i<=n; i++)
        for (j=1; j<=m; j++)
            {
                if (f[i-1][j]<f[i][j-1])
                f[i][j]=f[i][j]+f[i-1][j];
            else f[i][j]=f[i][j]+f[i][j-1];
            f[1][1]=temp;
            }
        if (f[n][m]<ans) ans=f[n][m];
       // cout<<f[n][n]<<"--"<<ave<<endl;
        ave=ave+Exp;
       // cout<<ave<<endl;
    }
    cout<<"Case #"<<c<<": "<<ans*(n+m-1)<<endl;
   }
   return 0;
}