【殺雞焉用牛刀?即便可以殺也要在乎雞的感受!選取合適的方法可以減少出錯率】(這就是爲什麼我要嗶嗶三種方法)
1:Pascal公式打表
const int N=3000;
long long C[N][N];
///組合數打表模板,適用於N<=3000
///c[i][j]表示從i箇中選j個的選法。
void get_C(int maxn)
{
C[0][0] = 1;
for(int i=1;i<=maxn;i++)
{
C[i][0] = 1;
for(int j=1;j<=i;j++)
C[i][j] = C[i-1][j]+C[i-1][j-1];
//C[i][j] = (C[i-1][j]+C[i-1][j-1])%MOD;
}
}
2:逆元求取組合數
注意:1:初始化caljc()函數
2:nm都不可以大於mod
typedef long long LL;
const LL maxn(1000005), mod(1e6 + 3);
LL Jc[maxn];
void calJc()
{
Jc[0] = Jc[1] = 1;
for(LL i = 2; i < maxn; i++)
Jc[i] = Jc[i - 1] * i % mod;
}
LL pow(LL a, LL n, LL p)
{
LL ans = 1;
while(n)
{
if(n & 1) ans = ans * a % p;
a = a * a % p;
n >>= 1;
}
return ans;
}
LL niYuan(LL a, LL b)
{
return pow(a, b - 2, b);
}
LL C(LL a, LL b)///在主函數中求C(ab)
{
return Jc[a] * niYuan(Jc[b], mod) % mod
* niYuan(Jc[a - b], mod) % mod;
}
盧卡斯(Lucas)定理(mn任何一個大於mod)
LL Lucas(LL a, LL b)
{
if(a < mod && b < mod)
return C(a, b);
return
C(a % mod, b % mod) * Lucas(a / mod, b / mod);
}
///只需要加這個到上面的逆元求組合數裏面(在主函數裏面用這個公式求C(ab))