标签：log int 197 init 分解 阶乘 pi 质数

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

输入格式
一个整数 N。

输出格式
N! 分解质因数后的结果，共若干行，每行一对 pi,ci，表示含有 pi^ci 项。按照 pi 从小到大的顺序输出。

数据范围
3≤N≤10^6

输入样例：
5

输出样例：
2 3
3 1
5 1

样例解释
5!=120=2^3∗3∗5

题意：

任给一个整数n，让我们将n!分解质因数。将每个质因子及其对应的次数求出。（n<=1e6）

本题思路：（两步骤）

时间复杂度分析：O(n)

1~n中共有n/logn个素数，每一个质数p求取一下它对应的次数,我们按照下面的方法求：

每一个质数算次数都需要 log{p,n} 次，算出来总运算次数是：

log{2,n} + log{3,n} + log{5,n} + ...（每一项都是以质数为底的对数，根据公式一共有n / log{2,n}项）

我们放缩一下，上式 <= n / log{2,n} * log{2,n} = n（n / log{2,n}表示1~n中质数的个数），

我们最后再加上线筛预处理的O(n)，因此这个算法时间复杂度最终是 O(n) 的，总次数大概在2e6，属于合法范围内。

两份代码：一份while，一份for

#include<bits/stdc++.h>

using namespace std;

const int N = 1e6+10;
int primes[N],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;
    init(n);
    for(int i=0;i<cnt;++i)
    {
        int t = n;
        int p = primes[i];
        int s = 0;
        while(t>=1) s += t/p, t /= p;
        cout<<primes[i]<<' '<<s<<endl;
    }
    
    return 0;
}

#include<bits/stdc++.h>

using namespace std;

const int N = 1e6+10;
int primes[N],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;
    init(n);
    for(int i=0;i<cnt;++i)
    {
        int p = primes[i];
        int s = 0;
        for(int j=n; j>=1; j/=p) s+=j/p;
        cout<<primes[i]<<' '<<s<<endl;
    }
    
    return 0;
}

标签：log,int,197,init,分解,阶乘,pi,质数
来源： https://blog.csdn.net/Jacob0824/article/details/122847894

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