Left-Right-Win

题目链接

神奇的题解

看了两天也没看懂上面链接的大佬的公式，直到我动手写了遍模拟。。

本文就当是上面题解的读后感吧，请先看完链接的题解，如果也像我一样一脸懵圈，再来看这一篇。下面主要分析五个公式的由来（我猜的）

公式一二先画图再解释

对不起，画完图我也晕了，不过我解释一下，一共6个人，赢的概率为w，左边概率为l，右边概率为r，最内圈的意思是第一次就能中奖概率，第二圈的意思是第二次中奖的概率，依次类推（我写晕了，再外圈就没写）箭头是如果上一次没中奖那么可以转移到这一次，那么总体中奖概率就是把1/6圆中所有的值加起来，这里先设状态为dp[i][j]表示i在第j次可能中奖的概率，显然第一层中dp[0][1]=w,其余全为0，其他层需要由[j-1]层转移而来，由图中1和5（其他的不明显）的中奖转移可以看出

dp[1][2] = dp[2][1]*r+dp[0][1]*l    dp[1][4] = dp[2][3]*r+dp[0][3]*l   ...，

原因就是如果上一次没中奖，那么本次中奖概率就是

(dp[(i+1)%n][j-1]/w*r+dp[(n+i-1)%n][j-1]/w*l)*w，

化简便得到了dp[(i+1)%n][j-1]*r+dp[(n+i-1)%n][j-1]*l，

而dp[j][i]也可以实现降维得到dp[(i+1)%n]*r+dp[(n+i-1)%n]*l。原因就是所有的转移都是通过上一次左右不中奖来的（太晚了，写的不太明白，好好理解一下）本题限制了最长天数天也就是说第天的转移不需要再乘以w但是个l或r相乘之后就很小了*1和*w不会影响结果（也是我猜的），所以忽略了，至此，公式一推出。而由于所有的dp[i]都是通过dp[0]得到的，而dp[0]比着其他人多了w的概率中奖（因为第一次就是他整的，如果第一次没中奖才会有以后的dp[0]），至此，公式二也出来了。

公式三我是想不到的，下次再写吧，我困死了

以下是根据dp[i][j]来做的，注意maxn不要取太大了，由于样例上给的n最大15，所以我一开始全部给的1<<15，一直WA，但是上界只有1<<n的话样例过不了，所以我写了1<<10就A了，原因就是如果层数过大的话，dp[i][j]会很小，不会影响结果，以上，公式三四五先不推了（因为我真的想不明白，或许这就是大佬吧）

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1<<10;
double dp[20][maxn+3];
double solve(int n, int k, double gl, double gr){
	dp[0][1] = 1-gl-gr;
	for (int i = 2;  i <= maxn; ++i) {
		for (int j = 0;  j < n; ++j)
			dp[j][i] = gr*dp[(j+1)%n][i-1] + gl*dp[(n+j-1)%n][i-1];
	}
	double sum = 0;
	for (int i = 1; i <= maxn; ++i) {
		sum += dp[k][i];
	}
	return sum;
}
int main () {
	int t,k,n,no;
	double gl,gr,gw;
	scanf("%d",&t);
	while(t--) {
		scanf("%d%d%d%lf%lf",&no,&n,&k,&gl,&gr);
		printf("%d %.2lf\n", no, 100* solve(n,k,gl,gr));
	}
	//system("pause");
	return 0;
}

经验：

1.动手写写画画，思路才可能出来

2.除法精度会有问题，类似分子分母也都是分数的情况下（特别是分母是分数）类似a/(b/c)的形式，先化成a*c/b再计算（本题就会遇到这样的坑）

3.如果遇到过多小数相乘，注意当个数到达一定值时就不会影响精度了