数论模板

组合数

const int MOD = 1e9 + 7;
const int M = 1e6 + 5;
int F[M], Finv[M], inv[M];//F是阶乘，Finv是逆元的阶乘 
void init()
{
    inv[1] = 1;
    for(int i = 2; i < M; i ++)
    {
        inv[i] = (MOD - MOD / i) * 1ll * inv[MOD % i] % MOD;
    }
    F[0] = Finv[0] = 1;
    for(int i = 1; i < M; i ++)
    {
        F[i] = F[i-1] * 1ll * i % MOD;
        Finv[i] = Finv[i-1] * 1ll * inv[i] % MOD;
    }
}
int C(int n, int m) // m <= n
{
    if(m < 0 || m > n) return 0;
    return F[n] * 1ll * Finv[n - m] % MOD * Finv[m] % MOD;
}

• 别忘了在主函数离调用打表函数啊QAQ

gcd

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

快速幂

LL pow_mod(LL a, LL b, LL p) //a的b次方求余p 
{
    LL ret = 1;
    while(b)
    {
        if(b & 1) ret = (ret * a) % p;
        a = (a * a) % p;
        b >>= 1;
    }
    return ret;
}

快速乘

LL mul(LL a, LL b, LL p) //快速乘，计算a*b%p 
{
    LL ret = 0;
    while(b)
    {
        if(b & 1) ret = (ret + a) % p;
        a = (a + a) % p;
        b >>= 1;
    }
    return ret;
}

（不断更新。。。感谢仓鼠学姐的博客orz）
http://www.cnblogs.com/linyujun/