2020牛客寒假算法基础集训营2 C

链接：https://ac.nowcoder.com/acm/contest/3003/C
来源：牛客网

牛牛刚刚考完了期末，尽管 牛牛 做答了所有 n 道题目，但他不知道有多少题是正确的。
不过，牛牛 知道第i 道题的正确率是 pi
牛牛 想知道这 n 题里恰好有 0,1…,n 题正确的概率分别是多少，对 10^9+7 取模。
对1e9+7取模的含义是：对于一个b!=0的不可约分数a/b，存在q使得b*qmod（1e9+7）=a，q即为a/b对1e9+7取模的结果。

输入描述
第一行，一个正整数 n 。
第二行n个整数 ，在模 10^9+7 意义下给出。
保证 1<=n<=2000

输出描述
输出一行n+1个用空格隔开的整数表示答案（对1e9+7取模）。

输入样例
1
500000004
输出样例
500000004 500000004

#include<iostream>
#include<bits/stdc++.h>
using namespace std;
const int N=1e9+7;
int n,i,j;
long long p[2005],f[2005][2005];//f[i][j]代表前i道题中有j道做对的概率
int main()
{
    scanf("%d",&n);
    for (i=1; i<=n; i++)
    {
        scanf("%lld",&p[i]);
    }
    f[0][0]=1;
    for (i=1; i<=n; i++)
    {
        f[i][0]=f[i-1][0]*(N+1-p[i])%N;//全部做错的概率
        for (j=1; j<=i; j++)
        {
            f[i][j]=(f[i-1][j]*(N+1-p[i])+f[i-1][j-1]*p[i])%N;//前i题有j题做对的概率
        }//      前i-1题有j题做对*第i题做错+前i-1题有j-1题做对*第i题做对
    }
    for (int i=0; i<=n; i++)
    {
        printf("%lld ",f[n][i]);
    }
    return 0;
}
//对于不可约分数a/b 使b=1 即q%1e9+7=a a为每道题做对的概率 q=k*1e9+7+a（k>=0）
//取k=1 q=1e9+7+a  a=q-1e9+7 每题错误的概率为1-a即 1e9+7+1-q
//取k=0 a=q 每道题正确的概率即为q

如有不正确的地方请指出XD