[素筛] hdu 5317 RGCDQ_hdu5317-CSDN博客

本文链接：https://blog.csdn.net/gdymind/article/details/52094107

博客介绍了如何利用素筛算法解决HDU 5317 RGCDQ问题，重点在于计算质因数种类数并在过程中求得前缀和，从而找出最大质因数种类数。分析了题目需求和代码实现细节，指出关键的限制条件和优化策略。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目分析

    令f[i]为i的质因数的种类数，则f[i]可以在求素筛的过程中顺便求出（见代码25~30行），注意for循环中约束条件是“i < maxn”而非“i*i < maxn”，因为后面的质数的f值要初始化为1。
    注意到2*3*5*7*11*13*17=510510<1e6,而2*3*5*7*11*13*17*19=9699690>1e6,故 $f_{i,max}=7$
    1~7这七个数字是互质的，因此
$gcd(f_i,f_j)= \begin{cases} 1& \text{$f_i \not= f_j$},\\ f_i& \text{$f_i = f_j$}. \end{cases}$
    然后用前缀和求一下[l,r]间f[i]各个值的个数，找最大的f[i]即可。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<map>
#include<algorithm>
#define lson l,mid,cur<<1
#define rson mid+1,r,cur<<1|1
typedef long long ll;
using namespace std;
const int maxn = 1e6+7;
int f[maxn];
int sum[maxn][8];
int l,r;
bool check(int x)
{
    if(sum[r][x]-sum[l-1][x]>=2) return true;
    return false;
}
int main()
{
    memset(f,0,sizeof f);
    memset(sum,0,sizeof sum);
    for(int i=2;i<maxn;i++)
    {
        if(f[i]) continue;
        f[i]=1;
        for(int j=(i<<1);j<maxn;j+=i) f[j]++;
    }
//    for(int i=2;i<=199;i++) cout<<i<<':'<<f[i]<<endl;
    for(int i=2;i<maxn;i++)
    {
        for(int j=1;j<=7;j++) sum[i][j] = sum[i-1][j];
        sum[i][f[i]]++;
    }

    int T;
    cin>>T;
    int i;
    while(T--)
    {
        scanf("%d%d",&l,&r);
        for(i=7;i>=1;i--)
            if(check(i)) break;
        printf("%d\n",i);
    }
    return 0;
}