Технокубок 2018 - Отборочный Раунд 3 E. Maximum Element DP

Description
求n的排列有多少满足：从左到右扫求最大值，当最大值不再变化k次后的最大值不为n。

---

Sample Input
5 2

---

Sample Output
22

---

这题氪了一波金。。。
其实有很多解决方法，在这里我提供两种。
首先定义合法为满足出现一个最大值它是假的。
第一种：
设$f[i]$为填到第$i$个数都是不合法的。
那么可以得到转移：
$f[i]={\sum }_{j=0}^{min(k-1,i-1)}f[i-j-1]\ast C(i-1,j)\ast jc[j]$
你可以化一下，变成：$f[i]={\sum }_{j=0}^{min(k-1,i-1)}f[i-j-1]\ast jc[i-j-1]\ast jc[i-1]$
这个$f[i-j-1]\ast jc[i-j-1]$你可以用一个前缀和来优化它。
最后统计答案时，你枚举n出现的位置。
那么你要保证n前面是合法的，后面就随便填。
那不合法的方案就是：${\sum }_{i=1}^{n}f[i-1]\ast C(n-1,i-1)\ast jc[n-i]$
用总数减掉即可。

第二种：
设$f[i]$为填到第$i$个数是合法的。
设$g[i]$为填到第$i$个数是不合法的。
那么你有两种情况的转移，第一种是前面已经构成合法方案：
$f[i]=f[i-1]\ast C(i,i-1)$ => $f[i]=f[i-1]\ast i$
第二种是当前这里构成了一个新和方案：
$f[i]=g[i-k-1]\ast C(i-1,k)\ast jc[k]$
对于$g$的转移直接用总数减去$f$即可。

我写了第一种。

---

#include <cstdio>
#include <cstring>

using namespace std;
typedef long long LL;
int _max(int x, int y) {return x > y ? x : y;}
const LL mod = 1e9 + 7;
int read() {
	int s = 0, f = 1; char ch = getchar();
	while(ch < '0' || ch > '9') {if(ch == '-') f = -1; ch = getchar();}
	while(ch >= '0' && ch <= '9') s = s * 10 + ch - '0', ch = getchar();
	return s * f;
}

LL inv[1000010], jc[1000010];
LL f[1000010], sum[1000010];

LL pow_mod(LL a, LL k) {
	LL ans = 1;
	while(k) {
		if(k & 1) (ans *= a) %= mod;
		(a *= a) %= mod; k /= 2;
	} return ans;
}

LL C(LL n, LL m) {return jc[n] * inv[m] % mod * inv[n - m] % mod;}

int main() {
	int n = read(), k = read();
	jc[0] = 1; for(int i = 1; i <= n; i++) jc[i] = (LL)jc[i - 1] * i % mod;
	inv[0] = 1; inv[n] = pow_mod(jc[n], mod - 2);
	for(int i = n - 1; i >= 1; i--) inv[i] = (LL)inv[i + 1] * (i + 1) % mod;
	sum[1] = 1; f[0] = 1;
	for(int i = 1; i <= n; i++) {
		f[i] = (sum[i] - sum[_max(0, i - k)]) % mod * jc[i - 1] % mod;
		(sum[i + 1] = sum[i] + f[i] * inv[i] % mod) % mod;
	} LL ans = 0;
	for(int i = 1; i <= n; i++) (ans += f[i - 1] * C(n - 1, i - 1) % mod * jc[n - i] % mod) % mod;
	printf("%lld\n", ((jc[n] - ans) % mod + mod) % mod);
	return 0;
}