Visit the Park

最新推荐文章于 2024-09-30 17:15:14 发布

梦中醉卧惊起

最新推荐文章于 2024-09-30 17:15:14 发布

阅读量171

点赞数 1

分类专栏：省赛文章标签：图论算法 c++

原文链接：https://blog.csdn.net/m0_55682843/article/details/123510947?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522166057295316781667815283%2522%252C%2522scm%2522%253A%252220140713.130102334..%2522%257D&request_id=166057295316781667815283&biz_id=0&utm_me

版权

省赛专栏收录该内容

20 篇文章 0 订阅

订阅专栏

该博客探讨了一种无向图中基于权值和行进序列的概率计算问题。通过给出的样例解析，解释了如何计算从一个节点到另一个节点的期望值，并利用数学期望的性质来求解总期望。AC代码展示了如何实现这个算法，包括路径数量的计算和概率的求和。

摘要由CSDN通过智能技术生成

Visit the Park

题目链接
题意 :
有n个点，m条边的无向图（多路径），每条边有一个权值，后给定一个行进序列，问概率是多少？期望值的求法：字符串拼接的样子将每段路的期望值合并

思路 :

在这里插入图片描述
样例解析 :
1 – > 2 有三条路,值为1 2 2
2 – > 3 有两条路,值为1 4
1 x 1 x 10 + 1 = 11
1 x 1 x 10 + 4 = 14
1 x 2 x 10 + 1 = 21
1 x 2 x 10 + 4 = 24

数学期望等于 = 每一种结果 x 概率求和
AC代码 :


#include <iostream>
#include <map>
#include <cstring>
#include <algorithm>
#include <math.h>

#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
using namespace std;
#define ll long long
#define mod 998244853
map<int, map<int, map<int, int>>> mp;
int arr[300005];

ll ksm(ll a, ll b,ll p) {
    ll res = 1;
    while (b) {
        if (b & 1)
            res = res * a % p;
        a = a * a % p;
        b >>= 1;
    }
    return res % p;
}

ll inv(ll a, ll p){
	return ksm(a, p - 2, p);
}
 
int main(){
	IOS;
	int n, m, k;
	while(cin >> n >> m >> k){
		mp.clear();
		
		for(int i = 0; i < m; i++){
			int u, v, w;
			cin >> u >> v >> w;
			mp[u][v][w]++;
			mp[v][u][w]++;
		}
		
		for(int i = 1; i <= k; i++){
			cin >> arr[i];
		}
		ll sum = 0; //求总的预算
		ll cnt = 1; //求总路径数量 
		bool f = 0;  
		for(int i = 2; i <= k; i++){
			if(mp[arr[i - 1]][arr[i]].size() == 0){
				f = 1;
				break;
			}
			ll ans = 0;
			for(auto l : mp[arr[i - 1]][arr[i]]){
				ans += l.second;
				ans %= mod;
			}
			cnt *= ans;
			cnt %= mod; 
		}
		if(f){
			cout << "Stupid Msacywy!" << endl;
		}else{
			sum = 0;
			for(int i = 2; i <= k; i++){
				ll ans = 0;
				sum *= 10;
				sum %= mod;
				for(auto l : mp[arr[i - 1]][arr[i]]){
					ans += l.second;
					ans %= mod;
				}
				for(auto l : mp[arr[i - 1]][arr[i]]){
					sum += (((l.first * l. second % mod) * cnt % mod)* inv(ans, mod) % mod) % mod;
					sum %= mod;
				}
			}
			cnt = inv(cnt, mod);
			cout << sum * cnt % mod << endl;
		}
	}
	return 0;
}