题面
分析
首先得到一个DP方程
高消
走人。
不能直接高消,我就想按照CF963E Circles of Waiting的做法,由于有很多系数为 ,那就标记一下不为零的地方,乱搞一下把高消变成 的。经过超过2小时的尝试,我证明了自己还是太菜太naive。
由于这个环面上可以随意“平移”,不妨把 放到 的地方,那么我们的问题就变成了求 到 的步数期望。
我们把第 行和第 列的 值作为主元,根据 式,我们可以用这些主元由下至上,由右至左表示出其他的 值,然后再由 式将第 行、第 列、第 行、第 列之间的关系表示成一个 元 次方程组,用高斯消元解这个方程组,时间复杂度 ,进而可以知道所有 值。
代码
点 不要处理,否则会算重。如果访问到 说明起点终点重合,答案是 。
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <vector>
const int MAXN = 100;
struct Node {
double A[MAXN * 2 + 5], B;
}P[MAXN + 5][MAXN + 5];
double Abs(double x) {
return x > 0 ? x : -x;
}
double A[MAXN * 2 + 5][MAXN * 2 + 5], Res[MAXN * 2 + 5];
void Gauss(int N) {
for (int i = 0; i < N; i++) {
double tmp = A[i][i];
for (int j = i; j <= N; j++)
A[i][j] /= tmp;
for (int j = i + 1; j < N; j++) {
tmp = A[j][i];
for (int k = i; k <= N; k++)
A[j][k] -= tmp * A[i][k];
}
}
for (int i = N - 1; i >= 0; i--) {
Res[i] = A[i][N];
for (int j = i + 1; j < N; j++)
Res[i] -= A[i][j] * Res[j];
}
}
class TorusSailing {
public:
double expectedTime(int N, int M, int X, int Y) {
for (int i = 0; i < N - 1; i++)
P[i][M - 1].A[i] = 1;
for (int i = 0; i < M - 1; i++)
P[N - 1][i].A[i + N - 1] = 1;
for (int i = N - 2; i >= 0; i--)
for (int j = M - 2; j >= 0; j--) {
P[i][j].B = 0.5 * (P[i][j + 1].B + P[i + 1][j].B) + 1;
for (int k = 0; k < N + M - 2; k++)
P[i][j].A[k] = 0.5 * (P[i][j + 1].A[k] + P[i + 1][j].A[k]);
}
for (int i = 0; i < N - 1; i++) {
A[i][i] = 1, A[i][N + M - 2] = 0.5 * P[i][0].B + 1;
for (int k = 0; k < N + M - 2; k++)
A[i][k] -= 0.5 * (P[i][0].A[k] + P[i + 1][M - 1].A[k]);
}
for (int i = 0; i < M - 1; i++) {
int j = i + N - 1;
A[j][j] = 1, A[j][N + M - 2] = 0.5 * P[0][i].B + 1;
for (int k = 0; k < N + M - 2; k++)
A[j][k] -= 0.5 * (P[0][i].A[k] + P[N - 1][i + 1].A[k]);
}
Gauss(N + M - 2);
X = N - 1 - X, Y = M - 1 - Y;
double Ans = P[X][Y].B;
for (int i = 0; i < N + M - 2; i++)
Ans += Res[i] * P[X][Y].A[i];
return Ans;
}
};