牛客寒假算法基础集训营4-J.区间合数的最小公倍数

标签：数论

思路

首先这个题不是遍历每个合数挨着求最小公倍数不会真有人这么做吧

前置知识：要求多个数的最小公倍数，只需要统计每个数的素因子的次数，然后将次数最高的素因子乘起来（指素因子连同次数）。例如：求4，8，10的最小公倍数，4 = 2 * 2，8 = 2 * 2 * 2，10 = 2 * 5，素因子中2的最高次为3，5的最高次为1，那么最小公倍数就是2³ * 5¹ = 40

由于本题数据范围较小，所以用朴素的素数筛法也能过，剩下的就是暴力枚举l到r区间的数，是合数就统计其素因子的次数并更新最大值，最后全部乘起来即可

代码实现

#include<stdio.h>
#include<iostream>
#include<algorithm>
using namespace std;
const int mod=1e9+7;
const int N=3e5+5;
int l,r;
int num[N]; 
int prime[N];
void getprime(){       //埃氏筛
	for(int i=1;i<=30005;i++) prime[i]=1;
	for(int i=2;i<=30005;i++){
		for(int j=2;i*j<=30005;j++){
			prime[i*j]=0;
		}
	}
}
long long qpow(long long x,long long y){    //快速幂
	long long c=1;
	while(y){
		if(y&1) c=c*x%mod;
		x=x*x%mod;
		y>>=1;
	}
	return c;
}
int main(){
    cin>>l>>r;
    int book=0;
    getprime();
    for(int i=l;i<=r;i++){
    	if(!prime[i]){
    		book=1;           //记录区间有没有出现合数，判断是否输出-1
    		for(int j=2;j<=i/2;j++){   //枚举因子
    			int p=i,cnt=0;
    			if(prime[j]){
    				while(p%j==0){
    					cnt++;
    					p/=j;
					}
					num[j]=max(num[j],cnt);	//更新因子的最大次数
				}	
			}
		}
	}
	long long ans=1;
	for(int i=2;i<=r/2;i++) ans=ans*qpow(i,num[i])%mod;
	if(book==0) cout<<-1;
	else cout<<ans;
	return 0;
}