CSP-S 模拟 七十和十七 （期望DP）（找规律）

传送门
瞪一下样例发现，每次排好序都是先把 $n-1$ 位排好序，再排第 $n$ 位
如果我们令 $f_{n-1}$ 表示到 $n-1$ 的所有排列的计数次数，然后算出第 $n$ 位的新增贡献就可以递推了
手玩一下发现如果 $a_n=n$ 不需要归位，否则需要 $2^{a_i-1}$ 次
然后就可以枚举最后一个填什么转移，发现是个等比数列
$f_i=i\ast f_{i-1}+(2^{i-1}-1)\ast (i-1)!$
除以总方案就是期望

#include<bits/stdc++.h>
#define cs const
using namespace std;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
	while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
	return cnt * f;
}
cs int Mod = 1e9 + 7;
int add(int a, int b){ return a + b >= Mod ? a + b - Mod : a + b; }
int mul(int a, int b){ return 1ll * a * b % Mod; }
int ksm(int a, int b){ int ans = 1; for(;b;b>>=1, a=mul(a,a)) if(b&1) ans=mul(ans, a); return ans;}
cs int N = 1e5 + 5;
int n, E[N], pw[N], fac[N];
int main(){
	n = read();
	pw[0] = fac[0] = 1;
	for(int i = 1; i <= n; i++) fac[i] = mul(fac[i-1], i), pw[i] = add(pw[i-1], pw[i-1]);
	for(int i = 1; i <= n; i++) pw[i] = add(pw[i], pw[i - 1]);
	E[1] = 0;
	for(int i = 2; i <= n; i++){
		E[i] = add(mul(E[i-1], i), mul(fac[i - 1], pw[i - 2]));
	} cout << mul(ksm(fac[n], Mod-2), E[n]); return 0;
}