N!分解素因子及若干问题【转载】

这里写的非常好http://www.cnblogs.com/openorz/archive/2011/11/14/2248992.html,感谢博主,我这里就直接用了。

将N!表示成

N! = p1^t1*p2^t2*…pi^ti…*pk^tk(其中p1,p2……pk是素数,1<N<= 10^6)

显然很容易通过素数筛选求出pi,因为1<pi<=N,关键是如何快速地求出ti。

我们先来看一下对于2这个素因子,把N!分成两部分,即奇偶两部分

假设N是偶数

N!

=1*2*3*4*5……N

=(2*4*6……) * (1*3*5……)

因为有N/2个偶数,所以偶数部分可以提出N/2个2,

=2^(N/2) * (1*2*3*……N/2) * (1*3*5*……)

=2^(N/2) * (N/2)! * (1*3*5*……)

看到了吗!神奇的事情发生了,N规模的问题转化成了N/2的问题了。上面假设了N是偶数,当然N是奇数时也是一样的,只要规定这里的除法是取整就可以了

于是有递推公式 f(n,2) = f(n/2,2) + n/2,表示n!中2的个数。

用同样的方法可以推出 f(n,p) = f(n/p) + n/p,表示n!中素数p的个数。

于是有代码

1 int f(int n,int p)
2 {
3     if(n==0) return 0;
4     return f(n/p) + n/p;
5 }

将问题推广一下:

问题1:N!的末尾有几个0?

因为 10 = 2*5,所以只要知道N!有多少个2和多少个5,问题就解决了。Min(f(n,2),f(n,5)) 显然f(n,2)>f(n,5),所以问题就转化成了求f(n,5)。

问题2:N!的转化成12进制之后,末尾有几个0?

和问题一样,12=2*2*3,所以只要求Min(f(n,2)/2,f(n,3)),就可以了。

问题3: 求组合数C(n,m)(mod p)

C(n,m)=n!/(m!*(n-m)!) ,只要对分子和分母分别分解素因子,然后因为C(n,m)肯定是整数,所以C(n,m)肯定可以表示成p1^t1*p2^t2*......pi^ti的形式,只要拿分子素因子的幂减去分母对应的素因子的幂即可。好了,后面就简单了,二分快速幂取模......

时间: 2024-11-08 09:42:58

N!分解素因子及若干问题【转载】的相关文章

【小方法--最简】C - 分解素因子

C - 分解素因子Crawling in process... Crawling failed Time Limit:2000MS     Memory Limit:10240KB     64bit IO Format:%lld & %llu Submit Status    Description 假设x是一个正整数,它的值不超过65535(即1< x < =65535),请编写一个程序, 将x分解为若干个素数的乘积.    Input 输入的第一行含一个正整数k (1<=k

Foj 1075 分解素因子

题目链接:http://acm.fzu.edu.cn/problem.php?pid=1075 思路:分解素因子 #include <iostream> #include <cstring> #include <cstdio> using namespace std; const int maxn=65535+5; typedef long long LL; bool prime[maxn]; int p[maxn/10]; int k; void isprime()

FZU 1075 分解素因子【数论/唯一分解定理/分解素因子裸模板】

[唯一分解定理]:https://www.cnblogs.com/mjtcn/p/6743624.html 假设x是一个正整数,它的值不超过65535(即1<x<=65535),请编写一个程序,将x分解为若干个素数的乘积. Input 输入的第一行含一个正整数k (1<=k<=10),表示测试例的个数,后面紧接着k行,每行对应一个测试例,包含一个正整数x. Output 每个测试例对应一行输出,输出x的素数乘积表示式,式中的素数从小到大排列,两个素数之间用“*”表示乘法. Samp

LightOJ 1340 - Story of Tomisu Ghost 阶乘分解素因子

http://www.lightoj.com/volume_showproblem.php?problem=1340 题意:问n!在b进制下至少有t个后缀零,求最大的b. 思路:很容易想到一个数通过分解素因子可以得到最大的指数.那么问题关键在于求得n!的素因子的指数,找到指数大于t的所有素因子,再将那些指数除去t,剩下的数就是最大的b了.分解阶乘时,对n不断除素数p,直到n为0时,此时商的和即该素因子的指数. /** @Date : 2016-11-30-19.35 * @Author : Lw

Xcode 6 beta3 中 Swift 的若干更新(转载)

Xcode 6 beta3 对Swift 自从公布以来的第一次比较大的更新,尤其是对数组的更新,主要体现在几个方面:1.数组被重写,现在声明为 let 的数组为不可变数组,声明为var 的是可变的2.声明数组和字典的语法糖也作了修改, 原来声明数组:Int[] 改为 [Int] ,同样的, 字典声明方式变为  [Key: Value]:3.区间运算符 .. 不再使用,而用  ..< 替代,使语意更明确. Xcode 6 beta3 中 Swift 的若干更新(转载)

hdu 1299 数论 分解素因子

题意:求1.1/x+1/y=1/z 给定z 求x,y的个数zsd:1: 要知道这个素数的因子的范围 范围为2--sqrt(n);2: 带入方程得:x = n * n / k + n ; 现在就变成了求x的值.又因为n一定大于k,所以转换成求n * n的分解数: 因为n = p1 ^ e1 * p2 ^ e2 *..........*pn ^ en sum ( n)= ( 1 + e1 ) * ( 1 +e2 ) * .........* ( 1 +en ); sum (n * n) = ( 1

hdu-1299 Diophantus of Alexandria(分解素因子)

思路: 因为x,y必须要大与n,那么将y设为(n+k);那么根据等式可求的x=(n2)/k+n;因为y为整数所以k要整除n*n; 那么符合上面等式的x,y的个数就变为求能被n*n整除的数k的个数,且k>=n; 那么k就在n到n*n之间,不过这样不则么好求. 因为x,y具有对称性,从y=(n+k)和x=(n2)/k+n;也可以看出(n*n/k,和k等价因为k*(n*n/k)=n*n): 由于任何数都可以拆成素数的乘积,形式为a1^x1*a2^x2*a3*^x3......; 那么因数的个数为(1+

fzu 1075 分解素因子

代码: #include<cstdio> #include<cstring> #include<iostream> using namespace std; int prime(int n) { if(n==2||n==3) { return 1; } if(n%6!=1&&n%6!=5) return 0; for(int i=5; i*i<=n; i+=6) { if(n%i==0||n%(i+2)==0) return 0; } return

Java hashCode() 和 equals()的若干问题解答&lt;转载自skywang12345&gt;

本章的内容主要解决下面几个问题: 1 equals() 的作用是什么? 2 equals() 与 == 的区别是什么? 3 hashCode() 的作用是什么? 4 hashCode() 和 equals() 之间有什么联系? 第1部分 equals() 的作用 equals() 的作用是 用来判断两个对象是否相等. equals() 定义在JDK的Object.java中.通过判断两个对象的地址是否相等(即,是否是同一个对象)来区分它们是否相等.源码如下: public boolean equ