ZJU 3582 期望DP

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#define LL long long
#define eps (1e-9)
using namespace std;

double e[55][55],p,pow1[200],pow2[200];
LL C[55][55];
bool vis[55][55];
int n,m;
void init()
{
    int i,j;
    C[0][0]=1;
    for(i=1;i<=50;i++)
    {
        C[i][0]=1;C[i][i]=1;
        for(j=1;j<i;j++) C[i][j]=C[i-1][j-1]+C[i-1][j];
    }
}
double E(int x,int y)
{
    if(x>=m && y>=m) return e[x][y]=0.0;
    if(vis[x][y]) return e[x][y];
    if(vis[y][x]) return e[y][x];
    int i,j;
    double ans=1.0;
    for(i=0;i+x<=n;i++) for(j=0;j+y<=n;j++) if(i || j) ans+=E(x+i,y+j)*C[n-x][i]*C[n-y][j]*pow1[i+j]*pow2[n-x-i+n-y-j];
    ans=ans/(1-pow2[n-x+n-y]);
    vis[x][y]=1;
    return e[x][y]=ans;
}

int main()
{
    init();
    int i,j;
    while(scanf("%d %d %lf",&n,&m,&p)==3)
    {
        if(n==0 && m==0 && fabs(p)<eps) return 0;
        pow1[0]=pow2[0]=1.0;
        for(i=1;i<200;i++)
        {
            pow1[i]=p*pow1[i-1];
            pow2[i]=(1-p)*pow2[i-1];
        }
        memset(vis,0,sizeof(vis));
        printf("%.7lf\n",E(0,0));
    }
    return 0;
}