HDU 6069 - Counting Divisors - (素数筛答案)

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

参考

题意:首先d(n)表示n的约数的个数,每组样例给出l,r,k让求:

                                                                  \biggl(\sum_{i=l}^rd(i^k)\biggr)\mod998244353

其中:l,r,k(1\leq l\leq r\leq 10^{12},r-l\leq 10^6,1\leq k\leq 10^7)

解析:首先要知道约数个数定理:

对于一个大于1正整数n可以分解质因数:

https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D231/sign=e19b382db04543a9f11bfdcf2f168a7b/a8ec8a13632762d0a5814266a7ec08fa503dc6a8.jpg

则n的正约数的个数就是:

https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D316/sign=ca4b36d4d243ad4ba22e40c1b4035a89/8ad4b31c8701a18b94aed93e9d2f07082838fe7e.jpg

那么d(i^k)=(ka_1+1)(ka_2+1)...(ka_m+1),考虑到r-l+1 \leq 1e^6 ,尝试枚举l-r间的每个数,然后对每个数分解质因子是超时的。转变思路想到枚举1-1e^6之间的所有素数prime[i],再枚举区间[l,r]中所有prime[i]的倍数x,对所有prime[i]的倍数x累计prime[i]对其的贡献,然后将x除以prime[i]表示己经算过prime[i]的贡献,最后区间中每一个数除完剩下值不为1的部分就是超过枚举的质数的剩余质数,只可能是0个或1个,对于有剩余的再统计一下这个最后的质数的贡献即可。具体看代码吧,代码还是很好理解的。

复杂度O(\sqrt{r}+(r-l+1)\log\log(r-l+1))

代码(2355ms)

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll MAXN=1e6+500;
const ll MOD=998244353;

ll len,prime[MAXN+50];
bool vis[MAXN+50];
ll l,r,k;
ll d[MAXN+50],arr[MAXN+50];

void init()
{
    ll tmp;
    int i,j;len=0;
    for(i=2;i<MAXN;i++)
    {
        if(!vis[i])
            prime[len++]=i;
        for(j=0;j<len;j++)
        {
            tmp=i*prime[j];
            if(tmp>MAXN) break;
            vis[tmp]=true;
            if(i%prime[j]==0)
                break;
        }
    }
}

ll work(ll l,ll r,ll k)
{
    for(ll i=l;i<=r;i++)//d[i-l]统计区间中i的答案,arr[i-l]是它的本身数值
    {
        d[i-l]=1;
        arr[i-l]=i;
    }
    for(ll i=0;i<len;i++)//枚举每一个素数
    {
        ll pos=(l/prime[i])*prime[i]+(l%prime[i]!=0)*prime[i];//pos是素数prime[i]区间中的第一个倍数
        while(pos<=r)//pos枚举每个倍数
        {
            ll cnt=0;
            while(arr[pos-l]%prime[i]==0)//prime[i]对pos的答案的贡献为cnt*k+1
            {
                cnt++;
                arr[pos-l]=arr[pos-l]/prime[i];
            }
            d[pos-l]=d[pos-l]*(cnt*k+1)%MOD;
            pos+=prime[i];
        }
    }
    ll ans=0;
    for(ll i=l;i<=r;i++)
    {
        if(arr[i-l]==1)//没有剩余素数
            ans=(ans+d[i-l])%MOD;
        else {
            ans=(ans+d[i-l]*(k+1))%MOD;
        }
    }
    return ans;
}
int main()
{
    init();
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lld%lld%lld",&l,&r,&k);
        printf("%lld\n",work(l,r,k));
    }
    return 0;
}

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值