【题目】HDU P6595（C849E）Everything Is Generated In Equal Probability

题目大意

给出一个整数$N$和一段程序（包括三个函数）：

• $\text{SUBSEQUENCE(}Array\text{)}$：随机返回$Array$的一个子序列（可为空）。
• $\text{CNT-INVERSION-PAIRS(}Array\text{)}$：返回$Array$的逆序对数。
• $\text{CALCULATE(}Array\text{)}$：若$Array$为空，返回$0$；否则返回$\text{CNT-INVERSION-PAIRS(}Array\text{)}+\text{CALCULATE(SUBSEQUENCE(}Array\text{))}$。

先随机生成一个整数$n\in [1,N]$，再随机生成一个长度为$n$的排列，求这个排列的$\text{CALCULATE()}$函数的期望值。
有多组测试数据，每组数据$N⩽3000$。

---

思路

先求出长度为$n$的排列的$\text{CNT-INVERSION-PAIRS()}$函数的期望值，有两种思路：

思路1： 长度为$n$的排列中一共有$C_n^2$对数对，每一对逆序的概率为$\frac{1}{2}$，累加起来就是
$$\frac{C_n^2}{2}=\frac{n!}{2(2!)(n-2)!}=\frac{n(n-1)}{4}$$

思路2： 设该期望值对应的状态为$d{p}_1[n]$，显然$d{p}_1[1]=0$，$d{p}_1[2]=\frac{1}{2}$。考虑将数字$n$插入$1$到$n-1$的排列中，一共有$n$个位置可插入，增加的逆序对数分别为$0,1,2,\dots ,n-1$，于是得到
$$\begin{array}{rl}d{p}_1[n]& =d{p}_1[n-1]+\frac{0+1+2+\cdots +(n-1)}{n}\\ & =d{p}_1[n-1]+\frac{n-1}{2}\\ d{p}_1[n]& =\frac{n(n-1)}{4}\end{array}$$

再求长度为$n$的排列的$\text{CALCULATE()}$函数的期望值。设该期望值对应的状态为$d{p}_2[n]$，其中$d{p}_2[1]=0$。由题意不难得到
$$\begin{array}{rl}d{p}_2[n]& =d{p}_1[n]+\frac{{\sum }_{i=0}^n{C}_n^i\cdot d{p}_2[i]}{{\sum }_{i=0}^n{C}_n^i}\\ & =\frac{n(n-1)}{4}+\frac{{\sum }_{i=0}^{n-1}{C}_n^i\cdot d{p}_2[i]}{2^n}+\frac{d{p}_2[n]}{2^n}\\ d{p}_2[n]& =\frac{2^{n-2}n(n-1)+{\sum }_{i=0}^{n-1}{C}_n^i\cdot d{p}_2[i]}{2^n-1}\end{array}$$

于是可以递推求出$d{p}_2[n]$，时间复杂度$O(n^2)$。
计算出$d{p}_2[n]$后，由于$n\in [1,N]$，显然
$$ans=\frac{{\sum }_{n=1}^{N}d{p}_2[n]}{N}$$

预处理$d{p}_2[n]$的前缀和后就可以$O({\log }_{2}n)$求解。如果预处理了$N$的逆元，求解可以达到$O(1)$。