hdu6069 Counting Divisors

根据约数个数定理：n=p1^a1×p2^a2×p3^a3*…*pk^ak,n的约数的个数就是(a1+1)(a2+1)(a3+1)…(ak+1).

若i=p1^a1×p2^a2×p3^a3*…*pk^ak,则i^K=p1^(a1*K)×p2^(a2*K)×p3^(a3*K)*…*pk^(ak*K),i^K的约数的个数就是(a1*K+1)(a2*K+1)(a3*K+1)…(ak*K+1)

题目重点转换为L~R的质因数分解，两次塞选，第一次塞出1e6以内的所有质数，第二次枚举区间[L,R]中之前塞选得到质数的倍数。

最好全部开ll，中间会爆。

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
typedef long long ll;
using namespace std;
const int maxn=1e6+10;
const int mod=998244353;
int tot,t;
ll l,r,k,ans,cnt[maxn],q[maxn],primes[maxn];
bool vis[maxn];
void init(){
    for(int i=2;i<maxn;i++){
        if(!vis[i])
            primes[tot++]=i;
        for(int j=0;j<tot;j++){
            int k=i*primes[j];
            if(k>maxn)break;
            vis[k]=1;
        }
    }
}
int main(){
    scanf("%d",&t);
    init();
    while(t--){
        scanf("%lld%lld%lld",&l,&r,&k);
        ans=0;
        if(l==1) ans++,l++;
        for(ll i=0;i<=r-l;i++) cnt[i]=1,q[i]=l+i;
        for(ll i=0;primes[i]*primes[i]<=r;i++){
            ll j=l/primes[i]+(l%primes[i]!=0);
            for(j=j*primes[i];j<=r;j+=primes[i]){
                ll tmp=0;
                while(q[j-l]%primes[i]==0) q[j-l]/=primes[i],tmp++;
                cnt[j-l]*=(tmp*k)%mod+1,cnt[j-l]%=mod;
            }
        }
        for(ll i=0;i<=r-l;i++){
            if(q[i]!=1) ans+=((k+1)*cnt[i])%mod;
            else        ans+=cnt[i];
            ans%=mod;
        }
        printf("%lld\n",ans);
    }
}