数学知识(一)

一、数论

1.1质数

**定义:**在所有大于1的自然数,如果只包含1和本身这两个约数,就被称为质数(素数).

**质数的判断:**试除法

cpp 复制代码
bool is_prime(int n)
{
    if(n < 2) return false;
    for(int i = 2;i <= n / i;i ++)
    {
        if(n % i == 0)
            return false;
    }
    return true;
}

分解质因数

基本思路:从小到大枚举n的所有约数。

cpp 复制代码
void divide(int n)
{
    for(int i = 2;i <= n;i ++)
    {
        if(n % i == 0)
        {
            int s = 0;
            while(n % i == 0)
            {
                n /= i;
                s ++;
            }
            printf("%d %d\n",i,s)
        }
    }
}

埃氏筛法

cpp 复制代码
#include<iostream>
#include<algorithm>

using namespace std;

int primes[N],cnt;
bool st[N];

void get_prime(int n)
{
    for(int i = 2;i <= n;i ++)
    {
        if(!st[i])
        {
            primes[cnt ++] = n;
            for(int j = i + i;j <= n;j += i) st[j] = true;
        }
    }
}

线性筛法

保证了每个数都是被他的最小质因子所筛掉

cpp 复制代码
#include<iostream>
#include<algorithm>

using namespace std;

int primes[N],cnt;
bool st[N];

void get_primes(int n)
{
    for(int i = 2;i<= n;i ++)
    {
        if(!st[i]) primes[cnt ++] = i;
        for(int j = 0;primes[j] <= n / i;j ++)
        {
            st[primes[j] * i = true;
            if(i % primes[j] == 0) break;  //primes[j]一定是i的最小质因子
         }
    }
}

1.2约数

试除法求约数

cpp 复制代码
vector<int> get_divisors(int n)
{
    vector<int> res;
    
    for(int i = 1;i <= n / i;i ++)
        if(n % i == 0)
        {
            res.push_back(i);
            if(i != n / i) res.push_back(n / i);
        }
    sort(res.begin(),res.end());
    return res;
}

约数个数与 约数之和

欧几里得算法

cpp 复制代码
//返回a与b的最大公约数
int gcd(int a,int b)
{
    return b ? gcd(b,(a % b) : a);
}

二、欧拉定理

2.1欧拉函数

定义:X(n)表示1~n中与n互质的个数

算法步骤:

  1. 从1~N中去掉p1,p2,...,pk的倍数
  2. 加上所有Pi * Pj的倍数
  3. 减去所有Pi * Pj * Pk的倍数
cpp 复制代码
#include<iostream>
#include<cstring>
#include<algorithm>

using namespace std;

int main()
{
    int n;
    cin>>n;

    while(n --)
    {
        int a;
        cin>>a;

        int res = a;
        //分解质因数
        for(int i = 2 ;i <= a;i ++)
        {
            if(a % i == 0)
            {
                res = res / i * (i - 1);
                while (a % i == 0) a /= i;
            }
        }
        if(a > 1) res = res / a * (a - 1);

        cout<<res<<endl;
    }
    return 0;
}

快速幂

cpp 复制代码
#include<iostream>
#include<algorithm>

using namespace std;

typedef long long LL;

//求a^k % p
int main(int a,int k,int p)
{
    int res = 1;
    while(k)
    {
        //如果当前k的末位为1,则
        if( k & 1) res = (LL)res * a % p;
        //删除k的末位
        k >>= 1;
        //把a平方
        a = (LL) a * a % p;
    }
    return res;
}
int main()
{
    int n;
    scanf("%d",&n);
    
    while(n --)
    {
        int a,k,p;
        scanf("%d%d%d",&a,&k,&p);
        
        printf("%d\n",qmi(a,k,p));
    }
    return 0;
}
相关推荐
sky_ph2 分钟前
JAVA-GC浅析(二)G1(Garbage First)回收器
java·后端
开发者工具分享6 分钟前
如何应对敏捷转型中的团队阻力
开发语言
gregmankiw12 分钟前
C#调用Rust动态链接库DLL的案例
开发语言·rust·c#
IDRSolutions_CN24 分钟前
PDF 转 HTML5 —— HTML5 填充图形不支持 Even-Odd 奇偶规则?(第二部分)
java·经验分享·pdf·软件工程·团队开发
hello早上好27 分钟前
Spring不同类型的ApplicationContext的创建方式
java·后端·架构
roman_日积跬步-终至千里27 分钟前
【Go语言基础【20】】Go的包与工程
开发语言·后端·golang
music&movie43 分钟前
算法工程师认知水平要求总结
人工智能·算法
秦少游在淮海1 小时前
C++ - string 的使用 #auto #范围for #访问及遍历操作 #容量操作 #修改操作 #其他操作 #非成员函数
开发语言·c++·stl·string·范围for·auto·string 的使用
const5441 小时前
cpp自学 day2(—>运算符)
开发语言·c++
心扬1 小时前
python生成器
开发语言·python