【XSY2786】Mythological VI 数学

题目描述

　　有\(1\sim n\)一共\(n\)个数。保证\(n\)为偶数。

　　你要把这\(2n\)个数两两配对，一共配成\(n\)对。每一对的权值是他们两个数的和。

　　你想要知道这\(n\)对里最大的权值的期望是多少。

　　请输出答案对\(1000000007\)取模的值。

　　\(n\leq 500000\)

题解

　　枚举\(v\)，计算最大权值\(\leq v\)的概率。

　　从大到小枚举\(> \frac{v}{2}\)的数，这些数每次都有\(v-n\)种选择，方案数为
\[ {(v-n)}^{n-\frac{v}{2}} \]
　　\(\leq \frac{v}{2}\)的数可以随便匹配。

　　记\(f(x)\)为\(x\)个数随便匹配的方案数，那么
\[ f(x)=\frac{x!}{2^{\frac{x}{2}}(\frac{x}{2})!} \]
　　（考虑\(x\)的全排列，\(a_{2i-1}\)和\(a_{2i}\)匹配。）

　　时间复杂度：\(O(n\log n)\)

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll p=1000000007;
ll fp(ll a,ll b)
{
    ll s=1;
    for(;b;b>>=1,a=a*a%p)
        if(b&1)
            s=s*a%p;
    return s;
}
ll f[500010];
ll g[1000010];
int main()
{
#ifndef ONLINE_JUDGE
    freopen("b.in","r",stdin);
    freopen("b.out","w",stdout);
#endif
    int n;
    scanf("%d",&n);
    ll ans=0;
    f[0]=1;
    for(int i=2;i<=n;i+=2)
        f[i]=f[i-2]*(i-1)%p;
    for(int i=n+1;i<=2*n;i++)
        g[i]=fp(i-n,n-i/2)*f[n-2*(n-i/2)]%p;
    for(int i=2*n;i>=n+1;i--)
        g[i]=(g[i]-g[i-1])%p;
    for(int i=n+1;i<=2*n;i++)
        ans=(ans+g[i]*i)%p;
    ans=ans*fp(f[n],p-2)%p;
    ans=(ans+p)%p;
    printf("%lld\n",ans);
    return 0;
}

转载于:https://www.cnblogs.com/ywwyww/p/8611973.html