cf546 Round #304 Div2-D【埃氏筛拓展】

Date：2022.02.04
题意：给定一个$n=\frac{a!}{b!}$，每轮游戏选定一个$x$让$n=n/x$，求最大游戏次数。

思路①：阶乘分解。
看见阶乘我就想阶乘分解！$a!$和$b!$阶乘分解为$O(N)$，但每次都要知道$a$和$b$才能分解，因此总复杂度为$O(T\ast N)$铁t。那么再来考虑预处理一下$x\in [0!,N!]$，对所有$x$阶乘分解，显然是$O(N^2+T)$，一样糟糕。

思路②：打一个质数表，预处理一下。约$O(\frac{N}{logN})$个质数，对$N$个数分别枚举质数，看有多少个能成为它的质因数，直到枚举的质数不可能成为当前数的质因子。复杂度大概$O(\frac{N^2}{logN})$再乘上一个什么，显然也会t。还是交了一发，试了试预处理到$N==1e4$就寄了。
代码如下：

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 5e6+10;
typedef long long LL;
LL sum[N],cnt,prime[N];
LL t,n,m;
bool st[N];
void init()
{
    for(int i=2;i<N;i++)
    {
        if(!st[i]) {prime[cnt++]=i;}
        for(int j=0;prime[j]*i<N;j++)
        {
            st[prime[j]*i]=true;
            if(i%prime[j]==0) break; 
        }
    }
}
int main()
{
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    init();
    for(int i=1;i<N;i++)
    {
        LL res=0,x=i;
        for(int j=0;j<cnt&&x>=prime[j];j++)
            while(x%prime[j]==0) {res++;x/=prime[j];}
        sum[i]=sum[i-1]+res;
    }
    cin>>t;
    while(t--)
    {
        cin>>n>>m;LL res=0;
        cout<<sum[n]-sum[m]<<endl;
    }
    return 0;
}

思路③：我们只需要求一个数的所有质因子的次幂之和，不需要直到它每个数是多少等等信息。因此可以考虑埃氏筛，因为它有个性质：每个数只会被它的所有质因子筛掉。注意，含有多个质因子的数可能被筛掉多次，而不是一次。这也正是为什么它的复杂度大概是$O(N\ast loglogN)$而和线性筛法有所不同。因此我们可在埃氏筛过程中求出每个数的所有质因子的次幂之和，之后前缀和处理。
代码如下：

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 5e6;
int sum[N+10],d[N+10];//d[j]：j的所有质因子的次幂之和
int t,n,m;
void init()
{
    for(int i=2;i<=N;i++)
    {
        if(!d[i])//i可作为质因子
        {
            d[i]++;
            for(int j=i+i;j<=N;j+=i)//质因子i的倍数
            {
                int x=j;
                while(x%i==0) {d[j]++;x/=i;}
            }
        }
        sum[i]=sum[i-1]+d[i];
    }
}
int main()
{
    init();
    scanf("%d",&t);
    while(t--)
    {
        int a,b;scanf("%d%d",&a,&b);
        printf("%d\n",sum[a]-sum[b]);
    }
    return 0;
}