[CSP-S模拟测试]:排列组合(数学 or 找规律)
Posted wzc521
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[CSP-S模拟测试]:排列组合(数学 or 找规律)相关的知识,希望对你有一定的参考价值。
题目描述
$T$组数据,每次给定$n$,请求出下式的值,对$10^9+7$取模:
$$C_n^0 imes C_n^0+C_n^1 imes C_n^1+C_n^2 imes C_n^2+...+C_n^n imes C_n^n$$
输入格式
第一行一个整数$T$,表示数据组数。
接下来$T$行,每一行包含一个整数$n$,含义如题所示。
输出格式
输出$T$行,每行包含一个整数,表示对$10^9+7$取模后的答案。
样例
样例输入:
2
1
2
样例输出:
2
6
数据范围与提示
对于$30\\%$的数据,$Tleqslant 500,nleqslant 10,000$。
对于$100\\%$的数据,$Tleqslant 100,000,nleqslant 1,000,000$。
题解
打表找规律可以发现答案其实就是$C_{2n}^n$,那么我们现在来讲一下这到底是为什么。
先来讲个故事:
我家门前有两棵树,一棵是枣树,另一棵也是枣树。——鲁迅
我家门前有两棵枣树,一棵枣树上有$n$颗枣,另一棵枣树上也有$n$颗枣。——$HEOI-$动动
麻麻让我在这两棵枣树上摘$n$棵枣,有多少种方案数呢?
我们假设在第一棵树上摘$i$颗枣,那么另一棵树上要摘$n-i$颗枣,方案数就是:$C_n^i imes C_n^{n-i}$,那么总的方案数就是$sum limits_{i=0}^n C_n^i imes C_n^{n-i}$。
然而,枣是一样的,树也是一样的,那么方案数也可以写成:$C_{2n}^n$。
也就是说$C_{2n}^n=sum limits_{i=0}^n C_n^i imes C_n^{n-i}$。
我们还知道$C_n^i=C_n{n-i}$,于是上式就变成了:$C_{2n}^n=sum limits_{i=0}^n {(C_n^i)}^2$。
那么$C_n^0 imes C_n^0+C_n^1 imes C_n^1+C_n^2 imes C_n^2+...+C_n^n imes C_n^n=C_{2n}^n$。
故事讲完啦,小朋友们都懂了嘛~
时间复杂度:$Theta(T imes log_{mod}n)$。
期望得分:$100$分。
实际得分:$100$分。
代码时刻
#include<bits/stdc++.h>
using namespace std;
const int mod=1000000007;
int n;
long long jc[2000001],inv[2000001];
long long qpow(long long x,long long y)
{
long long res=1;
while(y)
{
if(y%2)res=res*x%mod;
y>>=1;
x=x*x%mod;
}
return res;
}
void pre_work()
{
jc[0]=1;
for(long long i=1;i<=2000000;i++)
jc[i]=jc[i-1]*i%mod;
inv[2000000]=qpow(jc[2000000],1000000005)%mod;
for(long long i=2000000;i>0;i--)
inv[i-1]=inv[i]*i%mod;
}
long long get_C(int x,int y){return jc[x]*inv[y]%mod*inv[x-y]%mod;}
int lucas(int x,int y)
{
if(!y)return 1;
return get_C(x%mod,y%mod)*lucas(x/mod,y/mod)%mod;
}
int main()
{
pre_work();
int T;scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
printf("%d
",lucas(n<<1,n));
}
return 0;
}
rp++
以上是关于[CSP-S模拟测试]:排列组合(数学 or 找规律)的主要内容,如果未能解决你的问题,请参考以下文章