高斯消元+期望dp——light1151

高斯消元弄了半天没弄对。。

#include<bits/stdc++.h>
using namespace std;
#define maxn 205
#define eps 1e-8
double A[maxn][maxn],x[maxn],ans[maxn];
int nxt[maxn],n;
#define a A
void guess(int n){   //行，列
    for(int i=;i<n;i++){
        if(a[i][i]==){//求主元的时候不能直接用swap进行交换
            int id=;
            for(int j=i+;j<=n;j++)
                if(a[j][i]!=)
                    id=j;
            for(int j=i;j<=n+;j++)
                swap(a[i][j],a[id][j]);
        }
        for(int j=i+;j<=n;j++){//消下三角
            double t=a[j][i]/a[i][i];
            for(int k=i;k<=n+;k++)
                a[j][k]-=(a[i][k]*t);
        }
    }

    for(int i=n;i>=;i--)
    {
        for(int j=i+;j<=n;j++)
            a[i][n+]-=ans[j]*a[i][j];
        ans[i]=a[i][n+]/a[i][i];
    }
}
int main(){
    int t;cin>>t;
    for(int tt=;tt<=t;tt++){
        int n;cin>>n;
        memset(nxt,,sizeof nxt);
        memset(A,,sizeof A);
        for(int i=;i<n;i++){int x,y;cin>>x>>y;nxt[x]=y;}
        for(int i=;i<;i++){//建立矩阵
            if(nxt[i]){
                A[i][i]=;
                A[i][]=;
                A[i][nxt[i]]=-;
            }
            else {
                int k=;
                for(int j=;i+j<= && j<=;j++){
                    k++;
                    A[i][i+j]=-;
                }
                A[i][i]=k,A[i][]=;
            }
        }        

        A[][]=;A[][]=;
        guess();
        printf("Case %d: %.10lf\n",tt,ans[]);
    }
    return ;
}