hdu6069 Counting Divisors 晒区间素数

/**
题目：hdu6069 Counting Divisors
链接：http://acm.hdu.edu.cn/showproblem.php?pid=6069
题意：求[l,r]内所有数的k次方的约数个数之和。

思路：
用(1+e1)*(1+e2)*...*(1+en)的公式计算约数个数。
素数筛出[l,r]内的素因子，然后直接计算结果。（一开始我用vector存起来，之后再处理，结果超时，
时间卡的很紧的时候，vector也会很占用时间。）

*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<map>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
using namespace std;
typedef pair<int,int> P;
typedef long long LL;
const int INF = 0x3f3f3f3f;
const int mod = ;
const int maxn = 1e6+;
bool is_prime_small[maxn];
vector<int> v[maxn];
int prime[maxn], z;
LL value[maxn], cnt[maxn];
void init()
{
    for(int i = ; i<=; i++) is_prime_small[i] = true;
    z = ;
    for(int i = ; i<=; i++){
        if(is_prime_small[i]){
            z++;
            prime[z] = i;
        }
        for(int j = ; j <= z; j++){
            if((LL)i*prime[j]>=){
                break;
            }
            is_prime_small[i*prime[j]] = false;
            if(i%prime[j]==) break;
        }
    }
}
void segment_sieve(LL a,LL b,LL k)  //[a,b)
{

    for(int i = ; i < b-a; i++) value[i] = i+a, cnt[i] = ;
    int num;
    for(int i = ; (LL)prime[i]*prime[i]<b&&i<=z; i++){
        for(LL j = max((LL)prime[i],(a+prime[i]-)/prime[i])*prime[i]; j < b; j+=prime[i]){
            num = ;
            while(value[j-a]%prime[i]==){
                num++; value[j-a]/=prime[i];
            }
            cnt[j-a] = cnt[j-a]*(num*k+)%mod;
        }
    }
}
LL solve(LL a,LL b,LL k)
{

    LL ans = ;
    for(int i = ; i < b-a; i++){
        if(value[i]>){
          cnt[i] = cnt[i] * (k+) % mod;
        }
        ans = (ans + cnt[i])%mod;
    }

    return ans;
}
int main()
{
    int T;
    cin>>T;
    init();
    LL l, r, k;
    while(T--){
        scanf("%lld%lld%lld",&l,&r,&k);
        if(l==r&&l==){
            printf("1\n"); continue;
        }
        int flag = ;
        if(l==){
            l++;
            flag = ;
        }
        segment_sieve(l,r+,k);
        printf("%lld\n",(solve(l,r+,k)+flag)%mod);
    }
    return ;
}