hdu6069

根据约数个数定理：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，中间会爆。

[cpp] view plain copy

print ?

1. #include<iostream>  
2. #include<stdio.h>  
3. #include<algorithm>  
4. #include<string.h>  
5. typedef long long ll;  
6. using namespace std;  
7. const int maxn=1e6+10;  
8. const int mod=998244353;  
9. int tot,t;  
10. ll l,r,k,ans,cnt[maxn],q[maxn],primes[maxn];  
11. bool vis[maxn];  
12. void init(){  
13.     for(int i=2;i<maxn;i++){  
14.         if(!vis[i])  
15.             primes[tot++]=i;  
16.         for(int j=0;j<tot;j++){  
17.             int k=i*primes[j];  
18.             if(k>maxn)break;  
19.             vis[k]=1;  
20.         }  
21.     }  
22. }  
23. int main(){  
24.     scanf("%d",&t);  
25.     init();  
26.     while(t--){  
27.         scanf("%lld%lld%lld",&l,&r,&k);  
28.         ans=0;  
29.         if(l==1) ans++,l++;  
30.         for(ll i=0;i<=r-l;i++) cnt[i]=1,q[i]=l+i;  
31.         for(ll i=0;primes[i]*primes[i]<=r;i++){  
32.             ll j=l/primes[i]+(l%primes[i]!=0);  
33.             for(j=j*primes[i];j<=r;j+=primes[i]){  
34.                 ll tmp=0;  
35.                 while(q[j-l]%primes[i]==0) q[j-l]/=primes[i],tmp++;  
36.                 cnt[j-l]*=(tmp*k)%mod+1,cnt[j-l]%=mod;  
37.             }  
38.         }  
39.         for(ll i=0;i<=r-l;i++){  
40.             if(q[i]!=1) ans+=((k+1)*cnt[i])%mod;  
41.             else        ans+=cnt[i];  
42.             ans%=mod;  
43.         }  
44.         printf("%lld\n",ans);  
45.     }  
46. }  

#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);
    }
}