POJ 2096 概率dp

题意：给一个程序，甲每天能发现一个bug，bug分为n种，可以发生在s个位置，每次发现bug的位置是随机的且每次每种bug发生的概率都是1/n*s，问发现了n种和s种位置的bug都发现了的时候的概率

思路：啦啦啦~百度的题解~用dp[i][j]表示发现了i种bug，发生在j种位置时，到达完成是天数的期望。

那么dp[i][j]可以有四种转移方向，dp[i][j],dp[i+1][j],dp[i][j+1],dp[i+1][j+1]。

每种的概率也很好求，于是dp[i,j] = p1*dp[i,j] + p2*dp[i+1,j] + p3*dp[i,j+1] + p4*dp[i+1,j+1] + 1~

这里的+1是因为又过了一天，所以应该是发生四种事件的期望天数+1。于是整理出来dp[i][j]就可以递推了

代码：///随便换个语言看看是啥样的=。=

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
const int N = 130;
const int mod = 1e9+7;
#define ll long long
const int maxn = 1010;
double dp[maxn][maxn];
int n,s;
int main(){
    while(cin>>n>>s){
        dp[n][s] = 0.0;
        for(int i = n;i >= 0;i --){
            for(int j = s;j >= 0;j --){
                if(i == n&&j == s)continue;
                dp[i][j] = ((n-i)*j*dp[i+1][j]+i*(s-j)*dp[i][j+1]+(n-i)*(s-j)*dp[i+1][j+1]+n*s)/(n*s-i*j);
            }
        }
        printf("%.4f\n",dp[0][0]);
    }
    return 0;
}