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 ;

 }

巴特西

HDU4089（概率dp）

最新文章

热门文章