L - Lucky Draw概率dp

题目链接：https://vjudge.ppsucxtt.cn/contest/450938#problem/L
题目大意：n个人，每人k滴血,一个人一轮p的概率不掉一滴血,求他们一块死的概率;

#include<bits/stdc++.h>
using namespace std;
const int N=10005;
const int INF=0x3f3f3f3f;
const int m=1e+3;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PII;
double f[N][N];//f[i][j]表示一个人在k轮后j滴血的概率
int main()
{
    int n,k;
    double p;
    cin>>n>>k>>p;
    f[0][k]=1;
    for(int i=1;i<=m;++i)//轮数越大,对答案影响越小
    {
        for(int j=k;j>0;--j)
        {
            f[i][j]=f[i-1][j]*p+f[i-1][j+1]*(1-p);
          	 		//上一局就是j滴血*不掉血概率,上一局j+1滴血*掉血概率
        }
        f[i][0]=f[i-1][0]+(1-p)*f[i-1][1];//特判0血(上一局无法继续掉血了)
    }

    double ans=0;//ans表示一个人获胜的概率;
    for(int i=1;i<m;++i)//枚举每一局获胜的概率
    {
        ans+=(f[i][0]-f[i-1][0])*(pow(f[i-1][0],n-1));
        //在第i局获胜的概率=(正好在第i局死去)*(其他人在第i-1局之前死去)
        //正好在第i局死去=在第i局死亡-在第i局死亡；
        //其他人在第i-1局死亡=一个人在第i-1局之前死亡概率*其他人数;
    }
    printf("%.10f\n",1-n*ans);

}