bzoj3270 博物馆(概率dp+高斯消元)

设dp[i][j]表示第一个人在i,第二个人在j的概率。
我们有转移:
dp[i][j]=dp[i][j]p[i]p[j]+kmp[k][i]p[j](1p[k])/du[k]+kmp[k][j]p[i](1p[k])/du[k]+k1k2mp[k1][i]mp[k2][j](1p[k1])(1p[k2])/du[k1]/du[k2]
特别的,我们还有dp[a][b]=…+1.
转移带环,高斯消元来解这 n2 个方程。
复杂度 O(n6)

注意转移的时候,dp[i][i]是不能再转移出去的!因为这已经是结束状态了!所以正着写大概会好写一些。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 25
#define eps 1e-9
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,aa,b,mp[N][N],du[N];
double p[N],a[N*N][N*N];
inline void Gauss(){
    for(int i=1;i<=m;++i){
        int r=i;
        for(int j=i+1;j<=m;++j) if(abs(a[j][i])>abs(a[r][i])) r=j;
        if(r!=i) for(int j=i;j<=m+1;++j) swap(a[i][j],a[r][j]);
        for(int j=i+1;j<=m;++j){
            double t=a[j][i]/a[i][i];
            for(int k=i;k<=m+1;++k) a[j][k]-=a[i][k]*t;
        }
    }for(int i=m;i>=1;--i){
        for(int j=i+1;j<=m;++j) a[i][m+1]-=a[i][j]*a[j][m+1];
        a[i][m+1]/=a[i][i];
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();aa=read();b=read();
    while(m--){
        int x=read(),y=read();
        mp[x][y]=1;mp[y][x]=1;du[x]++;du[y]++;
    }m=n*n;for(int i=1;i<=n;++i) scanf("%lf",&p[i]);
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j){
            int x=(i-1)*n+j;a[x][x]=-1;
            if(i==j) continue;a[x][x]+=p[i]*p[j];
            for(int k=1;k<=n;++k){
                if(mp[i][k]) a[(k-1)*n+j][x]+=p[j]*(1-p[i])/du[i];
                if(mp[j][k]) a[(i-1)*n+k][x]+=p[i]*(1-p[j])/du[j];
            }for(int k1=1;k1<=n;++k1){
                if(!mp[i][k1]) continue;
                for(int k2=1;k2<=n;++k2){
                    if(!mp[j][k2]) continue;
                    a[(k1-1)*n+k2][x]=(1-p[i])*(1-p[j])/du[i]/du[j];
                }
            }
        }a[(aa-1)*n+b][m+1]=-1;Gauss();
    for(int i=1;i<=n;++i) printf("%.6lf ",a[(i-1)*n+i][m+1]);
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值