排列最小值

也许更好的阅读体验

Description\mathcal{Description}
problem
多组询问
1n,Q1051 \leq n,Q\leq 10^5

Solution\mathcal{Solution}

5050分解法

考虑DPDP,感觉上直接算答案不好算,所以考虑算长度为nn的所有排列改变的次数的排列的个数算出来
之后再把个数乘以次数的平方即可
fi,jf_{i,j}表示长度为ii的排列的改变次数jj排列个数
考虑把最大的ii插入到长度为i1i-1的改变次数为jj的排列中,则有ii个空可以插,其中最前面的空插进去会使改变次数加一,所以有i1i-1个空使改变次数不变
考虑插在最前面则原改变次数应为j1j-1
所以fi,j=fi1,j(i1)+fi1,j1f_{i,j}=f_{i- 1,j}*(i-1)+f_{i-1,j-1}
其实就是第一类斯特林数的递推公式
提前预处理一下即可做到 n2n^{2}

100100分解法

仍然是DPDP,原来的是把问题拆开,好做但是不能直接算
不妨大胆一点,设fif_i表示长度为ii答案
仍然考虑从i1i-1转移过来
还是上面那句
考虑把最大的ii插入到长度为i1i-1的改变次数为jj的排列中,则有ii个空可以插,其中最前面的空插进去会使改变次数加一,所以有i1i-1个空使改变次数不变
fi=fi1(i1)+f_i=f_{i-1}*(i-1)+插到最前面的贡献
考虑插在最前面,原本的改变次数为xx的都变为x+1x+1
也就是x2>(x+1)2=x2+2x+1x^2 -> (x+1)^2=x^2+2x+1
我们设gig_i表示原本算贡献时是按照xx来算的长度为ii的答案
对于那个11,因为插在最前面后面的i1i-1个数共有(i1)!(i-1)!种组合
所以fi=fi1(i1)+fi1+2gi1+(i1)!=fi1i+fi1+(i1)!f_i=f_{i-1}*(i-1)+f_{i-1}+2g_{i-1}+(i-1)!=f_{i-1}*i+f_{i-1}+(i-1)!
gig_i同样这么推
gi=gi1i+(i1)!g_i=g_{i-1}*i+(i-1)!

Code\mathcal{Code}

/*******************************
Author:Morning_Glory
LANG:C++
Created Time:2019年09月28日 星期六 09时34分38秒
*******************************/
#include <cstdio>
#include <fstream>
using namespace std;
const int maxn = 100005;
const int mod = 998244353;
//{{{cin
struct IO{
	template<typename T>
	IO & operator>>(T&res){
		res=0;
		bool flag=false;
		char ch;
		while((ch=getchar())>'9'||ch<'0')	flag|=ch=='-';
		while(ch>='0'&&ch<='9')	res=(res<<1)+(res<<3)+(ch^'0'),ch=getchar();
		if (flag)	res=~res+1;
		return *this;
	}
}cin;
//}}}
int T,n;
int f[maxn],g[maxn],fac[maxn];
int main()
{
	fac[0]=1;
	for (int i=1;i<=100000;++i){
		fac[i]=1ll*fac[i-1]*i%mod;
		f[i]=((1ll*f[i-1]*i%mod+2ll*g[i-1]%mod)%mod+fac[i-1])%mod;
		g[i]=(1ll*g[i-1]*i%mod+fac[i-1])%mod;
	}
	cin>>T;
	while (T--){
		cin>>n;
		printf("%d\n",f[n]);
	}
	return 0;
}

如有哪里讲得不是很明白或是有错误,欢迎指正
如您喜欢的话不妨点个赞收藏一下吧

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章