Codeforces 1265E Beautiful Mirrors（概率DP，逆元）

Yuhan の Blog

于 2019-12-09 15:57:59 发布

阅读量279

点赞数 2

分类专栏： Codeforces # 数学 # 动态规划文章标签： dp cf 逆元

本文链接：https://blog.csdn.net/qq_45228537/article/details/103459328

版权

Codeforces 同时被 3 个专栏收录

104 篇文章

订阅专栏

数学

34 篇文章

订阅专栏

动态规划

31 篇文章

订阅专栏

本文深入解析了使用动态规划(DP)算法解决概率问题的思路，特别是如何利用逆元处理取模运算，确保计算精度。通过具体案例，详细介绍了DP状态转移方程的构建过程，并提供了C++代码实现，展示了如何高效地求解通过一系列镜子的期望天数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

思路：

1.记dp[i]为通过第i面镜子的期望天数，pi为通过第i面镜子的概率，则我们可以得到dp[i]=dp[i-1]+pi*1+(1-pi)*(dp[i]+1)，其中的dp[i-1]为走到第i-1的期望天数，此时我们有pi概率只需要再走一天，有1-pi概率需要再走1+dp[i]天（往前走一天，失败，然后再走通过i的期望天数)；
2.除法取模需要用到逆元；

代码：

#include<bits/stdc++.h>
using namespace std;
#define pt(a) cerr<<a<<"---\n"
typedef long long LL;
void extgcd(LL a,LL b,LL& x,LL& y){
    if(b==0){x=1,y=0;return;}
    extgcd(b,a%b,x,y);
    LL t=x;x=y;y=t-(a/b)*y;
}
LL mod_inv(LL a,LL m){LL x,y;extgcd(a,m,x,y);return (m+x%m)%m;}
const int maxn=2e5+99;
const LL M=998244353;
int n,p[maxn];
LL dp[maxn],inv[105];
void solve(){
    for(int i=1;i<=n;i++){
        if(inv[p[i]]==0)inv[p[i]]=mod_inv(p[i],M);
        dp[i]=(dp[i-1]+1)*100ll%M*inv[p[i]]%M;
    }
    cout<<dp[n];
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(NULL);
	cin>>n;
	for(int i=1;i<=n;i++)cin>>p[i];
	solve();
	return 0;
}