I - Counting Divisors（区间素数筛）

这个题以前我就遇见过，结果不会写，现在终于搞懂了。


这道题整体思路就是首先我们需要知道一个约数定理：

这里的p1,p2…pn为质数，d(x)表示约数个数。知道这个定理之后我们可以推到出关于幂的结论：

多个x相乘 好理解吧，指数运算高中的；
d(x^k)=k+1如果x为素数很好理解带入公式(1不是素数)
这样就解决了一半了，下面一半我们就应该找每个的d(x^k)了，这里就有问题了，如果直接用素数去枚举每个数那肯定TLE（超时）。所以有人给我提示了区间素数筛，New konwledge！！！ 但是怎么筛呢？这点就需要一个数组num来存区间l–r的值然后在筛的时候就可以用来计算，那么每个的d(x^k)了，这也不可能一个一个算，所以这点用筛法 优越性就体现出来了；
先定素数，然后把每一个都扫一遍，变化后的值放在num中并且同时要记录他的d()值
代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll const mod = 998244353;
#define N 1000005
bool isPrime[N];
int prime[N];//用来存素数 
int tot;//表示素数的总数 
 
void get_prime()//得到素数表prime[素数]=1; 
{
   for(int i=2;i<N;i++){
    if(!isPrime[i]){
        prime[tot++] = i;
    }
    for(int j=0;j<tot;j++){
        if(prime[j] * i >= N) break;
        isPrime[prime[j] * i] = 1;
        if(i%prime[j] == 0) break;
    }
   }
}
ll num[N]; //num[i] 表示 i+l 被所有素数约分后剩下的数
ll ans[N]; //ans[i] 表示 i 这一位的答案，最后需要累加
int main()
{
    get_prime();
    int t;
    scanf("%d",&t);
    while(t--){
        ll l,r,k;
        ll Ans = 0;
        scanf("%lld%lld%lld",&l,&r,&k);
       if(l==1) Ans++,l++;  // l为1,则l从2开始 因为1不是素数
       for(ll i=l;i<=r;i++) ans[i-l] = 1,num[i-l] = i;
        for(int i=0;i<tot;i++){
            if(prime[i] > r) break;
            ll w = l % prime[i];//为了能从第一个能被i所对应的素数整除开始   降低时间复杂度
            if(w) w = prime[i] - w;
           for(ll j=l+w;j<=r;j+=prime[i]){  //枚举 l ~ r 之间所有能整除prime[i]的数
            ll cot = 0;
            while(num[j-l]%prime[i]==0){    //查询 j 有多少个 素因子 prime[i]
                cot++;
                num[j-l] /= prime[i];
            }
            if(cot) ans[j-l] *= (cot * k + 1) % mod ; //ans[j] 的答案 乘以 cot 个 prime[i]的贡献
            ans[j-l] %= mod;
           }
        }
 
        for(ll i=l;i<=r;i++){
            if(num[i-l]==1) Ans = (Ans + ans[i-l]) % mod;  
            else Ans = (Ans + ans[i-l] * (k+1) % mod) % mod; //i还剩下一个大于10^6的素数
        }
        printf("%lld\n",Ans);
    }
}