NOI模拟七十和十七【递推找规律】

最新推荐文章于 2022-07-28 20:39:24 发布

gigo_64

最新推荐文章于 2022-07-28 20:39:24 发布

阅读量272

点赞数

分类专栏： ----------数论-------------

本文链接：https://blog.csdn.net/qq_42835815/article/details/102944200

版权

----------数论------------- 专栏收录该内容

36 篇文章 0 订阅

订阅专栏

传送门

本质还是求按照规则sort好的次数。

我们把问题划分。我们发现，如果前i-1位没弄好，那第i位不会动。

所以考虑i-1位都弄好了，整第i位。

手玩一下发现，如果第i位值为k，需要花费 $\sum_{i=1}^{k-1}2^{i-1}-1$

不管阶乘，令 f(i) 表示长度为i的方案数。

考虑新加入一个数。不论这个数是多少，总能使前i-1个数选哪些数的方案从(i-1)!变成i!/1=i!。所以多乘了一个i。

然后第i个数就是2^{i-1}*(i-1)!。阶乘是因为前面可以乱排。

所以 $f(i)=i*f(i-1)+(2^{i-1}-1)(i-1)!$

最后f(n)乘阶乘。

但是我发现还有一种写法：

答案是 $\frac{f(n)}{n!}$

我们将f(n)强行展开，发现得到的系数刚好和n!的第一项抵消，以此类推。

所以 $Ans=\sum_{i=1}^n\frac{2^{i-1}-1}{i}$

这样更好写

#include<bits/stdc++.h>
using namespace std;
#define in read()
#define int long long
int in{
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){
		ch=getchar();if(ch=='-')f=-1;
	}
	while(isdigit(ch)){
		cnt=cnt*10+ch-48;
		ch=getchar();
	}return cnt*f;
}
int n;
int mod=1e9+7;
int ksm(int a,int b){
	int sum=1;
	while(b){
		if(b&1)sum=sum*a%mod;a=a*a%mod;b>>=1;
	}return sum;
}
signed main(){
	n=in;
	int ans=0;
	for(int i=1;i<=n;i++){
		ans=(ans+(ksm(2,i-1)-1)*ksm(i,mod-2)%mod)%mod;
	}cout<<ans;
	return 0;
}