【JXOI 2018】游戏

传送门

---

problem

可怜公司有 $n$ 个办公室，办公室编号是 $l$ 到 $r$（$r-l+1=n$），可怜会事先制定一个顺序，按照这个顺序依次检查办公室。一开始的时候，所有办公室的员工都在偷懒，当她检查完编号是 $i$ 的办公室时候，这个办公室的员工会认真工作，并且这个办公室的员工通知所有办公室编号是 $i$ 的倍数的办公室，通知他们老板来了，让他们认真工作。因此，可怜检查完第 $i$ 个办公室的时候，所有编号是 $i$ 的倍数（包括 $i$）的办公室的员工会认真工作。

可怜发现了员工们通风报信的行为，她发现，对于每种不同的顺序 $p$，都存在一个最小的 $t(p)$，使得可怜按照这个顺序检查完前 $t(p)$ 个办公室之后，所有的办公室都会开始认真工作。她把这个 $t(p)$ 定义为 $p$ 的检查时间。

可怜想知道所有 $t(p)$ 的和。

但是这个结果可能很大，她想知道和对 $10^9+7$ 取模后的结果。

数据范围：$1\le l\le r\le 10^7$。

---

solution

比较妙的一道题。

如果一个数 $x$ 没有因子（除自己本身）在 $[l,r]$ 内，那我们称 $x$ 为关键数。显然，$t(p)$ 表示 $p$ 这种顺序中最后一个关键数的位置。

那我们可以用线性筛筛出每个数 $x$ 的最小质因子 $mi{n}_x$，那如果 $\frac{x}{mi{n}_x}<l$，$x$ 就是一个关键数。

接下来考虑如何计算答案。

我们把 $[l,r]$ 中的关键数数量统计出来（设为 $sum$），那么有：

$$ans=\sum _{i=sum}^{n}i\times \left(\genfrac{}{}{0ex}{}{i-1}{sum-1}\right)\times sum!\times (n-sum)!$$

这相当于就是我们枚举最后一个关键数的位置。在位置 $i$ 首先有 $i$ 的贡献，$\left(\genfrac{}{}{0ex}{}{i-1}{sum-1}\right)$ 表示在前 $i-1$ 个位置放入 $sum-1$ 个关键数，$sum!$ 表示关键数是无序的，要乘上一个排列，$(n-sum)!$ 同理。

注意要特判 $l=1$ 的情况。

---

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 10000005
#define P 1000000007
using namespace std;
const int inv2=500000004;
int l,r,fac[N],inv[N];
int add(int x,int y)  {return x+y>=P?x+y-P:x+y;}
int dec(int x,int y)  {return x-y< 0?x-y+P:x-y;}
int mul(int x,int y)  {return 1ll*x*y%P;}
int power(int a,int b,int ans=1){
	for(;b;b>>=1,a=mul(a,a))
		if(b&1)  ans=mul(ans,a);
	return ans;
}
int tot,sum,prime[N],mark[N],Min[N];
void linear_sieves(){
	for(int i=2;i<=r;++i){
		if(!mark[i])  prime[++tot]=i,Min[i]=i;
		for(int j=1;j<=tot&&i*prime[j]<=r;++j){
			mark[i*prime[j]]=1;
			Min[i*prime[j]]=prime[j];
			if(i%prime[j]==0)  break;
		}
	}
}
void prework(){
	fac[0]=fac[1]=1;
	for(int i=2;i<=r;++i)  fac[i]=mul(fac[i-1],i);
	inv[r]=power(fac[r],P-2);
	for(int i=r-1;~i;--i)  inv[i]=mul(inv[i+1],i+1);
}
int C(int n,int m)  {return mul(fac[n],mul(inv[m],inv[n-m]));}
int main(){
	scanf("%d%d",&l,&r);
	linear_sieves(),prework();
	int n=r-l+1,ans=0;
	if(l==1)  return printf("%d\n",mul(mul(mul(n,n+1),inv2),fac[n-1])),0;
	for(int i=l;i<=r;++i)  if(i/Min[i]<l)  sum++;
	for(int i=sum;i<=n;++i)  ans=add(ans,mul(i,mul(C(i-1,sum-1),mul(fac[sum],fac[n-sum]))));
	printf("%d\n",ans);
	return 0;
}