##概率dp##

参考http://www.cnblogs.com/kuangbin/archive/2012/10/02/2710606.html

参考http://www.cnblogs.com/SolarWings/archive/2013/04/24/3040784.html

概率dp

两种求法

（以hdu4050为例子）

求法1：初始状态（dp[0]）概率为1，求dp[j]时会用到dp[i](j>i), ans=所有概率的累加（except初始的1）

求法2：终状态（dp[n]）期望为0，求dp[j]时会用到dp[i](j<i)，往往有+1,ans=dp[0]

PS:1、Dp[i]中转移式有本身的需要移项&除法

   2、如果不能递推，可能涉及到迭代法求根O(n) 或高斯消元求根O(n^3)

 

//有equ个方程，var个变元。增广矩阵行数为equ,分别为0到equ-1,列数为var+1,分别为a[][0到var-1]及x[].
void doit() //采用高斯消元模版2
{
    //step 1 读入,特判(如起点已是终点)及预处理
    //……
    //step 2 bfs搜索所有到得了的状态并重编号 判断目标状态是否可达
    memset(num,255,sizeof(num));
    equ=var=bfs(初状态); //num[o++]=new; return o;
    if (num[末状态]==-1) {printf("Impossible !\n"); return;}
    //step 3
    memset(x,0,sizeof(x));
    memset(a,0,sizeof(a));
    for (int i=0;i<var;i++)
        a[i][i]=1;
    for (int i=0;i<m;i++)//m为原始状态数
    if(num[i]!=-1)
        {
            if (i==0)
                {  x[num[i]]=0;
                    continue;
                }
            int kk; double pp;
            //对所有kk pp进行操作
            kk=……; pp=……;
            if (num[kk]!=-1) {
            a[num[i]][num[kk]]-=pp;
            x[num[i]]+=pp;}
        }
    if (!gauss())   {printf("INF\n"); return;}
    printf("%.2lf\n",x[num[初状态]]+eps);
}

Hdu 3853

方程略去。

除了右下角有没有一个点停留在原地的概率为1（这样除法就会跪），

如果有，答案应该是无穷大。

还有一种可能，虽然有，但这点根本不可达。（实践测试确实如此）

 

Zoj 3329

dp[i]=dp[0]*p0+∑dp[i+k]*P[k]+1;　

设dp[i]=fa[i]*dp[0]+fb[i];  //关键！维护系数，树状数组/线段树也有这样的题

dp[i]=(∑P[k]*fa[i+k]+p0)*dp[0]+∑P[k]*fb[i+k]+1;

得：fa[i]=∑P[K]*fa[i+k]+p0; (0<=i<=n)

　　fb[i]=∑P[K]*fb[i+k]+1;  (0<=i<=n)

ans=dp[0]=fb[0]/(1-fa[0])

 

 Hdu 4652不一定都需要向上题一样，结合高中# 数列 #知识灵活求dp[0]

 Hdu 4035树形概率dp维护系数

 Hdu 4326高斯消元

 POJ 3071  CF148D  ZOJ3380 ZOJ 3640的dp方程还是值得思考的