BZOJ5323 & 洛谷4562：[JXOI2018]游戏——题解

https://www.luogu.org/problemnew/show/P4562

https://www.lydsy.com/JudgeOnline/problem.php?id=5323

（BZOJ有点卡常数过不去。）时限已经开大。

实际上我们只需要求出l～r区间内有多少数是满足不存在l～r内的数a使得i*a=这个数。

我们欧拉筛实际上就是一个数可以分解成的最大的两个数（其中一个是最大质数）的乘积，于是我们判断那个合数是否<l且这个数是否在l~r的区间内，如果满足则这个数就是我们要求的。

然后数学公式导一波就行了。

#include<cmath>
#include<queue>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1e7+;
const int p=1e9+;
int l,r,tot,su[N],cnt;
ll jc[N],inv[N];
bool he[N];
int qpow(int k,int n){
    int res=;
    while(n){
    if(n&)res=(ll)res*k%p;
    k=(ll)k*k%p;n>>=;
    }
    return res;
}
void Euler(int n){
    for(int i=;i<=n;++i){
    if(!he[i]){
        su[++tot]=i;
        if(l<=i&&i<=r)++cnt;
    }
    bool ok=;
    for(int j=;j<=tot;++j){
        int q=i*su[j];
        if(q>n)break;
        he[q]=;
        if(l<=q&&q<=r&&i<l)++cnt;
        if(i%su[j]==)break;
    }
    }
}
ll C(int n,int m){
    return jc[n]*inv[m]%p*inv[n-m]%p;
}
void init(int n){
    jc[]=;
    for(int i=;i<=n;++i)jc[i]=jc[i-]*i%p;
    inv[n]=qpow(jc[n],p-);
    for(int i=n-;i;i--)inv[i]=inv[i+]*(i+)%p;
    inv[]=;
}
int main(){
    scanf("%d%d",&l,&r);
    init(r);
    ll ans=;
    if(l==){
    for(int i=;i<=r;++i){
        ans=(ans+i*jc[r-])%p;
    }
    printf("%lld\n",ans);
    }else{
    Euler(r);
    int n=r-l+;
    for(int i=cnt;i<=n;++i){
        ans=(ans+i*C(n-cnt,i-cnt)%p*jc[i-]%p*jc[n-i]%p*cnt%p)%p;
    }
    printf("%lld\n",ans);
    }
    return ;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++