Rating

题目链接

• 题意：
  起始状态是（0。0），每次转移的时候都是对两个数中的较小的数操作。

  1）以概率p转向（min（a + 50，1000）。b）    2）以概率1-p转向（max（a-100，0），b）

• 分析：
  首先发现状态转移的时候都是以50为单位，所以事实上就是除以50之后。即加1或者减2。达到20就可以
• 注意：
  题目精度要求比較高。eps至少要到1e-10

const double eps = 1e-10;
double a[22 * 22][22 * 22], x[22 * 22]; //方程的左边的矩阵和等式右边的值，求解之后x存的就是结果
int equ, var;                            //方程数和未知数个数
int Gauss()
{int i, j, k, col, max_r;for (k = 0, col = 0; k < equ && col < var; k++, col++){max_r = k;for (i = k + 1; i < equ; i++)if (fabs(a[i][col]) > fabs(a[max_r][col]))max_r = i;if (fabs(a[max_r][col]) < eps) return 0;if (k != max_r){for (j = col; j < var; j++)swap(a[k][j], a[max_r][j]);swap(x[k], x[max_r]);}x[k] /= a[k][col];for (j = col + 1; j < var; j++) a[k][j] /= a[k][col];a[k][col] = 1;for (i = 0; i < equ; i++)if (i != k){x[i] -= x[k] * a[i][k];for (j = col + 1; j < var; j++) a[i][j] -= a[k][j] * a[i][col];a[i][col] = 0;}}return 1;
}double P;
int s[22 * 22][22 * 22];
void build()
{CLR(a, 0); CLR(x, 0);FE(i, 0, 20) FE(j, i, 20){int cur = s[i][j];if (~cur){a[cur][cur] = 1;if (i == 20 || j == 20)x[cur] = 0;else{int tx = min(i + 1, 20), ty = j;if (tx > ty) swap(tx, ty);int nxt = s[tx][ty];a[cur][nxt] -= P;tx = max(i - 2, 0); ty = j;if (tx > ty) swap(tx, ty);nxt = s[tx][ty];a[cur][nxt] -= 1 - P;x[cur] = 1;}}}
}
void bfs()
{CLR(s, -1);queue<int> qx, qy;qx.push(0); qy.push(0);int cnt = 0;s[0][0] = cnt++;while (!qx.empty()){int x = qx.front(); qx.pop();int y = qy.front(); qy.pop();int tx = min(x + 1, 20), ty = y;if (tx > ty) swap(tx, ty);if (!~s[tx][ty]){s[tx][ty] = cnt++;qx.push(tx); qy.push(ty);}tx = max(x - 2, 0), ty = y;if (!~s[tx][ty]){s[tx][ty] = cnt++;qx.push(tx); qy.push(ty);}}equ = var = cnt;
}int main()
{bfs();while (cin >> P){build();Gauss();printf("%.6lf\n", x[s[0][0]]);}return 0;
}