素因子数GCD

题目描述
本题首先我们来定义素因子数函数F(x)，它表示x的素因子有多少种。例如F(2)=1, 因为2的素因子只有2，F(10)=2, 因为10的素因子有2和5，F(12)=2，因为12的素因子也是2和5。

给定区间[L, R]，我们可以任意找两个数i, j (L <= i < j <= R)，我们可以得到F(i)和F(j)。对于F(i)和F(j)，我们还可以求他们的最大公因数，即GCD( F(i), F(j) )。那么问题来了，区间[L, R]中，最大的GCD( F(i), F(j) )是多少？

例如区间[6, 15]，答案为2，可以选择F(6)和F(15)，其中F(6)=2, F(15)=2。

输入
首先输入测试组数T

对于每组测试数据，输入L, R

输出
对于每组测试数据，输出maxGCD( F(i), F(j) )

样例输入
3
2 3
3 5
6 15
样例输出
1
1
2
提示
【数据规模和约定】

1<=T<=1000000

2<= L < R <=1000000

题解
发现一百万以内，一个数的质因子个数不会超过7个，利用线性筛求出一百万以内每个数的质因子个数，求前缀和O（1）求解即可、

代码

#include<bits/stdc++.h>
#define sqr(x) (x*x) 
#define eps 1e-8
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
bool flag[1000005];
int tot,p[1000005],f[1000005],ans[8],sum[1000005][8];
void pre()
{
    for (int i=2;i<=1000000;i++)
    {
        if (!flag[i]){p[++tot]=i;f[i]=1;}
        for (int j=1;j<=tot&&p[j]*i<=1000000;j++)
        {
            flag[i*p[j]]=1;
            f[i*p[j]]=f[p[j]]+f[i];
            if (i%p[j]==0){f[i*p[j]]-=f[p[j]];break;}
        }
    }
    for (int i=2;i<=1000000;i++)
    {
        for (int j=1;j<=7;j++) sum[i][j]=sum[i-1][j];
        sum[i][f[i]]++;
    }
}
int main()
{
    pre();
    int Case=read();
    while (Case--)
    {
        int L=read(),R=read();ans[1]=2;
        for (int i=2;i<=7;i++) ans[i]=sum[R][i]-sum[L-1][i];
        ans[2]=ans[2]+ans[4]+ans[6];
        ans[3]=ans[3]+ans[6];
        for (int i=7;i;i--) if (ans[i]>=2){printf("%d\n",i);break;}
    }
    return 0;
}