hdu 5317 RGCDQ 2015 Multi-University Training Contest 3

开始看到100W的数据量吓cry了，后来才知道2,3,5,7,11,13,17乘起来就已经900W了，所以f[i]<=7

预处理中先筛选出每个合数有几种质数因子。枚举每一个质数(a[i]!=0)将其倍数对应的合数统计量+1。因为质数分解不包括1，所以要从i=2开始枚举。

令f[i][j]表示数组中前i个数中f[i]=j的数的个数，即为前缀和，所以一个区间内f[x]=j的个数即为sub=f[R][j]-f[L-1][j]。如果sub>1，那么可以取区间中的两个数x，y使得f[x]=j，f[y]=j，对应的GCD=j。如果sub==1，需要和其他sub==1的f[x]比较，最后暴力特判GCD。

#include<iostream>
#include<stdio.h>
#include<cstdio>
#include<stdlib.h>
#include<vector>
#include<string>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<stack>
#include<queue>
#include<ctype.h>
#include<map>
#include<time.h>
#include<bitset>
#include<set>
#include<list>
using namespace std;
//hdu 5317
const int maxn=1000010;
const int N=8;
int T;
int L;
int R;
int ans=0;
int a[maxn];
int f[maxn][N];//f[i][j] the times that number j appears in the first i numbers
void initial()
{
    for(int i=2;i<maxn;i++)//i starts from 2 since there is no factor 1 when decomposition
    {
        if(a[i]==0)//之前没有被赋值，所以i一定是质数
        {
            a[i]=1;//每个i可以被分解成多少个质数
            for(int j=2*i;j<maxn;j+=i)
            {
                a[j]++;
            }
        }
    }
    for(int i=1;i<maxn;i++)
    {
        for(int j=1;j<N;j++)
        {
            f[i][j]=f[i-1][j];
        }
        f[i][a[i]]++;
    }
}
int quickgcd(int a, int b)
{
    if (a == 0) return b;
    if (b == 0) return a;
    if (!(a & 1) && !(b & 1))
        return quickgcd(a>>1, b>>1)<<1;
    else if (!(b & 1))
    return quickgcd(a, b>>1);
    else if (!(a & 1)) return quickgcd(a>>1, b);
    else return quickgcd(abs(a - b), min(a, b));
}
int main()
{
    freopen("input.txt","r",stdin);
    initial();
    scanf("%d",&T);
    for(int ca=1;ca<=T;ca++)
    {
        scanf("%d %d",&L,&R);
        ans=0;
        int cnt=0;
        for(int i=1;i<N;i++)
        {
            int sub=f[R][i]-f[L-1][i];
            if(sub>1)
            {
                ans=max(ans,i);
            }
            else if(sub==1)
            {
                a[cnt++]=i;
            }
        }
        for(int i=0;i<cnt;i++)
        {
            for(int j=i+1;j<cnt;j++)
            {
                ans=max(ans,quickgcd(a[i],a[j]));
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}