算法（组合数）

给定n组询问，每组询问给定两个整数a，b，请你输出Cba mod (109+7)的值。

输入格式
第一行包含整数n。

接下来n行，每行包含一组a和b。

输出格式
共n行，每行输出一个询问的解。

数据范围
1≤n≤10000,
第一个范围1≤b≤a≤2000
输入样例：
3
3 1
5 3
2 2
输出样例：
3
10
1

公式如下 （递归） N^2

#include<iostream>
#include<cstdio>
using namespace std;
const int mod=1e9+7;
const int N=2010;
int n;
int c[N][N];
//从1-2000进行初始化
void calculate()
{
    
    for(int i=0;i<N;i++)
    for(int j=0;j<=i;j++)//j本身就要小于i,C 34 是错误的
    {
        if(!j)//j=0时
    c[i][j]=1;
     else c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
    
    }
    
    
}

int main()
{calculate();
    int a,b;
    cin>>n;
    while(n--)
    {
        scanf("%d%d",&a,&b);
       cout<<c[a][b]<<endl;
    }
    return 0;
}

第二个范围1≤b≤a≤10000
输入样例：
3
3 1
5 3
2 2
输出样例：
3
10
1

NlogN

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
using namespace std;
typedef long long LL;
//const int 、N=100010; // 会显示segmentation fault！！！改成1e+5即可
const int N=1e5+10,mod=1e9+7;
int fact[N],infact[N];//代表阶乘的逆元
//因为除法得到的数据不准确，改换 乘法
//快速幂求逆元
 int qmi(int a,int k,int p)

 {
     int res=1;
     while(k)
     {
         if(k&1)
            res=(LL)res*a%p;
         a=(LL)a*a%p;
         k>>=1;
     }
     return res;
 }

int main()
{
//预处理
   fact[0]=infact[0]=1;
   for(int i=1;i<N;i++)
   {
       fact[i]=(LL)fact[i-1]*i%mod;
       infact[i]=(LL)infact[i-1]*qmi(i,mod-2,mod)%mod;
   }
   int n;
   scanf("%d",&n);
   while(n--)
   {
       int a,b;
       cin>>a>>b;
       printf("%d\n",(LL)fact[a]*infact[b]%mod*infact[a-b]%mod);//及时取模
       //因为公式C ab，需要把b,a-b的阶乘放在分子上，为了避免除法出现的错误
       //使用乘法， 需要借助逆元
   }
    return 0;
}

第三个范围
1≤n≤20 ,
1≤b≤a≤10^18,
1≤p≤10^5,
输入样例
3
5 3 7
3 1 5
6 4 13
输出样例：
3
3
2

#include<iostream>
using namespace std;
typedef long long LL;
int p;//定义为全局变量，减少传值次数
LL qmi(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1)
        res=(LL)res*a%p;
        b>>=1;
        a=(LL)a*a%p;
        
    }
    return  res;
}
LL C(int a,int b)
{
    int res=1;
    for(int i=1,j=a;i<=b;i++,j--)
    {
        res=(LL)res*j%p;//a会受到i<=b的限制，不是a的阶乘
        res=(LL)res*qmi(i,p-2)%p;//p是全局变量，不需要总是传值
    }
    return res;
}
LL lucas(LL a,LL b)
{
    if(a<p&&b<p)
    return C(a,b);
    return C(a%p,b%p)*lucas(a/p,b/p)%p;
}

int main()
    {
        int n;
        cin>>n;
        while(n--)
        {
            LL a,b;
            cin>>a>>b>>p;
            cout<<lucas(a,b)<<endl;
            
        }
        return 0;
    }
    

第四种
输入a,b，求Cba的值。

注意结果可能很大，需要使用高精度计算。

输入格式
共一行，包含两个整数a和b。

输出格式
共一行，输出Cba的值。

数据范围
1≤b≤a≤5000
输入样例：
5 3
输出样例：
10

#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
const int N=5010;
int primes[N],cnt;
int sum[N];
bool st[N];
//筛素数,求质数 
void getprimes(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!st[i])primes[cnt++]=i;
        for(int j=0;primes[j]<=n/i;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)
            break;
        }
    }
}
 int get (int n,int p)
 {
     int res=0;
     while(n)
     {
         res+=n/p;
         n/=p;
     }
     return res;
 }
 //高精度乘法吧所有的质因子乘在一起
 vector<int>mul(vector<int> a,int b)
 {
     vector<int>c;
     int t=0;
     for(int i=0;i<a.size();i++)
     {
         t+=a[i]*b;
         c.push_back(t%10);
         t/=10;
     }
     while(t)
     {
         c.push_back(t%10);
         t/=10;
     }
     return c;
 }
 
 int main()
 {
     int a,b;
     cin>>a>>b;
     getprimes(a);
     for(int i=0;i<cnt;i++)
     {
         int p=primes[i];
         sum[i]=get(a,p)-get(b,p)-get(a-b,p);
     }
     vector<int> res;
     res.push_back(1);
     for(int i=0;i<cnt;i++)
     for(int j=0;j<sum[i];j++)
     res=mul(res,primes[i]);
     for(int i=res.size()-1;i>=0;i--)
     printf("%d",res[i]);
     puts("");
     return 0;
 }