1. 程式人生 > >復習1 數論

復習1 數論

class sqrt http pre bsp sin out png com

1.歐幾裏得

#include<cstdio>
using namespace std;
int a,b;
int gcd(int x,int y){
    if(y==0)return x;
    return gcd(y,x%y);
}
int main(){
    scanf("%d%d",&a,&b);
    printf("%d",gcd(a,b));
}

2.擴展歐幾裏得

#include<iostream>
#include<cstdio>
#include<cstdlib>
using
namespace std; #define ll long long ll a,b,x,y; ll exgcd(ll a,ll b){ if(b==0){ x=1;y=0; return a; } exgcd(b,a%b); ll tmp=x; x=y;y=tmp-(a/b)*y; } int main(){ scanf("%lld%lld",&a,&b); exgcd(a,b); while(x<0)x+=b; cout<<x; }

3.歐拉篩

void prepare(){
    for(int i=2;i<=n;i++){
        if(!vis[i])p[++cnt]=i;
        for(int j=1;j<=cnt&&i*p[j]<=n;j++){
            vis[i*p[j]]=1;
            if(i%p[j]==0)break;
        }
    }
}

4.歐拉函數

技術分享

int euler_phi(int n){//單個歐拉函數 
    int m=(int)sqrt(n+0.5
); int ans=n; for(int i=2;i<=m;i++)if(n%i==0){ ans=ans/i*(i-1); while(n%i==0)n/=i; } if(n>1)ans=ans/n*(n-1); }

5.中國剩余定理

6.卡特蘭數

7.逆元(求組合數)

8.盧卡斯定理(求組合數取模)

復習1 數論