【算法】常用的基础数论

时间:2021-01-12 01:20:05

作者:指针不指南吗
专栏:算法篇

????或许会很慢,但是不可以停下????

1.GCD&LCM

最大公约数&最小共倍数

欧几里得算法——高效

//最大公约数
int gcd(int x,int y)
{
	return y?gcd(y,x%y):x;
}

//最小公倍数
int lcm(int x,int y)
{
	return x*y/gcd(x,y);
}

2.判断素数(质数)

孪生素数法——高效

  • 分析过程如下

首先看一个关于质数分布的规律:大于等于5的质数一定和6的倍数相邻。例如5和7,11和13,17和19等等;

证明:令x≥1,将大于等于5的自然数表示如下:
······ 6x-1,6x,6x+1,6x+2,6x+3,6x+4,6x+5,6(x+1),6(x+1)+1 ······
可以看到,不在6的倍数两侧,即6x两侧的数为6x+2,6x+3,6x+4,由于2(3x+1),3(2x+1),2(3x+2),所以它们一定不是素数,再除去6x本身,显然,素数要出现只可能出现在6x的相邻两侧。这里要注意的一点是,在6的倍数相邻两侧并不是一定就是质数

此时用一个循环判断质数,可以以6为单元快进,i+=6,加快判断速度。

原因是,假如要判定的数为n,则n必定是6x-1或6x+1的形式,对于循环中6i-1,6i,6i+1,6i+2,6i+3,6i+4,其中如果n能被 6i,6i+2,6i+4整除,则n至少得是一个偶数,但是6x-1或6x+1的形式明显是一个奇数,故不成立;另外,如果n能被6i+3整除,则n至少能被3整除,但是6x能被3整除,故6x-1或6x+1(即n)不可能被3整除,故不成立。

综上,循环中只需要考虑6i-1和6i+1的情况,即循环的步长可以定为6,每次判断循环变量k和k+2的情况即可

bool isprime(int n)
{
	if(n<=3)   //特判几个较小的数
		return n>1;
		
	if(n%6!=1&&n%6!=5)  //不在6的倍数的两侧,一定不是素数
		return 0;
		
	for(int i=5;i<=sqrt(n);i+=6)  //判断在6的倍数的两侧的数,是不是素数
		if(n%i==0||n%(i+2)==0)
			return 0;
			
	return 1;
}

3.分解质因子

直接看例题吧

  • 题目

    链接: AcWing 867. 分解质因数- AcWing

    给定 n 个正整数 ai,将每个数分解质因数,并按照质因数从小到大的顺序输出每个质因数的底数和指数。

    输入格式

    第一行包含整数 n。

    接下来 n 行,每行包含一个正整数 ai。

    输出格式

    对于每个正整数 ai,按照从小到大的顺序输出其分解质因数后,每个质因数的底数和指数,每个底数和指数占一行。

    每个正整数的质因数全部输出完毕后,输出一个空行。

    数据范围

    1≤n≤100,
    2≤ai≤2× 1 0 9 10^9 109

    输入样例:

    2
    6
    8
    

    输出样例:

    2 1
    3 1
    
    2 3
    
  • 分析

    • x 的质因子最多只包含一个大于 根号x 的质数。如果有两个,这两个因子的乘积就会大于 x,矛盾。
    • i 从 2 遍历到 根号x。 用 x / i,如果余数为 0,则 i 是一个质因子。
    • s 表示质因子 i 的指数,x /= i 为 0,则 s++, x = x / i 。
    • 最后检查是否有大于 根号x 的质因子,如果有,输出。
  • 代码实现

    #include<bits/stdc++.h>
    using namespace std;
    
    void divide(int x)
    {
        for (int i = 2; i <= x / i; i ++ )//i <= x / i:防止越界,速度大于 i < sqrt(x)
            if (x % i == 0)//i为底数
            {
                int s = 0;//s为指数
                while (x % i == 0) x /= i, s ++ ;
                cout << i << ' ' << s << endl;//输出
            }
        if (x > 1) cout << x << ' ' << 1 << endl;//如果x还有剩余,单独处理
        cout << endl;
    }
    
    int main()
    {
        int n;
        cin >> n;
        while (n -- )
        {
            int x;
            cin >> x;
            divide(x);
        }
    
        return 0;
    }
    

【算法】常用的基础数论