算法中的数学知识

质数

判断质数的方法

模板

思路还是很简单的，质数的特点是只能被一和自身整除，那么从小到大枚举是否能被整除即可。

#include <iostream>
using namespace std;
void judge(int n) {
    if(n == 1) {
        printf("No\n");
            return;
    }
    for(int i = 2; i <= n / i; i ++) {
        if(n % i == 0) {
            printf("No\n");
            return;
        }
    }
    printf("Yes\n");
    
}
int main() {
    int m;
    cin >> m;
    while(m --) {
        int n;
        scanf("%d", &n);
        judge(n);
    }
    return 0;
}

分解质因数

首先补一个概念：质数是只有一和自身作为约数的数，而合数一定可以被几个质数作为约数，也就是一个合数一定可以由若干个质数相乘得来。

如果给定一个合数，要求给出所有可能出现的质因数搭配，应该怎么做？

其实还是简单的，从2开始枚举，每次枚举一直把当前数除到不能整除为止，然后加一，直到原数被除为1，代码如下：

#include <iostream>
using namespace std;
int main(){
    int n, m, s;
    cin >> n;
    while(n --) {
        scanf("%d", &m);
        for(int i = 2; i <= m / i; i ++) {
            if(m % i == 0) {
                s = 0;
                while(m % i == 0) {
                    s ++;
                    m /= i;
                }
                printf("%d %d\n", i, s);
            }        
        }
        if (m > 1) cout << m << ' ' << 1 << endl;
        printf("\n");
    }
    return 0;
}

判断2~n间有多少质数模板

根据上面，一个合数一定可以由若干个质数相乘得来，那么，如果一个数无法由质数相乘得来，就一定是质数。

如下，先把二的倍数去除：

再去除3：

五：

用一个数组来记录当前数字是否被枚举过倍数，如果没有，说明这是质数，当前质数的数量加一，并枚举他的所有倍数，2~n枚举完毕后，即可得出质数数量。

模板

#include <iostream>
#include <algorithm>

using namespace std;

const int N= 1000010;

int primes[N], cnt;
bool st[N];

void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ ) {
        if (st[i]) continue;
        primes[cnt ++ ] = i;
        for (int j = i + i; j <= n; j += i)
            st[j] = true;
    }
}
int main() {
    int n;
    cin >> n;
    get_primes(n);
    cout << cnt << endl;
    return 0;
}

约数

最大公约数

对于此类经典题型，有一个一行代码模板，使用的是欧几里得算法，也叫辗转相除法，记住可以事半功倍。

那么，结论在前，原理在后，先上代码：

#include <iostream>
#include <algorithm>

using namespace std;

int gcd(int a, int b)
{
    return b ? gcd(b, a % b) : a;
}

int main()
{
    int n;
    cin >> n;
    while (n -- )
    {
        int a, b;
        scanf("%d%d", &a, &b);
        printf("%d\n", gcd(a, b));
    }

    return 0;
}

核心原理在于：

如果 $\frac{d}{a}$ 并且 $\frac{d}{b}$ ，那么 $\frac{d}{ax+by}$。

若 $d$ 整除于 $a$，也整除于$b$，那么也整除于 $a$ 的若干倍数加 $b$ 的若干倍数。

注：整除不是整除以

整除：2能整除6
整除以：6整除以2（6能被2整除）

---

其他相信大家多少了解过，网络上也有很多版本的讲解了，本人也没有什么创新的思路，本着能原创就不抄袭的原则，如果有疑问可以看这篇大佬的博客：指路

快速幂

思路

快速幂是指快速算出一个数的幂运算，如 $2^{32}$，是 $2\ast 2\ast 2\ast 2\ast \dots \ast 2$，共乘了32次，但是很显然 $2^{32}$ 是可以由 $2^{16}$ 平方得来，而 $2^{16}$ 也可以由 $2^8$ 平方得来，以此类推。

例如 $3^{29}=3^1+3^4+3^8+3^{16}$

但这是如何得来的呢？

因为29的二进制表示为11101，一所在的位是幂需要累乘的，以下为该思路代码：

模板

#include <cstdio>
using namespace std;

// 求快速幂
int QuickPower(int x,int n){
    int answer = 1;
    while (n != 0) {
        if (n % 2 == 1) {
            answer *= x;
        }
        n /= 2;
        x *= x;
    }
    return answer;
}
int main(){
    int a,b;
    printf("%d\n", QuickPower(a, b));
    return 0;
}