【区间筛】2017多校训练四 HDU6069 Counting Divisors

http://acm.hdu.edu.cn/showproblem.php?pid=6069

【题意】

给定l,r,k，求

d(n)是n的因子个数

【思路】

【Accepted】

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<string>
 #include<cmath>
 #include<algorithm>
 #include<queue>
 #include<stack>
 #include<map>
 #include<vector>
 #include<set>
 using namespace std;
 typedef long long ll;
 const int maxn=1e6+;
 const ll mod=;
 ll l,r,k;
 bool isprime[maxn];
 int prime[maxn];
 int cnt;
 ll tot[maxn];
 vector<int> s[maxn];
 ll a[maxn];
 void init()
 {
     memset(isprime,true,sizeof(isprime));
     isprime[]=false;
     for(int i=;i<maxn;i++)
     {
         if(isprime[i])
         {
             for(int k=*i;k<maxn;k+=i)
             {
                 isprime[k]=false;
             }
         }
     }
     cnt=;
     for(int i=;i<maxn;i++)
     {
         if(isprime[i])
         {
             prime[cnt++]=i;
         }
     }
 }
 ll fac(ll& num,ll p)
 {
     ll cou=;
     while(num%p==)
     {
         cou++;
         num/=p;
     }
     return cou;
 }

 int main()
 {
     init();
     int T;
     scanf("%d",&T);
     while(T--)
     {
         scanf("%lld%lld%lld",&l,&r,&k);
         for(int i=;i<=r-l;i++)
         {
             a[i]=(ll)i+l;
             tot[i]=;
         }
         for(int i=;i<cnt;i++)
         {
             if(prime[i]>r) break;
             ll cou=l/(ll)prime[i];
             if(l%(ll)prime[i]) cou++;
             for(ll j=cou*(ll)prime[i];j<=r;j+=(ll)prime[i])
             {
                 int num=(int)(j-l);
                 if(a[num]==) continue;
                 ll sum=fac(a[num],prime[i]);
                 tot[num]=tot[num]*((k*sum)%mod+)%mod;
             }
         }
         ll ans;
         if(l==) ans=;
         else ans=;
         for(int i=;i<=(int)(r-l);i++)
         {
             if(a[i]!=) tot[i]=tot[i]*(k+)%mod;
             if(tot[i]!=)
             ans=(ans+tot[i])%mod;
         }
         cout<<ans<<endl;
     }
     return ;
 }