[JXOI2018]游戏 线性筛

题面
我们定义不存在 d d 使得d|x,Ld,dx的数 x x 为好数。
那么对于一个排列p t(p) t ( p ) 显然是最后一个好数的出现位置。
线性筛求出好数个数 m m
枚举最后一个好数的出现位置i,那么满足的排列的个数有 (i1m1)m!(nm)! ( i − 1 m − 1 ) m ! ( n − m ) ! 个,其中 n=rl+1 n = r − l + 1
代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 10000010
#define ll long long
using namespace std;
const int mod=1000000007;
int L,R,pri[N],num,minp[N],cnt;
ll fac[N],ifac[N],ans;
ll ksm(ll a,int b){ll r=1;for(;b;b>>=1){if(b&1)r=r*a%mod;a=a*a%mod;}return r;}
void getpri(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!minp[i]) {pri[++num]=i;minp[i]=i;}
        for(int j=1;j<=num&&i*pri[j]<=n;j++)
        {
            minp[i*pri[j]]=pri[j];
            if(i%pri[j]==0) break;
        }
    }
}
int main()
{
    scanf("%d%d",&L,&R);
    fac[0]=1;
    for(int i=1;i<=R;i++)
        fac[i]=fac[i-1]*i%mod;
    ifac[R]=ksm(fac[R],mod-2);
    for(int i=R-1;i>=0;i--)
        ifac[i]=ifac[i+1]*(i+1)%mod;
    getpri(R);
    cnt=0;
    for(int i=L;i<=R;i++)
        if(minp[i]==0||i/minp[i]<L) cnt++;
    int len=R-L+1;ans=0;
    for(int i=cnt;i<=len;i++)
        ans=(ans+fac[len-cnt]*ifac[i-cnt]%mod*cnt%mod*fac[i-1]%mod*i)%mod;
    printf("%lld",ans); 
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值