正题

题目链接:https://www.luogu.com.cn/problem/CF446D

题目大意

给出nnn个点mmm条边的一张无向图，一些点有陷阱，走到时会损失一条生命，总共有kkk条生命，求从111出发随机游走到nnn没有死亡且到终点时仅剩一条命的概率。

1≤n≤500,1≤m≤105,2≤k≤1091\leq n\leq 500,1\leq m\leq 10^5,2\leq k\leq 10^91≤n≤500,1≤m≤105,2≤k≤109

陷阱点个数不超过100100100。

解题思路

这个kkk很大，这个陷阱点个数又很少，我们可以考虑矩阵乘法，预处理ai,ja_{i,j}ai,j表示陷阱点iii走到陷阱点jjj且中间没有走陷阱点的概率，然后矩阵乘法转移即可。

但是现在的问题是我们如何快速预处理出ai,ja_{i,j}ai,j，可以考虑枚举终点xxx那么有fx=1f_x=1fx=1，然后其他的陷阱点处fx=0f_x=0fx=0，一般的点处fx=1degx∑x→yfyf_x=\frac{1}{deg_x}\sum_{x\rightarrow y}f_{y}fx=degx1∑x→yfy，这样我们就能对于每个起点预处理出gi,jg_{i,j}gi,j表示从无陷阱的节点iii走到陷阱点jjj且中间没有其他陷阱点的概率。

之后我们枚举起点陷阱点的出边就可以预处理出aaa了，因为上面的过程要用到高斯消元，所以这样的复杂度是O(n4)O(n^4)O(n4)的，无法通过本题。

不难注意到上面的消元中，我们只有陷阱点处的常数（且陷阱点处仅有常数）发生了变化，所以我们可以直接高斯消出每个非陷阱点和所有陷阱点的关系式，然后直接带入常数即可。

记陷阱点个数为rrr，时间复杂度O((n+r)3+r3log⁡k)O((n+r)^3+r^3\log k)O((n+r)3+r3logk)

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
const int N=610,S=110;
const double eps=1e-8;
struct matrix{double a[S][S];
}ans,m,c;
matrix operator*(const matrix &a,const matrix &b){memset(c.a,0,sizeof(c.a));for(int i=0;i<S;i++)for(int j=0;j<S;j++)for(int k=0;k<S;k++)c.a[i][j]+=a.a[i][k]*b.a[k][j];return c;
}
int n,h,k,deg[N],a[N][N];
double f[N][N];bool v[N];
vector<int> q;
int main()
{scanf("%d%d%d",&n,&h,&k);for(int i=1;i<=n;i++){scanf("%d",&v[i]);if(v[i]){q.push_back(i);f[i][n+q.size()]=f[i][i]=1;}}int r=n+q.size();for(int i=1,x,y;i<=h;i++){scanf("%d%d",&x,&y);a[x][y]++;a[y][x]++;deg[x]++;deg[y]++;}for(int i=1;i<=n;i++){if(v[i])continue;for(int j=1;j<=n;j++)f[i][j]=-1.0*a[i][j]/(double)deg[i];f[i][i]=1;}for(int i=1;i<=n;i++){for(int j=i;j<=n;j++)if(fabs(f[i][j])>eps){swap(f[i],f[j]);break;}double d=f[i][i];for(int j=i;j<=r;j++)f[i][j]=f[i][j]/d;for(int j=1;j<=n;j++){if(i==j)continue;double rate=-f[j][i]/f[i][i];for(int k=i;k<=r;k++)f[j][k]+=rate*f[i][k];}}for(int i=0;i<q.size();i++)ans.a[0][i]=f[1][n+i+1];for(int i=0;i<q.size();i++){for(int j=0;j<q.size();j++){for(int k=1;k<=n;k++)m.a[i][j]+=a[q[i]][k]*f[k][n+j+1];m.a[i][j]/=(double)deg[q[i]];}}k-=2;while(k){if(k&1)ans=ans*m;m=m*m;k>>=1;}printf("%.10lf\n",ans.a[0][q.size()-1]);return 0;
}

CF446D-DZY Loves Games【高斯消元,矩阵乘法】相关推荐

CF446D. DZY Loves Games
题目大意: 有一个连通的迷宫,你可以理解为一个 $n$ 个点 $m$ 条边的无向连通图.有些点可能有陷阱,保证 $1$ 号点没有陷阱而 $n$ 号点存在陷阱.一开始你有 $k$ 条命,每进入一个有 ...
BZOJ_1778_[Usaco2010_Hol]_Dotp_驱逐猪猡_(期望动态规划+高斯消元+矩阵)
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1778 炸弹从1出发,有$\frac{P}{Q}$的概率爆炸,如果不爆炸,等概率移动到连通的 ...
Covering（dfs打表+高斯消元+矩阵快速幂）
传送门这是个递推题,可以看下维基的各种解释:传送门假设递推式是一个四元方程,至于为什么也不大清楚(可能是操场是4*n的缘故) 然后得到递推式f(n)=f(n-1)+5*f(n-2)+f(n-3)- ...
矩阵树 Matrix-Tree 定理实现模板(高斯消元求解行列式)
大佬1博客:https://www.cnblogs.com/zj75211/p/8039443.html 大佬2博客:https://www.cnblogs.com/yangsongyi/p/1069 ...
计算矩阵的逆和行列式的值(高斯消元+LU分解)
计算矩阵的逆选主元的高斯消元法朴素的高斯消元法是将矩阵A和单位矩阵放在一起,通过行操作(或者列操作)将A变为单位矩阵,这个时候单位矩阵就是矩阵A的逆矩阵.从上到下将A变为上三角矩阵的复杂度为O(n ...
bzoj 3503: [Cqoi2014]和谐矩阵（高斯消元）
3503: [Cqoi2014]和谐矩阵 Time Limit: 10 Sec Memory Limit: 128 MBSec Special Judge Submit: 1101 Solved ...
【bzoj1778】[Usaco2010 Hol]Dotp 驱逐猪猡矩阵乘法+概率dp+高斯消元
题目描述奶牛们建立了一个随机化的臭气炸弹来驱逐猪猡.猪猡的文明包含1到N (2 <= N <= 300)一共N个猪城.这些城市由M (1 <= M <= 44,850)条由两 ...
HDU 4305 Lightning (高斯消元解kirchhoff矩阵+逆元)
题意是:给一些坐标点,如果两点之间的距离小于R,并且两点之间没有其他点,则这两个点保持连通,这样构成了一个图.问这个图中生成树的个数. 因为数据量并不大,O(N^3)的建图没有问题. 建好图以后就可以 ...
POJ 1830 开关问题【01矩阵高斯消元】
任意门:http://poj.org/problem?id=1830 开关问题 Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 1 ...

CF446D-DZY Loves Games【高斯消元,矩阵乘法】

正题

题目大意

解题思路

code

CF446D-DZY Loves Games【高斯消元,矩阵乘法】相关推荐

最新文章

热门文章