HDU 6069 Counting Divisors（数论）

Description

Input
第一行一整数T表示用例组数，每组用例输入三个整数l,r,k
(1<=T<=15,1<=l<=r<=1e12,r-l<=1e6,1<=k<=1e7)
Output
对于每组用例，输出答案
Sample Input
3
1 5 1
1 10 2
1 100 3
Sample Output
10
48
2302
Solution

Code

#include<cstdio>#include<cstring>using namespace std;typedef long long ll;const int maxn=1000011,mod=998244353;int mark[maxn],prime[maxn],res=0;void get_prime(int n=1000000){    for(int i=2;i<=n;i++)    {        if(!mark[i])mark[i]=prime[res++]=i;        for(int j=0;j<res&&prime[j]*i<=n;j++)        {            mark[i*prime[j]]=prime[j];            if(i%prime[j]==0) break;        }    }} int T,k,n,ans[maxn];ll l,r,a[maxn];int Solve(){    for(int i=0;i<res&&prime[i]*prime[i]<=r;i++)    {        int p=prime[i];        ll t=((l+p-1)/p)*p;        for(int j=t-l+1;j<=n;j+=p)        {            int num=0;            while(a[j]%p==0)a[j]/=p,num++;            ans[j]=(ll)ans[j]*((ll)k*num%mod+1)%mod;        }    }    int sum=0;    for(int i=1;i<=n;i++)    {        if(a[i]>1)ans[i]=(ll)ans[i]*(k+1)%mod;        sum=(sum+ans[i])%mod;     }    return sum;}int main(){    get_prime();    scanf("%d",&T);    while(T--)    {        scanf("%I64d%I64d%d",&l,&r,&k);        n=r-l+1;            for(int i=1;i<=n;i++)a[i]=l+i-1,ans[i]=1;        printf("%d\n",Solve());    }    return 0;}