http://poj.org/problem?id=2096
G++ WA C++ AC
题解看这个吧:我没看懂题目,所以找的题解翻译:http://blog.csdn.net/morgan_xww/article/details/6774708
可是题解里感觉有一点写的不好---初始化应该有两部分:1、dp[s][n]=0.0。以及dp[s+1][j]=0.0 上面那个博客里没有初始化dp[s+1][j]=0.0 由于全局变量
我的dp[i][j]意思是,前i个子系统。前j种bug的天数期望。
事实上这道概率DP入门题让我学到的是一种高中学概率无用的的思维方式----由后一阶段推前一阶段的概率/期望
#include <cstdio>
#include <iostream>
using namespace std;
#define rep(i,s,e) for(int i=s;i<e;i++)
#define reped(i,s,e) for(int i=s;i>=e;i--)
const int MAXN = 1000+100;
int n,s;
double dp[MAXN][MAXN];
double solve()
{
dp[s][n]=0.0;
rep(i,0,n+1)dp[s+1][i]=0.0;
reped(i,s,0)
{
reped(j,n,0)
{
if(i==s && j==n)continue;
dp[i][j]=(n*s+(s-i)*j*dp[i+1][j]+(s-i)*(n-j)*dp[i+1][j+1]+i*(n-j)*dp[i][j+1])/(n*s-i*j);
}
}
return dp[0][0];
}
int main()
{
while(~scanf("%d%d",&n,&s))
{
printf("%.4lf\n",solve());
}
return 0;
}