NOIP2016

数论

gcd与lcm

#include<algorithm>
using namespace std;
int a=12,b=20;
int main(){
    cout<<__gcd(a,b)<<endl;//求gcd 位于algorithm库 注意加两个下划线
    cout<<a/__gcd(a,b)*b<<endl;//求lcm 注意先除再乘防止溢出    
}

扩展欧几里得

#include<iostream>
using namespace std;

void exgcd(int a,int b,int &d,int &x,int &y){
    if(!b) {d=a;x=1;y=0;}
    else{exgcd(b,a%b,d,y,x);y-=x*(a/b);}
}
int main(){
    int x,y,d,a,b;
    a=6,b=15;
    exgcd(a,b,d,x,y);
    cout<<x<<' '<<y;
}

素数

//素数判断 s1单个 s2筛表

int n;
bool vis[10000000],flag;
bool s1(int n){
    int m=sqrt(n+0.5);
    for(int i=2;i<=m;i++)if(!vis[i]){
        if(!(n%i))return 0;
        for(int j=i*i;j<=n;j+=i)vis[j]=1;   
    }
    return 1;

}
void s2(int n){
    int m=sqrt(n+0.5);
    memset(vis,0,sizeof(vis));
    for(int i=2;i<=m;i++)if(!vis[i])
      for(int j=i*i;j<=n;j+=i) vis[j]=1;
}

int main(){
    cin>>n;
    cout << s1(n);
    s2(n);for(int i=1;i<=n;i++)if(!vis[i])cout<<i<<' ';

}

素数

埃氏筛法

int m=sqrt(n+0.5);
memset(vis,0,sizeod(vis);
for(int i=2;i<=m;i++)if(!vis[i])
    for(int j=i*i;j<=n;j+=i)vis[j]=1;

快速幂

int pow(int a,int b,int p)
{
    int ans=1;
    //a%=k;
    for(;b;b>>=1,a=a*a,a%=p)
    {
        if(b&1) ans=ans*a;
        ans%=p;
    }
    return ans;
}

int main()
{
    int a,b,p;
    cin>>a>>b>>p;
    cout<<pow(a,b,p);
    return 0;
}

排列组合

排列

STL 中有排列算法。头文件为algorithm
next_permutation(序列第一项的地址, 序列最后一项的地址+1)：产生下一排列。
prev_permutation(序列第一项的地址, 序列最后一项的地址+1)：产生上一排列。

组合

int C[1001][1001];//根据实际需要开数组，必要时采用高精度类型
memset(C,0,sizeof(C));
for (int i=0;i<=n;i++)
{
C[i][0]=1;
for (int j=0;j<=i;j++)
    C[i][j]=C[i-1][j-1] + C[i-1][j];
}