你在一栋楼房下面,楼房一共有n层,第i层每秒有pi的概率会扔下一个东西并砸到你
求第一秒内你被砸到的概率

输入描述:

第一行一个整数n之后有n行,第i+1行有两个整数a

i

,b

i

,表示

输出描述:

设答案为

,你只需要找到一个最小的非负整数T,使得

输出这个T就行了
示例1

输入

复制

2
1 2
1 2

输出

复制

750000006

说明

一共只有如下状态:

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

2. 第一层扔了下来

3. 第二层扔了下来

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

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

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

因此被砸到的概率是 3/4,这个值在模1e9+7意义下就是750000006

备注:

数据范围0 ≤ n ≤ 10

5

1 ≤ a

i

 ≤ b

i

 ≤ 10

5

思路:题意很好懂,主要是求逆元
对逆元不熟悉
例如:4关于1模7的乘法逆元为多少?
4X≡1 mod 7
这个方程等价于求一个X和K,满足
4X=7K+1
其中X和K都是整数。
若ax≡1 mod f, 则称a关于1模f的乘法逆元为x。也可表示为ax≡1(mod f)。
当a与f互素时,a关于模f的乘法逆元有解。如果不互素,则无解。如果f为素数,则从1到f-1的任意数都与f互素,即在1到f-1之间都恰好有一个关于模f的乘法逆元。
求此算法还可以使用费马小定理
只不过局限性比较大,要求模数是素数
a^(p-1)

1(mod p)

p要求是素数
那么a^(p-2)就是a的乘法逆元

首先是逆元的模板:
 1 // 扩展欧几里得做法;
 2 #include<iostream>
 3 #include<cstdio>
 4 #include<algorithm>
 5 #include<cstring>
 6 #define ll long long
 7 using namespace std;
 8
 9 ll ex_gcd(ll a,ll b,ll &x,ll &y)//扩展欧几里得(扩展gcd)
10 {
11     if (a==0&&b==0) return -1;
12     if (b==0){x=1;y=0;return a;}
13     ll d=ex_gcd(b,a%b,y,x);
14     y-=a/b*x;
15     return d;
16 }
17
18 ll mod_inverse(ll a,ll mod)//乘法逆元
19 {
20     ll x,y;
21     ll d = ex_gcd(a,mod,x,y);
22     return (x%mod+mod)%mod;
23 }
24 int low_bit(int x){return x&(-x);}
25 int main()
26 {
27     for(int i=0;i<=16;i++)
28         cout<<i<<' '<<low_bit(i)<<endl;
29     return 0;
30 }

解题代码如下:

 1 #include<iostream>
 2 using namespace std;
 3 #define ll long long
 4 ll ex_gcd(ll a,ll b,ll &x,ll &y)//扩展欧几里得(扩展gcd)
 5 {
 6     if (a==0&&b==0) return -1;
 7     if (b==0){x=1;y=0;return a;}
 8     ll d=ex_gcd(b,a%b,y,x);
 9     y-=a/b*x;
10     return d;
11 }
12
13 ll mod_inverse(ll a,ll mod)//乘法逆元
14 {
15     ll x,y;
16     ll d = ex_gcd(a,mod,x,y);
17     return (x%mod+mod)%mod;
18 }
19 long long gcd(long long a,long long b)
20 {
21     long long c=1;
22     while(c)
23     {
24         c=a%b;
25         a=b;
26         b=c;
27     }
28     return a;
29 }
30 int main()
31 {
32     int n;
33     cin>>n;
34     long long a,b,c,sum_a=1,sum_b=1;
35     for(int i=0;i<n;i++)
36     {
37         cin>>a>>b;
38         c=gcd(a,b);
39         a/=c;
40         b/=c;
41         sum_a*=b-a;
42         sum_b*=b;
43         sum_b%=1000000007;
44         sum_a%=1000000007;
45     }
46     c=gcd(sum_a,sum_b);
47     sum_a/=c;
48     sum_b/=c;
49     sum_a=(sum_b-sum_a+1000000007)%10000000007;
50     c=mod_inverse(sum_b,1000000007);
51     cout<<c*sum_a%1000000007;
52 }

转载于:https://www.cnblogs.com/lu1nacy/p/10013853.html

牛客小白月赛9 A签到(乘法逆元)相关推荐

  1. 牛客小白月赛16 小石的签到题(博弈)

    牛客小白月赛16 小石的签到题 链接:https://ac.nowcoder.com/acm/contest/949/A来源:牛客网 题目描述 输入描述: 共一行,输入一个数 nnn . 输出描述: ...

  2. 2021牛客小白月赛39

    2021牛客小白月赛39(A+B+C+E+G+H) A.憧憬 B.欢欣 C.奋发 D.绝望 E.迷惘 F.孤独 G.冷静 H.终别 A.憧憬 题目链接:https://ac.nowcoder.com/ ...

  3. F.孤独(牛客小白月赛39)

    F.孤独(牛客小白月赛39) 题意: 给定一棵树,寻找一个路径,将断掉所有与这个路径上的点相连的边,使得剩下的最大连通块的大小最小 题解: 这题有点印象,感觉做过,至少这个方法肯定遇到过 设dp[u] ...

  4. 牛客小白月赛58 B(暴力)C(思维)D(dp滚动数组优化)

    牛客小白月赛58 感觉没什么意思,这场月赛出题人有点问题. B 题意:给定一定的数据填充顺序和每一个填充层的名称,只有填充满了上一层才可以填充下一层. 但是每一层又属于某一个大层,一个大层中包含了若干 ...

  5. 牛客小白月赛65个人题解A-E

    1. 牛客小白月赛65 A. 牛牛去购物 题意:给定n元,购买价格为a元的篮球和价格为b的篮球,数量不定,要使得花掉的钱最多,也就是剩余的钱数最少,求这个值 (1 <= n, a, b < ...

  6. 牛客小白月赛24 J.建设道路

    牛客小白月赛24 J.建设道路 题目链接 题目描述 牛牛国有 nnn 个城市,编号为 1-n,第 iii 个城市有一个价值 aia_iai​ ,牛国的国王牛阔落特别喜欢在牛牛国旅游,并且他不想每次旅游 ...

  7. 牛客小白月赛4 D.郊区春游

    牛客小白月赛4 D.郊区春游 题目链接 题目描述 今天春天铁子的班上组织了一场春游,在铁子的城市里有 nnn 个郊区和 mmm 条无向道路,第 iii 条道路连接郊区 AiA_iAi​ 和 BiB_i ...

  8. 牛客小白月赛25 C.白魔法师

    牛客小白月赛25 C.白魔法师 题目链接 题目描述 你是一个白魔法师. 现在你拿到了一棵树,树上有 个点,每个点被染成了黑色或白色. 你可以释放一次魔法,将某个点染成白色.(该点不一定是黑色点,也可以 ...

  9. 牛客小白月赛60(A~C)题解

    原题地址:牛客小白月赛60_ACM/NOI/CSP/CCPC/ICPC算法编程高难度练习赛_牛客竞赛OJ 题目A:小竹与妈妈 题意:输出(x-b)/a即可,没什么好说 代码: #include< ...

最新文章

  1. 黑龙江科技大学计算机类分数线,2019年黑龙江科技大学优势专业排名及分数线...
  2. 多线程的操作与数据绑定
  3. flink报错:Error: Static methods in interface require -target:jvm-1.8 已解决
  4. MySQL复制常用拓扑结构详解
  5. 05.full_text match_bool_prefix match_phrase_prefix 查询
  6. 在Webstorm中配置Compass Watcher
  7. xamarin android 通知,在 Xamarin.Android 中使用 Notification.Builder 构建通知
  8. C语言 数组排序 – 冒泡法排序 - C语言零基础入门教程
  9. PlantUML in a nutshell(官方文档)
  10. Vue打包并发布项目
  11. jQuery文档就绪事件
  12. 嫡权法赋权法_三种权重赋权法的比较分析
  13. 十三水牌型 图片_十三水三分天下是什么 十三水三分天下牌型详解
  14. 可恶的零宽空格—ZWSP
  15. H5如何获取内网IP和公网IP
  16. 小米忙着营销,麻烦带上技术!
  17. 捣鼓车间 | 学生获奖作品:戒烟帽
  18. x64dbg 实现插件Socket反向通信
  19. JavaScript中文与阿拉伯数字互相转换
  20. 2021-“新“的开源项目之handpose_x(手势识别交互)

热门文章

  1. Microsoft Blazor——快速开发与SQL Forms开源平台Platz.SqlForms
  2. 边缘设备上的实时AI人员检测:使用预先训练的SSD模型检测人员
  3. 使用Newtonsoft.Json接受部分资源
  4. 沉船会有什么_世界六大沉船,最后一艘有40吨财宝!8吨黄金!你眼红了吗?
  5. 用科学计算机log21 p x,山东理工大学数据结构期末 试题及答案
  6. quartz 每30秒执行一次_Net Core 使用 Quartz
  7. xbox手柄 查看您的控制器保修_Xbox的“异形手柄”开始让残疾的老兵重新享受游戏的乐趣...
  8. java executor 异步_Java并发编程11-异步执行框架Executor
  9. mysql utf8mb4 bin_MySQL中utf8字符集、排序规则及utf8mb4_bin列大小写不敏感方法
  10. android 设置功能,NDK 使用入门  |  Android NDK  |  Android Developers