hdu 6069 区间筛

题意：求出i满足l到r区间所有 i^k次方的因子的数量的和。
思路：对10^6的区间进行素数筛， 求出其质因数，再根据每个质因数的(1+number*k)的乘积求出因子和。

#include<cstdio>#include<iostream>#include<algorithm>#include<cstring>using namespace std;const int MAXN = 1e6+1e3;  bool vis[MAXN];typedef long long LL;const int mod = 998244353;LL num[MAXN];LL prime[MAXN];LL a[MAXN];int cnt;void seg_vis(LL l,LL r, int k)   //区间筛法{    LL n=r-l;    for(LL i=0; i<=n; i++) num[i]=1, a[i]=l+i;    for (LL i=0; prime[i]<=n+1&&i<cnt; i++)    {        LL j = 0;        if (l % prime[i] != 0)//第一个需要筛掉的数(j+a) % isprime[i] == 0            j = prime[i] - l % prime[i];        for (; j<=n; j+=prime[i])        {            LL Count = 0;            while(a[j]%prime[i]==0)            {                Count++;                a[j]/=prime[i];            }            num[j]=(num[j]*(k*Count+1))%mod;        }    }    LL sum=0;    for(LL i=0; i<=n; ++i)    {        if(a[i]!=1)            num[i]=num[i]*(k+1)%mod;        sum=(sum+num[i])%mod;    }    printf("%lld\n", sum);    return ;}int main(){    cnt=0;    memset(vis, true, sizeof(vis));    for(int i=2; i<MAXN; i++) //预处理        if(vis[i])        {            prime[cnt++]=i;            for(int j=2*i; j<MAXN; j+=i)                vis[j]=false;        }    LL l, r;    int t;    int k;    scanf("%d", &t);    while(t--)    {        scanf("%lld %lld %d", &l, &r, &k);        seg_vis(l, r, k);    }    return 0;}