《算法笔记》5.5 质因子分解

5.5 质因子分解

所谓质因子分解是指将一个正整数n写成一个或多个质数的乘积的形式，例如6=2x3，8=222，180=22335，或者我们也可以写成指数的形式，例如6=2^!*31，8=2^3，180 =22335’。显然，由于最后都要归结到若干不同质数的乘积，因此不妨先把素数表打印出来。而打印素数表的方法上面已经阐述，下面我们主要就质因子分解本身进行讲解。注意：由于1本身不是素数，因此它没有质因子，下面的讲解是针对大于1的正整数来说的，而如果有些题目中要求对1进行处理，那么视题目条件而定来进行特判处理。
由于每个质因子都可以不止出现一次，因此不妨定义结构体factor，用来存放质因子及其个数，如下所示：

struct factor{
    int x ; //x为质因子
    int cnt ; //cnt为其个数
}fac[10] ;

对一个int型范围的数来说，fac数组的大小只需要开到10就可以了。

前面提到过，对一个正整数n来说，如果它存在1和本身之外的因子，那么一定是在sqrt（n）
的左右成对出现。而这里把这个结论用在“质因子”上面，会得到一个强化结论：对一个正整数n来说，如果它存在[2，n]范围内的质因子，要么这些质因子全部小于等于sqt（n），要么只存在一个大于sqrtmn）的质因子，而其余质因子全部小于等于sqrt（n），这就给进行质因子分解提供了一个很好的思路：

1. 枚举1~sart（n）范围内的所有质因子p，判断p是否是n的因子。

• 如果p是n的因子，那么给fac数组中增加质因子p，并初始化其不数为0。然后，只要p还是n的因子，就让n不断除以p，每次操作令p的个数加1，直到p不再是n的因子为止。

if( n % prime[i] == 0 ){
    fac[num].x = prime[i] ;
    fac[num].cnt = 0 ;
    while( n % prime[i] == 0 ){     //计算出质因子prime[i]的个数
        fac[num].cnt ++ ;
        n /= prime[i] ;
    }
    num ++  //不同质因子个数加1
    
}

• 如果p不是n的因子，就直接跳过。

2. 如果在上面步骤结束后n仍然大于1，说明n有且仅有一个大于sqrt（n）的质因子（有可能是n本身），这时需要把这个质因子加入fac数组，并令其个数为1.

if( n != 1 ){  //如果无法被根号n以内的质因数除尽
    fac[num].x = n ;    //那么一定有一个大于根号n的质因子
    fac[num++].cnt = 1 ;    
    
}

时间复杂度为O(sqrt(n))。