阶乘分解（数论）

最新推荐文章于 2022-10-17 08:08:18 发布

林苏泽

最新推荐文章于 2022-10-17 08:08:18 发布

阅读量713

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/weixin_45911397/article/details/110304106

版权

数论专栏收录该内容

72 篇文章 5 订阅

订阅专栏

题目链接
题意：
给定整数 $N(1≤N≤10^6)$ ，试把阶乘 $N!$ 分解质因数，按照算术基本定理的形式输出分解结果中的 $p_i$ 和 $c_i$ 即可。
输入一个整数 $n$
输出：
N! 分解质因数后的结果，共若干行，每行一对 $p_i$ , $c_i$ ，表示含有 $p_i^{c_i}$ 项。按照 $p_i$ 从小到大的顺序输出。

思路: 这个题用唯一分解会超时，因为外层 $n$ 内层 $l o g n$ 枚举的质数，所以复杂度就是 $n l o g n$ 的复杂度。所以不行
我们用到整除分块的思想，就是每个数记录它最小质因数，那么我们就可以一层一层的往下找他的质因数了。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=2e6+7;
bool isprime[maxn];
ll p[maxn],v[maxn],ans[maxn];
ll cnt;

int main(){
    int n;
    cin>>n;
    v[1]=1;
    for(int i=2;i<=n;i++){
        if(!v[i]) p[++cnt]=v[i]=i;///素数最小质因子就是它本身
        for(int j=1;j<=cnt&&i*p[j]<=n&&p[j]<=v[i];j++){
            v[i*p[j]]=p[j];///存放最小质因子　，每次除去最小质因子
        }
    }
    memset(ans,0,sizeof ans);
    for(int i=2;i<=n;i++){
        for(int j=i;j>1;j/=v[j]){
            ans[v[j]]++;///是ｖj的倍数　每次除去最小质因子
        }
    }
    for(int i=1;i<=cnt;i++) printf("%lld %lld\n",p[i],ans[p[i]]);
    return 0;
}

唯一分解的代码也放上吧：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn=1e6+7;
bool isprime[maxn];
ll prime[maxn],p[maxn];
ll cnt;
ll dnum[maxn],znum[maxn];
void Prime(int x) {
    for(int i=2; i<=x; i++) {
        if(!isprime[i]) prime[++cnt]=i,dnum[cnt]=i;
        for(int j=1; j<=cnt&&i*prime[j]<=x; j++) {
            isprime[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
void fe(int x){
    for(int i=1;i<=cnt;i++){
        while(x%prime[i]==0){
            znum[i]++;
            x=x/prime[i];
        }
        if(x<prime[i]) break;
    }
}
int main() {
    ll n;
    cin>>n;
    Prime(n);
    for(int i=1;i<=n;i++){
        fe(i);
    }
    for(int i=1;i<=cnt;i++){
        if(znum[i]){
            printf("%lld %lld\n",dnum[i],znum[i]);
        }
    }
    return 0;
}