本题我用一个更直观的方法来求解,与网上的做法不同。
提示:
1. 在每个城市中爆炸的概率和另外一个量有关,也就是说我们尝试间接求这个概率,容易许多。
2. 这个量我们经常在概率问题中碰到,可以考虑猜猜我们经常求的那些……
注意: 输出带上
fabs()
, 因为可能出现
−0.00000
代码后我将详细说明:
#include <bits/stdc++.h>
using namespace std;
const int maxn = 310;
int n , m ;
double p , q;
int g[maxn][maxn] , e[maxn];
double a[maxn][maxn];
void gaussJohdan()
{
for(int i=1;i<=n;i++)
{
int r = i;
for(int j=i+1;j<=n;j++) if(fabs(a[j][i]) > fabs(a[r][i])) r = j;
if(r != i) for(int j=1;j<=n+1;j++) swap(a[i][j] , a[r][j]);
for(int k=1;k<=n;k++) if(k != i) for(int j=n+1;j>=i;j--) a[k][j] -= a[k][i] / a[i][i] * a[i][j];
}
}
int main(int argc, char *argv[]) {
cin>>n>>m>>p>>q;
p /= q; p = 1-p;
while(m--)
{
int x , y;
scanf("%d%d" , &x , &y);
g[x][y] = g[y][x] = 1;
e[x]++; e[y]++;
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++) if(g[i][j]) a[i][j] = p / e[j];
a[i][i] --;
}
a[1][n+1]--;
gaussJohdan();
double all = 0;
for(int i=1;i<=n;i++) all += a[i][n+1]/=a[i][i];
for(int i=1;i<=n;i++) printf("%.9lf\n" , fabs(a[i][n+1]/all));
return 0;
}
我们求出每个城市的期望经过次数
Si
, 那么第
i
个城市的答案就是:
这个我想并不难理解,因为每经过一次爆炸的可能都是相同的,所以每个城市爆炸的概率就与每个城市的期望经过次数正相关。然而求期望经过次数我想难不到读者吧。