数学专题：同余定理、gcd与lcm、中位数定理、筛法、快速幂等

数学专题：同余定理、gcd与lcm、中位数定理、筛法、快速幂等

Java中数学相关的库函数

  Math类提供了一系列数学方法，这里列了两个不常用的

方法	描述
Math.sin(弧度)	求sin值，参数为弧度而不是角度
Math.PI	PI值

同余定理

原理

(a±b)%mod = (a%mod±b%mod)%mod
(a*b)%mod = (a%mod*b%mod)%mod
  当计算结果太大，计算机无法表示时，题目会让我们对这个数取模再输出。但如果在计算这个数的过程中，中间结果就已经非常大，计算机无法表示时，这时该怎么办？此时同余定理就登场了，只要最终结果是由加减或乘法运算得到的，那么在每一步计算中我们都可以进行取模，最终的结果不变。

例题

解题思路

  对于无法用long类型表示的非常大的数，我们一般都是用字符串表示，对于该题，最终我们还是需要将字符串变为数字，这时我们需要知道将字符串变为数字的方法
  一般的想法是将字符串从后向前遍历，将数字的每位乘以其权重，如"12345" = 5 * 10 ^ 0 + 4 * 10 ^ 1 + 3 * 10 ^ 2 ...。但实际上我们从前向后遍历更快，也更易编写程序，我们使用ans表示转换后得到的数字，其初值为0，每遍历到一个字符，我们将ans * 10 + 当前字符表示的数字。例如"12345"的处理过程为：

s=12345

ans = 0 * 10 + 1         1 %  p

ans = 1 * 10 + 2        12 % p

ans = 12 * 10 + 3        123 % p

...

  这样，最终的数字是由乘法和加法运算得到的，利用同余定理，我们可以在每一小步都进行取余操作，这样中间结果也不会溢出，最终将得到的数字取余即可

代码

import java.util.*;

public class Main {
   
    public static void main(String[] args) {
   
        String N;
        long p, num;
        int i;
        Scanner scanner = new Scanner(System.in);
        while(scanner.hasNext()) {
   
            N = scanner.next();
            p = scanner.nextLong();
            num = 0;
            for (i = 0;i < N.length();i++) {
   
                num = ((num * 10 % p) + (N.charAt(i) - '0')) % p;
            }
            System.out.println(num);
        }
    }
}

gcd与lcm

原理

  gcd表示最大公约数，lcm表示最小公倍数，a和b的最小公倍数等于a * b / gcd(a,b)，因此只要能求最大公约数，最小公倍数也能求得
  求最大公约数的方法为辗转相除法，具体做法为：用较大数除以较小数，再用较小数（除数）除以出现的余数(第一余数)，再用第一余数除以出现的余数(第二余数)，如此反复，直到最后余数是0为止。如果是求两个数的最大公约数，那么最后的除数就是这两个数的最大公约数。
  比如说:计算15和36的最大公约数。用辗转相除法是这么计算的：

36÷15=2······6 (第一余数)
15÷6=2······3 (第二余数)
6÷3=2······0(第三余数)

最终结果为3

  用符号来表示，即：

a ÷ b = a / b ... a % b
接下来令 a = b，b = a % b，再执行上述公式，以此类推，
直到a % b == 0，此时b即为最大公约数

  用递归式来表示，即：
gcd(a,b) = gcd(b, a % b)
  注意a,b的大小顺序是不用考虑的，若a小b大，执行一次后又会换回来，例如15 ÷ 36 = 0 ... 15，第二步变成36 ÷ 15 = 2 ... 6，即又换了回来
  gcd的算法程序如下：

int gcd(int a, int b) {
   
	int c;
	while (b != 0) {
   
		c = a % b;
		a = b;
		b = c;
	}
	return a;
}

例题

解题思路

  多个数的最小公倍数，就相当于不断求两个数间的最小公倍数

LCM(a1,a2,a3) = LCM(LCM(a1,a2),a3)
…
LCM(a1,a2,…,an) = LCM(LCM(LCM(a1,a2),a3)…)

  lcm(a,b) = a * b / gcd(a, b)

代码

import java.util.*;
 
public class Main {
   
    public static void main(String[] args) {
   
        long num, t, n, i, j, g = 1;
        long m = 1;
        Scanner scanner = new Scanner(System.in);
        t = scanner.nextLong();
        for (i = 0;i < t;i++) {
   
            m = 1;
            n = scanner.nextLong();
            for (j = 0;j < n;j++) {
   
                num = scanner.nextInt();
                m = m * num / (gcd(m, num));
            }
            System.out.println