洛谷：P3556 [POI2013]MOR-Tales of seafaring（拆点 + 奇偶最短路）

P3556 [POI2013]MOR-Tales of seafaring

如何评价不知道简单路径是什么，切割！

简单路径：路径上的所有点至多只访问了一次（2走到3、3走回2 就算重复）

同时此题边权全部都为 1 （显然引导我们可以用 dfs 求最短路优化时间）

一看这数据范围，询问次数那么多，长度范围又大，暴力 dfs 肯定没法处理

所以考虑下性质相关，长度为 d 的路径，不一定是 （不重复）简单路径，那么我们可以重复（来回）地走路径上的一条边，使总长度 +2+2+2… ，从而使路径长度增大；

但来回走代表了奇偶性的严格限制，若 a、b 两点间只有两条边，显然怎么都走不出长度为 奇数 的路径 （来回走长度 +2，奇偶性不变），反之偶数情况亦然

因此此题的做法就确定了：

有关拆点

• 维护一个奇偶最短路，将原来的每个点拆点，分成奇点和偶点，一个奇点显然只能由偶点走多 1步 得来，反之亦然；
• 边权都为1，用 dfs 优化时间，预处理出每个点（奇偶）到另一个点（奇偶）的 dist
• 如果询问某两点间 a、b 间是否有 偶数 长度的路径，那就判断下 偶点a 到 偶点b 的 dist 是否存在；否则判断 偶点a 到 奇点b 的 dist 是否存在；还要注意下路径的长度是否小于 d

代码：

#include<bits/stdc++.h>
#include<unordered_set>
#include<unordered_map>
#define mem(a,b) memset(a,b,sizeof a)
#define cinios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))
#define sca scanf
#define pri printf
#define ul (u << 1)
#define ur (u << 1 | 1)
#define fx first
#define fy second
//#pragma GCC optimize(2)
//[博客地址](https://blog.csdn.net/weixin_51797626?t=1) 
using namespace std;

typedef long long ll;
typedef pair<int, int> PII;
typedef pair <ll, PII> PI;

const int N = 10010, M = 40010, MM = N;
int INF = 0x3f3f3f3f, mod = 100003;
ll LNF = 0x3f3f3f3f3f3f3f3f;
int n, m, k, T, S, D;
int h[N], e[M], ne[M], idx;
int dist[N][N];

void add(int a, int b) {
	e[idx] = b, ne[idx] = h[a], h[a] = idx++;
}

void bfs(int x) {
	queue<int> q;
	q.push(x);

	while (q.size())
	{
		int t = q.front();
		q.pop();

		for (int i = h[t]; ~i; i = ne[i]) {
			int j = e[i];
			if (!dist[x][j]) { //第一次更新到一定是最短距离
				dist[x][j] = dist[x][t] + 1;//+1
				q.push(j);
			}
		}
	}
}

int main() {
	cinios;

	cin >> n >> m >> k;
	mem(h, -1);

	for (int i = 0; i < m; i++) {
		int a, b;
		cin >> a >> b;
		add(a, b + n), add(b + n, a);//偶点的编号在 1 ~ n 范围
		add(a + n, b), add(b, a + n);//奇点在 n+1 ~ 2*n
		//无向图，给出两点，让偶点连向奇点，反之亦然
		//方便之后维护长度
	}

	for (int i = 1; i <= n; i++)//每一个点到其他点
		bfs(i);

	while (k--)
	{
		int a, b, x;
		cin >> a >> b >> x;//判断长度的奇偶性
		if (dist[a][b + n * (x & 1)] <= x && dist[a][b + n * (x & 1)])
			cout << "TAK";//存在
		else cout << "NIE";
		cout << '\n';
	}

	return 0;
}
/*
*/