最大公约数与最小公倍数
最大公约数
辗转相除法:
int gcd(int a,int b){
if(b==0) return a;
else return gcd(b,a%b);
}
最小公倍数
在求出a b的最大公约数d之后,可以马上得到a b的最小公倍数是ab/d
素数
素数的判断
判断一个数字从2到其根号是否能被它整除
bool isPrime(int n) {
if (n <= 1)return false;//注意1是不算质数的
int sqr = (int)sqrt(1.0*n);
for (int i = 2; i <= sqr; i++) {
if (n%i == 0)return false;
}
return true;
}
可以基于这个配合枚举列出1-100的所有素数:
#include<iostream>
using namespace std;
const int maxn = 101;
int prime[maxn], pNum = 0;//素数表和素数个数
//bool p[101] = { 0 };
bool isPrime(int n) {
if (n <= 1)return false;//注意1不是素数
int sqr = sqrt(n);
for (int i = 2; i <= sqr; i++) {
if (n%i == 0)return false;
}
return true;
}
void find_prime() {
for (int i = 1; i < maxn; i++) {
if (isPrime(i)) {
prime[pNum++] = i;
//p[i] = true;
}
}
}
int main() {
find_prime();
for (int i = 0; i < pNum; i++) {
cout << prime[i] << " ";
}
return 0;
}
时间复杂度为O(n根号n)在10^5之内可以接收。更高效的算法:筛法 O(nloglogn)
例子:
筛选法
#include<iostream>
using namespace std;
const int maxn = 101;
int prime[maxn], pNum = 0;//素数表和素数个数
bool p[maxn];//默认是素数
void find_prime() {
for (int i = 2; i < maxn; i++) {
if (p[i] == true) {
//是素数
prime[pNum++] = i;
for (int j = i + i; j < maxn; j += i) {
p[j] = false;//筛掉i的倍数
}
}
}
}
int main() {
fill(p, p + maxn, true);
find_prime();
for (int i = 0; i < pNum; i++) {
cout << prime[i] << " ";
}
return 0;
}
质因子分解
将一个整数n写成一个或者多个质数的乘积
struct factor{
int x,cnt;
}fac[10];
为什么十个就够了 因为2 3 5 7 11 13 17 19 23 29这些数乘起来就已经超过int范围了。
结论:对一个正整数n,如果存在1和本身之外的因子,那么一定是在sqrt(n)的左右成对出现。如果它在[2,n]内存在质因子,要么这些质因子全部小于sqrt(n),要么只存在一个大于sqrt(n)的质因子,其余质因子全部小于等于sqrt(n)
if(n%prime[i]==0){
//如果prime[i]是n的因子
fac[num].x=prime[i];
fac[num].cnt=0;
while(n%prime[i]==0){
fac[num].cnt++;
n/=prime[i];
}
num++;
}
if(n!=1){
//如果无法被根号n以内的质因子除尽
fac[num].x=n;//一定有一个大于根号n的质因子
fac[num++].cnt=1;
}
综合应用:给一个长整型(64位)求其分解的因式
#include<vector>
#include<iostream>
using namespace std;
const int maxn = 2001;//判断上限
int prime[1000];
int primeNum = 0;
bool p[maxn];
struct factor {
int x, cnt;
};
vector<factor> facList;
//筛选法求质数表,注意筛选法是每个数都要看
void getPrime() {//生成2000以内质数
fill(p, p + maxn, true);//默认都是质数
for (int i = 2; i < maxn; i++) {
if (p[i] == true) {
prime[primeNum++] = i;
for (int j = i + i; j < maxn; j += i) {
p[j] = false;//i的倍数都不是质数
}
}
}
}
//求因式
void get_prime_re(long long n) {
factor tempFac;
if (n <= 2) {
tempFac.x = n;
tempFac.cnt = 1;
facList.push_back(tempFac);
return;
}
for (int i = 0; i < primeNum; i++) {
if (n%prime[i] == 0) {//如果可以呗prime[i]整除
tempFac.x = prime[i];//指定质因子种类
tempFac.cnt = 0;//初始化个数
while (n%prime[i] == 0) {//算这个质因子的个数
tempFac.cnt++;
n /= prime[i];
}
facList.push_back(tempFac);//加到总表
}
if (n == 0)break;
}
}
int main() {
//输入
long long inputNum;
cin >> inputNum;
getPrime();
get_prime_re(inputNum);
cout << inputNum << "=";
for (int i = 0; i < facList.size(); i++) {
if (facList[i].cnt != 0) {
if (facList[i].cnt == 1) {
cout << facList[i].x;
}
else {
cout << facList[i].x << "^" << facList[i].cnt;
}
if (i != facList.size() - 1)cout << "*";
}
}
return 0;
}