51NOD 1434 区间LCM(素数筛)

传送门

1434 区间LCM

题目来源: TopCoder

基准时间限制:1 秒 空间限制:131072 KB 分值: 40 难度:4级算法题 收藏 关注

一个整数序列S的LCM(最小公倍数)是指最小的正整数X使得它是序列S中所有元素的倍数,那么LCM(S)=X。

例如,LCM(2)=2,LCM(4,6)=12,LCM(1,2,3,4,5)=60。

现在给定一个整数N(1<=N<=1000000),需要找到一个整数M,满足M>N,同时LCM(1,2,3,4,…,N-1,N) 整除 LCM(N+1,N+2,….,M-1,M),即LCM(N+1,N+2,….,M-1,M)是LCM(1,2,3,4,…,N-1,N) 的倍数.求最小的M值。

Input

多组测试数据,第一行一个整数T,表示测试数据数量,1<=T<=5

每组测试数据有相同的结构构成:

每组数据一行一个整数N,1<=N<=1000000。

Output

每组数据一行输出,即M的最小值。

Input示例

3

1

2

3

Output示例

2

4

6

解题思路:

假设有质数K,可以求出t,使得K的t次方刚好小于等于m(K^t<=m)那么lcm(1,2,…,m)分解质因数中一定而且最多有t个质数K连乘,其实我们可以就是求的质数的最高次幂,这样就可以很快地吧lcm(1,2,…,m)分解质因数那么怎么把 lcm(n+1,n+2,…,m)分解质因数呢仍然假设质数K,可以求出最大的t,以及一个常数c(1<=c < K),使得 n+1<=c*K^t<=m

那么lcm(n+1,n+2,…,m)分解质因数中一定而且最多有t个质数K连乘。比如说质数3,n=16,m=22,可以求的c=2,t=2,即17<=2*3^2=18<=22,这样lcm(17,18,19,20,21,22)中最多有2个质数3连乘既然知道怎么求lcm(n+1,n+2,..,m)和lcm(1,2,..,m)了来探讨一下怎么求最小的m吧我们想让这两个lcm分解质因数后完全一样,也就是说连乘的质数个数也完全相等。也就是说,对于每个质数K都可以满足,存在c和最大的t 使得n+1<=c*K^t<=m对于大于n小于m的质数,我们假设是P,那么一定n+1<=P<=m,一定可以满足条件所以我们就只看小于等于n的质数就可以了因为要使每个小于N的质数K,都存在c和最大的t 使得n+1<=c*K^t<=m,我们枚举每一个质数,并求得c和t,使得刚好c*K^t>=n,

答案就取最大的c*K^t,即 max( c*K^t )这样lcm(1,2,…,m)和lcm(n+1,n+2,…,m)的分解质因数后均至少有t个质数K。如果最终m有 k^(t+1)<=m,那么这个K^(t+1)>n一定成立,故仍满足条件

代码:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
using namespace std;
typedef long long LL;
const int MAXN = 1e6+5;
LL p[MAXN];
bool prime[MAXN];
int k;
void isprime()
{
    k = 0;
    memset(prime, 0, sizeof(prime));
    prime[0] = prime[1] = 1;
    for(LL i=2; i<MAXN; i++)
    {
        if(!prime[i])
        {
            p[k++] = i;
            for(LL j=i*i; j<MAXN; j+=i)
                prime[j] = 1;
        }
    }
}
int main()
{
    isprime();
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int n;
        scanf("%d",&n);
        int ans = 2;
        for(int i=2; i<=n; i++)
        {
            if(!prime[i])///素数
            {
                int sum = i;
                while(sum <= n/i)
                    sum *= i;
                ans = max(ans, (n/sum+1)*sum);
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}
时间: 2024-08-04 17:24:57

51NOD 1434 区间LCM(素数筛)的相关文章

51nod 1434 区间LCM 素数

1434 区间LCM 题目来源: TopCoder 基准时间限制:1 秒 空间限制:131072 KB 分值: 40 难度:4级算法题 一个整数序列S的LCM(最小公倍数)是指最小的正整数X使得它是序列S中所有元素的倍数,那么LCM(S)=X. 例如,LCM(2)=2,LCM(4,6)=12,LCM(1,2,3,4,5)=60. 现 在给定一个整数N(1<=N<=1000000),需要找到一个整数M,满足M>N,同时LCM(1,2,3,4,...,N- 1,N) 整除 LCM(N+1,N

51nod 1434 区间LCM (质因数分解)

分析:考虑从1到n所有数的质因数分解,记录每个质数的最高次数,同理从n+1循环到2n,如果循环到m时每个质因子的次数都不低于所记录的,则跳出循环,结果即为m.先预处理质数,复杂度为O(nlongn). 1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<algorithm> 5 using namespace std; 6 const int maxn=1000005;

light_oj 1197 区间素数筛

light_oj 1197 区间素数筛 M - Help Hanzo Time Limit:2000MS     Memory Limit:32768KB     64bit IO Format:%lld & %llu Submit Status Practice LightOJ 1197 Description Amakusa, the evil spiritual leader has captured the beautiful princess Nakururu. The reason

区间素数筛模版

区间素数筛模版 筛出区间[a,b]的素数.(b-a<=10000,1<=a<=b<=2^31) 存在P中,素数个数即为P的size(). ll a,b; bool isprime[maxn]; vector<ll> prime; bool isP[maxn]; vector<ll> P; void play_prime() { memset(isprime,1,sizeof(isprime)); isprime[1]=0; for(int i=2;i<

Help Hanzo (素数筛+区间枚举)

Help Hanzo 题意:求a~b间素数个数(1 ≤ a ≤ b < 231, b - a ≤ 100000).     (全题在文末) 题解: a~b枚举必定TLE,普通打表MLE,真是头疼.. b - a ≤ 100000 是关键. 类似素数筛的方法: 1.初始化vis[]=0 ; 2.素数的倍数vis[]=1; 3.  b较小时,素数筛解决   b很大时,素数筛的vis[]会MLE,此时用vis2[i-a]保存vis[i]就不会MLE 了.. #include<iostream>

Light oj 1197 - Help Hanzo (素数筛技巧)

题目链接:http://lightoj.com/volume_showproblem.php?problem=1197 给你a和b求a到b之间的素数个数. 先在小区间素数筛,大区间就用类似素数筛的想法,把a到b之间不是素数的标记出来.因为b-a最多1e5的大小,所以每组数据的时间复杂度最多就o(1e5 log1e5). 1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 using names

HDU 4497 素数筛,合数分解

题意:给你两个数,G和L ,它们是一组数的GCD(最大公约数)和LCM(最小公倍数),求出满足条件的组合数,每个组合有三个数,排序不同,也算不同组合. L : p1^t1 * p2^t2 ... * pi^ti G: q1^s1 * q2^s2... * qi^si (pi和qii都是素数ii) G(最大公约数)里出现的所有数都要在L(最小公倍数)里出现,或者说L mod G=0,有人要问为什么了?其实仔细想想就清楚了,如果a,b,c的最小公倍数肯定可以被a,b,c任意一个整除,而他们的最大公约

【模版】素数筛, 最大公约数(辗转相除法),最小公倍数

素数筛: //数除了{2,3,5}为素数,其他的数可以写成6N,6N+1,6N+2,6N+3,6N+4,6N+5 N>=1 可以表示全部的数 //6N,6N+2,6N+4都为偶数,不是素数,6N+3 == 3(2N+1) 不是素数,那么就只筛6N+1和6N+5就可以了 int prime[1000000]={2,3,5}; void is_prime() { int i,j; int flag=0;//标记 int gcd=2; int k=3;//素数的下标,因为已有三个,所以下一个出现的宿舍

[luoguP1835] 素数密度_NOI导刊2011提高(04)(素数筛)

传送门 数据辣么大,怎么搞?(L≤R≤2147483647) 注意到R-L≤1000000 所以可以直接筛R-L区间内的数, 但是需要用已知的小的素数筛, R-L区间内的大部分数肯定能用较小的素数筛去,但是还有一些较大的数,可能等于两个大质数的乘积,没法被筛去. 但是又注意到,数据最大才10位,也就是说我们只需要用位数<=5的素数筛就可以了,所以先预处理出来,直接筛就ok了. #include <cstdio> #define N 1000001 #define max(x, y) ((