求1~n的最小公倍数(LCM)

给你一个整数n,求出[1,n]所有数的最小公倍数,并对1000000007取模(1<=n<=10000000)

要求1~n的最小公倍数,首先要知道,互质的两个数的最小公倍数就是两个相乘,之后在考虑合数。

我们先使用筛法筛出质数,合数中又可以分为两类,一类是合数等于一个质数的k次方,也就是说他和其中一个质数的最小公倍数等于他本身,另一类则可以用质数和质数相乘表示,化简可以变成1。

举个例子,1-10的最小公倍数,lcm(1 2 3 4 5 6 7 8 9 10) 可以化为 lcm(1 1 1 1 5 1 7 8 9 1)相乘。

首先6是2和3的最小公倍数,所以2和3同时存在的时候,6就可以去掉,保留2和3即可。

10也同理,2和5就可以i表示10。

然后再来某个数字是某个素数的k次方的,当该数字存在时,他的那个素数因子就可以被去掉,例如4,8。因为8是2的3次方,4是2的2次方,当4存在的时候可以直接把2删去,当8存在的时候可以把4删去。

故而我们得到lcm[1,10]的最小公倍数为5*7*8*9 = 2520。

因此我们只需要用n以内所有质数的小于n的k次方那个数来替换质数即可。

#include<bits/stdc++.h>
using namespace std;
#define maxnum 10000051
bool Isprim[maxnum];
int  primes[maxnum];

int main(){
    unsigned long long nums = 1, n = 0,count_prime = 0,temp = 0;
    scanf("%lld",&n);
    memset(Isprim,1,sizeof(Isprim));
    Isprim[1] = false;
    for(unsigned long long i = 2;i<=n;i++){
        if(Isprim[i]){
            primes[count_prime] = i;
            count_prime++;
        }
        for(unsigned long long j = 0;j<count_prime&&(temp = primes[j]*i)<=n;j++){
            Isprim[temp] = false;
            if(i%primes[j]==0){
                break;
            }
        }
    }
 
    for(unsigned long long i = 0;i<count_prime;i++){
        unsigned long long emm = 1;
        while(emm<=n){//通过循环让emm从1开始不断乘以质数直到大于输入的数字n
            emm *= primes[i];//这时候让emm除以该质数就可以得到小于等于n的数
        }
        emm /= primes[i];
        primes[i] = emm;
    }
    for(unsigned long long i = 0;i<count_prime;i++){
        nums = nums*primes[i]%1000000007;
    }
    printf("%llu\n",nums);
    return 0;
}

这是萌新第一次写博客,不好请见谅,谢谢。

  • 7
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
1-n的最小公倍数可以使用辗转相除法和穷举法两种方法。 使用辗转相除法的步骤如下: 1. 初始化最小公倍数为1。 2. 从2开始遍历到n,对每个数执行以下操作: a. 判断当前数与最小公倍数的最大公约数是否为1,如果是,则将最小公倍数乘以当前数。 b. 如果最大公约数不是1,则将最小公倍数除以最大公约数,再乘以当前数。 3. 返回最小公倍数的值。 使用穷举法的步骤如下: 1. 初始化最小公倍数为n。 2. 从n-1开始递减到1,对每个数执行以下操作: a. 如果当前数与最小公倍数取模不为0,将最小公倍数加上n。 3. 返回最小公倍数的值。 两种方法都能够正确得1-n的最小公倍数,选择哪种方法主要取决于具体的情况和n的大小。 举个例子,如果我们要1-10的最小公倍数,可以使用辗转相除法: 1. 初始化最小公倍数为1。 2. 遍历2到10的每个数: a. 对于2,判断1和2的最大公约数是否为1,是的话,将最小公倍数乘以2,得到2。 b. 对于3,判断2和3的最大公约数是否为1,是的话,将最小公倍数乘以3,得到6。 c. 对于4,判断6和4的最大公约数是否为1,不是的话,将最小公倍数除以2,再乘以4,得到12。 d. 对于5,判断12和5的最大公约数是否为1,是的话,将最小公倍数乘以5,得到60。 e. 对于6,判断60和6的最大公约数是否为1,不是的话,将最小公倍数除以6,再乘以6,得到60。 f. 对于7、8、9、10,都与60的最大公约数为1,所以将最小公倍数分别乘以7、8、9、10,得到最终结果2520。 4. 返回最小公倍数2520。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* [1~n的最小公倍数LCM)](https://blog.csdn.net/coldwind902/article/details/124548995)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *3* [最大公约数最小公倍数n-s盒图](https://download.csdn.net/download/d52370/11026280)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值