景驰无人驾驶 1024 编程邀请赛 E题 （被大佬智商碾压的）DP

对于每一个合法串，都有很多种方案，我们不妨把这些方案记为，f1，f2，f3……fx

p = （f1+f2+f3+……+fk）*（f1+f2+f3+……+fk）；

把他展开是不是很简单！！！但是我怎么就想不到！！！

展开之后 再来看，是两个 fi * fj 的和

这两个fi和fj的方案相同，那么 可以把他看成两个人，两个人分别走，并且他们序列相同的概率。

dp[i][j][k] 第i秒，第一个人走在j位置，第二个人走在k位置。

#include <bits/stdc++.h>
const long long mod = 1e9+7;
const double ex = 1e-10;
#define inf 0x3f3f3f3f
using namespace std;
int ran[60];
double p[60][60];
double dp[35][60][60];
int main()
{
    int T;
    scanf("%d",&T);
    while (T--){
        int n,m;
        scanf("%d%d",&n,&m);
        for (int i = 1; i<=n; i++){
            for (int j = 1; j <=n; j++){
                cin >> p[i][j];
            }
        }
        int z;
        scanf("%d",&z);
        for (int i = 1; i<=n;i++){
            scanf("%d",&ran[i]);
        }
        dp[1][1][1] = 1;
        for (int i = 2; i<=m ;i++){
            for (int j = 1; j <=n ;j++){
                for (int k = 1; k <=n ;k++){
                    dp[i][j][k] = 0;
                    if (ran[j]!=ran[k]) continue;
                    for (int frj = 1; frj <=n ; frj++){
                        if (abs(ran[j] - ran[frj]) > 2) continue;
                        for (int frk = 1; frk <=n ; frk++){
                            dp[i][j][k] += dp[i-1][frj][frk]*p[frj][j]*p[frk][k];
                        }
                    }
                }
            }
        }
        double ans = 0;
        for (int i = 1; i<=n; i++)
            for (int j =  1; j <=n; j++)
                ans += dp[m][i][j];
        printf("%.18f
",ans);
    }
    return 0;
}