其他分享
首页 > 其他分享> > 2019西安联训B层 Day 5 test T2 排列组合

2019西安联训B层 Day 5 test T2 排列组合

作者:互联网

做法:组合数取模

其实我们观察那么长一个式子,其实它有一个结论,在数学竞赛中也常提到,结果就是C(2n,n),相当于从2n个数选n个数出来。因为阶乘的结果太大,所以我们还需要用到逆元。

对于逆元我也会写博客来讲解,求逆元有多种方法,扩欧,费马小定理+快速幂,递推打表,递归(会爆栈)

扩欧求逆元

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
const int maxn=1e7+7;
int T,n;
long long f[maxn];
long long ans;
void init(){
    f[0]=1;
    for(int i=1;i<=2000010;i++){
        f[i]=f[i-1]*i%mod;
    }
}
long long Exgcd(long long a,long long &x,long long b,long long &y){
    if(b==0){
        x=1;
        y=0;
        return a;
    }
    long long GCD=Exgcd(b,x,a%b,y);
    long long tmp=x;
    x=y;
    y=tmp-(a/b)*y;
    return GCD;
}
long long inv(long long a,long long b){
    long long x,y;
    if(Exgcd(a,x,b,y)!=1) return -1;
    else return (x%b+b)%b;
}
int main()
{
    scanf("%d",&T);
    init();
    while(T--){
        scanf("%d",&n);
        printf("%lld\n",f[2*n]*inv(f[n],mod)%mod*inv(f[n],mod)%mod);
    }
    return 0;
}

快速幂求逆元

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
const int maxn=1e7+7;
long long f[maxn];
void jc(){
    f[0]=1;
    for(int i=1;i<=2000010;i++){
        f[i]=f[i-1]*i%mod;
    }
}
long long ksm(long long a,long long b){
    long long base=1;
    while(b){
        if(b&1) base=base*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return base;
}
int T,n;
int main()
{
    scanf("%d",&T);
    jc();
    while(T--){
        scanf("%d",&n);
        printf("%lld\n",f[2*n]*ksm(f[n],mod-2)%mod*ksm(f[n],mod-2)%mod); 
    }
    return 0;
}

结果

 

标签:const,int,逆元,T2,long,联训,2019,maxn,mod
来源: https://www.cnblogs.com/LJB666/p/11005902.html