2020百度之星复赛 Battle for Wosneth(期望+逆元)附:快速幂求逆元

题目链接: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)

分析题意:

  1. A打中B的期望为p,则A有m/p次命中B
  2. 完整的回合(A–>B,B–>A)共 m / p − 1 m/p-1 m/p1次。(最后一次B已经牺牲,无法攻击A,构不成一个完整的回合)
  3. 每一个完整的回合A的生命期望为 p − q p-q pq,则 m / p − 1 m/p-1 m/p1个完整的回合,A的生命期望为 ( m / p − 1 ) ∗ ( p − q ) (m/p-1)*(p-q) (m/p1)(pq)
  4. 最后一次A打死B,A的生命期望为p

综上分析:
A的生命总期望为: ( m / p − 1 ) ∗ ( p − q ) + p (m/p-1)*(p-q)+p (m/p1)(pq)+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=pinv(100)q=qinv(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) (mpinv(100))(pq)inv(p)+pinv(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;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

zaiyang遇见

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值