#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; }