简单数论

算术基本定理：
任何一个自然数N可以分解成有限个质数的乘积 N=P1^a1*P2^a2*P3^a3…*Pn^an
质数：
定义：在大于1的整数当中，只包含1和它本身能被该数字整除。
1。试除法：(从定义入手解决该问题)

#include<iostream>
using namespace std;
int x;
bool check_primes(int x){
    if(x<2) return false;
    for(int i=2;i<x/i;i++){  //判断一半就行
        if(x%i==0) return false;
    }
    return true;
}
int main(){
    cin>>x;
    if(check_primes(x)==true){
        puts("Yes\n");
    }
    else puts("No\n");
    return 0;
}

分解质因数
其中质数n分解出来的质因数比根号n大的应该只有一个因为如果有两个比根号n大，则不成立。

#include<iostream>
using namespace std;
const int N=110;
void divide(int x ){
    for(int i=2;i<=x/i;i++){
        if(x%i==0){
        int ans=0;
        while(x%i==0){
            x=x/i;
            ans++;
        }
        cout<<i<<" "<<ans<<endl;
    }
    }
    if(x>1) cout<<x<<" "<<"1"<<endl;
    cout<<endl;
}
int main(){
    int n;
    int x;
    cin>>n;
    while(n--){
        cin>>x;
        divide(x);
    }
    return 0;
}

筛法求质数：

/朴素筛法
#include<iostream>
using namespace std;
const int N=1000010;
int cnt;
bool ck[N];
void primes(int n){
    for(int i=2;i<=n;i++){
        if(ck[i]==false) cnt++;
        for(int j=i+i;j<=n;j=j+i){
            ck[j]=true;
        }
    }
}
int main(){
    int n;
    cin>>n;
    primes(n);
    cout<<cnt<<endl;
}
//埃氏筛法
#include<iostream>
using namespace std;
const int N=1000010;
int cnt;
bool ck[N];
void primes(int n){
    for(int i=2;i<=n;i++){
        if(ck[i]==false){
        cnt++;
        for(int j=i+i;j<=n;j=j+i){
            ck[j]=true;
        }
        }
    }
}
int main(){
    int n;
    cin>>n;
    primes(n);
    cout<<cnt<<endl;
}
//线性筛法

#include<iostream>
using namespace std;
const int N=1000010;
int cnt;
int pri[N];
bool ck[N];
void primes(int n){
    for(int i=2;i<=n;i++){
        if(!ck[i]) pri[cnt++]=i;
        for(int j=0;pri[j]*i<=n;j++){
            ck[ pri[j]*i ]=true;
            if(i%pri[j]==0) break;
        }
    }
}
int main(){
    int n;
    cin>>n;
    primes(n);
    cout<<cnt<<endl;
}

约数个数定理:
(a1+1)(a2+1)(a3+1)…(an+1)
约数个数：

#include<iostream>
#include<cstdio>
#include<unordered_map>
using namespace std;
const int N=1e9+7;
int n;
unordered_map<int,int>m;
void answer(int x){
    for(int i=2;i<=x/i;i++){  //优化算法遍历到一半，每个n只有一个分解出来的质因数是比根号n大的
        if(x%i==0){
            while(x%i==0){
                x=x/i;
                m[i]++;
            }
        }
    }
    if(x>1) m[x]++;
}
int main(){

    int n;
    int x;
    cin>>n;
    while(n--){
        cin>>x;
        answer(x);  //分解质因数
    }
    long long ans=1;
    for(auto t:m){
        ans=ans*(t.second+1)%N;
    }
    cout<<ans<<endl;
}  

约数之和定理：
（1+p1¹+p1²…p1^a1）(1+p2¹+p2²+…+p2^a2)(1+pn¹+pn²+…+pn^an)






约数之和~：

#include<iostream>
#include<cstdio>
#include<unordered_map>
using namespace std;
const int N=1e9+7;
int n;
unordered_map<int,int>m;
void answer(int x){
    for(int i=2;i<=x/i;i++){  //优化算法遍历到一半，每个n只有一个分解出来的质因数是比根号n大的
        if(x%i==0){
            while(x%i==0){
                x=x/i;
                m[i]++;
            }
        }
    }
    if(x>1) m[x]++;
}
int main(){

    int n;
    int x;
    cin>>n;
    while(n--){
        cin>>x;
        answer(x);  //分解质因数
    }
    long long ans=1;
    for(auto t:m){
        long long a=t.first;
        long long b=t.second;
        long long p=1;
        while(b--){
            p=(p*a+1)%N;
        }
        ans=ans*p%N;
    }
    cout<<ans<<endl;
}  

求两个数的最大公约数~

#include<iostream>
#include<cstdio>
using namespace std;
int gcb(int a , int b){
    return b?gcb(b,a%b):a;
}
int main(){
    int a,b;
    cin>>a>>b;
    cout<<gcb(a,b)<<endl;
}

筛法求一定范围内的质数和一个数的最小质因子：

#include<cstdio>
#include<iostream>
using namespace std;
const int N=100010;
bool st[N];   //当前数有没有被筛过
int primes[N],cnt; //存所有质数
int minp[N]; //存最小质因子5
int get_primes(int n){
    for(int i=2;i<=n;i++){
        if(!st[i])minp[i]=i,primes[cnt++]=i; //如果没有被筛过表示他是一个质数
        for(int j=0;primes[j]*i<=n;j++){
            st[primes[j]*i]=true;
            //筛掉的一定是和数
            minp[primes[j]*i]=primes[j];
            if(i%primes[j]==0)break;
        }
    }
}
int main(){
    get_primes(100);
        for(int i=0;i<100;i++){
        printf("%d ",primes[i]);
    }
}

本文章的思想及代码均根据Acwing y总的讲解总结的hhh