家人们,你们想要的代码和讲解来啦,内容摘自:浚浚的二师兄
[题目通道]([CERC2015] Looping Labyrinth - 洛谷)
这条城墙把平面分成三个部分:上面,下面 和 里面。不难发现 11 路径如果和 00 路径无交,那就必定完全存在于这三部分之一。不难证明它不可能在里面,要在就在上面或下面,不妨设为上面。
那么 22 路径在 11 路径的上面,33 路径在 22 路径的的上面,……𝜆λ 路径在 𝜆−1λ−1 路径的上面,自然也就在 00 路径的上面。
——等等!𝜆λ 路径难道不就是 00 路径自己吗?!一条路径怎么可能在自己的上面呢?这引出了矛盾,原假设不成立。
有了以上知识储备,我们可以发现,可达的 (∗,∗,0,0)(∗,∗,0,0) 构成的集合的形态是极度有限的。
{(0,0,0,0)}{(0,0,0,0)}。最简单的形态。
{(𝑢,𝑣,0,0)},gcd(𝑢,𝑣)=1{(u,v,0,0)},gcd(u,v)=1。即找到了一组可达 (𝑢,𝑣,0,0)(u,v,0,0)。
如果找到了两组 (𝑢1,𝑣1,0,0),(𝑢2,𝑣2,0,0)(u1,v1,0,0),(u2,v2,0,0),那么不难发现任何一个 (𝑥,𝑦,0,0)(x,y,0,0) 都能被表示(具体证明可以见鱼大的博客),即形态为 {(𝑥,𝑦,0,0),𝑥,𝑦∈𝑍}{(x,y,0,0),x,y∈Z}。
考虑求解这个"块可达性形态"。我们从 (0,0,0,0)(0,0,0,0) 开始 BFS,限制每个块内坐标至多访问一次,若试图访问多次便说明我们找到了一组 (𝑢,𝑣)(u,v)。这样复杂度是 𝑂(𝑛𝑚)O(nm) 的,还顺利完成了形态求解的任务。
然后又发现这个形态好像只和 BFS 出的集合有关,换句话说:除了那些根本到不了的点,其他所有点的"块可达性形态"是一样的。直接用相同的做法套就完事了。
代码如下:
#include<bits/stdc++.h>
using namespace std;
int n, m;
bool H[105][105];
struct node {
int xB, yB, xI, yI;
};
node calc(int x, int y) {
node ans;
ans.xB = x >= 0 ? (x / n) : ((x + 1) / n - 1);
ans.yB = y >= 0 ? (y / m) : ((y + 1) / m - 1);
ans.xI = x - ans.xB * n; assert(ans.xI >= 0 && ans.xI < n);
ans.yI = y - ans.yB * m;
return ans;
}
int U = 0, V = 0;
void insert(int nU, int nV) {
if (U == -998244353) return; // case 2
if (nU == 0 && nV == 0) return; // meaningless
if (nU < 0) nU = -nU, nV = -nV;
else if (nU == 0 && nV < 0) nV = -nV;
int g = __gcd(nU, abs(nV)); nU /= g, nV /= g;
if (nU == U && nV == V) return;
if (U || V) U = V = -998244353; // case 1 -> case 2
else U = nU, V = nV; // case 0 -> case 1
}
bool vis[105][105];
pair<int, int> pB[105][105];
bool isvalid(node u) {
if (H[u.xI][u.yI]) return 0;
if (!vis[u.xI][u.yI]) {
vis[u.xI][u.yI] = 1;
pB[u.xI][u.yI] = make_pair(u.xB, u.yB);
return 1;
}
insert(u.xB - pB[u.xI][u.yI].first, u.yB - pB[u.xI][u.yI].second);
return 0;
}
int d[4][2] = {{1, 0}, {-1, 0}, {0, 1}, {0, -1}};
int main() {
scanf("%d%d", &n, &m);
for (int i = 0; i < n; i++) {
char c = getchar();
while (c != '.' && c != '#') c = getchar();
for (int j = 0; j < m; j++)
H[i][j] = c == '#', c = getchar();
}
queue<pair<int, int> > Q;
isvalid((node){0, 0, 0, 0}); Q.push(make_pair(0, 0));
while (!Q.empty()) {
pair<int, int> u = Q.front(); Q.pop();
for (int i = 0; i < 4; i++) {
int vx = u.first + d[i][0], vy = u.second + d[i][1];
if (isvalid(calc(vx, vy))) Q.push(make_pair(vx, vy));
}
}
int T; scanf("%d", &T);
while (T--) {
int x, y; scanf("%d%d", &x, &y);
node qaq = calc(x, y);
if (!vis[qaq.xI][qaq.yI]) { printf("no\n"); continue; }
if (U == 0 && V == 0) {
if (qaq.xB != pB[qaq.xI][qaq.yI].first || qaq.yB != pB[qaq.xI][qaq.yI].second)
printf("no\n");
else printf("yes\n");
continue;
}
if (U == -998244353 && V == -998244353) {
printf("yes\n"); continue;
}
int nU = qaq.xB - pB[qaq.xI][qaq.yI].first,
nV = qaq.yB - pB[qaq.xI][qaq.yI].second;
if (nU == 0 && nV == 0) { printf("yes\n"); continue; }
if (nU < 0) nU = -nU, nV = -nV;
else if (nU == 0 && nV < 0) nV = -nV;
int g = __gcd(nU, abs(nV)); nU /= g, nV /= g;
if (nU == U && nV == V) { printf("yes\n"); continue; }
printf("no\n");
}
}
————————————————
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
原文链接:https://blog.csdn.net/2302_78855882/article/details/141112254