sincerit 签到(分数取模，逆元，费马小定理)

链接：https://ac.nowcoder.com/acm/contest/275/A?&headNav=acm
来源：牛客网

题目描述
你在一栋楼房下面，楼房一共有n层，第i层每秒有pi的概率会扔下一个东西并砸到你
求第一秒内你被砸到的概率
输入描述:
第一行一个整数n
之后有n行，第i+1行有两个整数ai,bi，表示
输出描述:
设答案为，你只需要找到一个最小的非负整数T，使得
输出这个T就行了
示例1
输入
复制
2
1 2
1 2
输出
复制
750000006
说明
一共只有如下状态：

1. 第一层和第二层都扔了下来

2. 第一层扔了下来

3. 第二层扔了下来

4. 第一层和第二层都没有扔下来

以上四种都是等概率发生的

除了第四种情况外，都会被砸到

因此被砸到的概率是 3/4，这个值在模1e9+7意义下就是750000006
备注:
数据范围
0 ≤ n ≤ 105
1 ≤ ai ≤ bi ≤ 105

概率题目从正向去做会很麻烦，可以考虑反向即考虑不被砸中的概率s， 然后p = 1 - s;
这里的难点是分数取模

#include <stdio.h>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const ll MOD  = 1e9 + 7;
ll spow(ll s, ll t) {
  ll ans = 1;
  while (t) {
    if (t&1) ans = ans * s % MOD;
    t >>= 1;
    s = (s % MOD * s) % MOD;
  }
  return ans;
}
// 求b的逆元 == spow(b, MOD-2); 
ll inv(ll b) {
  return spow(b, MOD-2);
}

int main() {
  ll n, a, b, ans = 1;
  scanf("%lld", &n);
  for (int i = 0; i < n; i++) {
    scanf("%lld%lld", &a, &b);
    ll k = b - a;
    ans = (ans * k % MOD * inv(b)) % MOD;
  }
  // ans求得的是未被砸中的答案
  // 被砸得概率 1-ans+MOD
  printf("%lld\n", (1-ans+MOD) % MOD);
  return 0;
}

求x关于模m的逆元 逆元就是 x^m-2

费马小定理解释求逆元
费马小定理： 假如p是质数，且a,p互质，那么 a的(p-1)次方除以p的余数恒等于1
则等式就是： a^(p-1) = 1 (mod p) ----a^(p-1)除以p的余数为1
假设x为a的逆元(a的倒数) 则有 ax = 1 (mod p) -----一个数乘以它的倒数为1
所以 a^(p-1) = ax , x = a^(p-2)—快速幂可求

当求解公式：(a/b)%m 时，因b可能会过大，会出现爆精度的情况，所以需变除法为乘法：
设c是b的逆元，则有bc≡1(mod m)；
则(a/b)%m = (a/b)1%m = (a/b)bc%m = ac(mod m);
即a/b的模等于ab的逆元的模；
逆元就是这样应用的；