HZOJ 大佬(kat)

及其水水水的假期望（然而我已經被期望嚇怕了……）。

數據范圍及其沙雕導致丟掉5分……

因為其實每天的期望是一樣的，考慮分開。

f[i][j]表示做k道題，難度最大為j的概率。

則f[i][j]=(f[i-1][j])*(j-1)*tem+q[j]*tem;q為前綴和，tem為分母。

ps.if(k>n){puts("0");return 0;}

?

#include<iostream>
#include<cstring>
#include<cstdio>
#define MAXN 510
#define LL long long
#define ma(x,y) memset(x,y,sizeof(x))
using namespace std;
const int mod=1e9+7;
int n,m,k;
LL wt[MAXN];
LL sum;
LL f[MAXN][MAXN];
LL q[MAXN];
LL poww(LL a,int b,int mod)
{
    LL ans=1;
    while(b)
    {
        if(b&1)ans=ans*a%mod;
        a=a*a%mod;
        b=b>>1;
    }
    return ans;
}
inline int read();
signed main()
{
//    freopen("in.txt","r",stdin);
//    freopen("1.out","w",stdout);

    n=read();m=read();k=read();
    if(k>n){puts("0");return 0;}
    for(int i=1;i<=m;i++)wt[i]=read(),sum+=wt[i];
    if(k==1)
    {
        LL tem=sum%mod*poww(m,mod-2,mod)%mod;
        printf("%lld\n",tem*n%mod);
        return 0;
    }
    f[0][0]=1;
    LL tem=poww(m,mod-2,mod);
    for(int j=0;j<=m;j++)q[j]=1;
    for(int i=1;i<=k;i++)
    {
        for(int j=1;j<=m;j++)
            f[i][j]=(f[i-1][j]*(j-1)%mod*tem%mod+q[j]*tem%mod)%mod;
        q[0]=0;
        for(int j=1;j<=m;j++)
            q[j]=(q[j-1]+f[i][j])%mod;
    }
    LL ans=0;
    for(int i=1;i<=m;i++)ans=(ans+f[k][i]*wt[i]%mod)%mod;
    printf("%lld\n",ans*(n-k+1)%mod);
}
inline int read()
{
    int s=0;char a=getchar();
    while(a<'0'||a>'9')a=getchar();    
    while(a>='0'&&a<='9'){s=s*10+a-'0';a=getchar();}
    return s;
}

?