Subpermutation-2021中国大学生程序设计竞赛（CCPC）- 网络选拔赛（重赛）

题意：
把所有$n$元排列按字典序从小到大排成一列，求其中有多少子串为$m$元排列,$t$组询问。
题解：
(1)考虑单个$n$元组中的$m$元组个数。
将$m$元组看成一个元素，排列组合有$m!$的情况，现在共有$(n-m+1)$个元素，排列组合为$(n-m+1)!$。共有$m!(n-m+1)!$的$m$元组。
(2)考虑跨两个$n$元组的$m$元组个数。
假设$k$为某个排列中的最后一个使得$p_k<p_{k+1}>pk+2$的位置。
即$p_1,p_2,p_3...p_k<p_{k+1}>p_{k+2}>p_{k+3}...>p_n$
那么下一个字典序的排列为:
$p1,p2,p3...p_j<p_n<...<p_{j+1}<p_k<p_{j-1}<....<p_{k+1}(p_j>p_k)$
假设$i$为我们所要的$m$元组开头的位置。

1.当$i<=k$时：

第一个排列中存在$n-i+1$个数，剩下的数在第二排列中，因为$m<n$,所以有$m-(n-i+1)<k$。
因为$i<=k$所以在$[i,n]$内必须存在一个位置$k$使得$p_k<p_{k+1}$。方案数转化为全部的方案减去$p_i>p_{i+1}>...>p_n$的方案数。
也就是$m!(全部的方案)-{\complement }_m^{n-i+1}(m-(n-i+1))!(m个数选出n-i+1个数在第一个排列中,且为降序排列的方案)$
剩下的$n-m$个数随便排列。
于是此时的方案数为$(n-m)!(m!-{\complement }_m^{n-i+1}(m-(n-i+1))!)$

2.当$i>=k+1$时：

此时$p_k$不在$m$元组中,$p_k>m$，而$p_j>p_k>m$，所以$p_j$也不再$m$元组中。
同样：(标黄的是m元组出现的范围)
$p_1,p_2,p_3...p_k<p_{k+1}>$$p_{k+2}>p_{k+3}...>p_n$
那么下一个字典序的排列为:
$p1,p2,p3$$...p_j<p_n<...<p_{j+1}<p_k<p_{j-1}<....<p_{k+1}(p_j>p_k)$
这就意味着在第二个排列中，最后一个元素的下标为$i+m-1-n<j$
前$n-i+1$个数是递减的,剩下$(m-(n-i+1))$个数没有限制，方案数为$(m-(n-i+1))!$。对于$m$元组以外的$m-n$个数，必须存在一个$k$使得$p_k<p_{k+1}$,方案数为$((n-m)!-1)$

$i$的范围:
第一个排列中的数要小于$m$.
$(n-i+1)>=m-1,i>=n+2-m$

最后答案就是$m!(n-m+1)!+{\sum }_{i=n+2-m}^n((n-m)!(m!-{\complement }_m^{n-i+1}(m-(n-i+1))!)+{\complement }_m^{n-i+1}(m-(n-i+1))!((n-m)!-1))$
化简得：$m!(n-m+1)!+(m-1)(n-m)!m!-m!{\sum }_{i=1}^{m-1}\frac{1}{i!}$

#include <cstdio>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <map>
#include <string>
#include <stack>
#include <cctype>
#include <vector>
#include <queue>
#include <set>
#include <utility>
#include <cassert>
#include <iomanip>
#include <deque>
#include <time.h>
#include <bitset>
using namespace std;
#define ll long long
#define maxn 1000005
#define mod 1000000007
#define MOD 998244353
#define Mod 1000000009
#define eps 1e-10
const ll inf=0x3f3f3f3f3f3f3f3f;
const ll INF=0x3f3f3f3f;
const ll mod1=1e9+7;
const ll mod2=1e9+9;
template <typename T>
inline void read(T& X) {X = 0; int w = 0; char ch = 0;while (!isdigit(ch)) { w |= ch == '-'; ch = getchar(); }while (isdigit(ch)) X = (X << 3) + (X << 1) + (ch ^ 48), ch = getchar();if (w) X = -X;}
char F[200];inline void write(int x){if(x == 0){putchar('0');return;}int tmp = x > 0 ? x : -x;int cnt = 0;if(x < 0)putchar( '-' );while(tmp > 0){F[cnt++] = tmp % 10 + '0';tmp /= 10;}while(cnt > 0)putchar(F[--cnt]) ;}
template<typename T> void print(T x){if(x>9) print(x/10);putchar(x%10+'0');}
ll q_pow(ll x,ll y,ll M){ll ans=1;while(y){if(y%2){y--;ans=ans*x%M;}else {y/=2;x=x*x%M;}}return ans;}
ll a[maxn],b[maxn],f[maxn];
void init(){
	a[0]=1;
	for(ll i=1;i<=1000000;i++)a[i]=a[i-1]*i%mod;
	b[1000000]=q_pow(a[1000000],mod-2,mod);
	for(ll i=999999;i>=0;i--)b[i]=b[i+1]*(i+1ll)%mod;
	for(int i=1;i<=1000000;i++)f[i]=(f[i-1]+b[i])%mod;
}
ll C(ll n,ll m){
	return a[n]*b[m]%mod*b[n-m]%mod;
}
int main() 
{
	ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
	int t;
	ll n,m;
	init();
	cin>>t;
	while(t--){
		cin>>n>>m;
		ll ans=a[m]*a[n-m+1]%mod+(m-1)*a[n-m]%mod*a[m]%mod-a[m]*f[m-1]%mod;
		ans=(ans%mod+mod)%mod;
		cout<<ans<<"\n";
	}
    return 0;
}