HDU 5317 RGCDQ

题意
就是给你一个函数F(x),代表x的质因子数,让你求区间[L,R]范围內GDC(F(i),F(j))的最大值(L <= i < j <= R).
(题目链接)

分析

先考虑F(x)这个函数的特性.即对于很大的x,其分解质因数后得到的F(x)将很小.题目中x最大为1000000,则F(x)最大值为7(当F(x)为8时对应的x最小为2*3*5*7*11*13**17*19,其值超过1000000).
因此我们可以先预处理,求出[1,1000000]范围所有的F(x)值,用count[i][j]代表从1到i范围内F(x)的值为j的数的个数,这样只需要做差就能获得[L,R]范围内F(x)值的统计情况.预处理完成后,每次询问时先做差得出[L,R]区间内各种F(x)值的统计,然后计算所有i,j组合情况(不超过8*8种)对应的GDC(F(i),F(j))的值,取最大值即可.

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <vector>
#include <cmath>

using namespace std;

#define MAXN    1000000+10

#define PSIZE       100000
#define MAXX        9

int plist[PSIZE],pcount=0;

int gcd(int a,int b)
{
    if (a<b)
        return gcd(b,a);
    if (a%b==0)
        return b;
    else
        return gcd(b,a%b);
}

int prime(int n)
{
    int i;
    if ((n!=2 && n%2==0) ||(n!=3 && n%3==0) ||(n!=5 && n%5==0) ||(n!=7 && n%7==0) 
        ||(n!=11 && n%11==0) ||(n!=13 && n%13==0))
        return 0;
    for (i=0;plist[i]*plist[i]<=n;i++)
            if (!(n%plist[i]))
                return 0;
    return n>1;
}

void initprime()
{
    int i;
    for (plist[pcount++]=2,i=3;i<PSIZE;++i)
        if (prime(i))
            plist[pcount++]=i;
}

int prime_factor(int n)
{
    int cnt=0;
    int n2=sqrt((double)n);
    for (int i=0;n>1&& plist[i]<=n2;++i)
    {
        if (n%plist[i]==0)
        {
            cnt++;
            while(n%plist[i]==0)
                n/=plist[i];
        }
    }
    if (n>1)
        cnt++;
    return cnt;
}

int count[MAXN][MAXX+1];

int main()
{
    int T;
    int L,R;
    initprime();
    int k;
    for (int j=0;j<=MAXX;j++)
    {
            count[1][j]=0;
    }
    for (int i=2;i<MAXN;++i)
    {        
        for (int j=0;j<=MAXX;j++)
        {
            count[i][j]=count[i-1][j];
        }
         k=prime_factor(i);
        count[i][k]++;
    }
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d %d",&L,&R);
        int ans=0;
        int a,b;
        int g;
        for (int i=1;i<=MAXX;++i)
        {
            for  (int j=1;j<=MAXX;++j)
            {
                 a=count[R][i]-count[L-1][i];
                 b=count[R][j]-count[L-1][j];
                 if (i!=j)
                 {
                    if ((a>0)&&(b>0))
                    {
                        g=gcd(i,j);
                        ans=g>ans?g:ans;
                    }
                }
                else
                {
                    if (a>=2)
                    {
                        g=i;
                        ans=g>ans?g:ans;
                    }
                }
            }
        }             
        printf("%d\n",ans);   
    }
    return 0;
}