P4562 [JXOI2018]游戏(思维+组合数)

这题真是米奇妙妙屋呀~妙到家了

披着期望的外观,其实是个思维+组合数

每个数有两种消除方式

$Ⅰ.被自己的因子消除$

$Ⅱ.被自己消除$

如果一个数在$[l,r]$没有因子,这个数就必须要选择,我们称之为关键数

而且当你选择了所有关键数,其他所有数字选不选都不重要了,已经完成了目的

所以发现$tp(p)$的定义就是排列$p$中关键数字的最后一个的位置

假设有$ans$个关键数字,$len=r-l+1$

枚举最后一个关键数字放在$pos$位置$(pos\in [1,len])$

那么$[pos+1,len]不允许存在关键数字$

所以是$C_{len-ans}^{len-pos}\ast (len-pos)!\ast (pos-1)!$

$又因为最后一个关键数有ans种可能,再乘以ans$

枚举计算答案即可

求逆元我用的是线性递推

#include <bits/stdc++.h>
using namespace std;
#define int long long
const int mod=1e9+7;
const int maxn=1e7+10;
int l,r,ans,len,fac[maxn],ifac[maxn];
bool vis[maxn];
int quick(int x,int n)
{
	int ans=1;
	while( n )
	{
		if( n&1 )	ans = ans*x%mod;
		x = x*x%mod;
		n>>=1; 
	}
	return ans;
}
inline int C(int n,int m){return n<m?0:1ll*fac[n]*ifac[m]%mod*ifac[n-m]%mod;}
signed main()
{
	fac[0]=1;
	for(int i=1;i<=10000000;i++)	fac[i]=fac[i-1]*i%mod;
	cin >> l >> r;
	len=r-l+1;
    ifac[len]=quick(fac[len],mod-2);
    for(int i=len-1;i>=0;i--)
        ifac[i] = ifac[i+1]*(i+1)%mod;
	for(int i=l;i<=r;i++)
		if( !vis[i] )
		{
			ans++;
			for(int j=i+i;j<=r;j+=i)	vis[j]=1;
		}
	int he=0;
	for(int i=ans;i<=len;i++)//最后一个数放在哪里 
	{
		he += i*fac[len-i]%mod*fac[i-1]%mod*C(len-ans,len-i)%mod*ans%mod;
		he%=mod;
	}
	cout << he;
}