当前位置:首页 > 编程笔记 > 正文
已解决

组合数 2.1 2.2

来自网友在路上 155855提问 提问时间:2023-09-23 22:25:56阅读次数: 55

最佳答案 问答题库558位专家为你答疑解惑

 O(nlogn)预处理, O(1)查询

#include<bits/stdc++.h>
#define IOS ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
#define endl '\n'using namespace std;typedef pair<int, int> PII;
typedef long long ll;
typedef long double ld;const int N = 100010, mod = 1e9 + 7;int n;
int fact[N], infact[N];int qmi(int a, int k)
{int res = 1;while(k){if(k & 1)res = (ll)res * a % mod;a = (ll)a * a % mod;k >>= 1;}return res;
}int C(int a, int b)
{return (ll)fact[a] * infact[b] % mod * infact[a - b] % mod;
}int main()
{IOScin >> n;fact[0] = infact[0] = 1;for(int i = 1; i < N; i ++){fact[i] = ((ll)fact[i - 1] * i) % mod;infact[i] = qmi(fact[i], mod - 2);}while(n --){int a, b;cin >> a >> b;cout << C(a, b) << endl;}return 0;
}

O(n)预处理,O(1)查询

#include<bits/stdc++.h>
#define IOS ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
#define endl '\n'using namespace std;typedef pair<int, int> PII;
typedef long long ll;
typedef long double ld;const int N = 100010, mod = 1e9 + 7;int n;
int fact[N], infact[N], inv[N];int C(int a, int b)
{return (ll)fact[a] * infact[b] % mod * infact[a - b] % mod;
}int main()
{IOScin >> n;fact[0] = fact[1]= infact[0] = infact[1] = inv[1] = 1;for(int i = 2; i < N; i ++){fact[i] = ((ll)fact[i - 1] * i) % mod;inv[i] = (ll)(mod - mod / i) * inv[mod % i] % mod;infact[i] = (ll)infact[i - 1] * inv[i] % mod;}while(n --){int a, b;cin >> a >> b;cout << C(a, b) << endl;}return 0;
}

查看全文

99%的人还看了

猜你感兴趣

版权申明

本文"组合数 2.1 2.2":http://eshow365.cn/6-12344-0.html 内容来自互联网,请自行判断内容的正确性。如有侵权请联系我们,立即删除!