hdu6069 Counting Divisors(质数分解+技巧)

d(n)里n的取值最大1e12，那么要判断它是不是素数只需计算到sqrt(n),所以打个1e6的素数表，如果它是个合数，就可以用打好的表里的素数进行分解，如果是个素数，
   那么它对最后结果的贡献就是k+1。
遍历和计算中的技巧：
数组开到1e6，从l到r拿0到r-l对应。
中间遍历时先遍历素数表，假设现在遍历到prime[i],先找到l~r第一个能整除prime[i]的数a[j],以此开始，每次下标加prime[i],在分解a[j]时一直除以prime[i]直到除不尽，这时已经消除了prime[i]因子，把分解的幂值加进sum[j]。

具体看代码（有注释）

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int mod = 998244353;
const int maxn=1000009;
bool Isprime[maxn];
LL prime[maxn];
LL a[maxn],sum[maxn];
int cnt;
LL l,r,k;
void Prime()
{
    cnt = 0;
    memset(Isprime, true, sizeof(Isprime));
    Isprime[1] = false;
    for(int i = 2 ; i < maxn ; i++)
    {
        if(Isprime[i])
        {
            prime[cnt++] = i;
            for(int j = 2 ; i * j < maxn ;j++)
                Isprime[i * j] = false;
        }
    }
}

int main()
{
    Prime();
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld%lld%lld",&l,&r,&k);
        for(int i=0;i<=r-l;i++)
        {
            sum[i]=1;///sum是要做乘法的，初值为1
            a[i]=i+l;
        }
        for(int i=0;i<cnt;i++)
        {
            LL fir=(l/prime[i]+(l%prime[i]?1:0))*prime[i];///fir是从l开始第一个能整除prime[i]的
            for(LL j=fir;j<=r;j+=prime[i])///每个j都能整除prime[i]
            {
                LL res=0;
                while(a[j-l]%prime[i]==0){res++;a[j-l]/=prime[i];}///一直除下去,让数变小，找到prime[i]的幂值
                sum[j-l]=(sum[j-l]*((res*k+1)%mod))%mod;///sum[j-l]乘一下，取模
            }
        }
        LL ans=0;
        for(int i=0;i<=r-l;i++)
        {
            if(a[i]!=1)sum[i]=(sum[i]*(k+1))%mod;///经过一系列除法没除到1，说明1e6内的素数分解不了，那么它就是一个大素数
            ans=(ans+sum[i])%mod;
        }
        printf("%lld\n",ans);
    }

    return 0;
}