UPC 6887 游戏(组合计数+线性筛)

注意到我们只关心不能被其他数筛掉的数 , 也就是不存在取值在 [l,r][l,r] 的约数的数.
最后一个满足这个条件的数出现的位置就是游戏的轮数.
我们只关心满足这个条件的数的数量 , 设为 cc , 那么枚举最后一个出现的位置 ii 来算贡献.

我们知道,如果一个数在区间内没有除了它自己以外的因数,那么这个数就必须要取,不然不能满足所有数都被取过。再然后,如果一个数在区间内有除了它自己以外的因数,那它就没有必要取了,也就是说取这个数是多余的,因为它的某一个因数是一定要取的,那么在它被取的时候这个数就已经被取了。
由这两个结论,我们得出一个排列p的次数t(p)就等于那些在区间内没有除了自己以外因数的数,在排列p中最后一个出现的位置。
现在有两个问题,
一个是如何求一个数是不是在区间内没有除自己以外的因数,
另一个是如何统计所有排列p的t(p)之和。
对于第一个问题,我们显然要找到该数除了自己以外的最大的因数,如果这个因数不在区间中,就表示这个数在区间中没有除自己以外的因数。显然这个因数等于原数除以它最小的质因子,而区间[l,r]中所有数的最小的质因子可以通过线性筛来O(n)算出,于是我们解决了这一个问题。
对于第二个问题,令所有在区间中没有除自己以外因数的数的个数为tot,考虑求满足t(p)=i(tot≤i≤r−l+1)的排列数目,显然此时满足条件的排列数为Ci−1tot−1⋅tot!⋅(n−tot)!,可以化简为tot⋅(n−tot)!⋅(i−1)!(i−tot)!。于是有:
ans=tot⋅(n−tot)!⋅∑r−l+1i=toti⋅(i−1)!(i−tot)!
O(n)预处理阶乘和阶乘的逆元后即可O(n)算出这个式子,那么我们就解决了这一题。
学了一个求阶乘和阶乘逆元的常数很小的方法先简单地求出阶乘,再快速幂算出n!的逆元,然后倒着推出阶乘的逆元,即 inv( (i−1) !)=inv( i ! )⋅i mod mod,这样可以极大地减小取模的常数,在这种卡常题中特别有用。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=1000000007;
int l,r,prime[2000010],low[10000010];
ll fac[10000010],inv[10000010],ans=0;

void calc_prime()//线性筛
{
    prime[0]=0;
    low[1]=1;
    for(int i=2;i<=r;i++)//从二开始
    {
        if (!low[i])//low里面存的是零的时候执行
        {
        prime[++prime[0]]=i;//存的是下标+1,也可以理解成编号+1
            low[i]=i;//里面存的是下标
        }
        for(int j=1;j<=prime[0]&&i*prime[j]<=r;j++)
        {
            low[i*prime[j]]=prime[j];
            if (!i%prime[j]) break;
        }
    }
}

ll power(ll a,ll b)
{
    ll s=1,ss=a;
    while(b)
    {
        if (b&1) s=s*ss%mod;
        ss=ss*ss%mod;b>>=1;
    }
    return s;
}

int main()
{
    scanf("%d%d",&l,&r);//标号范围
    calc_prime();//线性筛出所有质数

    fac[0]=fac[1]=1;
    inv[0]=inv[1]=1;
    for(ll i=2;i<=r;i++)
        fac[i]=fac[i-1]*i%mod;//求阶乘
    inv[r]=power(fac[r],mod-2);
    for(ll i=r;i>=1;i--)
        inv[i-1]=inv[i]*i%mod;//求阶乘逆元

    int tot=0;
    for(int i=l;i<=r;i++)
        if (low[i]==0||i/low[i]<l) tot++;
   if (l==1) tot++;//必须有这一行
    for(ll i=tot;i<=r-l+1;i++)
        ans=(ans+fac[i]*inv[i-tot]%mod)%mod;
    printf("%lld",ans*tot%mod*fac[r-l+1-tot]%mod);
   /*  for(int i=tot;i<=r-l+1;i++)
        ans=(ans+fac[r-l+1-tot]*i*fac[i-tot]%mod*tot%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、付费专栏及课程。

余额充值