hdu6069 Counting Divisors

题意：l，r之间所有i^k（l≤i≤r）的因子数之和。

想法：

对于每一个素数，去筛这段区间[l,r]。对于第i个素数pri[i]，它只会每隔pri[i]个数出现一次，所以找到起点后，遍历所有有素因子pri[i]的数即可。

#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<deque>
#include<vector>
#include<functional>
using namespace std;
#define LL long long
#define mm(a,b) memset(a,b,sizeof(a))
const double eps=1.0e-6;
const double PI=acos(-1.0);
template<typename T>T gcd(T a,T b){return b==0?a:gcd(b,a%b);}
template<typename T>T lcm(T a,T b){return a/gcd(a,b)*b;}
template<typename T>T _abs(T a){return a>0?a:-a;}
const int maxn=1000000;
const LL mod=998244353;
int vis[maxn+1];
int pri[80000];
int spri;
LL a[maxn+1],s[maxn+1];
void sieve(int n)
{
	spri=0;
	mm(vis,0);
	for(int i=2;i<=sqrt(n+0.5);i++)
		if(!vis[i])
			for(int j=i*i;j<=n;j+=i)
				vis[j]=1;
	for(int i=2;i<=n;i++)
		if(!vis[i])
			pri[spri++]=i;
}
int main()
{
    sieve(maxn);
    int t;
    LL l,r,k;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld%lld%lld",&l,&r,&k);
        LL len=r-l+1;
        for(int i=0;i<len;i++)
        {
            a[i]=i+l;
            s[i]=1;
        }
        for(int i=0;i<spri;i++)
        {
            int st=(l+pri[i]-1)/pri[i]*pri[i]-l;//寻找起点
            for(int j=st;j<len;j+=pri[i])//从起点开始每个pri[i]个数会有素因子pri[i]出现
            {
                int cnt=0;
                while(a[j]%pri[i]==0)
                {
                    cnt++;
                    a[j]/=pri[i];
                }
                s[j]=s[j]*(1+k*cnt%mod)%mod;
            }
        }
        LL ans=0;
        for(int i=0;i<len;i++)
        {
            if(a[i]>1)
                s[i]=s[i]*(1+k)%mod;
            ans=(ans+s[i])%mod;
        }
        printf("%lld\n",ans);
    }
    return 0;
}