1294. 樱花（求约数）

最新推荐文章于 2024-10-02 23:40:03 发布

seez

最新推荐文章于 2024-10-02 23:40:03 发布

阅读量187

点赞数

分类专栏：数论结论文章标签：算法数据结构动态规划

本文链接：https://blog.csdn.net/qq12323qweeqwe/article/details/122888129

版权

数论同时被 2 个专栏收录

17 篇文章 0 订阅

订阅专栏

结论

5 篇文章 0 订阅

订阅专栏

本文探讨了如何利用n!的阶乘分解特性和约数计数原理，解决AcWing题库中涉及数对问题的优化方法。通过转化模型，将问题转化为求n!^2的约数个数，进而应用阶乘分解和约数模板求解。关键步骤包括筛选质因子和计算每个质因子的倍数分布。

摘要由CSDN通过智能技术生成

原题链接 - AcWing

分析：思考怎么进行优化，就是找出一个特殊的性质,由于看到了n!,可以想到197. 阶乘分解 - AcWing题库这道题，并且看到了数对数量，可能是和约数的个数是有关的

限制条件：

x,y必须是正整数
满足 1/x + 1/y =1/n!

那么可以通过这两个限制条件来转化一下模型，过程如下

由于x>n!,x-n!>0

那么问题就转化成为了

n!^2 的约数个数

n!^2的约数个数，应该先算出n!的每个质因子和个数，n!^2的每个质因子的个数=n!的每个质因子的个数*2，那么套用求约数个数的模板就可以了

由于是求n!,那么n!的各个质因子个数就套用阶乘分解就可以了，要求约数个数就套用约数个数的模板就可以了

#include <iostream>
#include <algorithm>
using namespace std;
const int mod=1e9+7,N=1e6+10;
int primes[N];
int cnt;
bool st[N];
void init(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!st[i]) primes[cnt++]=i;
        for(int j=0;primes[j]*i<=n;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0) break;
        }
    }
}
int main()
{
    int n;
    cin>>n;
    //由于求n!,筛一遍质数
    init(n);
    int res=1;
    for(int i=0;i<cnt;i++)
    {
        int p=primes[i];
        int s=0;
        for(int j=n;j;j/=p)
            s+=j/p;
        res=1ll*res*(2*s+1)%mod;
    }
    cout<<res;
    return 0;
    
}