计算组合数

有两种计算方式，直接计算或者递推计算：

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
#define ll long long
const int MOD = 1000007;
const int maxk = 500;
ll C[maxk + 2][maxk + 2];
// 线性算法，可以加取模
void get_C(){
    memset(C,0,sizeof C);
    C[0][0] = 1;
    for(int i = 0 ; i <= maxk ; i++){
        C[i][0] = C[i][i] = 1;
        for(int j = 1 ; j < i ; j++)
            C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]);
    }
}
// 直接计算，不要随便取模，计算量过大时会有误差
long long cal_C(long long n,long long m){
    double ans = 1;
    for(int i = 0 ; i < m ; i++) ans *= n - i;
    for(int i = 0 ; i < m ; i++) ans /= i + 1;
    return (long long)(ans + 0.5);
}
int main(){
    get_C();
    for(int i = 0 ; i < 225 ; i++)
        for(int j = 0 ; j <= i && j < 10 ; j++)if(C[i][j] != cal_C(i,j))
            printf("%d %d :%lld,%lld
",i,j,C[i][j],cal_C(i,j));
    return 0;
}