[BZOJ5323][线性筛][组合数学]JXOI2018：游戏

BZOJ5323

很简单，我们类似线性筛处理出有多少个数只能被它自己筛掉，然后剩下的组合数随便算一下就完了
枚举最后一个只能被自己筛掉的数的位置，然后其他的这样的数只能放它前面，再乘两个全排列就是答案

Code：

#include<bits/stdc++.h>
#define mod 1000000007
using namespace std;
inline int read(){
	int res=0,f=1;char ch=getchar();
	while(!isdigit(ch)) {if(ch=='-') f=-f;ch=getchar();}
	while(isdigit(ch)) {res=(res<<1)+(res<<3)+(ch^48);ch=getchar();}
	return res*f;
}
const int N=1e7+5;
inline int add(int x,int y){x+=y;if(x>=mod) x-=mod;return x;}
inline int dec(int x,int y){x-=y;if(x<0) x+=mod;return x;}
inline int mul(int x,int y){return 1ll*x*y%mod;}
inline void inc(int &x,int y){x+=y;if(x>=mod) x-=mod;}
inline void Dec(int &x,int y){x-=y;if(x<0) x+=mod;}
inline void Mul(int &x,int y){x=1ll*x*y%mod;}
inline int ksm(int a,int b){int res=1;for(;b;b>>=1,a=mul(a,a)) if(b&1) res=mul(res,a);return res;}
int fac[N],ifac[N],pri[N],cnt=0;
bool pt[N];
inline void init(int s,int n){
	if(s==1) for(int i=2;i<=n;i++) pt[i]=1;
	else{
		for(int i=s;i<=n;i++) if(!pt[i]){
			++cnt;
			for(int j=i<<1;j<=n;j+=i) pt[j]=1;
		}
	}
	fac[0]=ifac[0]=1;
	for(int i=1;i<=n;i++) fac[i]=mul(fac[i-1],i);
	ifac[n]=ksm(fac[n],mod-2);
	for(int i=n-1;i;i--) ifac[i]=mul(ifac[i+1],i+1);
}
inline int C(int n,int m){if(n<0 || m<0 || n<m) return 0;return mul(fac[n],mul(ifac[m],ifac[n-m]));}
int n,m=0;
int main(){
	int l=read(),r=read();n=r-l+1;
	init(l,r);m=cnt;
	if(l==1) m=1;
	int ans=0;
	for(int i=m;i<=n;i++)
		inc(ans,mul(mul(C(i-1,m-1),fac[m]),mul(fac[n-m],i)));
	cout<<ans;
	return 0;
}