因子和与因子个数 (乘性函数)

乘性函数的详解:http://blog.csdn.net/luyuncheng/article/details/8017016#t3

在非数论的领域,积性函数指所有对于任何a,b都有性质f(ab)=f(a)f(b)的函数。  

在数论中的积性函数:对于正整数n的一个算术函数 f(n),若f(1)=1,且当a,b互质时f(ab)=f(a)f(b),在数论上就称它为积性函数。

若对于某积性函数 f(n),就算a, b不互质,也有f(ab)=f(a)f(b),则称它为完全积性的。[1]

因子和函数与因子个数函数都是乘性函数。

因子和:

s(6)=s(2)*s(3)=3*4=12;

s(20)=s(4)*s(5)=7*6=42;

再看 s(50)= 1+2+5+10+25+50=93=3*31=s(2)*s(25),s(25)=1+5+25=31.

这在数论中叫积性函数,gcd(a,b)=1时 s(a*b)=s(a)*s(b);

如果p是素数 //特别重要的性质

s(p^n)=1+p+p^2+...+p^n= (p^(n+1)-1) /(p-1) (1)

T(p^n)=n+1;//因子的个数

定理:设正整数n的所有素因子分解n=p1^a1*p2^a2*p3^a3****Ps^as,那么

s(n)=[(p1^a1+1 -1)/(p1-1)]*[(p2^a2+1 -1)/(p2-1)]*[(p3^a3+1 -1)/(p3-1)]***[(ps^as+1 -1)/(ps-1)];

T(n)=(a1+1)*(a2+1)*(a3+1)***(an+1);

 

个数:

1. N的因子个数

条件:给定任意一个一个正整数N

要求:求其因子的个数

首先给出结论:对于任意的整型N,分解质因数得到N= P1^x1 * P2^x2* …… * Pn^xn;

则N的因子个数M为 M=(x1+1) * (x2+1) * …… *(xn+1);

证明过程:

首先 举个例子吧

24 = 2^3 * 3^1;

其质因子有:为2和3  指数为 3和1

那么对于2 有0 1 2 3四种指数选择,对于3 有0 1两种指数选择

所以 就是4 * 2 = 8 个因子个数

如果还是不懂,那么我们就列举出来吧

2 3

2^0*3^0=1             2^0*3^1=3

2^1*3^0=2             2^1*3^1=6

2^2*3^0=4             2^2*3^1=12

2^3*3^0=8             2^3*3^1=24

结果很清晰了吧??其实这里用到了数学的排列组合的知识

也就是说每一个质因子的不同指数幂与其它质因子相乘,得到的结果一定不会重复

因此能够将所有的因子都列举出来。

所以N的因子数M,我们可以用M=(x1+1) * (x2+1) * …… *(xn+1)表示

 2.求n!中含有某个因子个数的方法

2. N!的因子个数

有上面的结论,这个问题就变得明朗多了吧?嘿嘿,不要着急,这里面还有许多细节问题需要我们考虑。

a.  最大的质因子一定不会大于N

b.  N的质因子并不完全包含N!所有的质因子

至于原因是什么,自己想想吧,嘿嘿

那我们就直接说思路了:

首先,我们可以把所有的N以内的质数给打表求出来

然后,求每一个质因子的指数个数,这里用到了一个公式,:

     ei=[N/pi^1]+ [N/pi^2]+ …… + [N/pi^n]  其中[]为取整

       附:这一步最近又想到了一个更好的方法  int ei=0;while(N)  ei+=(N/=pi);   怎么样?? 

           (想一想为什么,实在想不通你就举个例子试一下)

最后,就是套公式计算了,M=(e1+1)*(e2+1)*……*(en+1)

算了,还是举个例子吧

比如5!

质因子2的指数是 2+1=3;

质因子3的指数是 1;

质因子5的指数是 1;

所以因子个数为 4 * 2 * 2 = 16

5!=120 因子有1 2 3 4 5 6 8 10 12 15 20 24 30 4060 120  刚好16

转载于:http://www.cnblogs.com/dolphin0520/archive/2011/04/11/2012891.html

求n的阶乘某个因子a的个数,如果n比较小,可以直接算出来,但是如果n很大,此时n!超出了数据的表示范围,这种直接求的方法肯定行不通。其实n!可以表示成统一的方式。

n!=(k^m)*(m!)*a   其中k是该因子,m=n/k,a是不含因子k的数的乘积

下面推导这个公式

n!=n*(n-1)*(n-2)*......3*2*1

   =(k*2k*3k.....*mk)*a      a是不含因子k的数的乘积,显然m=n/k;

   =(k^m)*(1*2*3...*m)*a

   =k^m*m!*a

接下来按照相同的方法可以求出m!中含有因子k的个数。

因此就可以求出n!中因子k的个数

int count(int n,int k)
{
    int num=0;
    while(n)
    {
        num+=n/k;
        n/=k;
    }
    return num;
}

 

3.给定数列的乘积因子个数

其实这个也是基于第一个结论得到的。

给定 a1 a2 a3 …… an;

我们可以找到最大的一个元素Max(a);

把Max以内的素数打表

然后把质因子清零,进行如下循环,就可以找到各个质因子的个数:

for(a=1;a<=n;a++)

for(p=1;p<_;p++ )

     if(__) e(p)++;

这样质因子的质数个数就求出来了,然后就可以根据公式M=(e1+1)*(e2+1)*……*(en+1)求出因子个数

posted @ 2015-01-22 12:32  人艰不拆_zmc  阅读(2178)  评论(0编辑  收藏  举报