sdut2878 环形依赖的DP(高斯消元，剪枝后的模板

这题的状态是循环依赖的有环。。

之前一道概率DP，类似有环。。但是它是可以消掉的

比如dp[i]=0.3*dp[i+1]+0.2*dp[i+2]+0.5*dp[i];

完全可以变成，0.5*dp[i]=0.3*dp[i+1]+0.2*dp[i+2]

然后把系数除过去就好了，

然而这个题是，dp[i]=0.5*dp[i+1]+0.5*dp[i-1]+1;

这个+1是什么意思呢，dp[i]->要到i+1,i-1任意两个状态之一，一定要付出1步的代价！

想一想背包问题。。类似的，

然后你会发现dp[i]还没递归完。。dp[i-1]和dp[i+1]又跑回来找它了。。。这不可能DP得出来

所以这个题正确的消环方法是，高斯消元

高斯消元我还不太会写。。待我学习一波回来写一写！

（学了高斯消元回来补辣！

我博客里的模板正常用是会TLE的！

需要剪枝！

#include <iostream>
#include <cstdio>
#include <vector>
#include <cmath>
#include <cstring>
using namespace std;
const double EPS=1e-8;
int n,p;double A[1005][1005];
double x[1005];
int Gauss(){
    for(int i=0;i<n;++i){
        int pivot=i;
        for(int j=i+1;j<n;++j)if(abs(A[j][i]>abs(A[pivot][i]))) pivot=j;
       if(pivot!=i) for(int k=0;k<n;++k) swap(A[i][k],A[pivot][k]);
        if(abs(A[i][i])<EPS) return 0;
        for(int j=i+1;j<=n;++j) A[i][j]/=A[i][i];
        for(int j=i+1;j<n;++j)
            if(j!=i) {
                if(abs(A[j][i])<EPS) continue;
                for(int k=i+1;k<=n;++k) A[j][k]-=A[j][i]*A[i][k];
            }
    }
    for(int i=0;i<n;++i) x[i]=A[i][n];  double ans;
    for(int i=n-1;i>=0;--i){
        ans=x[i];
        for(int j=i+1;j<n;++j) ans-=A[i][j]*x[j];x[i]=ans;
    }
    return 1;
}
int main(){
    int T;scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&p);
        //构造系数矩阵
        memset(A,0,sizeof(A));
        for(int i=0;i<n;++i){
            if(i==p) {
                A[i][i]=1;A[i][n]=0;
                continue;
            }
            A[i][i]=1;A[i][(i-1+n)%n]=-0.5;
            A[i][(i+1)%n]=-0.5;A[i][n]=1;
        }
        Gauss();
        printf("%.4f
",x[0]);
    }
    return 0;
}

14行和19行是剪枝的地方，不减的话这题最高1000^3还是会T的