程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
 程式師世界 >> 編程語言 >> C語言 >> C++ >> C++入門知識 >> UVA - 1498 Activation (DP+概率)

UVA - 1498 Activation (DP+概率)

編輯:C++入門知識

UVA - 1498 Activation (DP+概率)


Description

Download as PDF

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:


Activation failed: This happens with the probability of p1. The queue remains unchanged and the server will try to deal with the same request the next time.Connection failed: This happens with the probability of p2. Something just happenedand the first player in queue lost his connection with the server. The server will then remove his request from the queue. After that, the player will immediately connect to the server againand starts queuing at the tail of the queue.Activation succeeded: This happens with the probability of p3. Congratulations, theplayer will leave the queue and enjoy the game himself.Service unavailable: This happens with the probability of p4. Something just happened and the server is down. The website must shutdown the server at once. All the requests that are still in the queue will never be dealt.

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$ \le$M$ \le$N$ \le$2000), K (K$ \ge$1), p1, p2, p3, p4(0$ \le$p1, p2, p3, p4$ \le$1, 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;
kn 遞推求解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 
#include 
#include 
#include 
#include 
#include 
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;
}


  1. 上一頁:
  2. 下一頁:
Copyright © 程式師世界 All Rights Reserved