题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6838
Problem Description
你在打游戏的时候碰到了如下问题:
有两个人记作Alice和Bob,Bob的生命值为m,Alice的生命值很高,所以可以认为是无限的。两个人的攻击命中率分别为p%,q%。两个人轮流攻击对方。从Alice开始攻击,每次攻击的时候,如果Alice命中,那么能让对方的生命值减低1,同时自己的生命值能恢复1,如果Bob命中,那么能让对方的生命值减低1,注意Bob不会自己回血。
直到Bob的血量变为0,游戏结束。Alice想知道,游戏结束的时候,自己期望生命值变化是多少,对998244353取模。
注意这里的变化量不是绝对值,也就是如果50%的概率加一,50%的概率减一,那么期望的变化量就是0。
对于一个分数a/b,其中gcd(a,b)=1,那么我们认为这个分数对998244353取模的值为一个数c(0≤c<998244353)满足bc≡a(mod998244353)。
Input
第一行一个正整数T(1≤T≤104)表示数据组数。
对于每组数据,第一行三个整数m,p,q(1≤m≤109,1≤p,q≤100)。
Output
每组测试数据,输出一个数,表示答案。
Sample Input
2
4 100 100
1 50 50
Sample Output
1
499122177
Hint
第一组数据中,每次都能命中,所以Alice能恢复 4 点生命值,减低 3 点生命值,变化量必定为 1。
第二组数据中,对应的分数为 1/2,在Alice命中Bob之前,Bob能期望命中Alice 1/2 次。
快速幂求逆元的思想:
a / b ≡ a * x (mod n)
两边同乘b可得 a ≡ a * b * x (mod n)
即:
b * x ≡ 1 (mod n)
由费马定理可知,当n为质数时b ^ (n - 1) ≡ 1 (mod n)
拆一个b出来可得 b * b ^ (n - 2) ≡ 1 (mod n)
故当n为质数时,b的乘法逆元 x = b ^ (n - 2)
分析题意:
- A打中B的期望为p,则A有m/p次命中B
- 完整的回合(A–>B,B–>A)共 m / p − 1 m/p-1 m/p−1次。(最后一次B已经牺牲,无法攻击A,构不成一个完整的回合)
- 每一个完整的回合A的生命期望为 p − q p-q p−q,则 m / p − 1 m/p-1 m/p−1个完整的回合,A的生命期望为 ( m / p − 1 ) ∗ ( p − q ) (m/p-1)*(p-q) (m/p−1)∗(p−q)
- 最后一次A打死B,A的生命期望为p
综上分析:
A的生命总期望为:
(
m
/
p
−
1
)
∗
(
p
−
q
)
+
p
(m/p-1)*(p-q)+p
(m/p−1)∗(p−q)+p
p,q为概率,题目中输入的是没有除以100。实际p=p/100,q=q/100。最后结果还有对mod取余。
将
p
=
p
∗
i
n
v
(
100
)
,
q
=
q
∗
i
n
v
(
100
)
p=p*inv(100),q=q*inv(100)
p=p∗inv(100),q=q∗inv(100)代入生命总期望的式子,得:
( m − p ∗ i n v ( 100 ) ) ∗ ( p − q ) ∗ i n v ( p ) + p ∗ i n v ( 100 ) (m-p*inv(100))*(p-q)*inv(p)+p*inv(100) (m−p∗inv(100))∗(p−q)∗inv(p)+p∗inv(100)
完整代码:
#include<cstdio>
#include<cstring>
#include<math.h>
#include<stdlib.h>
#include<iostream>
#include<queue>
#include<algorithm>
using namespace std;
typedef long long LL;
#define mod 998244353
LL quick_mi(LL a,LL n)
{
LL res=1;
while(n)
{
if(n&1)
res=res*a%mod;
a=a*a%mod;
n>>=1;
}
return res%mod;
}
LL m,p,q;
LL inv(LL x)
{
return quick_mi(x,mod-2);
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
scanf("%lld%lld%lld",&m,&p,&q);
LL ans;
ans=((((m-p*inv(100))*(p-q)%mod)*inv(p)%mod+p*inv(100)%mod)%mod+mod)%mod;
printf("%lld\n",ans);
}
return 0;
}