码迷,mamicode.com
首页 > 其他好文 > 详细

CH3101 阶乘分解

时间:2019-02-26 18:59:54      阅读:243      评论:0      收藏:0      [点我收藏+]

标签:namespace   bit   c++   text   char   阶乘   pac   顺序   mes   

题意

描述

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

输入格式

一个整数N。

输出格式

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

样例输入

5

样例输出

2 3
3 1
5 1

样例解释

5! = 120 = 2^3 * 3 * 5

分析

对1~N中每个质数统计即可。时间复杂度\(O(\frac{N}{\ln n}\log n)\)

代码

#include<bits/stdc++.h>
#define rg register
#define il inline
#define co const
template<class T>il T read(){
    rg T data=0,w=1;
    rg char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-') w=-1;
        ch=getchar();
    }
    while(isdigit(ch))
        data=data*10+ch-'0',ch=getchar();
    return data*w;
}
template<class T>il T read(rg T&x){
    return x=read<T>();
}
typedef long long ll;
using namespace std;

vector<int> prime;
bool v[1000010];
int n;
void primes(int n){
    for(int i=2;i<=n;++i){
        if(v[i]) continue;
        prime.push_back(i);
        for(int j=i;j<=n/i;++j) v[i*j]=1;
    }
}
int main(){
//  freopen(".in","r",stdin);
//  freopen(".out","w",stdout);
    primes(read(n));
    for(unsigned i=0;i<prime.size();++i){
        int p=prime[i],c=0;
        for(int j=n;j;j/=p) c+=j/p;
        printf("%d %d\n",p,c);
    }
    return 0;
}

CH3101 阶乘分解

标签:namespace   bit   c++   text   char   阶乘   pac   顺序   mes   

原文地址:https://www.cnblogs.com/autoint/p/10439156.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!