【模板】一坨数学算法

求GCD

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

线性筛求[1,n]的质数

bool isprime[1000];
int prime[100],tot;
void pri(int n)
{
    tot = 0;
    memset(isprime,true,sizeof(isprime));
    int i,j;
    for(i=2;i<=n;i++)
    {
        if(isprime[i]) prime[++tot] = i;
        for(j=1;j<=tot&&i*prime[j]<=n;j++)
        {
            isprime[i*prime[j]] = false;
            if(i%prime[j]==0) break;
        }
    }
}

扩展欧几里得算法

int exgcd(int a,int b,int &x,int &y)
{
    if(!b)
    {
        x = 1,y = 0;
        return a;
    }
    int ret = exgcd(b,a%b,x,y);
    int t = x;
    x = y;
    y = t - a/b * y;
    return ret;
}