其他分享
首页 > 其他分享> > P3317 [SDOI2014]重建

P3317 [SDOI2014]重建

作者:互联网

\(\color{#0066ff}{ 题目描述 }\)

T国有N个城市,用若干双向道路连接。一对城市之间至多存在一条道路。

在一次洪水之后,一些道路受损无法通行。虽然已经有人开始调查道路的损毁情况,但直到现在几乎没有消息传回。

幸运的是,此前T国政府调查过每条道路的强度,现在他们希望只利用这些信息估计灾情。具体地,给定每条道路在洪水后仍能通行的概率,请计算仍能通行的道路恰有N-1条,且能联通所有城市的概率。

\(\color{#0066ff}{输入格式}\)

输入的第一行包含整数N。

接下来N行,每行N个实数,第i+l行,列的数G[i][j]表示城市i与j之间仍有道路联通的概率。

输入保证G[i][j]=G[j][i],且G[i][i]=0;G[i][j]至多包含两位小数。

\(\color{#0066ff}{输出格式}\)

输出一个任意位数的实数表示答案。

你的答案与标准答案相对误差不超过10^(-4)即视为正确。

\(\color{#0066ff}{输入样例}\)

3
0 0.5 0.5
0.5 0 0.5
0.5 0.5 0

\(\color{#0066ff}{输出样例}\)

0.375

\(\color{#0066ff}{数据范围与提示}\)

1 < N < =50

数据保证答案非零时,答案不小于10^-4

\(\color{#0066ff}{题解}\)

根据题目,我们要求的就是

\[ ans=\sum_{E}\prod_{k\in E}P_k\prod_{k\notin E} (1-P_k) \]

如果没有后面那个东西,显然就是裸的矩阵树定理,但是后面的东西很不好处理,尤其是因为\(k\notin E\)

那么,考虑容斥一下, 把\(\notin换成\in\)

\[ ans=\sum_{E}\prod_{k\in E}P_k\frac{\prod_{k}(1-P_k)}{\prod_{k\in E} (1-P_k)} \]

然后把上面提出来,就成这样了

\[ ans=\prod_{k}(1-P_k)\sum_{E}\prod_{k\in E}\frac{P_k}{(1-P_k)} \]

这。。。。这是新的边权!!可以矩阵树直接做!

然后把前面的累乘处理一下即可

矩阵树第一题

有两点需要注意

答案是矩阵的余子式的值,也就是矩阵去掉任一行任一列的行列式的值

度数矩阵的变化

#include<bits/stdc++.h>
#define LL long long
LL in() {
    char ch; LL x = 0, f = 1;
    while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
    for(x = ch ^ 48; isdigit(ch = getchar()); x = (x << 1) + (x << 3) + (ch ^ 48));
    return x * f;
}
const int maxn = 55;
const double eps = 1e-8;
double ans = 1, mp[maxn][maxn];
int n;
void gauss() {
    for(int i = 1; i < n; i++) {
        for(int j = i + 1; j < n; j++) {
            double now = mp[j][i] / mp[i][i];
            for(int k = i; k < n; k++) mp[j][k] -= mp[i][k] * now;
        }
        ans *= mp[i][i];
    }
    ans = fabs(ans);
}
int main() {
    n = in();
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= n; j++) 
            scanf("%lf", &mp[i][j]);
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= n; j++) {
            if(fabs(mp[i][j]) <= eps) mp[i][j] = eps;
            if(fabs(1 - mp[i][j]) <= eps) mp[i][j] = 1 - eps;
            if(i < j) ans *= (1.0 - mp[i][j]);
            mp[i][j] = mp[i][j] / (1.0 - mp[i][j]);
        }
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= n; j++)
            if(i ^ j) mp[i][i] += mp[i][j], mp[i][j] = -mp[i][j];
    gauss();
    printf("%.5f", ans);
    return 0;
}

标签:ch,color,矩阵,0.5,0066ff,P3317,SDOI2014,prod,重建
来源: https://www.cnblogs.com/olinr/p/10422879.html