数论学习

最新推荐文章于 2022-06-07 16:18:11 发布
最新推荐文章于 2022-06-07 16:18:11 发布
阅读量711 收藏
点赞数 2
分类专栏： 数论 文章标签： c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Clarehehe/article/details/50539642
版权

1⃣️扩展欧几里得算法

sgu 141 扩展欧几里德＋枚举大法，不过sgu好像崩了。。。。


poj 2891

合并模方程式

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 100010
typedef long long ll;

ll K;
ll a[N],r[N];

inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void Exgcd(ll a,ll b,ll &x,ll &y,ll &d)
{
    if(b==0)
    {
        d=a;x=1;y=0;
        return;
    }
    Exgcd(b,a%b,y,x,d);
    y-=a/b*x;
}

ll work()
{
    ll Last_a=a[1],Last_r=r[1],x,y,d;
    for(int i=2;i<=K;i++)
    {
        Exgcd(Last_a,a[i],x,y,d);
        if((r[i]-Last_r)%d!=0)
            return -1;
        ll tmp=a[i]/d;
        x=((r[i]-Last_r)/d*x)%tmp;
        Last_r+=x*Last_a;
        Last_a=Last_a/d*a[i];
        Last_r%=Last_a;
    }
    return Last_r>0?Last_r:Last_r+Last_a;
}

int main()
{
    while(scanf("%lld",&K)!=EOF)
    {
        for(ll i=1;i<=K;i++)
        {
            a[i]=read();r[i]=read();
        }
        printf("%lld\n",work());
    }
    return 0;
}


2⃣️Miller-Rabbin 素数检测法

1.费马小定理

2.二次探测定理

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

typedef long long ll;
const int S=3;

int K,Ans;
ll A[4]={0,2,7,61};     //黑科技一波

inline ll read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

ll Quick_mul(ll a,ll b,ll n)
{
    ll base=a,r=0;
    while(b)
    {
        if(b&1)
            r=(r+base)%n;
        base=(base+base)%n;
        b>>=1;
    }
    return r;
}

ll Quick_pow(ll a,ll exp,ll n)
{
    ll base=a,r=1;
    while(exp)
    {
        if(exp&1)
            r=r*base%n;
        base=base*base%n;
        exp>>=1;
    }
    return r;
}

bool Miller_Rabbin(ll n)
{
    if(n==2)return true;
    int t=0;
    ll a,x,y,u=n-1;
    while(!(u&1))
    {
        t++;u>>=1;
    }
    for(int i=1;i<=S;i++)
    {
        a=A[i];
        x=Quick_pow(a,u,n);
        for(int j=1;j<=t;j++)
        {
            y=Quick_mul(x,x,n);
            if(y==1&&x!=1&&x!=n-1)
                return false;
            x=y;
        }
        if(x!=1)
            return false;
    }
    return true;
}

int main()
{
    while(scanf("%d",&K)!=EOF)
    {
        Ans=0;
        for(int i=1;i<=K;i++)
        {
            ll n=read();
            if(Miller_Rabbin(n))
                Ans++;
        }
        printf("%d\n",Ans);
    }
    return 0;
}


3⃣️.快速求逆元

4⃣️卢卡斯定理  －－－组合数取模

Lucas(n,m,p)=c(n%p,m%p)*Lucas(n/p,m/p,p)


1.hdu 3037 Saving Beans

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 100010
typedef long long ll;

int T;
ll n,m,p;
ll pre[N],S_pre[N];

ll Quick_mul(ll a,ll b)
{
    ll base=a,r=0;
    while(b)
    {
        if(b&1)
            r=(r+base)%p;
        base=(base+base)%p;
        b>>=1;
    }
    return r;
}

ll Quick_pow(ll a,ll exp)
{
    ll base=a,r=1;
    while(exp)
    {
        if(exp&1)
            r=Quick_mul(r,base);
        base=Quick_mul(base,base);
        exp>>=1;
    }
    return r;
}

void Prepare()
{
    pre[0]=1;
    for(int i=1;i<=p;i++)
        pre[i]=pre[i-1]*i%p;
    S_pre[p-1]=Quick_pow(pre[p-1],p-2);
    for(int i=(int)p-2;i>=1;i--)
        S_pre[i]=S_pre[i+1]*(i+1)%p;
    S_pre[0]=1;
}

int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lld%lld%lld",&n,&m,&p);
        n+=m;
        Prepare();
        ll Ans=1;
        while(n||m)
        {
            ll t1=n%p,t2=m%p;
            if(t1<t2)
            {
                Ans=0;break;
            }
            ll tmp=((pre[t1]*S_pre[t1-t2])%p*S_pre[t2])%p;
            Ans=Ans*tmp%p;
            n/=p;m/=p;
        }
        printf("%lld\n",Ans);
    }
    return 0;
}


5⃣️.中国剩余定理

网上看了好多都是在乱扯。。。强行把欧几里德题解说成中国剩余定理 如上poj 2891


1.hdu 1370  Biorhythms

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 10
typedef long long ll;

ll A[N],B[N],C[N];

inline ll read()
{
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void Exgcd(ll a,ll b,ll &x,ll &y,ll &d)
{
    if(b==0)
    {
        d=a;x=1;y=0;
        return;
    }
    Exgcd(b,a%b,y,x,d);
    y-=a/b*x;
}

int main()
{
    ll Zy=read();
    ll Lim_min,x,y,d;
    ll Mul=23*28*33;
    Exgcd(Mul/23,23,x,y,d);B[1]=x;C[1]=23;
    Exgcd(Mul/28,28,x,y,d);B[2]=x;C[2]=28;
    Exgcd(Mul/33,33,x,y,d);B[3]=x;C[3]=33;
    int Case=0;
    while(scanf("%lld%lld%lld%lld",&A[1],&A[2],&A[3],&Lim_min))
    {
        if(A[1]==-1&&A[2]==-1&&A[3]==-1&&Lim_min==-1)
            break;
        ll Ans=0;
        for(int i=1;i<=3;i++)
        {
            Ans=(Ans+Mul/C[i]*B[i]*A[i])%Mul;
        }
        if(Ans<=Lim_min)
        {
            ll qwer=(Lim_min-Ans)/Mul+1;
            Ans+=qwer*Mul;
        }
        printf("Case %d: the next triple peak occurs in %lld days.\n",++Case,Ans-Lim_min);
    }
    return 0;
}


6⃣️欧拉函数与莫比乌斯函数


1.hdu 2824 （线性筛求欧拉函数）

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 3000010

int a,b,tot;
int phi[N],Prime[N/3];
bool flag[N];

void Prepare()
{
    for(int i=2;i<=3000000;i++)
    {
        if(!flag[i])
        {
            Prime[++tot]=i;
            phi[i]=i-1;
        }
        for(int j=1;(long long)Prime[j]*i<=3000000;j++)
        {
            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()
{
    Prepare();
    while(scanf("%d%d",&a,&b)!=EOF)
    {
        long long Ans=0;
        for(int i=a;i<=b;i++)
            Ans+=phi[i];
        printf("%lld\n",Ans);
    }
    return 0;
}


2.求 ∑gcd(i,j)      1≤i≤n ,   1≤j≤m

Orz KZF


其过程就是不断替换求和公式的变量得到新的求和公式（瞎bb的）

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 100010

int n,m,tot;
int Prime[N],phi[N];
bool flag[N];

void Prepare()
{
    phi[1]=1;
    for(int i=2;i<=100000;i++)
    {
        if(!flag[i])
        {
            Prime[++tot]=i;
            phi[i]=i-1;
        }
        for(int j=1;(long long)Prime[j]*i<=100000;j++)
        {
            flag[i*Prime[j]]=true;
            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()
{
    Prepare();
    int a,b;scanf("%d%d",&a,&b);
    n=min(a,b);m=max(a,b);
    long long Ans=0;
    for(int i=1;i<=n;i++)
    {
        Ans+=(long long)(n/i)*(m/i)*phi[i];
    }
    cout<<Ans<<endl;
    return 0;
}

那么如果n是10^9数量级的呢？

让我们继续orz


复杂度s

是O(n^0.67)


3.poj 2480 （被网上坑了。。。看了好久。。暴力不就行了么，推个毛公式啊）

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 1000010

int n,tot;
int Prime[N];
bool flag[N];

void Prepare()
{
    for(int i=2;i<=1000000;i++)
    {
        if(!flag[i])
        {
            Prime[++tot]=i;
        }
        for(int j=1;(long long)Prime[j]*i<=1000000;j++)
        {
            flag[i*Prime[j]]=true;
            if(i%Prime[j]==0)
                break;
        }
    }
}

int Get_phi(int k)
{
    if(k==1)return 1;
    int num=k;
    for(int i=1;i<=tot;i++)
    {
        if(k%Prime[i]==0)
        {
            num=num/Prime[i]*(Prime[i]-1);
            while(k%Prime[i]==0)
                k/=Prime[i];
            if(k==1)
                break;
        }
    }
    if(k!=1)
        num=num/k*(k-1);
    return num;
}

int main()
{
    Prepare();
    while(scanf("%d",&n)!=EOF)
    {
        long long Ans=0;
        for(int i=1;i<=sqrt((double)n);i++)
        {
            if(n%i==0)
            {
                Ans+=(long long)i*Get_phi(n/i);
                if(i!=(n/i))
                {
                    int t=n/i;
                    Ans+=(long long)t*Get_phi(i);
                }
            }
        }
        printf("%lld\n",Ans);
    }
    return 0;
}


4.NYOJ 1066 CO—PRIME  

莫比乌斯反演

此题中，设F(d)表示n个数中gcd为d的倍数的数有多少对，f(d)表示n个数中gcd恰好为d的数有多少对，

则F(d)=∑f(n) (n % d == 0)

f(d)=∑mu[n / d] * F(n) (n %d == 0)   

--------->  f(1)=∑mu[n] * F(n)   

                                         .......（摘抄From   http://blog.csdn.net/lyhvoyage/article/details/38455415）




#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 100010

int n,tot,Max;
int A[N],Prime[N],mu[N],num[N],cnt[N];
bool flag[N];

inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void Prepare()
{
    mu[1]=1;
    for(int i=2;i<=100000;i++)
    {
        if(!flag[i])
        {
            Prime[++tot]=i;
            mu[i]=-1;
        }
        for(int j=1;(long long)i*Prime[j]<=100000;j++)
        {
            flag[i*Prime[j]]=true;
            if(i%Prime[j]==0)
            {
                mu[i*Prime[j]]=0;
                break;
            }
            else
                mu[i*Prime[j]]=mu[i]*(-1);
        }
    }
}

int main()
{
    Prepare();
    while(scanf("%d",&n)!=EOF)
    {
        Max=0;
        memset(num,0,sizeof(num));
        memset(cnt,0,sizeof(cnt));
        for(int i=1;i<=n;i++)
        {
            A[i]=read();
            Max=max(Max,A[i]);num[A[i]]++;
        }
        for(int i=1;i<=Max;i++)
        {
            for(int j=i;j<=Max;j+=i)
                cnt[i]+=num[j];
        }
        long long Ans=0;
        for(int i=1;i<=Max;i++)
        {
            Ans+=(long long)mu[i]*cnt[i]*(cnt[i]-1)/2;
        }
        printf("%lld\n",Ans);
    }
    return 0;
}


5.bzoj 2820 YY的GCD

莫比乌斯函数

orz Hzwer.......


在维护一个前缀和做。。。。。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 10000010

int n,m,tot,Max;
int A[N],Prime[N],mu[N],num[N],cnt[N],F[N];
long long sum[N];
bool flag[N];

inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void Prepare()
{
    mu[1]=1;
    for(int i=2;i<=10000000;i++)
    {
        if(!flag[i])
        {
            Prime[++tot]=i;
            mu[i]=-1;
        }
        for(int j=1;(long long)i*Prime[j]<=10000000;j++)
        {
            flag[i*Prime[j]]=true;
            if(i%Prime[j]==0)
            {
                mu[i*Prime[j]]=0;
                break;
            }
            else
                mu[i*Prime[j]]=mu[i]*(-1);
        }
    }
    for(int i=1;i<=tot;i++)
    {
        for(int j=Prime[i];j<=10000000;j+=Prime[i])
            F[j]+=mu[j/Prime[i]];
    }
    for(int i=1;i<=10000000;i++)
        sum[i]=sum[i-1]+F[i];
}

int main()
{
    Prepare();
    int T=read();
    while(T--)
    {
        n=read();m=read();
        if(n>m)swap(n,m);
        long long Ans=0;
        for(int i=1,j;i<=n;i=j+1)     //这段666
        {
            j=min(n/(n/i),m/(m/i));
            Ans+=(sum[j]-sum[i-1])*(n/i)*(m/i);
        }
        printf("%lld\n",Ans);
    }
    return 0;
}

6.bzoj 2154 Crash的数字表格 


#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 10000010
#define Mod 20101009

int n,m,tot;
int Prime[N/3],F[N],mu[N];
int sum[N];
bool flag[N];
long long Ans;

void Prepare(int n)
{
    mu[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!flag[i])
        {
            Prime[++tot]=i;
            mu[i]=-1;F[i]=-i;
        }
        for(int j=1;(long long)i*Prime[j]<=n;j++)
        {
            flag[i*Prime[j]]=true;
            if(i%Prime[j]==0)
            {
                mu[i*Prime[j]]=0;
                break;
            }
            mu[i*Prime[j]]=mu[i]*(-1);
        }
    }
    for(int i=1;i<=n;i++)
        sum[i]=(sum[i-1]+((long long)i*i%Mod*mu[i]))%Mod;
}

long long Cal(int x,int y)
{
    long long t1=(long long)(x+1)*x/2%Mod,t2=(long long)(y+1)*y/2;
    return t1*(t2%Mod)%Mod;
}

int Get(int x,int y)
{
    long long Sum=0;
    for(int i=1,j;i<=x;i=j+1)
    {
        j=min(x/(x/i),y/(y/i));
        long long qwer=Cal(x/i,y/i)*(sum[j]-sum[i-1])%Mod;
        Sum=(Sum+qwer)%Mod;
    }
    return (int)Sum;
}

int main()
{
    cin>>n>>m;
    if(n>m)swap(n,m);
    Prepare(m);
    for(int i=1,j;i<=n;i=j+1)
    {
        j=min(n/(n/i),m/(m/i));
        long long qwer=(long long)(i+j)*(j-i+1)/2%Mod*Get(n/i,m/i)%Mod;
        Ans=(Ans+qwer)%Mod;
    }
    Ans=(Ans+Mod)%Mod;
    cout<<Ans<<endl;
    return 0;
}


7⃣️扩展大步小步算法

bzoj 2480 Mod

真是醉了，写了个快速乘比普通的还慢。。tle了找了好久

/**************************************************************
    Problem: 2480
    User: Edward2173
    Language: C++
    Result: Accepted
    Time:6368 ms
    Memory:2604 kb
****************************************************************/
 
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <map>
#include <queue>
#include <vector>
using namespace std;
 
#define N 100010
#define INF 0x3f3f3f3f
typedef long long ll;
 
int a,b,p;
int len[N];
 
int GCD(int a,int b)
{
    return !b?a:GCD(b,a%b);
}
 
ll Quick_mul(ll a,ll b)
{
    ll base=a,r=0;
    while(b)
    {
        if(b&1)
            r=(r+base)%p;
        base=(base+base)%p;
        b>>=1;
    }
    return r;
}
 
ll Quick_pow(int a,int exp)
{
    ll base=a,r=1;
    while(exp)
    {
        if(exp&1)
            r=r*base%p;
        base=base*base%p;
        exp>>=1;
    }
    return r;
}
 
int work(int a,int b)
{
    a%=p;b%=p;
    if(!b)return 0;
    ll Add=0,tmp=1;
    while(1)
    {
        int qwer=GCD(a,p);
        if(b%qwer!=0)
            return INF;
        if(qwer==1)
            break;
        tmp=tmp*a/qwer%p;
        b/=qwer;p/=qwer;
        Add++;
        if(tmp==b)
            return (int)Add;
    }
    map<ll,int> Hash;
    int Unit=sqrt((double)p+1.0),Ans=INF;
    ll now=b;Hash[now]=0;
    for(int i=1;i<Unit;i++)
    {
        now=now*a%p;
        Hash[now]=i;
    }
    ll base=Quick_pow(a,Unit),S_now=tmp;
    for(int i=1;i<=Unit+1;i++)
    {
        S_now=S_now*base%p;
        if(Hash.count(S_now))
            return i*Unit-Hash[S_now]+(int)Add;
    }
    return Ans;
}
 
int main()
{
    while(scanf("%d%d%d",&a,&p,&b))
    {
        if(a==0&&b==0&&p==0)
            break;
        int qwer=work(a,b);
        if(qwer==INF)
            printf("No Solution\n");
        else printf("%d\n",qwer);
    }
    return 0;
}


8⃣️斐波那契数列

codeforces 446 C   推下公式就行了

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

#define N 300010
#define Mod 1000000009LL
typedef long long ll;

int n,m;
ll f[N],F[N],G[N],A[N];
ll a[N*4],b[N*4],sum[N*4];

inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void Prepare()
{
    G[1]=f[1]=F[0]=F[1]=1;
    for(int i=2;i<=300000;i++)
    {
        f[i]=(f[i-1]+f[i-2])%Mod;
        F[i]=(F[i-1]+f[i-1])%Mod;
        G[i]=(G[i-1]+f[i])%Mod;
    }
}

void Pushup(int i)
{
    sum[i]=(sum[i*2]+sum[i*2+1])%Mod;
}

void Change(int i,int len,ll v1,ll v2)
{
    a[i]=(a[i]+v1)%Mod;
    b[i]=(b[i]+v2)%Mod;
    sum[i]=(sum[i]+F[len-1]*v1%Mod+G[len-1]*v2%Mod)%Mod;
}

void Pushdown(int i,int l,int r)
{
    if(a[i]==0&&b[i]==0)
        return;
    int mid=(l+r)/2,len=mid-l+1;
    Change(i*2,len,a[i],b[i]);
    ll v1=(f[len-1]*a[i]+f[len]*b[i])%Mod,v2=(f[len]*a[i]+f[len+1]*b[i])%Mod;
    Change(i*2+1,r-mid,v1,v2);
    a[i]=b[i]=0;
}

void Build(int i,int l,int r)
{
    if(l==r)
    {
        a[i]=b[i]=0;sum[i]=A[l];
        return;
    }
    int mid=(l+r)/2;
    Build(i*2,l,mid);
    Build(i*2+1,mid+1,r);
    Pushup(i);
}

void Update(int i,int l,int r,int left,int right,ll v1,ll v2)
{
    if(l>=left&&r<=right)
    {
        Change(i,r-l+1,v1,v2);
        return;
    }
    Pushdown(i,l,r);
    int mid=(l+r)/2;
    if(left<=mid)
        Update(i*2,l,mid,left,right,v1,v2);
    if(right>mid)
    {
        ll t1=v1,t2=v2;
        if(left<=mid)
        {
            int len=mid-max(left,l)+1;
            t1=(f[len-1]*v1+f[len]*v2)%Mod;
            t2=(f[len]*v1+f[len+1]*v2)%Mod;
        }
        Update(i*2+1,mid+1,r,left,right,t1,t2);
    }
    Pushup(i);
}

ll Query(int i,int l,int r,int left,int right)
{
    if(l>=left&&r<=right)
    {
        return sum[i];
    }
    Pushdown(i,l,r);
    int mid=(l+r)/2;
    if(right<=mid)
        return Query(i*2,l,mid,left,right);
    else if(left>mid)
        return Query(i*2+1,mid+1,r,left,right);
    else
        return (Query(i*2,l,mid,left,mid)+Query(i*2+1,mid+1,r,mid+1,right))%Mod;
}

int main()
{
    Prepare();
    n=read();m=read();
    for(int i=1;i<=n;i++)
        A[i]=read();
    Build(1,1,n);
    for(int i=1;i<=m;i++)
    {
        int kind=read(),x=read(),y=read();
        if(kind==1)
            Update(1,1,n,x,y,1,1);
        else
            printf("%lld\n",Query(1,1,n,x,y));
    }
    return 0;
}




暂时更到这儿了，以后再继续。。。    2016.1.23

Mitsuha三叶
关注
  • 2
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
数论学习系列之取整函数.pdf
07-15
数论相关知识讲解，函数y=[x]称为取整函数，也称高斯函数。其中不超过实数x的最大整数称为x的整数部分，记作[x]或INT(x)。该函数被广泛应用于数论，函数绘图和计算机领域。
数论部分学习笔记-by BeiYu
02-21
声明:该篇博文大多摘抄自网络，zky学长的ppt，及神犇的博客(%%%)，个人整理，有所不好请见谅！
数论学习路线
m0_53373033的博客
09-25 372
数论学习大纲 一. 数论学习笔记1 整除 同余 最大公约数 扩展欧几里得 逆元 二. 数论学习笔记2 素数 中国剩余定理 扩展中国剩余定理 欧拉定理 三. 数论学习笔记3 扩展欧拉定理 大步小步算法（施工中） 扩展大步小步算法（施工中） 四. 数论学习笔记4 积性函数和各种筛法（施工中） 莫比乌斯反演（施工中） ...
学习数论的一些基础知识
To be a better man
06-07 567
自然数 自然数概念指用以计量事物的件数或表示事物件数的数，即用数码0，1，2，3，4，……所表示的数 。自然数由0开始 ， 一个接一个，组成一个无穷集体。 互质 互质是公约数只有1的两个整数，叫做互质整数。公约数只有1的两个自然数，叫做互质自然数，后者是前者的特殊情形。 a|b在数学中代表什么 “a|b是整除的意思。 若整数b除以非零整数a,商为整数,且余数为零, 我们就说b能被a整除(或说a能整除b),b为被除数,a为除数,即a|b(“|”是整除符号),读作“a整除b”或“b能被a整除”。”..
数论基础
weixin_43184669的博客
08-18 402
1.唯一分解定理 2.容斥定理：等号右边，每一项前面的符合取决于每一项集合的元素个数，奇数个的集合为正，偶数个的为负。
详解数论从入门到入土
Dawn的博客
10-24 9433
一篇写得非常认真的数论入门讲义啊 希望可以帮到你 如果你认为这篇博客写的很好可以给个赞嘛
1、[高中竞赛] 高中联赛数论学习总结(2020.04.26).pdf
05-20
1、[高中竞赛] 高中联赛数论学习总结(2020.04.26).pdf
数论学习课件
06-18
学习资源，课件，整体把握
数论导引 数论学习 华罗庚编写
01-06
这是关于数论导引的一本书 是rar格式 是学习数论的必备书目
数论学习之起步篇（一）
445822357
10-03 139
1. 除法表达式（gcd法求最大公约数） //给你一个除法表达式：X1/X2/X3/X4……/Xk //通过添加括号，问能否得到整数 //分析可得，只有X2是不能做分子的。题目就转化为求一个分数能否为整数 #include&lt;iostream&gt; #include&lt;cstdio&gt; #include&lt;cstring&gt; #include&lt;cmath&gt; ...
数论学习-初等数论基础总览
麒麟001号的博客
02-08 2897
文章目录初等数论基础（二）建议先看一、数论只会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,
Codeforces 86C Genetic engineering
三叶的专栏
01-11 612
orz题解的Dp，竟然没调就1A了。。。。。 #include #include #include #include #include #include #include using namespace std; #define N 1010 #define Mod 1000000009 int n,m; int Dp[N][110][11],id[3
数论学习总结
LjcoderDCC的博客
05-24 5863
目录 一、数论基本概念       1、整除性       2、素数            a.素数与合数            b.素数判定            c.素数定理            d.素数筛选法       3、因数分解            a.算术基本定理            b.素数拆分            c.因子个数
三叶青
weixin_30949361的博客
07-25 133
植物：三叶青 中文学名：三叶崖爬藤 拉丁学名：Tetrastigma hemsleyanum Diels et Gilg 别名： 三叶青、丝线吊金钟、蛇附子、土经丸、石老鼠、石猴子、小扁藤 分布范围： 产云南东南部。生于海拔1100-1300米山谷中。分布于印度（安达曼群岛）、缅甸、泰国、越南、马来西亚及印度尼西亚。 生长环境： 喜凉爽气候，适温在25℃左右生长健壮，冬季...
教你认识 阿迪达斯（ADI） 鞋标
热门推荐
lj89168的专栏
05-20 2万+
1. AD 标的第一部分为 鞋子 的产地栏。 迄今为止 AD 主要是在中国福建，越南，印尼以及印度地区设立了 OEM 工厂。每个地区 OEM 工厂制造鞋子的类型不同。与 NIKE 的标一样，只要你看得懂英文就知道产地是哪里了。不过现在香港 专卖 的鞋子都没有标注产地，也就是如下面这两个标一样，产地栏为空。 此时，很多筒子就会问了 “ 那我们如何知道鞋子的产地呢
基础数论_杜德利.pdf
最新发布
10-30
总的来说，这本《基础数论_杜德利.pdf》文件提供了一个很好的入门级的数论学习资源。通过学习这本文件，读者可以掌握基础数论中的重要概念、定理和算法，并且了解它们的应用。对于对数论感兴趣或需要应用数论知识的...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值