【hdu 6069】 Counting Divisors 【思维+数论】

Problem Description
In mathematics, the function d(n) denotes the number of divisors of positive integer n.

For example, d(12)=6 because 1,2,3,4,6,12 are all 12’s divisors.

In this problem, given l,r and k, your task is to calculate the following thing :

(∑i=lrd(ik))mod998244353

Input
The first line of the input contains an integer T(1≤T≤15), denoting the number of test cases.

In each test case, there are 3 integers l,r,k(1≤l≤r≤1012,r−l≤106,1≤k≤107).

Output
For each test case, print a single line containing an integer, denoting the answer.

Sample Input
3
1 5 1
1 10 2
1 100 3

Sample Output
10
48
2302

Source
2017 Multi-University Training Contest - Team 4
分析
1>首先明确 一个数的因数个数为多少 ：
因为 n=p1^r1*p2^r2….pn^rn ;所以 因子个数为（r1+1）(r2+1)(rn+1) ， 其k次方的因子个数为（r1*k+1）（r2*k+1)(rn*k+1) ，所以我们就要找素数来解。
2>剩下的不好描述看代码吧。

代码

#include<bits/stdc++.h>#define LL long longusing namespace std;const int MAXN = 1e6+500; const int MAXM = 1e5;const LL mod =998244353;bool su[MAXN];int prm[MAXN],sz=0;void init(){ // 这个素筛还是快一些      for(int i=2;i<MAXN;++i){        if(!su[i]) prm[sz++]=i;        for(int j=0;j<sz;++j){            int t=i*prm[j];            if(t>=MAXN) break;            su[t]=1;            if(i%prm[j]==0) break;         }     }}LL arr[MAXN],ans[MAXN];LL size;void solve(LL l,LL r,LL k){    for(LL i=0;i<sz;i++){        LL pos=(l+prm[i]-1)/prm[i]*prm[i]; // 这里通过计算可以求得 大于l的最小的可以被prm[i]整除的数字        pos-=l;//得到位置        while(pos<size){// 遍历区间，将prm[i]的倍数都筛掉             LL cnt=0;            while(arr[pos]%prm[i]==0) { cnt++; arr[pos]/=prm[i];}            ans[pos]=(ans[pos]*(cnt*k%mod+1))%mod;//更新个数            pos+=prm[i];//优化，不断的加prm[i]一定是其倍数        }    }//  puts("---------");    LL sum=0;    for(LL i=0;i<size;i++){        if(arr[i]==1) ;// 说明这个数字可以被完全分解为1-1e6之内的素数        //否则就说明 其剩下的一定是 一个素数，且次幂为1         else ans[i]=(ans[i]*(1*k%mod+1))%mod;        sum=(sum+ans[i])%mod;    }    printf("%lld\n",sum%mod);}int main(){    init();    int t;cin>>t;    while(t--){        LL l,r,k;        scanf("%lld%lld%lld",&l,&r,&k);         size=0;        for(LL i=l;i<=r;i++) {            ans[size]=1;            arr[size++]=i;        }        solve(l,r,k);     }    return 0;}