HDU6069---Counting Divisors（2017多校联赛：数学）

【题目来源】：http://acm.hdu.edu.cn/showproblem.php?pid=6069；
【题意】
求所示式子的答案。
【思路】
这道题对于我挺绕的，当时比赛的时候想了一二十分钟，想了出来，但是实现的时候由于代码水平有限，以至于今天打了一遍过了，依旧找不出昨天的错误，思路基本和题解一致，用到了几个数论小知识，比如：求一个数的所有因子的数目可以通过素数的次幂进行相乘求解，然后就是一个模拟了。
【代码】

//先统计素因子数目，记录，最后相加。#include<cstdio>#include<cstring>#include<algorithm>using namespace std;typedef long long LL;const int maxn=1e6+10;const int mod=998244353;LL prime[maxn];bool check[maxn];LL p[maxn],q[maxn];LL l,r,k;int tot;void pri(){    check[1]=1;    for(int i=2; i<maxn; i++)    {        if(!check[i])prime[tot++]=i;        for(int j=0; j<tot&&i*prime[j]<maxn; j++)        {            check[i*prime[j]]=1;            if(i%prime[j]==0)break;        }    }}int main(){    int T;    tot=0;    pri();    scanf("%d",&T);    while(T--)    {        scanf("%lld%lld%lld",&l,&r,&k);        for(int i=0;i<=r-l;i++)        {            p[i]=l+i;            q[i]=1;        }        for(int i=0; i<tot; i++)        {            if(prime[i]*prime[i]>r) break;            LL L=l/prime[i];            LL j=prime[i]*(L+(l%prime[i]!=0));            for(;j<=r;j+=prime[i])            {                int tmp=0;                while(p[j-l]%prime[i]==0)                {                    tmp++;                    p[j-l]/=prime[i];                }                q[j-l]*=(k*tmp+1);                q[j-l]%=mod;            }        }        LL sum=0;        for(LL i=l; i<=r; i++)        {            if(p[i-l]!=1)            {                q[i-l]*=(k+1);                if(q[i-l]>=mod)                    q[i-l]%=mod;            }            sum+=q[i-l];            if(sum>=mod)                sum%=mod;        }        printf("%lld\n",sum);    }}