JXOI 2018 游戏 题解

题目传送门

题目大意： 给出 $l,r$，当你消除 $i$ 时，它的倍数也会被消除，，对于一个 $l$ ~ $r$ 的排列，从左到右依次消除，全部消除完时，记消除次数为 $P$，求 $l$ ~ $r$ 组成的所有排列的 $P$ 之和。

题解

一开始想错方向了qwq……如果去考虑每个数在所有情况下的贡献，也就是考虑有多少种排列中它的前面没有它的约数，这样的话时间复杂度大概是 $O(n^{2}logn)$ 的……

于是考虑每种排列的贡献，可以发现，有一些数在 $l$ ~ $r$ 间没有约数，也就是说，它无法成为别人的约数从而在消除别人的时候消除它，所以这些数是一定要被主动消除掉的，设这些数有 $s$ 个。

显然，我们只需要从左到右消除到这 $s$ 个数都消除完，其它数也肯定消除完了。

于是我们可以枚举这些数中位置最后的那个的位置，设为 $i$，然后前面 $i-1$ 个位置中有剩下的 $s-1$ 个特殊的数，还剩下 $i-s$ 个空位，要从 $n-s$ 个数中选择然后放进去，方案数 $C_{n-s}^{i-s}$，然后这 $i-1$ 个数顺序是随意的，所以乘上 $(i-1)!$，以及剩下的在后面的 $n-i$ 个数的位置也随意，所以再乘上 $(n-i)!$，最后，这些方案的贡献都是 $i$，所以再乘 $i$ 即可。

代码如下：

#include <cstdio>
#define mod 1000000007
#define maxn 10000010

int l,r,s=0,ans=0;
bool v[maxn];
int fac[maxn],inv_fac[maxn];
int ksm(int x,int y)
{
	int re=1;
	while(y)
	{
		if(y&1)re=1ll*re*x%mod;
		x=1ll*x*x%mod;y>>=1;
	}
	return re;
}
#define inv(x) ksm(x,mod-2)
void work()
{
	fac[0]=1;
	for(int i=1;i<=maxn-10;i++)fac[i]=1ll*fac[i-1]*i%mod;
	inv_fac[maxn-10]=inv(fac[maxn-10]);
	for(int i=maxn-11;i>=0;i--)inv_fac[i]=1ll*inv_fac[i+1]*(i+1)%mod;
}
int C(int x,int y){return 1ll*fac[x]*inv_fac[y]%mod*inv_fac[x-y]%mod;}

int main()
{
	scanf("%d %d",&l,&r); work();
	#define n r-l+1
	for(int i=l;i<=r;i++)if(!v[i])
	{s++;for(int j=i;j<=r;j+=i)v[j]=true;}
	for(int i=s;i<=r-l+1;i++)
	ans=(ans+1ll*s*C(n-s,i-s)%mod*fac[i-1]%mod*fac[n-i]%mod*i%mod)%mod;
	printf("%d",ans);
}