[gym101981F]Frank

在本题中，每一步是独立的，因此即可以看作从$s$移动到$t$的期望步数（对于每一对$s$和$t$都求出答案）

令$f_{i,j}$表示当$s=i$且$t=j$时的答案，则有$f_{i,j}=\begin{cases}\sum_{(i,k)\in n}w_{(i,k)}f_{k,j}+1&(i\ne j)\\0&(i=j)\end{cases}$

事实上，对于每一个$j$，其仅有在$i=j$时的方程不同，我们需要利用这个公共部分

更具体的，令$g_{i}=\sum_{(i,j)\in n}w_{(i,j)}f_{j,i}+1$，即$s=t=i$且强制第一次不能结束的期望步数

将其简单变形，即有$f_{i,i}=\sum_{(i,j)\in n}w_{(i,j)}f_{j,i}+1-g_{i,i}=0$（即恰好为0）

构造矩阵$M_{i,j}=w_{(i,j)}$，$F_{i,j}=f_{i,j}$，$G_{i,j}=\begin{cases}g_{i}&(i=j)\\0&(i\ne j)\end{cases}$以及全1矩阵$J_{i,j}=1$，那么根据前面的递推式，就可以得到$MF+J-G=F$

进而，即$(M-I)F=G-J$，问题即变为求$G$以及$\frac{G-J}{M-I}$

$G$的计算

再构造一个$1\times |V|$的矩阵$\pi$，满足$\sum_{i=1}^{|V|}\pi_{1,i}=1$且$\pi M=\pi$

结论5：$\forall 1\le i\le n,g_{i}\pi_{1,i}=1$

考虑条件$MF+J-G=F$，两边同时左乘$\pi$即$\pi MF+\pi J-\pi G=\pi F$

注意到$\pi M=\pi$，化简后即有$\pi J=\pi G$

由此，即$(\pi G)_{i}=g_{i}\pi_{1,i}$，而$(\pi J)_{i}=\sum_{i=1}^{|V|}\pi_{1,i}=1$，即得证

对于$\pi$的计算，相当于解一个$|V|$元的方程，使用高斯消元，其中$\pi M=\pi$中的$|V|$个方程中是线性相关的（所有方程之和为0），删去其中任意一个即可

$\frac{G-J}{M-I}$的计算

这四个矩阵通过前面，都可以在$o(|V|^{3})$内求出，但$M-I$并不一定满秩，因此不能用前面矩阵求逆的手段来计算

结论6（矩阵树定理）：对于一张带权有向图$G$，令$M_{i,j}=\begin{cases}-e_{(i,j)}&(i\ne j)\\\sum_{k=1}^{n}e_{(i,k)}&(i=j)\end{cases}$，则矩阵$M$中去掉第$i$行和第$i$列，得到矩阵$M'$的行列式即所有以$i$为根的外向树边权乘积和

结论7：矩阵$M-I$的秩为$|V|-1$

将矩阵$M-I$的每一列求和，其结果都恰好为0，即$M-I$的秩小于$|V|$

显然$I-M$的秩与$M-I$相同，根据边权和为1的性质，其恰为$G$的矩阵$M$

根据矩阵树定理，去掉$M-I$的第$i$行和第$i$列，此时其行列式为以$i$为根的外向树边权乘积和，根据其强连通以及边权$\in R^{+}$，不难证明其大于0

因此，去掉第$i$行第$i$列后其满秩，因此$M-I$的秩大于等于$|V|-1$

综上，即证明了$M-I$的秩为$|V|-1$

令$A=G-J$，$B=M-I$，即求$AF=B$

考虑令$A$的第$i$行乘上一个非0数或减去$A$的第$j$行，此时不难发现对于$B$的影响（在$F$不变下）也恰恰是让$B$的第$i$行乘上一个非0数或减去第$j$行

由此进行高斯消元，可以使得$A$中第$i$行仅有第$i$和$|V|$个元素非0（第$|V|$行全为0）

注意到$F_{i,i}=0$，因此有$B_{i,j}=\sum_{k=1}^{|V|}A_{i,k}F_{k,j}$，对$i$分类讨论：

1.若$1\le i<|V|$，即$B_{i,j}=F_{i,j}+A_{i,n}F_{n,j}$

2.若$i=|V|$，由于$A$的第$|V|$行为0，因此$B_{i,j}=0$，且与$F$的值无关

通过$i=j$时，可以解出$F_{n,i}=\frac{B_{i,i}}{A_{i,n}}$，进而即可解出$F_{i,j}=B_{i,j}-A_{i,n}F_{n,j}$

综上，我们就得到了一个时间复杂度为$o(|V|^{3})$的做法

（实现中可能因为一些精度误差而挂掉了QAQ）

#include<bits/stdc++.h>
using namespace std;
#define N 405
#define eps 1e-8
int n,m,q,x,y,z,deg[N];
double ans,M[N][N],A[N][N],B[N][N],f[N][N];
void guess1(){
    for(int i=1;i<=n;i++){
        int k=-1;
        for(int j=i;j<=n;j++)
            if (fabs(A[j][i])>=eps){
                k=j;
                break;
            }
        if (k!=i){
            for(int j=i;j<=n+1;j++)swap(A[i][j],A[k][j]);
        }
        double x=A[i][i];
        for(int j=i;j<=n+1;j++)A[i][j]/=x;
        for(int j=i+1;j<=n;j++)
            if (fabs(A[j][i])>=eps){
                double x=A[j][i];
                for(int k=i;k<=n+1;k++)A[j][k]-=A[i][k]*x;
            }
    }
    for(int i=n;i;i--)
        for(int j=i+1;j<=n;j++){
            A[i][n+1]-=A[i][j]*A[j][n+1];
            A[i][j]=0;
        }
}
void guess2(){
    for(int i=1;i<n;i++){
        int k=-1;
        for(int j=i;j<=n;j++)
            if (fabs(A[j][i])>=eps){
                k=j;
                break;
            }
        if (k!=i){
            for(int j=1;j<=n;j++){
                swap(A[i][j],A[k][j]);
                swap(B[i][j],B[k][j]);
            }
        }
        double x=A[i][i];
        for(int j=1;j<=n;j++){
            A[i][j]/=x;
            B[i][j]/=x;
        }
        for(int j=i+1;j<=n;j++)
            if (fabs(A[j][i])>=eps){
                double x=A[j][i];
                for(int k=1;k<=n;k++){
                    A[j][k]-=A[i][k]*x;
                    B[j][k]-=B[i][k]*x;
                }
            }
    }
    for(int i=n-1;i;i--)
        for(int j=1;j<i;j++)
            if (fabs(A[j][i])>=eps){
                double x=A[j][i];
                for(int k=1;k<=n;k++){
                    A[j][k]-=A[i][k]*x;
                    B[j][k]-=B[i][k]*x;
                }
            }
}
int main(){
    scanf("%d%d%d",&n,&m,&q);
    for(int i=1;i<=m;i++){
        scanf("%d%d",&x,&y);
        x++,y++;
        deg[x]++;
        M[x][y]++;
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)M[i][j]=1.0*M[i][j]/deg[i];
    for(int i=1;i<n;i++)
        for(int j=1;j<=n;j++)A[i][j]=M[j][i]-(i==j);
    for(int i=1;i<=n;i++)A[n][i]=1;
    A[n][n+1]=1;
    guess1();
    for(int i=1;i<=n;i++)B[i][i]=1.0/A[i][n+1];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)B[i][j]--;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)A[i][j]=M[i][j]-(i==j);
    guess2();
    for(int i=1;i<=n;i++)
        if (fabs(A[i][n])>=eps)f[n][i]=B[i][i]/A[i][n];
    for(int i=1;i<n;i++)
        for(int j=1;j<=n;j++)f[i][j]=B[i][j]-A[i][n]*f[n][j];
    for(int i=1;i<=q;i++){
        scanf("%d%d",&x,&y);
        y++;
        ans=0;
        for(int j=1;j<x;j++){
            scanf("%d",&z);
            z++;
            ans+=f[y][z];
            y=z;
        }
        printf("%.9f\n",ans);
    }
}