BZOJ1778: [Usaco2010 Hol]Dotp 驅逐豬玀
阿新 • • 發佈:2018-11-06
題意
給你一張 個點 條邊的無向圖,最開始有一顆炸彈在一號節點,它有 的概率爆炸,如果沒有爆炸,它會等概率的移動到另一個與當前節點相連的點,問炸彈分別在每個點爆炸的概率;
題解
考慮一維的向量矩陣 ,第 位表示炸彈停留在第 個點的概率為 ,那麼初始時炸彈停在 號節點上 ;考慮轉移,設 表示 點的度數,那麼有 ( 為 能到達的點);那麼我們建立一個 的矩陣 , ,那麼最終的答案矩陣 ;用無窮等比數列求和公式得到 ( 為單位矩陣);推到這裡其實有個很巧的地方,那就是 和 都是一維的向量矩陣,有 , 能直接算,再模擬一下矩陣乘法的過程,發現 其實是可以用高斯消元解出來的,最後再注意一下精度問題就行了;
#include<bits/stdc++.h>
#define Fst first
#define Snd second
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef long long LL;
typedef unsigned int UI;
typedef unsigned long long ULL;
template<typename T> inline void read(T& x) {
char c = getchar();
bool f = false;
for (x = 0; !isdigit(c); c = getchar()) {
if (c == '-') {
f = true;
}
}
for (; isdigit(c); c = getchar()) {
x = x * 10 + c - '0';
}
if (f) {
x = -x;
}
}
template<typename T, typename... U> inline void read(T& x, U& ... y) {
read(x), read(y...);
}
const int N=3e2+10;
const double eps=1e-15;
int n,m,p;
int In[N],head[N];
double P,Q;
double A[N][N];
struct Edge {
int to,last;
Edge () {}
Edge (int a,int b) :to(a),last(b) {}
}edge[100010];
void ADD(int a,int b) {
edge[++p]=Edge(b,head[a]); head[a]=p;
edge[++p]=Edge(a,head[b]); head[b]=p;
}
void Gauss() {
for(int i=0;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=0;j<=n;++j) swap(A[r][j],A[i][j]);
for(int j=i+1;j<n;++j)
for(int k=n;k>=i;--k)
A[j][k]-=A[j][i]/A[i][i]*A[i][k];
}
for(int i=n-1;~i;--i) {
for(int j=i+1;j<n;++j)
A[i][n]-=A[i][j]*A[j][n];
A[i][n]/=A[i][i];
}
}
int main() {
read(n,m); scanf("%lf%lf",&P,&Q); P/=Q;
for(int i=1;i<=m;++i) {
int u,v; read(u,v); --u; --v;
++In[u]; ++In[v];
ADD(u,v);
}
A[0][n]=1;
for(int u=0;u<n;++u) {
A[u][u]=1;
for(int i=head[u];i;i=edge[i].last) {
int v=edge[i].to;
A[u][v]=-(1.0-P)/In[v];
}
}
Gauss();
for(int i=0;i<n;++i) printf("%.9lf\n",A[i][n]*P+eps);
return 0;
}