hdu 6069 Counting Divisors(筛法)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6069

题意：
i∈[l,r]，求d(i^k)在[l,r]的和，d(a)表示a的因子和(d(12)=6,d(9)=3)
(1≤l≤r≤10^12,r−l≤10^6,1≤k≤10^7)

解析:
求一个数的因子和s=p1^a1*p2^a2*p3^a3…..
d(s)=(a1+1)(a2+1)(a3+1)*……
d(s^k)=(a1*k+1)(a2*k+1)(a3*k+1)*……
打表先枚举10^6以内的素数(因为10^12以内的合数，因子不可能超过√10^6)
再之后枚举所有√r以内的素数，求这些素数p的倍数x，再对所计算出的倍数x*p分解该素数p，分解得到的指数,与记录该数x*p的之前因子乘积相乘。这样就得到每个数的因子和了
最后就是对剩下在[l,r]里面的质数和分解后得到的质数>√r的数求因子和
这样所有数的因子和就求出来了，最后再想加就好了。这有有点类似与筛法的性质，用一个数取不断的筛取待求数。

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<math.h>
#define MAXN 1000005
typedef long long ll;

const int MOD = 998244353;

ll prime[MAXN];
ll check[MAXN];
ll a[MAXN],b[MAXN];   //a表示在第i个位置上的现在的数字是几,b表示在第i个位置上的数的因子个数
ll l,r,k;
ll tot;

void calprime()  //求1-1e6以内的所有素数
{
    tot=0;
    memset(check, 0, sizeof(check));
    for (ll i = 2; i < MAXN; ++i)
    {
        if (!check[i])
        {
            prime[tot++] = i;
        }
        for (ll j = 0; j < tot&& i * prime[j]<MAXN; ++j)
        {
            if (i * prime[j] > MAXN)
            {
                break;
            }    
             check[i*prime[j]] = 1;
            if (i % prime[j] == 0)
            {
                break;
            }
        }
    }
}

int main()
{
    int t;
    ll ans=0;
    calprime();
    scanf("%d",&t);
    while(t--)
    {
        ans=0;
        scanf("%lld%lld%lld",&l,&r,&k);
        for(ll i=l;i<=r;i++) a[i-l]=i;  //预处理
        for(ll i=l;i<=r;i++) b[i-l]=1;
        for(ll i=0;i<tot&&prime[i]*prime[i]<=r;i++)
        {

            ll x=l/prime[i];  //素数prime[i]在[l,r]最小的倍数
            if(x*prime[i]<l)
                x++;

            for(;x*prime[i]<=r;x++)   //x表示倍数,对prime[i]的倍数x*prime[x]求解里面有多少个prime[i]
            {
                ll s=0;
                while(a[x*prime[i]-l]%prime[i]==0)   //在x*prime[i]中分解prime[i]的因子
                {
                    s++;
                    a[prime[i]*x-l]/=prime[i];   //将第i个位置的数减小，并把该减小的次数加到s中，为了将循环变量递减，直至将x*prime[i]中质因子prime[i]的指数完全读取出来
                }
                b[x*prime[i]-l]=(b[x*prime[i]-l]*(s*k+1))%MOD;   //将prime[i]的指数s与其他因子指数相乘记录与b中


            }
        }
        for(ll i=l;i<=r;i++)  //将所有底数i是质数或被分解质因子处理后的数(质数)>sqrt(r)的情况处理掉
        {
            if(a[i-l]>1)
                b[i-l]=(b[i-l]*(k+1))%MOD;

        }
        for(ll i=l;i<=r;i++)
        {
            ans=(ans+b[i-l])%MOD;
        }
        printf("%lld\n",ans);
    }
    return 0;
}