数论算法总结-素数

1. 1.如果给定一个数n,要你计算小于n的素数的个数，如果n足够大，那么

素数的个数等价于x/ln(x)；

代码如下：

#include<iostream>
#include<vector>
#include<math.h>
using namespace std;
/*
这是判断素数个数的位数
*/
void pr(int n)
{
	int k = 1;
	for (int i = 0; i < n; i++)
	{
		k *= 10;
	}
	cout << static_cast<int>(log10(k / log(k)))+1 << endl;

}
void main()
{
	int n;
	cin >> n;
	pr(n);
	system("pause");

}

1. 2。哥德巴赫猜想

哥德巴赫猜想的意思是每一个偶数都可以写成一对素数之和。

我们可以以3为开头，因为以2为开头，2+3=5，2+5=7，2+素数不为偶数，进行判断。

代码如下

#include<iostream>
#include<vector>
#include<math.h>
using namespace std;
/*
歌德巴赫猜想，我们假设第一个数数是3，那么排除偶数，我们每次加2
*/
bool judge(int n)//判断是否是素数
{
	int k = sqrt(n);
	for (int i = 2; i < k; i++)
	{
		if (k%i == 0)
			return 0;
	}
	return 1;
}
int pr(int n)
{
	if (n < 6)
	{
		cout << "error";
		return 0;
	}
	else
	{
		int i = 3;//每次加2，排除偶数
		int j = n - 3;
		int k = 0;
		while (i <= j)
		{
			if (judge(i) && judge(i))
			{
				k++;
				i += 2;
				j -= 2;
			}
			else
			{
				i += 2;
				j -= 2;
			}
		}
		return k;
	}
}
void main()
{
	int n;
	cin >> n;
	cout<<pr(n);
	system("pause");

}

1. 3.素数测试，埃拉托斯尼斯筛法

代码如下：

#include<iostream>
#include<vector>
#include<math.h>
using namespace std;
/*
先进行开方，找到其中的素数，然后对其他的素数进行筛选，剩下的就是
素数
*/
void pr(int n)
{
	int k = sqrt(n);
	vector<int> a;
	bool judge = true;
	for (int i = 2; i < k; i++)//寻找
	{
		for (int j = 2; j < i; j++)
		{
			if (i%j == 0)
			{
				judge = false;
				break;
			}
		}
		if (judge)
		{
			a.push_back(i);
			cout << i << endl;
		}
		judge = true;
	}
	//筛选
	for (int i = k; i < n; i++)
	{
		for (int j = 0; j < a.size(); j++)
		{
			if (i%a[j] == 0)
			{
				judge = false;
				break;
			}
		}
		if (judge)
			cout << i << endl;
		judge = true;
	}
}
void main()
{
	int n;
	cin >> n;
	pr(n);
	system("pause");

}

1. 4。6N+1法

任何一个数都可以写成6N+1 6N+2 6N+3 6N+4 6N+5 ,我们很容易就能知道

6N+1与6N+5必然是素数，那么我们可以使用这种方法来判断是否是素数。

代码如下：

#include<iostream>
#include<vector>
using namespace std;
/*
6n+1法，如果最后余数为1或者5，那么可以判断这是一个素数

  */
bool judge(int n)
{
	int r = n % 6;
	if (r == 1 || r == 5)
	{
		return 1;
	}
	else
		return 0;

}
void main()
{
	cout << judge(15) << endl;
	system("pause");

}

1. 5.分拆素数

把一个偶数分解成不同的素数和

代码如下：

#include<iostream>
#include<vector>
#include<math.h>
using namespace std;
/*
把偶数拆分成2个不同素数
与之前拆分代码不同
这段代码没有之前的代码好，时间复杂度提高了，主要是因为，我是从2开始，其实可以从3开始
以奇数判断，可以使算法更快。
*/
bool judge(int n)
{
	int k = sqrt(n);
	if (n == 2)
	{
		return 1;
	}
	for (int i = 2; i <= k; i++)
	{
		if (n%i== 0)
			return 0;
	}
	return 1;
}
int pr(int n)
{
	if (n % 2 != 0)
	{
		cout << "error" << endl;
		return 0;
	}
	int i = 2;
	int j = n - i;
	int k = 0;
	while (i < j)
	{
		if (judge(j) && judge(i))
		{
			k++;
			cout << i << " " << j << endl;
		}
		i += 1;
		j = n - i;
	}
	return k;
}

void main()
{
	cout << pr(26) << endl;
	system("pause");

}

1. 6.素数距离

判断m和n之间的素数的距离。

我们可以利用一个数组存储这之间的素数，并且如果数组的大小为1或者0，那么

返回error

代码如下：
 

#include<iostream>
#include<vector>
using namespace std;
/*
输入两个数 m和n 找出这两个数之间素数差最小的数子，和最大的数
*/
struct Max
{
	int max = 0;
	int i;
	int j;
};
struct Min
{
	int min = 99999;
	int i;
	int j;
};
bool judge(int n)//可以使用6N+1法判断更加快
{
	if (n == 1)
	{
		return 0;
	}
	if (n == 2)
	{
		return 1;
	}
	int k = sqrt(n);
	for (int i = 2; i <= k; i++)
	{
		if (n%i == 0)
		{
			return 0;
		}
	}
	return 1;
}
void pr(int m, int n)
{
	vector<int> a;
	for (int i = m; i <= n; i++)
	{
		if (judge(i))
		{
			a.push_back(i);
		}
		else
		{
			continue;
		}
	}
	if (a.size() == 0||a.size()==1)
	{
		cout << " there are no adjacent primes" << endl;
		return;
	}
	Max my_max;
	Min my_min;
	for (int i = 0; i < a.size()-1; i++)
	{
		int key = a[i + 1] -a[i];
		if (key>my_max.max)
		{
			my_max.max = key;
			my_max.i = a[i];
			my_max.j = a[i+1];
		}
		if (key<my_min.min)
		{
			my_min.min = key;
			my_min.i = a[i];
			my_min.j = a[i+1];
		}
	}
	cout << my_max.i << " " << my_max.j <<" "<< my_min.i << " " << my_min.j << endl;
}
void main()
{
	pr(14, 17);
	system("pause");

}

1. 7.判断n!的0的个数

判断这个东西需要使用算数基本定理中的一种定理

n!中素数p的个数（幂）为

[n/p]+[n/p^2]+[n/p^3]+....

注意[]是计算机中默认的向下取整也就是4.5与4.9在计算机int中都是4，而向上取整则是5；

代码如下

#include<iostream>
#include<vector>
using namespace std;
/*计算n!末尾0的个数
由于我们计算的是n!所以我们可以知道2的因子个数
肯定大于5，而一个2与一个5，可以构成一个0，所以我们是
求5的因子个数。
*/
int pr(int n)
{
	int sum = 0;
	int k = 5;
	int r = n / k;
	while (r)
	{
		sum += r;
		k *= 5;
		r = n / k;
	}
	return sum;
}

void main()
{
	cout << pr(1024) << endl;
	system("pause");

}

1.8 算数定理判断一个数的正因子的个数与所以因子的和

任何一个正数可以写成素数的积

n=p1^a1+p2^a2+...+pn^an;

设m为一个正数的正因子的个数

m=(a1+1)*(a2+1)*...*(an+1);

设n为正因子的总数和

n=(p1^(a1+1)-1)/(p1-1)*...*(pn&(an+1)-1)/pn-1;

代码如下：

#include<iostream>
#include<vector>
using namespace std;
/*
利用算数定理求正因子的个数 
*/
bool judge(int n)
{
	if (n == 1)
	{
		return 0;
	}
	if (n == 2)
	{
		return 1;
	}
	int k = sqrt(n);
	for (int i = 2; i <= k; i++)
	{
		if (n%i == 0)
		{
			return 0;
		}
	}
	return 1;
}
int pr(int n)
{
	int i = 2;
	vector<int> a;
	int k = sqrt(n);//注意相乘的素数应该在这里
	for (int i = 2; i <= k; i++)
	{
		if (judge(i))
		{
			a.push_back(i);
		}
	}
	int sum = 1;
	int r = n;
	int val = 0;
	for (int i = 0; i < a.size(); i++)
	{
		int m = r%a[i];//判断是否可以存在这个素数因子
		while (!m)
		{
			val++;
			r = r / a[i];
			m = r%a[i];
		}
		sum *= (val + 1);
		val = 0;
		if (r == 1)//最后一个素数因子
			break;
	}
	return sum;
}
void main()
{
	cout << pr(10780) << endl;
	system("pause");

}