POJ2096 Collecting Bugs

題目連接

題意

有n個bug和s個系統。天天會等機率的在某個系統中找一個某一個bug。問找到所有n個bug而且在每一個系統中都找到bug的指望天數。

思路

f[i][j]表示在j個系統中找到了i個bug的指望天數。

某天有四種可能：
1.在新的系統中找到新的bug，機率爲\(\frac{(n-i)\times(s-j)}{n\times s}\)
2.在舊的系統中找到新的bug，機率爲\(\frac{(n-i)\times j}{n\times s}\)
3.在新的系統中找到舊的bug，機率爲\(\frac{i \times (s-j)}{n\times s}\)
4.在舊的系統中找到舊的bug，機率爲\(\frac{i \times j}{n\times s}\)

即\[f[i][j]= \\ f[i+1][j+1] \times \frac{(n-i)\times(s-j)}{n\times s} + \\ f[i + 1][j] \times \frac{(n-i)\times j}{n\times s} + \\ f[i][j + 1] \times \frac{i \times (s-j)}{n\times s} + \\ f[i][j] \times \frac{i \times j}{n\times s} + 1 \]

左右都有\(f[i][j]\)，移項並同乘\(n\times s\)得

\[f[i][j] \times (n\times s -i \times j)= \\ f[i+1][j+1] \times (n-i)\times(s-j) + \\ f[i + 1][j] \times (n-i)\times j + \\ f[i][j + 1] \times i \times (s-j) + \\ n \times s \]

\[ f[i][j] =\frac{f[i+1][j+1] \times (n-i)\times(s-j) +f[i + 1][j] \times (n-i)\times j + f[i][j + 1] \times i \times (s-j) + n \times s}{n \times s} \]

轉移便可。lj POj的評測姬，在\(G++\)下輸出\(lf\)會\(wa\)

代碼

/*
* @Author: wxyww
* @Date:   2019-05-14 15:24:19
* @Last Modified time: 2019-05-14 15:59:55
*/
#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<ctime>
using namespace std;
typedef long long ll;
const int N = 1010;
ll read() {
    ll x=0,f=1;char c=getchar();
    while(c<'0'||c>'9') {
        if(c=='-') f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9') {
        x=x*10+c-'0';
        c=getchar();
    }
    return x*f;
}
double f[N][N];
int main() {
    int n,s;
    while(~scanf("%d%d",&n,&s)) {
        for(int i = n;i >= 0;--i) {
            for(int j = s;j >= 0;--j) {
                if(i == n && s == j) continue;
                f[i][j] = 
                    (f[i + 1][j + 1]* (n - i) * (s - j)  + //在新系統裏發現新漏洞
                    f[i + 1][j] * (n - i) * j  + //在舊系統裏發現新漏洞
                    f[i][j + 1]* i * (s - j)  + //在新系統裏發現舊漏洞
                    + n * s) //在舊系統發現舊漏洞
                    / (n * s - i * j );
            }
        }
        printf("%.4lf\n",f[0][0]);
    }
    return 0;
}