区间合数的最小公倍数

12 篇文章 0 订阅

题目链接
 

小红拿到了两个正整数 lll 和 rrr 。

她想知道 [l,r][l,r][l,r] 区间所有合数的最小公倍数是多少?由于答案可能过大,请对1000000007取模。

合数的定义:因子的个数不小于3的正整数。例如:9的因子有1、3、9这三个,所以9是合数。

输入描述:

两个正整数 lll 和 rrr ,用空格隔开。
数据范围:1≤l≤r≤300001\leq l \leq r \leq 300001≤l≤r≤30000

输出描述:

若区间 [l,r][l,r][l,r] 中没有任何合数,则输出-1。
否则输出区间 [l,r][l,r][l,r] 所有合数的最小公倍数,答案对1000000007取模。

示例1

输入

2 3

输出

-1

说明

区间 [2,3] 的两个数 2 和 3 都是素数,不是合数。

示例2

输入

4 6

输出

12

说明

区间[4,6]中有三个数:4、5、6,其中 5 是素数, 4 和 6 是合数,它们的最小公倍数是12。

示例3

输入

25000 30000

输出

187554966

思路:求最小公倍数其实就是求每个数的质数的幂的最大值,如:a=2^3 * 3 , b=2 * 3^3,最小公倍数就是:c=2^3 * 3^3。官方解释:区间所有合数的最小公倍数,也就是需要统计区间内每个素数出现的最高的幂。我们只需要统计出区间内的质数即可,计算它的最大幂,每个指数幂相乘即可。

完整代码:

#include <bits/stdc++.h>

using namespace std;

#define int long long
const int mod=1e9+7;

const int N=3e4+10;
int prem[N];
bool st[N];
int a[N];

int qmi(int a,int x)
{
    int res=1;
    while(x)
    {
        if(x&1)res=res*a%mod;
        a=a*a%mod;
        x>>=1;
    }
    return res;
}


signed main() {
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    int l,r;
    cin>>l>>r;
    int idx=0;
    st[1]=1;
    for(int i=2;i<=r;i++)
    {
        if(!st[i])prem[++idx]=i;
        for(int j=i+i;j<=r;j+=i)
        {
            st[j]=true;
        }
    }

    bool ok=false;
    for(int i=l;i<=r;i++)
    {
        if(!st[i])continue;
        ok=true;
        int pre=i;
        for(int j=1;j<=idx;j++)
        {
            int num=1;
            while(pre%prem[j]==0)pre/=prem[j],num=num*prem[j]%mod;
            a[j]=max(a[j],num);
        }
    }

    int ans=1;
    for(int i=1;i<=idx;i++)
    {
        //ans=ans*qmi(prem[i],a[i])%mod;
        ans=ans*a[i]%mod;
    }
    if(!ok)cout<<"-1"<<endl;
    else cout<<ans<<endl;
    return 0;
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值