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

关于阶乘

时间:2016-04-02 23:00:37      阅读:419      评论:0      收藏:0      [点我收藏+]

标签:

描述:给定两个数n,m,其中m是一个素数。

将n(0<=n<=2^31)的阶乘分解质因数,求其中有多少个m。

注:^为求幂符号。

输入:

  第一行是一个整数s(0<s<=100),表示测试数据的组数
  随后的s行, 每行有两个整数n,m。 

输出:

  输出m的个数

样例输入

  3

  100 5

  16 2

  1000000000 13

样例输出
  

24   15   83333329

  当n,m体量很小的时候,用这个代码就可以AC:

 

技术分享
 1  #include <iostream>
 2 using namespace std;
 3 
 4 //测试分解test中有多少个m 
 5 int num_m(int test,int m){
 6     int res = 0;
 7     while(test%m==0){
 8         res++;
 9         test/=m;
10     }
11     return res;
12 }
13 
14 int main(){
15     int s; 
16     cin>>s;
17     while(s--){
18         int n,m ;
19         cin>>n>>m ;
20         int res = 0;
21         for(int i=n;i>=1;i-- ){
22             res+= num_m(i,m);                    
23         }
24         cout<<res<<endl;
25     }
26     
27 }         
View Code

 

但是实际情况是这个数还是比较大的,已经达到了2^31量级,所以上面暴力的解法是超时的。

接下来我们考虑:比如n=100,m=5.的情况,这个时候100当中有5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100

这些数能够整除5;其中25,50,75,100又能整除25=5^2,所以可以分解出两个5,因为第一遍的时候已经加过一次了,所以最多只能分解出20+4=24个5.

按照这个思想,我们不妨设置这样一个算法:让n循环除以m,每次得到的商是本次循环中能够整除m的个数,然后n=n/m赋予n新值。n缩小为原来的1/m,继续循环直到n=0为止。这样以来,"数值上接近于n的那些数能够分解质因数m的个数"恰好等于"接近于m及m的小倍数的这些数的计算的次数",如此对称下去,从而巧妙的完成了计算。

技术分享
 1 #include<stdio.h>
 2 #include<string.h>
 3 #include<stdlib.h>
 4 #include<math.h>
 5 
 6 int main(){
 7     int N;
 8     scanf("%d",&N);
 9     while(N--){
10         int n,m,sum=0;
11         scanf("%d%d",&n,&m);
12         while(n){
13             sum+=n/m;//计算n中有多少个能整除m的数;  
14             n/=m;     // 计算n中有多少个能够整除m^2的数;
15             
16             //n = n/m  这样就把n的大小缩小成了原来的(1/m)  
17         } 
18         printf("%d\n",sum);
19     }
20     
21     return 0;
22 }
View Code

 

关于阶乘

标签:

原文地址:http://www.cnblogs.com/liugl7/p/5348442.html

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