hdu多校第二场 1005 (hdu6595) Everything Is Generated In Equal Probability

题意:

给定一个N,随机从[1,N]里产生一个n,然后随机产生一个n个数的全排列,求出n的逆序数对的数量,加到cnt里,然后随机地取出这个全排列中的一个非连续子序列(注意这个子序列可以是原序列),再求出这个子序列的逆序数对,加到cnt里,重复这个过程,直到最后取出的为空。

题解:

先不考虑第一步随机从[1,N]里产生一个n,只考虑n给定的情况,求出了f[n],那么最后的结果就是

$  ans[N]=frac{sum_{n=1}^N f[n]}{N} $

赛时和队友利用找规律法和暴力模拟法推出

$ f[i]=sum_{i=1}^{n-1}frac{2i}{3} $

下面给出证明:

因为任意一个长度为n的全排列,其所含的逆序对的期望为

$ inom{n}{2}/2 $ (不难理解,就是随便取两个点交换一下就出来一个逆序对)

而取出的那一个非连续子序列,我们把它里面的数字离散化以后也是一个全排列,所以

$ f[i]=inom{n}{2}/2+frac{1}{2^i}sum_{j=0}^iinom{i}{j}f[j] $

移项,得到递推公式

$ f[i]=frac{2^{i-1}}{2^i-1}inom{n}{2}/2+frac{1}{2^i-1}sum_{j=0}^{i-1}inom{i}{j}f[j] $

f[1]=0,不难推出通项公式

$ f[i]=sum_{i=1}^{n-1}frac{2i}{3} $

代码:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cctype>
#include <string>
#include <vector>
#include <map>
#include <set>
#include <vector>
#include <queue>
#include <stack>
const int MOD = 998244353;
#define rep(i,n,m) for(int i=n;i<=m;++i)
const double eps = 1e-16;
#define ll long long
using namespace std;
const int maxn = 10000000;
const int maxm = 2e5 + 10;
const int inf = 1 << 28;
typedef pair<int, int> P;
#define zero(a) fabs(a)<eps
inline int read()
{
    int x = 0, f = 1;
    char ch = getchar();
    while (ch < '0' || ch > '9')
    {
        if (ch == '-')
            f = -1;
        ch = getchar();
    }
    while (ch >= '0' && ch <= '9')
    {
        x = 10 * x + ch - '0';
        ch = getchar();
    }
    return x * f;
}


ll quick_mod(ll x, ll n) {
    ll res = 1;
    while (n) {
        if (n & 1) res = res * x % MOD;
        x = x * x % MOD;
        n = n >> 1;
    }
    return res;
}
ll a[3005];
void init()
{
    a[0] = a[1] = 0;
    a[2] = 2;
    for (int i = 3; i <= 3000; ++i)
        a[i] = a[i - 1] + (i - 1) * 2;
    for (int i = 1; i <= 3000; ++i)
        a[i] += a[i - 1];
}
int main()
{
    ll n;
    init();
    while (~scanf("%lld", &n))
    {
        if (n == 1) { cout << "0" << endl; continue; }
        ll a1 = a[n];
        ll b = 3 * n;
        ll x = quick_mod(b, MOD - 2);
        cout << a1 % MOD * x % MOD<<endl;
    }
}
原文地址:https://www.cnblogs.com/isakovsky/p/11247444.html