Advertisement
Mlxa

### Комбинаторика

Feb 19th, 2019
135
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.15 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. typedef long long ll;
  4. #define long ll
  5. #define all(x) begin(x), end(x)
  6.  
  7. const int MOD = 1e9 + 7;
  8. const int N = 1e6;
  9.  
  10. int add(int a, int b) {
  11.     return (a + b) % MOD;
  12. }
  13.  
  14. int sub(int a, int b) {
  15.     return (a + MOD - b) % MOD;
  16. }
  17.  
  18. int mul(int a, int b) {
  19.     return (int)((long)a * b % MOD);
  20. }
  21.  
  22. int power_mod(int a, int p) {
  23.     int x = 1;
  24.     for (; p > 0; p >>= 1) {
  25.         if (p & 1) {
  26.             x = mul(x, a);
  27.         }
  28.         a = mul(a, a);
  29.     }
  30.     return x;
  31. }
  32.  
  33. int divide(int a, int b) {
  34.     return mul(a, power_mod(b, MOD - 2));
  35. }
  36.  
  37. int fact[N], invfact[N];
  38.  
  39. int c_func(int n, int k) {
  40.     assert(0 <= k && k <= n);
  41.     return mul(fact[n], mul(invfact[k], invfact[n - k]));
  42. }
  43.  
  44. int a_func(int n, int k) {
  45.     assert(0 <= k && k <= n);
  46.     return mul(fact[n], invfact[n - k]);
  47. }
  48.  
  49. int main() {
  50. #ifdef LC
  51.     //assert(freopen("input.txt", "r", stdin));
  52. #endif
  53.     ios::sync_with_stdio(false);
  54.     cin.tie(nullptr);
  55.    
  56.     fact[0] = 1;
  57.     for (int i = 1; i < N; ++i) {
  58.         fact[i] = mul(fact[i - 1], i);
  59.     }
  60.     invfact[N - 1] = divide(1, fact[N - 1]);
  61.     for (int i = N - 2; i >= 0; --i) {
  62.         invfact[i] = mul(invfact[i + 1], i + 1);
  63.     }
  64.        
  65.     return 0;
  66. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement