首页 > 其他分享> 文章详细

阶乘分解

2021-10-20 16:00:34 阅读：205 来源： 互联网

标签：prime frac int 样例分解阶乘质因数

题目链接

acwing197. 阶乘分解

题目描述

给定整数 \(N\)，试把阶乘 \(N!\) 分解质因数，按照算术基本定理的形式输出分解结果中的 \(p_i\) 和 \(c_i\) 即可。

输入格式

一个整数 \(N\)。

输出格式

\(N!\) 分解质因数后的结果，共若干行，每行一对 \(p_i,c_i\)，表示含有 \({p_i}^{c_i}\) 项。按照 \(p_i\) 从小到大的顺序输出。

数据范围

\(1≤N≤10^6\)

输入样例：

输出样例：

2 3
3 1
5 1

样例解释

\(5!=120=23∗3∗5\)

解题思路

\(n!\) 中质因数 \(p\) 的个数至少有 \(\frac{n}{p}\) 个，\(p^2\) 的个数至少有 \(\frac{n}{p^2}\)个，\(\dots\)
其中，\(p^2\) 中质因子的个数在 \(p\) 中也算过一次，累加即可~
故 \(n!\) 中质因子的个数：

\[\frac{n}{p}+\frac{n}{p^2}+\dots+\frac{n}{p^{log_pn}} \]

时间复杂度：\(O(nlogn)\)

代码

#include<cstdio>
using namespace std;
int prime[1000000];
int v[1000005];
int m,n;
void primes(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!v[i])
        {
            v[i]=i;
            prime[++m]=i;
        }
        for(int j=1;j<=m;j++)
        {
            if(v[i]<prime[j]||i*prime[j]>n)break;
            v[i*prime[j]]=prime[j];
        }
    }
}
int main()
{
    scanf("%d",&n);
    primes(n);
    for(int i=1;i<=m;i++)
    {
        int p=prime[i],c=0,t=n;
        while(t>=p)c+=t/p,t/=p;
        printf("%d %d\n",p,c);
    }
    return 0;
}

标签：prime,frac,int,样例,分解,阶乘,质因数
来源： https://www.cnblogs.com/zyyun/p/15429189.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9