试除法及其应用

目录

原理

简单应用

1.判断 n 是否为质数

2.求约数

3.分解质因数

---

原理

试除法最简单的思路便是：对于整数 n，从 1 ~ n 分别判断是否为它的因数，判断 n 次。

对它进行稍加优化：因数都是成对出现的。

如果 2 是 n 的因数，那么 n / 2 肯定也是 n 的因数，且在 n / 2 ~ n 的范围内的因数可以通过 1 判断出来，因此判断范围可以缩小至 1 ~ n / 2；

同样，如果 3 是 n 的因数，那么 n / 3 肯定也是 n 的因数，且在 n / 3 ~ n 的范围内的因数可以通过 1 和 2 判断出来，因此判断范围可以缩小至 1 ~ n / 3；

以此类推，sqrt（n）~ n 的范围内的因数都可以通过 1 ~ sqrt（n） 判断出来。

所以优化后的最终判断范围为1 ~ sqrt（n ）。

简单应用

1.判断 n 是否为质数

判断范围从 1 ~ sqrt（n），如果其中存在除 1 外 n 的因数则跳出循环

通过最终循环变量 i 的值是否大于sqrt（n)判断是中途跳出循环（不是质数）还是整个范围内都没有找到 n 的因数（质数）。

#include <stdio.h>
#include <math.h>
int main()
{
    int n, i;
    scanf("%d", &n);
    for (i = 1; i <= sqrt(n); i++) {
        if (i != 1 && n % i == 0) break;
    }
    if (n != 1 && i > sqrt(n)) printf("%d是质数", n);
    else printf("%d不是质数", n); 
    return 0;
}

优化：除2外的偶数一定不是质数，且奇数的因数肯定不是偶数，因此可以先判断 n 的奇偶，奇数判断范围可以去掉其中的偶数部分

#include <stdio.h>
#include <math.h>
int main()
{
    int n, i;
    scanf("%d", &n);
    if (n == 2) printf("%d是质数", n);
    else if (n % 2 == 0) printf("%d不是质数", n);
    else {
        for (i = 1; i <= sqrt(n); i += 2) {
            if (i != 1 && n % i == 0) break;
        }
        if (n != 1 && i > sqrt(n)) printf("%d是质数", n);
        else printf("%d不是质数", n); 
    }
    
    return 0;
}

2.求约数

利用数组res[]来储存sqrt（a） ~  a 范围的约数，count用来记录个数（方便输出）

注意：当 i 为 a 的平方根时不要重复记录

#include <stdio.h>
#include <math.h>
int main()
{
    int n;
    scanf("%d", &n);
    while (n--) {
        int a, res[400], count = 0;
        scanf("%d", &a);
        for (int i = 1; i <= sqrt(a); i++) {
            if (a % i == 0) {
                if (i != a / i) res[count++] = a / i;
                printf("%d ", i);
            }
        }
        for (int i = count - 1; i >= 0; i--) {
            printf("%d ", res[i]);
        }
        printf("\n");
    }
    return 0;
}

优化：当 a 为奇数时其因数必不是偶数。

int temp = (1 & a) + 1;

利用位运算，当 a 为偶数时，temp = 1；

当 a 为奇数时，temp = 2，就将 1 ~ sqrt（a）范围内的偶数部分排除了

#include <stdio.h>
#include <math.h>
int main()
{
    int n;
    scanf("%d", &n);
    while (n--) {
        int a, res[400], count = 0;
        scanf("%d", &a);
        int temp = (1 & a) + 1;
        for (int i = 1; i <= sqrt(a); i += temp) {
            if (a % i == 0) {
                if (i != a / i) res[count++] = a / i;
                printf("%d ", i);
            }
        }
        for (int i = count - 1; i >= 0; i--) {
            printf("%d ", res[i]);
        }
        printf("\n");
    }
    return 0;
}

3.分解质因数

 思路：用数组res[x]来记录质因数x的个数，每次都要初始化为0，temp来记录最大的质因数

对拆分数 a 从 2 开始试除，找到它的最小质因数，该质因数对应的个数res[]++；

拆分数更新为 a / 该质因数；

立即判断新的拆分数是否为质数，如果是则拆分结束，若不是则继续上述步骤

#include <stdio.h>
#include <math.h>
#include <string.h>
//判断质数的函数
int zhishu(int x)
{
    int i;
    for (i = 1; i <= sqrt(x); i++) {
        if (i != 1 && x % i == 0) break;
    }
    if (x != 1 && i > sqrt(x)) return 1;
    else return 0;
}
int main()
{
    int n;
    scanf("%d", &n);
    while (n--) {
        int a, res[200000], temp = 0;
        memset(res, 0, sizeof(res));
        scanf("%d", &a);
        //特殊情况 1 单独列出
        if (a == 1) {
            printf("1 1\n\n");
            continue;
        }
        for (int i = 2; i <= a; i++) {
            if (a % i == 0 && zhishu(i)) {
                res[i]++;
                if (i > temp) temp = i;
                a /= i;
                if (zhishu(a)) {
                    res[a]++;
                    temp = a;
                    break;
                }
                //拆分后新的a从2开始重新试除
                //因为循环有i++，故i赋值为1
                i = 1;
            }
        }
        //输出质因数及个数
        for (int i = 0; i <= temp; i++) {
            if (res[i]) printf("%d %d\n", i, res[i]);
        }
        printf("\n");
    }
    return 0;
}