计蒜客 2017 ACM-ICPC 亚洲区(西安赛区)网络赛 B coin(求乘法逆元)

2024-05-16 03:40:18

Bob has a not even coin, every time he tosses the coin, the probability that the coin's front face up isqp(qp≤12)\frac{q}{p}(\frac{q}{p} \le \frac{1}{2})pq(pq≤21).

The question is, when Bob tosses the coin kkk times, what's the probability that the frequency of the coin facing up is even number.

If the answer is XY\frac{X}{Y}YX, because the answer could be extremely large, you only need to print (X∗Y−1)mod(109+7)(X * Y^{-1}) \mod (10^9+7)(X∗Y−1)mod(109+7).

Input Format

First line an integer TTT, indicates the number of test cases (T≤100T \le 100T≤100).

Then Each line has 333 integer p,q,k(1≤p,q,k≤107)p,q,k(1\le p,q,k \le 10^7)p,q,k(1≤p,q,k≤107) indicates the i-th test case.

Output Format

For each test case, print an integer in a single line indicates the answer.

样例输入

2
2 1 1
3 1 2

样例输出

500000004
555555560

题目来源

2017 ACM-ICPC 亚洲区（西安赛区）网络赛

题意:

计算扔k次硬币，正面朝上的次数为偶数的概率。再将答案mod(1e9+7)

解析:

正面朝上概率为a，反面朝上概率为b

a+b=1；

ans=C(k,0)*a^0*b^k+C(k,2)*a^2*b^(k-2)+...

(a+b)^k= C(k,0)*a^0 *b^k+ C(k,1)*a^1*b^(k-1) +C(k,2)*a^2*b^(k-2) .......+C(k,k)*a^k*b^0
(b-a)^k= C(k,0) *(-a)^0 *b^k +C(k,1)*(-a)* b^(k-1) .........+C(k,k)*(-a)^k*b^0

可得ans=((a+b)^k+(b-a)^k)/2 再将a=p/q,b=(1-p/q)带入

得(（p^k）+(p-2q)^k)/(2*p^k)

除法的模运算就是求乘法逆元

x*(x^-1)=1(mod p) 乘法逆元性质

d=(x/y)(mod p) => y*d=x(mod p) => y*y^-1*d=x*y^-1(mod p) => d=x*y^-1(mod p)

求乘法逆元的三种方法

http://blog.csdn.net/rain722/article/details/53170288

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;const long long int N = 1e9+7;long long p,q,k;
long long zi,mu;long long quickmulti(long long a,long long b)
{long long ans=1;while(b){if(b&1) ans = ((ans%N)*(a%N)) % N;a = ((a%N)*(a%N)) % N;b=b>>1;}return ans;
}long long inv2(long long b)
{return quickmulti(b,N-2);
}int main()
{int t;long long x,y,d,x0;scanf("%d",&t);while(t--){scanf("%lld%lld%lld",&p,&q,&k);zi=(quickmulti(p,k)+quickmulti(p-2*q,k))%N;mu=(2*quickmulti(p,k))%N;long long  ans=(zi*inv2(mu))%N;  printf("%lld\n",ans);}return 0;
}

计蒜客 2017 ACM-ICPC 亚洲区(西安赛区)网络赛 B coin(求乘法逆元)相关推荐

Maximum Flow（2017 ACM-ICPC 亚洲区(西安赛区)网络赛 E）
Problem Description Given a directed graph with nn nodes, labeled 0,1,⋯,n−1. For each <i, j> s ...
2017 ACM-ICPC 亚洲区(西安赛区)网络赛 E Maximum Flow
原题解链接:http://blog.csdn.net/kkkkahlua/article/details/78009087 他用的最小割的来求解最大流.认为只要讨论每一个点到0和n-1那个更小哪条边就 ...
计蒜客2020蓝桥杯大学A组模拟赛题解
计蒜客2020蓝桥杯大学A组模拟赛题解蓝桥杯的话,去年拿了C++组的国二.今年报名了新成立的Python组,不知道能不能摸到国一的鱼模拟赛链接如下: https://www.jisuanke.co ...
计蒜客 2020 蓝桥杯大学 A 组省赛模拟赛（一）题目及解析
新手,欢迎交流,剩下D题不知如何下手,望大佬指教. 计蒜客 2020 蓝桥杯大学 B 组省赛模拟赛 (一)题目及解析文章目录 A. 结果填空:计算周长 B. 结果填空:七巧板 C. 结果填空:苹果 ...
2016 ACM/ICPC亚洲区青岛站现场赛（部分题解）
摘要本文主要列举并求解了2016 ACM/ICPC亚洲区青岛站现场赛的部分真题,着重介绍了各个题目的解题思路,结合详细的AC代码,意在熟悉青岛赛区的出题策略,以备战2018青岛站现场赛. HDU 5 ...
计蒜客 2020 蓝桥杯大学 B 组省赛模拟赛（一）题目及解析
新手,欢迎交流,更新完毕. 计蒜客 2020 蓝桥杯大学 A 组省赛模拟赛 (一)题目及解析文章目录 A. 结果填空:有趣的数字 B. 结果填空:爬楼梯 C. 结果填空:七巧板 D. 结果填空:苹果 ...
【计蒜客模拟赛系列】-计蒜客2021年8月普及组模拟赛
提前:本文中部分代码和思路有借鉴或摘抄计蒜客官方题解赛后总结本次模拟赛的难度总算正常了些个人战绩: 220/400,排名61 ,太弱了,一大堆AK爷题目质量评价: 题目相比CSP-J还是简单了 ...
计蒜客-2017 ACM-ICPC 亚洲区(乌鲁木齐赛区)网络赛H题Skiing(拓扑序求DAG最长路)
题意: 给定一个有向无环图,求该图的最长路. 思路: 由于是有向无环图,所以最长路肯定是一个入度为0到出度为0的路径,拓扑序在确定当前点之前能够考虑到所有到它的情况,所以最后取个最值即可. 代码: # ...
2017 ACM/ICPC（西安）赛后总结
早上8:00的高铁,所以不得不6点前起床,向火车站赶--到达西安后已经是中午,西工大距离西安北站大概3小时车程的距离,只好先解决午饭再赶路了--下午3.30的热身赛,一行人在3.35左右赶到了赛场,坐 ...

最新文章

热门文章