HDU-6069 Counting Divisors

最新推荐文章于 2021-08-09 20:41:35 发布

wjhshuai

最新推荐文章于 2021-08-09 20:41:35 发布

阅读量206

点赞数

分类专栏：数论文章标签：数论求因子数

本文链接：https://blog.csdn.net/wjhshuai/article/details/77001205

版权

数论专栏收录该内容

9 篇文章 0 订阅

订阅专栏

Counting Divisors

这道题其实就是数论的题目（题目链接）

超时思路是就是从l到r所有的因子都枚举出来·····然后相加起来

当然是个TL

那么就要想办法减小时间复杂度

最好下手的地方就是在枚举因子的时候。

如何做到快速枚举因子？

其实是有个一个相关，可以用的上的结论

任意的一个数都可以表示为 n=p1^a1*p2^a2*p3^a3·····(p和a都是未知的并且p均为素数)

如 200= 2^3*5^2

那么这个数的因子数量为d(n)=（a1+1）*（a2+1）*（a3+1）·······

至于为什么？

（对于两个互素的数，这两个数的积的因子数等于这两个数的因子数相乘，可以写两个互素的数试一下，很容易想明白）p1^a1和p2^a2肯定是互素的·····

再拓展一下，n^k=p1^ka1*p2^ka2*p3^ka3······

那么d(n^k)=(k*a1+1)*(k*a2+1)*(k*a3+1)

由于我还未研究过数论，只能给出这么粗糙的 ‘猜想’

有了上面的结论，我们就可以对求因子数进行优化了

那么我就需要把一个数表示成上式的形式，求出因子数储存起来

那么首先我们要暴力打表，那么打多少的素数表？

打sqrt(n )的表就可以了，如果大于sqrt(n)，可能的素数只可能有0个或者1个

那么问题来了，为什么那？

由于计算因子的时候一般都是成对计算

一个数最多可以由一对有素数的数相乘得到（或者没有）

比如 68 因子为：1 2 4 17 34 68 只有 4 和17一对

而打表到sqrt(n)的时候，前一半的因子都打到了

而后一半的因子中，最多只有一个

然后代码：

#include<bits/stdc++.h>
using namespace std;

typedef long long LL;
const LL maxn=1e6+10;
const LL mod=998244353;
LL prime[maxn/10],num[maxn],sum[maxn];//sum[i]记录i+le的约数个数
bool vis[maxn];
LL k,tot,ri,le;

void Table()//1e6以内的素数打表
{
    tot=0;
    for(LL i=2; i<maxn; ++i)
    {
        if(!vis[i])
            prime[tot++]=i;
        for(LL j=0; j<tot&&i*prime[j]<maxn; ++j)
        {
            vis[i*prime[j]]=true;
            if(i%prime[j]==0)
                break;
        }
    }
}
LL solve()
{
    for(LL i=0; i<tot&&prime[i]*prime[i]<=ri; ++i)
    {
        LL tmp=le;
        if(tmp%prime[i])//将tmp变成>=le的第一个能整除prime[i]的数
            tmp=(tmp/prime[i]+1)*prime[i];
        for(LL j=tmp; j<=ri; j+=prime[i])//j都能整除prime[i]
        {
            LL cnt=0;
            while(num[j-le]%prime[i]==0)
            {
                ++cnt;
                num[j-le]/=prime[i];
            }
            sum[j-le]=(sum[j-le]*(cnt*k+1)%mod)%mod;
        }
    }
    LL ans=0;
    for(LL i=0; i<=ri-le; ++i)
    {
        if(num[i]>1)//对于有大于sqrt(r)的素因子的情况进行判断
            ans=(ans+sum[i]*(k+1)%mod)%mod;
        else
            ans=(ans+sum[i])%mod;
    }
    return ans;
}
int main()
{
    Table();
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld%lld%lld",&le,&ri,&k);
        for(LL i=0; i<=ri-le; ++i)
            num[i]=i+le,sum[i]=1;
        printf("%lld\n",solve());
    }
    return 0;
}

Alzh大佬的代码，方法很巧妙，尤其是找q和a的时候

看不懂的可以多看几遍