[素筛] hdu 5317 RGCDQ

题目分析

    令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,故 fi,max=7
    1~7这七个数字是互质的,因此
gcd(fi,fj)={1fififj,fi=fj.
    然后用前缀和求一下[l,r]间f[i]各个值的个数,找最大的f[i]即可。

代码

#include
#include
#include
#include
#include
#include
#include
#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;iif(f[i]) continue;
        f[i]=1;
        for(int j=(i<<1);j//    for(int i=2;i<=199;i++) cout<
    for(int i=2;ifor(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;
}

你可能感兴趣的:(acm,hdu,数论)