2022蓝桥杯C/C++ B组国赛 故障
题目描述
在软件或系统开发中,我们会遇到各种各样的故障。为了从故障现象反推故障原因,工程师们会总结一种叫做相关性矩阵的二维表格,来表示故障原因与故障现象之间的关系。比如:
其中每行表示一种故障原因,每一列表示一种故障现象。该矩阵表示故障原因 A A A 可能产生故障现象 2 2 2、 3 3 3、 4 4 4,故障原因 B B B 可能产生故障现象 1 1 1、 3 3 3。
在实际开发过程中,如果出现了故障原因,工程师就可以根据故障现象,去计算每种故障原因产生的概率,并按照概率大小对故障原因进行排查,以达到快速定位故障原因的目的。
现在,我们假设系统开发中同一时间只会出现一种故障原因, 并且故障原 因引起各故障现象是独立事件。举个例子来说:
假设系统现在发生了故障原因 A A A, 有 1 3 \frac{1}{3} 31 的概率出现故障现象 2 2 2,有 1 4 \frac{1}{4} 41 的概率出现故障现象 3 3 3,有 1 2 \frac{1}{2} 21 的概率出现故障现象 4 4 4。由于 3 3 3 种现象是独立发生的,因此有 1 2 × 3 × 4 \frac{1}{2 \times 3 \times 4} 2×3×41 的概率同时出现故障 2 2 2、 3 3 3、 4 4 4。
约定若相关性矩阵中没有 x
记号, 则表示该故障原因一定不会产生某故障现象,比如故障原因
A
A
A,一定不会产生故障现象
1
1
1。根据历史经验数据,我们统计得到了每一种故障原因出现的概率以及每一种故障原因对应的故障现象产生概率。
现在已知系统出现的故障现象,求问各个故障原因发生的概率。
输入格式
第 1 1 1 行: 2 2 2 个正整数 N , M , N N, M, N N,M,N 表示故障原因的个数(编号 1 … N ) , M 1 \ldots N),M 1…N),M 表示故障现象的个数(编号 1 … M 1 \ldots M 1…M)。
第 2 2 2 行: N N N 个整数,第 i i i 个数表示故障原因 i i i 产生的概率 P i P_{i} Pi。
第 3 … N + 2 3 \ldots N+2 3…N+2 行:每行 M + 1 M+1 M+1 个整数,第 i + 1 i+1 i+1 行第 j j j 个整数 P i j P_{i j} Pij 表示故障原因 i i i 出现故障现象 j j j 的概率(百分比)。
第 N + 3 N+3 N+3 行: 1 1 1 个正整数 K K K,表示目前出现的故障现象数量。
第 N + 4 N+4 N+4 行: K K K 个正整数,依次为当前出现的故障现象编号,不会重复。
输出格式
第 1 … N 1 \ldots N 1…N 行:按概率从高到低输出每种故障原因及其可能的概率,若出现 概率相同则优先输出编号小的故障原因。第 1 1 1 个数为故障原因编号, 第 2 2 2 个数为故障概率(百分比),保留 2 2 2 位小数。
样例 #1
样例输入 #1
3 5
30 20 50
0 50 33 25 0
30 0 35 0 0
0 0 0 25 60
1
3
样例输出 #1
2 56.89
1 43.11
3 0.00
提示
对于所有测试用例, 1 ≤ N ≤ 40 , 1 ≤ M ≤ 20 , 0 ≤ P i ≤ 100 , ∑ ( P i ) = 100 1 \leq N \leq 40,1 \leq M \leq 20,0 \leq P_{i} \leq 100, \sum\left(P_{i}\right)=100 1≤N≤40,1≤M≤20,0≤Pi≤100,∑(Pi)=100, 0 ≤ P i j ≤ 100 0 \leq P_{i j} \leq 100 0≤Pij≤100。
蓝桥杯 2022 国赛 B 组 G 题。
题解
设 A i A_i Ai 为故障原因 i ( 1 ≤ i ≤ N ) i(1≤i≤N) i(1≤i≤N), B j B_j Bj 为故障现象 j ( 1 ≤ j ≤ M ) j(1≤j≤M) j(1≤j≤M)。
利用贝叶斯公式:
P ( A i ∣ B ) = P ( A i B ) P ( B ) = P ( B ∣ A i ) × P ( A i ) ∑ ( P ( A i ∣ B ) ) P(A_i|B)=\frac{P(A_iB)}{P(B)}=\frac{P(B|A_i)×P(A_i)}{\sum\left(P(A_i|B)\right)} P(Ai∣B)=P(B)P(AiB)=∑(P(Ai∣B))P(B∣Ai)×P(Ai)
B1 | B2 | B3 | B4 | B5 | |
---|---|---|---|---|---|
A1 | 0 | 50 | 33 | 25 | 0 |
A2 | 30 | 0 | 35 | 0 | 0 |
A3 | 0 | 0 | 0 | 25 | 60 |
代码
#include <iostream>
#include <algorithm>
using namespace std;
const int N=47, M=27;
double p1[N], p2[N][M]; //p1[i]为故障原因i出现的概率,p2[i][j]为故障原因i出现的情况下故障现象j出现的概率
struct ANS { int id; double res; } ans[N]; //记录最终答案的
bool vis[M]; //vis[i]标记出现的故障现象i
bool cmp(ANS x, ANS y) { return x.res!=y.res?x.res>y.res:x.id<y.id; }
int main()
{
int n, m, k, x;
scanf("%d%d", &n, &m);
for(int i=1; i<=n; ++i)
scanf("%lf", &p1[i]);
for(int i=1; i<=n; ++i)
for(int j=1; j<=m; ++j)
scanf("%lf", &p2[i][j]);
scanf("%d", &k);
while(k--) scanf("%d", &x), vis[x]=1; //标记出现过的现象x
double sum=0; //sum作公式中的分母
for(int i=1; i<=n; ++i)
{
for(int j=1; j<=m; ++j)
vis[j]?p1[i]*=p2[i][j]/100:p1[i]=p1[i]*(100-p2[i][j])/100; //出现的故障现象概率计算使用p2[i][j],未出现的故障现象的概率计算使用(100-p2[i][j])
sum+=p1[i]; //计算分母
}
for(int i=1; i<=n; ++i)
ans[i].id=i, ans[i].res=p1[i]/sum*100; //每种故障原因记录对应的序号和概率结果(便于后续排序)
sort(ans+1, ans+1+n, cmp); //结构体排序
for(int i=1; i<=n; ++i)
printf("%d %.2lf\n", ans[i].id, ans[i].res);
return 0;
}