数论

最新推荐文章于 2022-02-08 23:41:20 发布
最新推荐文章于 2022-02-08 23:41:20 发布
阅读量178 收藏 1
点赞数 1
分类专栏： 算法基础
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_39953680/article/details/82113606
版权

一、欧几里得算法

int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}

若gcd (a, b)=m，则gcd (a/m, b/m)=1

最小公倍数 LCM(a,b)=a/gcd(a,b)*b

题目：POJ2773---Happy 2006

由题目可知k的范围明显比m大，想到往周期性靠拢

由gcd性质可得gcd(a, b)=gcd(b*t+a, b)，题目求的是与m互质的数即gcd(i,m)=1=gcd(m*t+i,m)，具有周期性，周期为m，注意k%num==0的情况

#include<cstdio>
using namespace std;
const int n_max=1e6+5;
int f[n_max];
int gcd(int a,int b)
{
    if(b==0) return a;
    return gcd(b,a%b);
}
int main()
{
    int m,k;
    while(~scanf("%d%d",&m,&k))
    {
        int num=0;
        for(int i=1; i<=m; ++i)
            if(gcd(i,m)==1)
                f[num++]=i;
        printf("%d\n",k%num==0?(k/num-1)*m+f[num-1]:k/num*m+f[k%num-1]);
    }
    return 0;
}

二、

矩阵快速幂运算

struct matr
{
    int m[n_max][n_max];
};
matr x,y;
const int mod=1e9+7;
for(int i=0; i<n; ++i)
    y.m[i][i]=1;
matr operator*(matr a,matr b)
{
    matr c;
    for(int i=0; i<n; ++i)
        for(int j=0; j<n; ++j)
        {
            c.m[i][j]=0;
            for(int k=0; k<n; ++k)
                c.m[i][j]=(c.m[i][j]+a.[i][k]*b.m[k][j])%mod;
        }
}
matr qpow(matr a,matr b,int k)
{
    while(k)
    {
        if(k&1) b=b*a;
        a=a*a;
        k>>1;
    }
    return b;
}
qpow(x,y,k);

快速积运算

int qmul(int a,int b)
{
    int ans=0;
    while(b)
    {
        if(b&1)
            ans+=a;
        a<<=1;
        b>>=1;
    }
    return ans;
}

快速幂运算

int qpow(int a,int b)
{
    int ans=1;
    while(b)
    {
        if(b&1)
            ans*=a;
        a*=a;
        b>>=1;
    }
    return ans;
}

快速积取模

int Modmul(int a,int b,int mod)
{
    int ans=0;
    while(b)
    {
        if(b&1)
            ans=(ans+a)%mod;
        a<<=1;
        a%=mod;
        b>>=1;
    }
    return ans;
}

快速幂取模

int Modpow(int a,int b,int mod)
{
    int ans=1;
    while(b)
    {
        if(b&1)
            ans=Modmul(ans,a,mod);
        a=Modmul(a,a,mod);
        b>>=1;
    }
    return ans;
}

素数筛

bool vis[n_max];
int prime[n_max],cnt;
void sieve()
{
    memset(vis,0,sizeof(vis));
    int m=sqrt(n_max);
    for(int i=2;i<=m;++i)
    {
        if(!vis[i])
        prime[cnt++]=i;
        for(int j=i;j<n_max;j+=i)
            vis[j]=true;
    }
}

三、质因数分解

题目1：POJ2992---Divisors

题目2：ACM-ICPC 2018 南京赛区网络预赛J---Sum

定义f(N)为N的分解组合数

①如果N本身是一个质数，它只能分解为1*N和N*1，f(N)=2

②如果N分解的两个乘数质数p和N/p=x，如果x可以整除p*p，很明显f(N)=0(N无论如何分解，其中一个乘数必定含有p的次方大于等于2)；如果x不能整除p，那么对于x的每一个分解，p都可以乘以其中一个乘数而不产生重复，所以f(N)=f(x)*2；如果x仅能整除p，x分解出来的每一个组合必定有其中一个乘数含有p，但f(N)!=f(x)，因为此时产生了重复(例如6*3)，但可以换一种思维来想，如果再将x的p因子剔除得到y，那么对于y的每一种分解，p都必须同时与两个乘数相乘(将两个p分开)，f(N)=f(y)=f(x/p)

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstring>
#include<string>
#include<cctype>
#include<cmath>
#include<map>
#include<vector>
#include<utility>
#define For(i,x,y) for(int i=x;i<=y;++i)
#define Fov(i,x,y) for(int i=x;i>=y;--i)
#define midf(a,b) ((a)+(b)>>1)
#define Num1(_) (_)<<1
#define Num2(_) ((_)<<1)|1
#define lowbit(x) x&-x
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
inline int read()
{
    char ch=getchar(); int x=0, f=1;
    while(ch<'0'||ch>'9') { if(ch=='-') f=-1; ch=getchar();}
    while('0'<=ch&&ch<='9') { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
const int inf=0x3f3f3f3f;
const double pi=acos(-1.0);
const int n_max=2e7+5;
ll f[n_max];
bool vis[n_max];
int prime[n_max];

int main()
{
	int cnt=0;
	f[1]=1;
	for(int i=2;i<n_max;++i)
	{
	    if(!vis[i])
        {
            prime[cnt++]=i;
            f[i]=2;
        }
        for(int j=0;j<cnt&&prime[j]*i<n_max;++j)
        {
            int x=i*prime[j];
            vis[x]=true;
            if(i%prime[j]==0)
            {
                if(i%(prime[j]*prime[j])==0)
                    f[x]=0;
                else
                    f[x]=f[i/prime[j]];
                break;
            }
            else
                f[x]=2*f[i];
        }
	}
	for(int i=2;i<n_max;++i)
        f[i]+=f[i-1];
	int T;
	scanf("%d",&T);
	while(T--)
    {
        int n;
        scanf("%d",&n);
        printf("%lld\n",f[n]);
    }
	return 0;
}

四、欧拉函数

欧拉函数φ(x):对于整数n，小于n的正整数中与n互质的数的数目

①求单个欧拉函数

int euler(int n)
{
    int ans=n;
    for(int i=2; i*i<=n; ++i)
        if(n%i==0)
        {
            ans=ans/i*(i-1);
            while(n%i==0) n/=i;
        }
    if(n>1) return ans=ans/n*(n-1);
    return ans;
}

②筛法求欧拉函数表

int euler[n_max];
for(int i=2; i<n_max; ++i)
    euler[i]=i;
for(int i=2; i<n_max; ++i)
    if(euler[i]==i)
        for(int j=i; j<n_max; j+=i)
            euler[j]=euler[j]/i*(i-1);

题目：POJ2480---Longge's problem

五、费马大定理

题目：HDU6441---Find Integer

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstring>
#include<string>
#include<cctype>
#include<cmath>
#include<map>
#include<vector>
#define For(i,x,y) for(int i=x;i<=y;++i)
#define Fov(i,x,y) for(int i=x;i>=y;--i)
#define midf(a,b) ((a)+(b)>>1)
#define Num1(_) (_)<<1
#define Num2(_) ((_)<<1)|1
using namespace std;
typedef long long ll;
inline int read()
{
    char ch=getchar(); int x=0, f=1;
    while(ch<'0'||ch>'9') { if(ch=='-') f=-1; ch=getchar();}
    while('0'<=ch&&ch<='9') { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
const int inf=0x3f3f3f3f;

int main()
{
    int n,a,T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&a);
        if(n==0||n>2)
            printf("-1 -1\n");
        else if(n==1)
            printf("%d %d\n",1,a+1);
        else
        {
            if(a&1)
            {
                int x=(a-1)>>1,b=2*x*x+2*x;
                printf("%d %d\n",b,b+1);
            }
            else
            {
                int x=a>>1,b=x*x-1;
                printf("%d %d\n",b,b+2);
            }
        }
    }
	return 0;
}

六、费马小定理

如果p为质数，且a与p互质，则有 ，可推出

取模运算法则

①若a≡b (% p)，则对于任意的c，都有(a + c) ≡ (b + c) (%p)和(a * c) ≡ (b * c) (%p)；

②若a≡b (% p)，c≡d (% p)，则 (a + c) ≡ (b + d) (%p)，(a - c) ≡ (b - d) (%p)，(a * c) ≡ (b * d) (%p)，(a / c) ≡ (b / d) (%p)；

题目1：HDU6440---Dream

对于任意比质数p小的数a，b， a、b、a+b与p肯定互质，因而根据费马小定理和取模运算法则有 ，因而可以定义a+b为(a+b)%p，a*b为(a*b)%p

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstring>
#include<string>
#include<cctype>
#include<cmath>
#include<map>
#include<vector>
#define For(i,x,y) for(int i=x;i<=y;++i)
#define Fov(i,x,y) for(int i=x;i>=y;--i)
#define midf(a,b) ((a)+(b)>>1)
#define Num1(_) (_)<<1
#define Num2(_) ((_)<<1)|1
using namespace std;
typedef long long ll;
inline int read()
{
    char ch=getchar();
    int x=0, f=1;
    while(ch<'0'||ch>'9')
    {
        if(ch=='-') f=-1;
        ch=getchar();
    }
    while('0'<=ch&&ch<='9')
    {
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
const int inf=0x3f3f3f3f;
const double pi=acos(-1.0);

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int p;
        scanf("%d",&p);
        For(i,0,p-1)
        For(j,0,p-1)
        printf("%d%c",(i+j)%p,j==p-1?'\n':' ');
        For(i,0,p-1)
        For(j,0,p-1)
        printf("%d%c",(i*j)%p,j==p-1?'\n':' ');
    }
    return 0;
}

欧拉定理

若正整数 a , n 互质，则  a^φ(n)≡1(mod n) ；推论：若正整数 a , n 互质，a^p≡a^(p mod φ(n)) (mod n)

题目2：ACM-ICPC 2018 焦作赛区网络预赛G---Give Candies

尝试几个数后得出规律ans=2^(n-1)；n的范围极大，题目很明显求的是超高次幂；根据欧拉定理的推论可实现降幂。

#include <bits/stdc++.h>
#define For(i,x,y) for(int i=(x);i<=(y);++i)
#define Fov(i,x,y) for(int i=(x);i>=(y);--i)
#define Fo(i,x,y) for(int i=(x);i<(y);++i)
#define midf(a,b) ((a)+(b)>>1)
#define L(_) (_)<<1
#define R(_) ((_)<<1)|1
#define fi first
#define se second
#define ss(_) scanf("%s",_)
#define si(_) scanf("%d",&_)
#define sii(x,y) scanf("%d%d",&x,&y)
#define siii(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define mem(x,y) memset(x,y,sizeof(x))
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
inline int read()
{
    char ch=getchar(); int x=0, f=1;
    while(ch<'0'||ch>'9') { if(ch=='-') f=-1; ch=getchar();}
    while('0'<=ch&&ch<='9') { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
const int inf=0x3f3f3f3f;
const double pi=acos(-1.0);
const int n_max=1e5+10;
const ll mod=1e9+7;
char s[n_max];
int main()
{
	//freopen("in.txt","r",stdin);
	int T; si(T);
	ll mod1=mod-1;
	while(T--)
    {
        ss(s);
        int len=strlen(s);
        ll n=0,ans=1,x=2;
        Fo(i,0,len) n=(n*10+s[i]-'0')%mod1;
        --n;
        while(n)
        {
            if(n&1) ans=ans*x%mod;
            x=x*x%mod;
            n>>=1;
        }
        printf("%lld\n",ans);
    }
	return 0;
}

七、莫比乌斯反演

整除分块：求

for(int l=1,r;l<=n;l=r+1)
{
    r=n/(n/l);
    ans+=(r-l+1)*(n/l);
}

 

King_Tom
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
数论专题(1)数论函数,整数分块
cyyyyds857的博客
12-26 1323
定义(数论函数)从今天起,我们将要开始数论的学习,是不是感觉很难?难的话就听我讲吧,讲了后就不难了(bushi)
初等数论
CYCKCN的博客
04-02 2889
初等数论 大佬们一个晚上讲完了n个知识点，蒟蒻只能事后默默的记下来了…… 一点点关于超级初等的数论的听课笔记，可能会有不充足或者顺序不对的地方，请见谅（毕竟我是蒟蒻嘛）…… 前言 １、求出所有因子 //求x的所有因子 for (int i = 1; i * i <= x; ++i) if (x % i == 0) { ans[++cnt] = i;
我要成为数论king！！！
Artist's Garden
12-07 306
唯一分解定理 牛牛的幂运算 /* 设a = p^k1 c = p^k2 （0）因为1的幂次永远是1,会使得代码错误,要特判 （1）k1==k2，此时a == c,那么必然有b == d。答案等于n*n （2）k1<k2, (2.1)p<=根号n 爆搜1~根号n 爆搜k1与k2 （2.1.1）k1与k2不互素,即gcd(k1,k2)!=1,那么在他们等于gcd(k1,k2)的时候 答案已经被贡献。 （2.1.1.1）答案怎么贡献？此时我们已经得到k1和k2了,就
数论-基础
king_fdy的博客
02-10 220
质数筛 #include <iostream> using namespace std; const int N = 1e6+10; int prime[N], cnt; bool st[N]; //朴素筛法-O(nlogn) void get_primes(int x) { for(int i = 2; i <= n; i++) { if(!st[i...
数学知识--错排公式
weixin_30634661的博客
04-03 261
错排公式 核心递推公式： D(n) = (n-1) [D(n-2) + D(n-1)] 特殊地，D(1) = 0, D(2) = 1. 问题: 十本不同的书放在书架上。现重新摆放，使每本书都不在原来放的位置。有几种摆法? 这个问题推广一下，就是错排问题，是组合数学中的问题之一。考虑一个有n个元素的排列，若一个排列中所有的元素都不在自己原来的位置上，那么这样...
数论约数
weixin_44203780的博客
03-27 192
数论，公约数 落谷p1414 根据题意，统计每个因数出现的个数，然后从大到小依次找出就好了； #include<bits/stdc++.h> using namespace std; #define fin(a,n) for(int i=a;i<=n;i++) #define ll long long int const int inf=1e6+10; const int max...
数论入门知识
热门推荐
徒手拆机甲的博客
08-01 1万+
数论入门 提到数论，可能很多人都感到很头疼，甚至很多时候遇到一些问题，看到成篇的证明都会感到恐惧，而且由于关于ACM方面的数论资料，网上资料都比较驳杂。有时候很容易出现知其然不知其所以然的情况。所以今天给大家介绍一些关于数论入门最基础的知识和算法，内容会尽量从0基础开始，所以内容会尽量详细。不过其中部分证明还是需要一些高中数学基础的。由于内容很多，我整理了一个目录，按顺序讲今天的内容. 目录 一.合数，质数，整除，互质，同余，取模等基础概念。 二.欧几里得算法 三.扩展欧几里得 四.费马小定理 五.欧拉函数
数论函数
Greninja的博客
12-10 1394
数论函数
数论学习-初等数论基础总览
麒麟001号的博客
02-08 2980
文章目录初等数论基础（二）建议先看一、数论只会gcd1.1 gcd(a,b) = (a,a+b) 的证明(a,b) = (b,a%b)的证明辗转相除法代码：辗转相除法时间复杂度1.2 exgcd1.2.1 裴蜀定理1.2.2 ax+by = gcd(a,b)的求解二、欧拉相关2.1 欧拉函数2.2 欧拉定理相关2.2.1 欧拉定理2.2.1 扩展欧拉定理2.2.2 光速幂 初等数论基础（二） 建议先看 初等数论基础（一）~ 一、数论只会gcd 1.1 gcd (a,b) = (a,a+b) 的证明 (a,
《初等数论》第七章习题答案.rar
01-19
初等数论是一门研究整数性质的数学分支，它主要关注整数的性质和问题，包括素数、同余类、除法定理等。在《初等数论》（第四版）一书中，作者闵嗣鹤和严士健深入浅出地介绍了这一领域的基本概念和定理。第七章聚焦...
《初等数论》第六章习题答案.rar
01-19
初等数论是数学的一个重要分支，主要研究整数的性质，特别是关于除法和同余类的问题。《初等数论》（第四版）由闵嗣鹤和严士健两位数学家编著，是一本深入浅出的经典教材。本章节我们将探讨第六章的核心概念——原根...
初等数论100例.pdf
03-10
初等数论100例
《初等数论》第九章习题答案.rar
01-19
初等数论是数学的一个重要分支，主要研究整数的性质，尤其是质数的分布和其他基本概念。在《初等数论》（第四版）这本经典教材中，第九章深入探讨了数论函数和质数分布这两个核心主题。本章节的习题答案涵盖了9.1可...
初等数论习题解答第三版
最新发布
10-14
"初等数论习题解答第三版" 初等数论是数学的基础课之一，对于学习数学的学生来说极为重要。初等数论习题解答第三版为广东石油化工学院编写的习题解答，涵盖了初等数论的主要内容。本资源摘要信息将对初等数论习题...
关于数论分块的证明
01-06
2.训练指南数论公因数部分 当 i∈[1,n]i\in[1,n]i∈[1,n] 时， ⌊ni⌋=⌊n⌊n⌊ni⌋⌋⌋\lfloor\frac{n}{i}\rfloor=\lfloor\frac{n}{\lfloor\frac{n}{\lfloor\frac{n}{i}\rfloor}\rfloor}\rfloor⌊in​⌋=⌊⌊⌊in​...
lecture8_数论.pdf
07-15
数论讲义， 数论是纯粹数学的分支之一，主要研究整数的性质。整数可以是方程式的解（丢番图方程）。有些解析函数（像黎曼ζ函数）中包括了一些整数、质数的性质，透过这些函数也可以了解一些数论的问题。透过数论也...
《初等数论》第二章习题答案.rar
12-16
初等数论是数学的一个重要分支，主要研究整数性质及其相互关系，是数学理论的基础之一。本压缩包文件“《初等数论》第二章习题答案.rar”提供了《初等数论》（第四版）一书第二章中关于不定方程的习题解答，由著名...
数论（常见数论知识点总结）.docx
07-20
"此文档是关于数论的常见知识点总结，涵盖了整除的性质、数论定理、模运算与余数、素数、莫比乌斯函数、逆序数、原根以及离散对数等多个核心概念。" 一. 整除的性质 整除性质是数论的基础，它们描述了整数之间除法...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值