19 hdu多校 Everything Is Generated In Equal Probability // 逆序对+期望+组合数

http://acm.hdu.edu.cn/showproblem.php?pid=6595

题意：给N，随机的从[1,N]中得到一个n，然后求CALCAULATE(n)。求这种操作的期望。

思路：

任意长度为n的数组得到逆序对的期望是：1/2*C(n,2)   (是否为逆序对的期望是相同的）

注意： C[i][j]=(C[i-1][j]+C[i-1][j-1])

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define FI first
#define SE second
#define MP make_pair
#define PII pair<int,int>
const int MXLOG = 30;
const LL mod = 998244353;
const int MX = 1e6+5;
LL ans[3005],f[3005],C[3005][3005],two[3005];
void exgcd(LL a, LL b, LL &x, LL &y, LL &d) {  //欧几里得函数
    if (!b) {d = a, x = 1, y = 0;}
    else{
        exgcd(b, a % b, y, x, d);
        y -= x * (a / b);
    }
}
LL inv(LL t, LL p){                   //返回t对p的逆元
    LL d, x, y;
    exgcd(t, p, x, y, d);
    return (x % p + p) % p;        //x可能为负，也可能过大
}

int main(){
    two[0]=1;
    for(LL i=1;i<=3000;i++){
        two[i]=two[i-1]*2LL%mod;
        C[i][i]=1,C[i][0]=1;
        for(LL j=1;j<i;j++){
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;

        }
    }
    //cout<<C[5][2]<<endl;
    for(int i=2;i<=3000;i++){
        f[i]=two[i-1]*C[i][2]%mod;
        for(int j=2;j<=i-1;j++){
            f[i]=(f[i]+(C[i][j]*f[j])%mod)%mod;
        }
        f[i]=f[i]*inv(two[i]-1LL,mod)%mod;
    }
    for(int i=2;i<=3000;i++) ans[i]=(ans[i-1]+f[i])%mod;
    int n;while(cin>>n){
        cout<<ans[n]*inv(n,mod)%mod<<'\n';
    }
    return 0;
}