标签：des style blog http os io ar for 2014
Description
After 4 years‘ waiting, the game "Chinese Paladin 5" finally comes out. Tomato is a crazy fan, and luckily he got the first release. Now he is at home, ready to begin his journey. But before starting the game, he must first activate the product on the official site. There aretoo many passionate fans that the activation server cannot deal with all the requests at the same time, so all the players must wait in queue. Each time, the server deals with the request of the first player in the queue, and the result may be one of the following, each has aprobability:
Tomato thinks it sucks if the server is down while he is still waiting in the queue and there are no more than K - 1 guys before him. And he wants to know the probability that this ugly thing happens.
To make it clear, we say three things may happen to Tomato: he succeeded activating the game; the server is down while he is in the queue and there are no more than
K - 1 guys before him; the server is down while he is in the queue and there are at least
K guys before him.Now you are to calculate the probability of the second thing.
Input
There are no more than 40 test cases. Each case in one line, contains three integers and four real numbers:
N, M (1
MN2000),
K (K
1),
p1, p2, p3, p4(0p1,
p2, p3, p41,
p1 + p2 + p3 + p4 = 1), indicating there are
N guys in the queue (the positions are numbered from 1 to
N), and at the beginning Tomato is at the
Mth position, with the probability p1,
p2, p3, p4 mentioned above.
Output
A real number in one line for each case, the probability that the ugly thing happens.
The answer should be rounded to 5 digits after the decimal point.
Sample Input
2 2 1 0.1 0.2 0.3 0.4 3 2 1 0.4 0.3 0.2 0.1 4 2 3 0.16 0.16 0.16 0.52
Sample Output
0.30427 0.23280 0.90343 题意：：有n个人排队等着在官网上激活游戏。Tomato排在第m个。对于队列中的第一个人。有一下情况： 1、激活失败，留在队列中等待下一次激活（概率为p1) 2、失去连接，出队列，然后排在队列的最后（概率为p2） 3、激活成功，离开队列（概率为p3） 4、服务器瘫痪，服务器停止激活，所有人都无法激活了。 求服务器瘫痪时Tomato在队列中的位置<=k的概率 思路：kuangbinGG写的很好，我就不说了 借鉴kuangbinGG的题解：概率DP； 设dp[i][j]表示i个人排队,Tomato排在第j个位置，达到目标状态的概率(j<=i) dp[n][m]就是所求 j==1: dp[i][1]=p1*dp[i][1]+p2*dp[i][i]+p4; 2<=j<=k: dp[i][j]=p1*dp[i][j]+p2*dp[i][j-1]+p3*dp[i-1][j-1]+p4; k<j<=i: dp[i][j]=p1*dp[i][j]+p2*dp[i][j-1]+p3*dp[i-1][j-1]; 化简： j==1: dp[i][1]=p*dp[i][i]+p41; 2<=j<=k: dp[i][j]=p*dp[i][j-1]+p31*dp[i-1][j-1]+p41; k<j<=i: dp[i][j]=p*dp[i][j-1]+p31*dp[i-1][j-1]; 其中: p=p2/(1-p1); p31=p3/(1-p1) p41=p4/(1-p1) 可以循环i=1->n 递推求解dp[i].在求解dp[i]的时候dp[i-1]就相当于常数了。 在求解dp[i][1~i]时等到下列i个方程 j==1: dp[i][1]=p*dp[i][i]+c[1]; 2<=j<=k:dp[i][j]=p*dp[i][j-1]+c[j]; k<j=i: dp[i][j]=p*dp[i][j]+c[j]; 其中c[j]都是常数了。上述方程可以解出dp[i]了。 首先是迭代得到 dp[i][i].然后再代入就可以得到所有的dp[i]了。 注意特判一种情况。就是p4<eps时候，就不会崩溃了，应该直接输出0#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> using namespace std; const int maxn = 2020; const double eps = 1e-5; double c[maxn], f[maxn], dp[maxn][maxn]; int main() { int n, m, k; double p1, p2, p3, p4; while (scanf("%d%d%d%lf%lf%lf%lf", &n, &m, &k, &p1, &p2, &p3, &p4) != EOF) { if (p4 < eps) { printf("0.00000\n"); continue; } double p = p2/(1-p1); double p41 = p4/(1-p1); double p31 = p3/(1-p1); f[0] = 1.0; for (int i = 1; i <= n; i++) f[i] = p * f[i-1]; dp[1][1] = p41 / (1-p); c[1] = p41; for (int i = 2; i <= n; i++) { for (int j = 2; j <= k; j++) c[j] = p31 * dp[i-1][j-1] + p41; for (int j = k+1; j <= i; j++) c[j] = p31 * dp[i-1][j-1]; double tmp = c[1] * f[i-1]; for (int j = 2; j <= i; j++) tmp += c[j] * f[i - j]; dp[i][i] = tmp / (1 - f[i]); dp[i][1] = p * dp[i][i] + c[1]; for (int j = 2; j < i; j++) dp[i][j] = p * dp[i][j-1] + c[j]; } printf("%.5lf\n", dp[n][m]); } return 0; }
标签：des style blog http os io ar for 2014
原文地址：http://blog.csdn.net/u011345136/article/details/39082379
