[Acwing] 197. 阶乘分解 筛法预处理

前言

传送门 :

思路

根据算数基本定理
N = p 1 a 1 ∗ p 2 a 2 . . . p n a n N = p1^{a1}*p2^{a2}...pn^{an} N=p1a1p2a2...pnan

N ! = N ∗ ( N − 1 ) ∗ ( N − 2 ) . . . 1 N! =N*(N-1)*(N-2)...1 N!=N(N1)(N2)...1

所以暴力解法我们可以知道

s u m ( N ! ) = p 1 a 1 + a 2 + . . . + a 3 + a k + p 2 a 1 + a 2 + . . + a k . . . sum(N!) = p1^{a_1+a_2+...+a_3+a_k}+p2^{a_1+a_2+..+a_k}... sum(N!)=p1a1+a2+...+a3+ak+p2a1+a2+..+ak...

所有的 p p p都在 N N N内,所以我们可以先处理出所有的 p p p

然后遍历 p p p,求出所有 p p p的指数即可

Mycode

b

ool st[N];
int primes[N],cnt;

void get_Prime(int x){
	for(int i=2;i<=x;i++){
		if(st[i]) continue;
		for(int j=2;j<=x/i;j++){
			st[i*j] = 1;
		}
	}
	
	for(int i=2;i<=x;i++){
		if(!st[i])
		primes[++cnt] = i ;
	}
}
void solve()
{
	int n;cin>>n;
	get_Prime(n);//获取 1到n的所有质数
	
	for(int i=1;i<=cnt;i++){
		int temp =  n;
		int u = primes[i];
		int ans =  0 ;
		
		while(temp){
			ans+=(temp/u);
			temp/=u;
		}
		
		cout<<u<<" "<<ans<<endl;
		
	}
	
	
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值