HDU4089（概率dp）

题解

要点：

1.转移方程分三段，这个……有点复杂但是还好吧……大概就是求啥设啥，然后只通过可行的状态过来。在纸上记一记。

2.每层里面必须先求dp[i][i]，简直就是我求我自己……用类似进制数那种方式解个方程。

3.居然还有eps特判这事情。果然去掉标准同步以后就不能用puts了，会WA～

 const int maxn = ;
 const db eps = 1e-;
 int n, m, k;
 db p1, p2, p3, p4;
 db dp[][maxn];
 
 int main() {
     ios_base::sync_with_stdio(false);
     cin.tie();
     while (cin >> n >> m >> k >> p1 >> p2 >> p3 >> p4) {
         if (p4 < eps) {
             cout << fixed << setprecision() << (db) << endl;
             continue;
         }
         init(dp, );
         dp[&][] = p4 / (p3 + p4);
         rep(i, , n) {
             db sum = , pp = ;
             irep(j, i, ) {
                 if (j == )    sum += p4 / ( - p1) * pp;
                 else if ( < j && j <= k)    sum += (p3 * dp[i-&][j - ] + p4) / ( - p1) * pp;
                 else    sum += p3 / ( - p1) * dp[i-&][j - ] * pp;
                 pp *= p2 / ( - p1);
             }
             dp[i&][i] = sum / ( - pp);
 
             rep(j, , i - ) {
                 if (j == )    dp[i&][j] = (p2 * dp[i&][i] + p4) / ( - p1);
                 else if ( < j && j <= k)    dp[i&][j] = (p2 * dp[i&][j - ] + p3 * dp[i-&][j - ] + p4) / ( - p1);
                 else    dp[i&][j] = (p2 * dp[i&][j - ] + p3 * dp[i-&][j - ]) / ( - p1);
             }
         }
         cout << fixed << setprecision() << dp[n&][m] << endl;
     }
     return ;
 }