GYM - 101147 K.Touristic Trip

题意：

　　一个人从城市0开始旅行。一共有N座城市，他每到一座城市都会寄一张明信片。给出从一座城市到另一座城市的概率和在每座城市寄出每张明信片的概率。给出长度为k的寄明信片的序列。问在该序列的条件下在第Z个城市寄出第Q张明信片的序列。

题解：

　　条件概率：P(A|B) = P(A,B)/P(B)

　　dp[i][j]代表第i个城市是城市j且寄出第a[i]张明信片的概率。

　　dp[i][j] += sum{dp[i-1][k]*p[k][j]*c[j][a[i]}(0<=j<n)

　　P(B) = sum{dp[k-1][j]}(0,=j<n)

　　令dpp[Q][Z] = dp[Q][Z]　dpp[Q][j] = 0(j!=Z)

　　后面dpp像dp一样算。

　　P(A,B) = sum{dpP[k-1][j]}(0,=j<n)

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int t;
int n, m, k, q, z;
double p[][], c[][];
double dp[][], dpp[][];
int a[];
int main() {
    freopen("trip.in","r",stdin);
    scanf("%d", &t);
    while(t--) {
        memset(dp, , sizeof(dp));
        memset(dpp, , sizeof(dpp));
        scanf("%d%d%d%d%d", &n, &m, &k, &q, &z);
        for(int i = ; i < n; i++)
        for(int j = ; j < n; j++) scanf("%lf", &p[i][j]);
        for(int i = ; i < n; i++)
        for(int j = ; j < m; j++) scanf("%lf", &c[i][j]);
        for(int i = ; i < k; i++) scanf("%d", &a[i]);
        dp[][] = c[][a[]];
        for(int i = ; i <= q; i++) {
            for(int j = ; j < n; j++)
            for(int k = ; k < n; k++)
            dp[i][j] += dp[i-][k]*p[k][j]*c[j][a[i]];
        }
        dpp[q][z] = dp[q][z];
        for(int i = q+; i < k; i++) {
            for(int j = ; j < n; j++)
            for(int k = ; k < n; k++) {
                dp[i][j] += dp[i-][k]*p[k][j]*c[j][a[i]];
                dpp[i][j] += dpp[i-][k]*p[k][j]*c[j][a[i]];
            }
        }
        for(int i = ; i < n; i++) {
            dp[k-][i] += dp[k-][i-];
            dpp[k-][i] += dpp[k-][i-];
        }
        printf("%.3lf\n", dpp[k-][n-]/dp[k-][n-]);
    }
}