组合数的求解
将n的阶乘和x ^ mod – 2(1 <= x <= 2000000)都初始化出来,大约可以减少一半以上的时间(具体不太清楚,应该会随着n的变化而变化)
代码:
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 2e6 + 5;
const LL mod = 1e9 + 7;
LL num[N];
LL ans[N];
//快速幂
LL qpow(LL a, LL b) {
LL res = 1;
while (b) {
if(b & 1) res = res * a % mod;
b >>= 1;
a = a * a % mod;
}
return res;
}
void init() {
num[1] = 1;
for (int i = 2; i <= 2000000; i ++ ) {
num[i] = (num[i - 1] * i) % mod;
}
ans[2000000] = qpow(num[2000000], mod - 2) % mod;
for (int i = 2000000; i >= 1; i -- ) {
ans[i - 1] = ans[i] * i % mod;
}
}
LL solve(LL m, LL n) {
return num[n] * ans[m] % mod * ans[n - m] % mod;
}
int main() {
int t;
cin >> t;
//将n的阶乘和x ^ mod - 2(1 <= x <= 2000000)都初始化出来,大约可以减少一半以上的时间
init();
while (t -- ) {
int n;
scanf("%d", &n);
if(n == 1) {
cout << "1 1" << endl;
continue;
}
printf("%lld %lld\n", solve(n - 1, 2 * n), (solve(n, 2 * n) - solve(n - 1, 2 * n) + mod) % mod);
}
return 0;
}
原文地址:http://www.cnblogs.com/luoyefenglin/p/16870770.html
1. 本站所有资源来源于用户上传和网络,如有侵权请邮件联系站长!
2. 分享目的仅供大家学习和交流,请务用于商业用途!
3. 如果你也有好源码或者教程,可以到用户中心发布,分享有积分奖励和额外收入!
4. 本站提供的源码、模板、插件等等其他资源,都不包含技术服务请大家谅解!
5. 如有链接无法下载、失效或广告,请联系管理员处理!
6. 本站资源售价只是赞助,收取费用仅维持本站的日常运营所需!
7. 如遇到加密压缩包,默认解压密码为"gltf",如遇到无法解压的请联系管理员!
8. 因为资源和程序源码均为可复制品,所以不支持任何理由的退款兑现,请斟酌后支付下载
声明:如果标题没有注明"已测试"或者"测试可用"等字样的资源源码均未经过站长测试.特别注意没有标注的源码不保证任何可用性