数论初步学习

最新推荐文章于 2024-07-21 21:42:58 发布

weixin_52109411

最新推荐文章于 2024-07-21 21:42:58 发布

阅读量116

点赞数 1

分类专栏：数论文章标签：线性代数算法

本文链接：https://blog.csdn.net/weixin_52109411/article/details/120806717

版权

数论专栏收录该内容

4 篇文章 0 订阅

订阅专栏

队伍需要，我开始学习数论
错过暑假，懊悔万分，尽快弥补

埃氏筛法与线性筛法

int pri[100010]={2},vis[100010],cnt;
for(int i=2;i<=100000;i++) {
	if(!vis[i]){
        for(int j=i;j<=100000;j+=i)  vis[j]=1;
        pri[++cnt]=i;
	}
}//cnt==9592

for(int i=2;i<100000;i++)
{
    if(!check[i])pri[tot++]=i;
    for(int j=0; j<tot,i*pri[j]<=100000;j++)
    {
        check[i*pri[j]]=1;
        if(i%pri[j]==0)break;
    }
}

质因数分解

任何一个大于1的自然数，如果N不为质数，都可以唯一分解成有限个质数P的乘积
在这里插入图片描述

	for(int i=1;i<=cnt;i++){
		while(n%pri[i]==0){
			n/=pri[i];
		}
	}

扩展欧几里德

关于x,y的方程 ax + by = gcd(a,b) 的所有整数解
求关于x的同余方程 ax≡1(mod b) 的所有整数解
上述两种叙述都可以用拓展欧几里得算法

int egcd(int a,int b){
    if(!b){
        x=1;
        y=0;
        return a;
    }
    int ans=egcd(b,a%b);
    int t=x;
    x=y;
    y=t-a/b*y;
    return ans;
}
//(x%b+b)%b   最小正整数解

乘法逆元

如果有ax≡1(mod n)，则称x是mod n意义下a的乘法逆元。
一个数有逆元<=>gcd(a,n)=1，此时逆元存在且唯一

费马小定理

在这里插入图片描述
p是素数，a是任意整数且a!=0(mod p)

欧拉定理

费马小定理的拓展
$欧拉函数φ ( n ) ( n ∈ N ∗ ) \varphi (n)(n\in N^*)φ(n)(n∈N ∗ )是小于等于 n nn 的正整数中与 n nn 互质的数的个数。$
在这里插入图片描述
性质：
1.对于质数p， φ( p) = p - 1。注意φ(1)=1
2. 若m,n互质，φ(mn)=φ(m)φ(n)
3. 当n为奇数时，φ(2n)=φ(n)
4.当n大于2时，所有φ(n)都是偶数
5.当n大于6时，所有φ(n)都是合数

中国剩余定理

在这里插入图片描述

莫比乌斯反演

莫比乌斯函数μ(d)：
如果d=1，μ(d)=1
如果d为互异质数p1,p2…pk的乘积，则μ(d)=(−1)^k，否则，μ(d)=0

常见的莫比乌斯反演有两种形式：
在这里插入图片描述

几道数论题
洛谷P1082 [NOIP2012 提高组] 同余方程
题目描述
求关于xx的同余方程 a x \equiv 1 \pmod {b}ax≡1(modb) 的最小正整数解。

输入格式
一行，包含两个正整数 a,ba,b，用一个空格隔开。

输出格式
一个正整数 x_0x
0

，即最小正整数解。输入数据保证一定有解。
拓展欧几里得算法

#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
int a,b,x,y;

int egcd(int a,int b){
    if(!b){
        x=1;
        y=0;
        return a;
    }
    int ans=egcd(b,a%b);
    int t=x;
    x=y;
    y=t-a/b*y;
    //cout<<x<<' '<<y<<' '<<a<<' '<<b<<endl;
    return ans;
}

int main(){
    ios::sync_with_stdio(0);
    cin>>a>>b;
    egcd(a,b);
    cout<<(x%b+b)%b;
	return 0;
}

洛谷P1214 [USACO1.4]等差数列 Arithmetic Progressions
一个等差数列是一个能表示成 a, a+b, a+2b, \dots ,a+nb\space (n \in \mathbb N)a,a+b,a+2b,…,a+nb (n∈N) 的数列。

在这个问题中 aa 是一个非负的整数，bb 是正整数。
写一个程序来找出在双平方数集合：

{ x | x = p^2 + q^2 \wedge p,q \in \mathbb N \cap [0,m]}
{x∣x=p
2
+q
2
∧p,q∈N∩[0,m]}

中长度为 nn 的等差数列。

输入格式
第一行一个正整数 nn，表示要找的数列长度。
第二行一个非负整数 mm，表示 p,qp,q 的上界。

输出格式
如果没有找到数列，输出 NONE。

如果找到了，输出一行或多行，每行由二个整数组成：a,ba,b。

这些行应该以 bb 为第一关键字，aa 为第二关键字升序排序。

所求的等差数列将不会多于10,000个。

思路：筛选出平方数，枚举首项和公差

#include<iostream>
#include<algorithm>
using namespace std;
#define inf 0x3f3f3f3f
const int N=2005;
int n,m,p[2*250*250+10],vis[2*250*250+10],cnt=0,k,f;
struct node{
    int x,y;
    bool friend operator <(const node&v,const node&u){
        if(v.y!=u.y)return v.y<u.y;
        else return v.x<u.x;
    }
}res[100010];

int main(){
    ios::sync_with_stdio(0);
    cin>>n>>m;
    for(int i=1;i<=m;i++){
        for(int j=0;j<=m;j++){
            p[++cnt]=i*i+j*j;
            vis[p[cnt]]=1;
        }
    }
    sort(p+1,p+1+cnt);
    cnt=unique(p+1,p+1+cnt)-p-1;
    for(int i=0;i<=cnt;i++){
        for(int j=i+1;j<=cnt;j++){
            int d=p[j]-p[i],a=p[i],t=p[j],tot=2;
			
            while(vis[t+d]){
				
                tot++;
                t+=d;
                if(tot==n){
                    res[++k].x=a;
                    res[k].y=d;
                    f=1;
                    break;
                }

            }

        }
    }
    sort(res+1,res+1+k);
    if(f)for(int i=1;i<=k;i++)cout<<res[i].x<<' '<<res[i].y<<endl;
    else cout<<"NONE";

	return 0;
}//

洛谷P1072 [NOIP2009 提高组] Hankson 的趣味题
题目描述
Hanks 博士是 BT（Bio-Tech，生物技术) 领域的知名专家，他的儿子名叫 Hankson。现在，刚刚放学回家的 Hankson 正在思考一个有趣的问题。

今天在课堂上，老师讲解了如何求两个正整数 c_1c
1

和 c_2c
2

的最大公约数和最小公倍数。现在 Hankson 认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数a_0,a_1,b_0,b_1a
0

,a
1

,b
0

,b
1

，设某未知正整数 xx 满足：

1． xx 和 a_0a
0

的最大公约数是 a_1a
1

；

2． xx 和 b_0b
0

的最小公倍数是 b_1b
1

。

Hankson 的“逆问题”就是求出满足条件的正整数 xx。但稍加思索之后，他发现这样的 xx 并不唯一，甚至可能不存在。因此他转而开始考虑如何求解满足条件的 xx 的个数。请你帮助他编程求解这个问题。

输入格式
第一行为一个正整数 nn，表示有 nn 组输入数据。接下来的nn 行每行一组输入数据，为四个正整数 a_0,a_1,b_0,b_1a
0

,a
1

,b
0

,b
1

，每两个整数之间用一个空格隔开。输入数据保证 a_0a
0

能被 a_1a
1

整除，b_1b
1

能被 b_0b
0

整除。

输出格式
共 nn 行。每组输入数据的输出结果占一行，为一个整数。

对于每组数据：若不存在这样的 xx，请输出 00，若存在这样的 xx，请输出满足条件的 xx 的个数；
思路：求解如下，从1到sqrt(b1)开始枚举，找到b1的因数判断即可
在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f

int main(){
    ios::sync_with_stdio(0);
    int n,a0,a1,b1,b0;
    cin>>n;
    while(n--){
            int cnt=0;
        cin>>a0>>a1>>b0>>b1;
        for(int i=1;i*i<=b1;i++){
            if(b1%i==0){
            if(i%a1==0&&__gcd(i/a1,a0/a1)==1&&__gcd(b1/i,b1/b0)==1){
                cnt++;
                //cout<<i<<endl;
            }
            int y=b1/i;
            if(i==y)continue;
            if(y%a1==0&&__gcd(y/a1,a0/a1)==1&&__gcd(b1/y,b1/b0)==1)cnt++;
            }
        }
        cout<<cnt<<endl;


    }


	return 0;
}

weixin_52109411

关注

1
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
数论初步学习

队伍需要，我开始学习数论错过暑假，懊悔万分，尽快弥补埃氏筛法与线性筛法int pri[100010]={2},vis[100010],cnt;for(int i=2;i<=100000;i++) { if(!vis[i]){ for(int j=i;j<=100000;j+=i) vis[j]=1; pri[++cnt]=i; }}//cnt==9592for(int i=2;i<100000;i++){ if(!check[i
复制链接

扫一扫