数论总结

一、知识梳理。

1、整除：把握好性质，另外特别注意几种常用的整除的情况，例如能11整除的数的特征等等。

https://blog.csdn.net/qq_41661919/article/details/81291452

2、同余：同余这个问题感觉还是不好理解的。。。注意，同余不满足同除性，一般要用逆元替换（快速幂）。

3、求最大公约数：

①辗转相除法求最大公约数就不说了：     %;

②进一步可以用二进制优化，不断去2，实现代码的话可以作为模板使用。

int GCD(int x,int y)
{
    int i,j;
    if(x==0)return y;
    if(y==0)return x;
    for(i=0;0==(x&1);++i)x>>=1;
    for(j=0;0==(y&1);++j)y>>=1;
    if(j<i)i=j;
    while(1)
    {
        if(x<y)x^=y,y^=x,x^=y;
        if(0==(x-=y))return y<<i;
        while(0==(x&1))x>>=1;
    }
}

③求最小公倍数：直接套公式。。。可以先除后乘。

④扩展欧几里得解方程：，主要利用递归边界中想，x，y的确定值在退栈时求出最终的x，y。

模板：

#include<bits/stdc++.h>
using namespace std;
int exgcd(int a,int b,int &x,int &y)
{
    int ret,tmp;
    if(!b)//注意递归边界
    {
        x=1;
        y=0;
        return a;
    }
    ret=exgcd(b,a%b,x,y);
    tmp=x;
    x=y;
    y=tmp-a/b*y;//对应关系自己推导一下就可以了
    return ret;
}
int main()
{
    
}

 

4、逆元：

求逆元的方法也不少，不过学长说最常用的还是费马小定理+快速幂，注意p一定要为质数。

其次，逆元的应用在于同余运算中没有除法运算，所以要用逆元开替换分母。

快速幂取模 模板（非递归实现）：

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

 

5、中国剩余定理：求解同余方程组。

原理以及模板：https://blog.csdn.net/qq_41661919/article/details/81322682

6、斐波那切数列：

求的方法：

①递推公式 ；

②待定系数法求的通项公式。（n较大时可以考虑）。结论也可以考虑记一下。

7、卡特兰数：

找规律，往公式上靠。。。

8、素数：

主要内容是怎样判断素数，以及素数的相关定理（了解）。

①穷举  ②线性筛 ③快速线性筛

重点是③的模板：

void is_prime()
{
    for(int i=2;i<n;++i)
    {
        if(flag[i]==0)prime[num++]=i;
        for(int j=0;j<num&&i*prime[j]<n;++j)
        {
            flag[i*prime[j]]=1;
            if(i%prime[j]==0)break;//保证无重复删除
        }
    }
}

④：M-R素数测试：模板。。。

9、欧拉函数的线性筛方法：求欧拉函数的值。

板子：

#include<bits/stdc++.h>
using namespace std;
int phi[10010],prime[10010],tot,ans;
bool flag[10010];
int n;
void getphi()
{
    int i,j;
    phi[1]=1;//**
    for(i=2;i<=n;++i)
    {
        if(flag[i]==0)
        {
            prime[++tot]=i;
            phi[i]=i-1;//**
        }
        for(j=1;j<=tot;++j)
        {
            if(i*prime[j]>n)break;
            flag[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];//**
                break;
            }
            else
            {
                phi[i*prime[j]]=phi[i]*(prime[j]-1);//**

            }
        }
    }
}
int main()
{

}