XOR和路径 (HYSBZ - 2337 ,高斯消元解后效性 DP)

2024-03-30 13:48

本文主要是介绍XOR和路径 (HYSBZ - 2337 ,高斯消元解后效性 DP),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一.题目链接:

HYSBZ-2337

二.题目大意:

给一张无向边权图,在每个节点都会等概率地选择一条边,求 1 ~ n 路径的权值异或和的期望值.

三.分析:

由于是异或,不妨按答案的二进制位逐位考虑.

假设当前考虑第 i 位

设 dp[u] 表示 u ~ n 路径的权值异或和二进制第 i 位的期望值.

设 v 是与顶点 u 相关联的顶点集合,de[u] 表示 u 的度, wi(u, v) 表示 u 与 v 之间边的二进制第 i 位.

可得:dp[u] = \frac{\sum_{v_{j} \in v , \; w_{i}(u, v_{j}) = 1}(1 - dp[v_{j}]) + \sum_{v_j \in v , \; w_{i}(u, v_{j}) = 0} dp[v_{j}]}{de[u]}

剩下的高斯消元套上模板就好啦.

最后答案等于 \sum 2^{i}dp[1]

四.代码实现:

#include <bits/stdc++.h>
using namespace std;const int M = (int)1e4;
const int N = (int)1e2;int cnt;
int head[N + 5];
struct node
{int v, w, nx;
}Edge[M * 2 + 5];int de[N + 5];
double a[N + 5][N + 5];void init(int n)
{cnt = 0;for(int i = 1; i <= n; ++i){head[i] = -1;}
}void add(int u, int v, int w)
{Edge[cnt].v = v;Edge[cnt].w = w;Edge[cnt].nx = head[u];head[u] = cnt++;
}void Gauss(int n)
{for(int i=1;i<=n;i++){int r=i;for(int j=i+1;j<=n;j++) if(fabs(a[r][i])<fabs(a[j][i])) r=i;if(r!=i) for(int j=1;j<=n+1;j++) std::swap(a[r][j],a[i][j]);double t=a[i][i];for(int j=i+1;j<=n+1;j++) a[i][j]/=t;for(int j=1;j<=n;j++)if(i!=j){double t=a[j][i];for(int k=1;k<=n+1;k++) a[j][k]-=t*a[i][k];}}
}double work(int n)
{double ans = 0.0;for(int i = 0; i < 30; ++i){memset(a, 0, sizeof(a));for(int u = 1; u < n; ++u){a[u][u] = 1.0;for(int j = head[u]; ~j; j = Edge[j].nx){int v = Edge[j].v;int w = Edge[j].w;if((w>>i) & 1)a[u][v] += 1.0 / de[u], a[u][n + 1] += 1.0 / de[u];elsea[u][v] -= 1.0 / de[u];}}a[n][n] = 1.0;Gauss(n);ans += (1<<i) * a[1][n + 1];}return ans;
}int main()
{int n, m;scanf("%d %d", &n, &m);init(n);for(int i = 0, u, v, w; i < m; ++i){scanf("%d %d %d", &u, &v, &w);add(u, v, w), de[u]++;if(u != v)  add(v, u, w), de[v]++;}printf("%.3f\n", work(n));return 0;
}

 

这篇关于XOR和路径 (HYSBZ - 2337 ,高斯消元解后效性 DP)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/861653

相关文章

hdu4826(三维DP)

这是一个百度之星的资格赛第四题 题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1004&cid=500 题意:从左上角的点到右上角的点,每个点只能走一遍,走的方向有三个:向上,向下,向右,求最大值。 咋一看像搜索题,先暴搜,TLE,然后剪枝,还是TLE.然后我就改方法,用DP来做,这题和普通dp相比,多个个向上

hdu1011(背包树形DP)

没有完全理解这题, m个人,攻打一个map,map的入口是1,在攻打某个结点之前要先攻打其他一个结点 dp[i][j]表示m个人攻打以第i个结点为根节点的子树得到的最优解 状态转移dp[i][ j ] = max(dp[i][j], dp[i][k]+dp[t][j-k]),其中t是i结点的子节点 代码如下: #include<iostream>#include<algorithm

hdu2544(单源最短路径)

模板题: //题意:求1到n的最短路径,模板题#include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#i

hdu4865(概率DP)

题意:已知前一天和今天的天气概率,某天的天气概率和叶子的潮湿程度的概率,n天叶子的湿度,求n天最有可能的天气情况。 思路:概率DP,dp[i][j]表示第i天天气为j的概率,状态转移如下:dp[i][j] = max(dp[i][j, dp[i-1][k]*table2[k][j]*table1[j][col] )  代码如下: #include <stdio.h>#include

usaco 1.1 Broken Necklace(DP)

直接上代码 接触的第一道dp ps.大概的思路就是 先从左往右用一个数组在每个点记下蓝或黑的个数 再从右到左算一遍 最后取出最大的即可 核心语句在于: 如果 str[i] = 'r'  ,   rl[i]=rl[i-1]+1, bl[i]=0 如果 str[i] = 'b' ,  bl[i]=bl[i-1]+1, rl[i]=0 如果 str[i] = 'w',  bl[i]=b

poj 1734 (floyd求最小环并打印路径)

题意: 求图中的一个最小环,并打印路径。 解析: ans 保存最小环长度。 一直wa,最后终于找到原因,inf开太大爆掉了。。。 虽然0x3f3f3f3f用memset好用,但是还是有局限性。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#incl

uva 10154 DP 叠乌龟

题意: 给你几只乌龟,每只乌龟有自身的重量和力量。 每只乌龟的力量可以承受自身体重和在其上的几只乌龟的体重和内。 问最多能叠放几只乌龟。 解析: 先将乌龟按力量从小到大排列。 然后dp的时候从前往后叠,状态转移方程: dp[i][j] = dp[i - 1][j];if (dp[i - 1][j - 1] != inf && dp[i - 1][j - 1] <= t[i]

uva 10118 dP

题意: 给4列篮子,每次从某一列开始无放回拿蜡烛放入篮子里,并且篮子最多只能放5支蜡烛,数字代表蜡烛的颜色。 当拿出当前颜色的蜡烛在篮子里存在时,猪脚可以把蜡烛带回家。 问最多拿多少只蜡烛。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cs

uva 10069 DP + 大数加法

代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <queue>#include <map>#include <cl

uva 10029 HASH + DP

题意: 给一个字典,里面有好多单词。单词可以由增加、删除、变换,变成另一个单词,问能变换的最长单词长度。 解析: HASH+dp 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc