【笔试相关】数论相关基础算法

本文深入讲解了数学算法中的关键概念,包括质因数分解、最大公约数、最小公倍数的计算方法,以及快速幂和素数判断算法。提供了详细的代码实现,如辗转相除法、辗转相减法和穷举法等。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

质因数分解

vector<int> dividenum(int n)
{
	vector<int> ret;
	int n_sqrt = sqrt(n);
	//循环截至条件到根号n就行
	//从2开始,遇到的第一个因数,一定是质数
	//比如24,遇到了2,就一直除,直到余数不为0
	//有人说,4也是24的因数啊,但是不是质数
	//这个在遇到2时,一直除,这里已经将4除掉了
	//24已经变成3了
	for (int i = 2;i <= n_sqrt;i++)
	{
		if (n%i == 0)//
		{
			while (n%i == 0)
			{
				n /= i;
				ret.push_back(i);
			}
		}
	}
	//像24这样的数 ,24 = 2*2*2*3,最后n变成了1
	//但是像22这样的数,22=2*11,但是11>sqrt(22),所以需要判断一下
	//又如11这样的数,11 = 11*1,这个判断也能处理
	if (n != 1) ret.push_back(n);
	return ret;
}

最大公约数

辗转相除法

辗转相除法, 又名欧几里德算法(Euclidean algorithm),是求最大公约数的一种方法。它的具体做法是:用较小数除较大数,再用出现的余数(第一余数)去除除数,再用出现的余数(第二余数)去除第一余数,如此反复,直到最后余数是0为止。如果是求两个数的最大公约数,那么最后的除数就是这两个数的最大公约数。

//辗转相除法 

int measure(int x, int y)
{

	int z = y;
	while (x%y != 0)
	{
		z = x%y;
		x = y;
		y = z;
	}
	return z;
}

辗转相减法

“较大数” 减 “较小数”,循环,当两数相同时,相同的数即为“最大公约数”

//辗转相减法 
int measure(int a, int b)
{
	while (a != b)
	{
		if (a>b)
		{
			a = a - b;
		}
		else
		{
			b = b - a;
		}
	}
	return a;
}

穷举法

int measure(int x, int y)
{
	int temp = 0;
	for (temp = x; ; temp--{
		if (x%temp == 0 && y%temp == 0)

			break;
	}
	return temp;
}

最小公倍数

通过最大公约数求解

int minmulti(int x, int y)
{
		int res = x*y;
		int z = y;
		while (x%y != 0)
		{
			z = x%y;
			x = y;
			y = z;
		}
		return res/z;
}

穷举法

void minmulti()
{
	int m;  //定义一个标记
	int a, b, c;  //定义三个正整数
	cout << "输入三个正整数" << endl;
	cin >> a >> b >> c;
	m = a>b ? (a>c ? a : c) : (b>c ? b : c);  //令m等于a,b,c的最大值
	while ((m%a != 0) || (m%b != 0) || (m%c != 0))
		m++;
	cout << "最小公倍数为:" << m << endl;
}

快速幂

过程描述

在这里插入图片描述

实现

int Pow(int a,int b){
    int ans = 1;
    int base = a;
    while(b){
        if(b & 1) ans *= base;
        base *= base;
        b >>= 1;
    }
    return ans;
}

素数的快速判断

bool is_Prime( int num )
{
    //快速判断一个数是否为素数

    //两个较小数另外处理
    if(num ==2|| num==3 )
        return 1 ;
    //不在6的倍数两侧的一定不是质数
    if(num %6!= 1&&num %6!= 5)
        return 0 ;
    int tmp =sqrt( num);
    //在6的倍数两侧的也可能不是质数
    for(int i= 5; i <=tmp; i+=6 )
        if(num %i== 0||num %(i+ 2)==0 )
            return 0 ;
    //排除所有,剩余的是质数
    return 1 ;
}

判断N以内有多少个质数

找出n以内的所有素数:做法:做法其实很简单,首先将2到n范围内的整数认定为素数,其中2是最小的素数。将表中所有的2的倍数划去,表中剩下的最小的数字就是3,他不能被更小的数整除,所以3是素数。再将表中所有的3的倍数划去……以此类推,如果表中剩余的最小的数是m,那么m就是素数。然后将表中所有m的倍数划去,像这样反复操作,就能依次枚举n以内的素数,这样的时间复杂度是O(nloglogn)。
const int MAXN=1000;
int prime[MAXN];//记录素数
bool is_pri[MAXN + 10]; //is_pri[i]表示i是否为素数
//返回n以内素数的个数
int sieve(int n)
{
    int p = 0;
    for(int i = 0; i <= n; i++)
        is_pri[i] = true; //初始化全部是素数
    is_pri[0] = is_pri[1] = false;//0和1不是素数
    for(int i = 2; i <= n; i++)
    {
        if(is_pri[i])
        {
            prime[++p] = i;
            for(int j = 2 * i; j <= n; j += i)
                if(is_pri[j])
                is_pri[j] = false;//i的倍数全部剔除
        }
    }
    return p;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值