HDU 6069 Counting Divisors 多校联合第四场

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

Titanzyh

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

阅读量323

点赞数

分类专栏：暑期多校联合

本文链接：https://blog.csdn.net/Titanzyh/article/details/76666072

版权

暑期多校联合专栏收录该内容

23 篇文章 0 订阅

订阅专栏

题目链接

用约数个数定理，题目转化为求[l,r]之间的数分解素因数求素因子指数，先枚举小于根下r的质数，因为所求大于根下r的质数的素因子的指数一定为1.

标程

#include<cstdio>
using namespace std;
typedef long long ll;
const int N=1000010,P=998244353;
int Case,i,j,k,p[N/10],tot,g[N],ans;
ll n,l,r,f[N];
bool v[N];
inline void work(ll p)
{
    for(ll i=l/p*p; i<=r; i+=p)
    {
        if(i>=l)
        {
            int o=0;
            while(f[i-l]%p==0)
                f[i-l]/=p,o++;//o 素因子的指数
            g[i-l]=1LL*g[i-l]*(o*k+1)%P;//除数个数
        }
    }
}
int main()
{
    ///记下小于N的素数
    for(i=2; i<N; i++)
    {
        if(!v[i])p[tot++]=i;
        for(j=0; j<tot&&i*p[j]<N; j++)
        {
            v[i*p[j]]=1;//i*p[j]两个素数相乘
            if(i%p[j]==0)
                break;
        }
    }

    scanf("%d",&Case);
    while(Case--)
    {
        scanf("%lld%lld%d",&l,&r,&k);
        n=r-l;//l与r之间的距离

        for(i=0; i<=n; i++)
            f[i]=i+l,g[i]=1;//f[i]指l至r之间的数
        ///tot 素数个数
        for(i=0; i<tot; i++)
        {
            if(1LL*p[i]*p[i]>r)
                break;
            work(p[i]);
        }
        ///处理大于根下r的质数
        for(ans=i=0; i<=n; i++)
        {
            if(f[i]>1)
                g[i]=1LL*g[i]*(k+1)%P;
            ans=(ans+g[i])%P;
        }
        printf("%d\n",ans);
    }
    return 0;
}