Advertisement
Mlxa

### Модульная арифметика

Feb 12th, 2019
125
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 0.89 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 N = 1.1e5;
  8. const int MOD = 1e9 + 7;
  9.  
  10. int add(int a, int b) {
  11.     return (a + b) % MOD;
  12. }
  13.  
  14. int mul(long a, long b) {
  15.     return (int)(a * b % MOD);
  16. }
  17.  
  18. int power_mod(int a, int p) {
  19.     int x = 1;
  20.     for (; p > 0; p /= 2) {
  21.         if (p % 2) {
  22.             x = mul(x, a);
  23.         }
  24.         a = mul(a, a);
  25.     }
  26.     return x;
  27. }
  28.  
  29. int divide(int a, int b) {
  30.     return mul(a, power_mod(b, MOD - 2));
  31. }
  32.  
  33. int fac[N], inv[N];
  34.  
  35. void precalc_factorials() {
  36.     fac[0] = 1;
  37.     for (int i = 1; i < N; ++i) {
  38.         fac[i] = mul(i, fac[i - 1]);
  39.     }
  40.     inv[N - 1] = divide(1, fac[N - 1]);
  41.     for (int i = N - 2; i >= 0; --i) {
  42.         inv[i] = mul(inv[i + 1], i + 1);
  43.     }
  44. }
  45.  
  46. int main() {
  47. #ifdef LC
  48.     assert(freopen("input.txt", "r", stdin));
  49. #endif
  50.     ios::sync_with_stdio(false);
  51.     cin.tie(nullptr);
  52.    
  53.     return 0;
  54. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement