[HNOI2019]校园旅行

题目

过于神仙啊，抄题解.jpg

首先\(n\)并不是很大啊，我们可以直接用\(f_{i,j}\)表示\(i\)到\(j\)是否存在一个回文路径

对于一条回文路径，如果在两端分别添加一个相同的字符，那么仍然是一个回文路径，于是我们可以利用这个来打一个暴力\(bfs\)

就像这样

while(!q[0].empty()) {int x=q[0].front(),y=q[1].front();q[0].pop(),q[1].pop();for(re int i=0;i<v[x].size();i++)for(re int j=0;j<v[y].size();j++) {int xx=v[x][i],yy=v[y][j];if(xx>yy) std::swap(xx,yy);if(S[xx]!=S[yy]||f[xx][yy]) continue;f[xx][yy]=dp[xx][yy]=1;q[0].push(xx),q[1].push(yy);}}

非常显然这个复杂度一点也不科学，卡满就是\(O(m^2)\)的，肯定是要\(T\)的

但是\(n\)却不是很大，能不能让边数减小一点呢

之后就不会啦，愉快地抄题解

首先我们注意到一个问题，就是我们只需要关注奇偶性就好了，因为我们可以来回走一条边使得长度增加，但是并不能改变奇偶性

之后我也不知道为什么通过这一点就想到了二分图

我们把边分成两类，一类是连接相同颜色点的边，一类是连接不同颜色点的边

我们考虑一个连通块，这个连通块仅由相同颜色的点构成，显然这个连通块里的边都是第一类边

如果这个连通块是一个二分图，由于二分图不存在奇环，对于任意两个点，所有连接这两个点的路径的奇偶性都是一样的

由于我们也只关注奇偶性，所以只需要对这个二分图搞出来一棵生成树即可

如果这个联通块不是一个二分图，那么就一定存在至少一个奇环，那么我们就可以通过走这个奇环使得路径的奇偶性改变

我们需要让边的数量尽量少，那么只需要求出来一棵有奇环的基环树即可

这些我们都能通过\(dfs\)来完成

对于第二类边，只考虑这些边的话整张图就是一个二分图，于是我们还是对这些边求一个生成树即可

最后我们发现我们这张图的边数已经和\(n\)同级了，于是我们直接上最开始的那个暴力，复杂度就是\(O(n^2)\)了

代码

#include<queue>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define re register
inline int read() {char c=getchar();int x=0;while(c<'0'||c>'9') c=getchar();while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
const int maxn=5e3+5;
std::vector<int> v[maxn];
std::queue<int> q[2];
struct E{int v,nxt;}e[1000005];
int n,m,Q,num,tot,flag;
char S[maxn];
int f[maxn][maxn],dp[maxn][maxn],vis[maxn],fa[maxn],sz[maxn],pre[maxn];
int head[maxn],a[500005],b[500005];
inline void add(int x,int y) {e[++num].v=y;e[num].nxt=head[x];head[x]=num;
}
void dfs(int x,int fa) {vis[x]=1;for(re int i=head[x];i;i=e[i].nxt) {if(S[e[i].v]!=S[x]||fa==e[i].v) continue;if(vis[e[i].v]&&S[e[i].v]==S[x]&&!flag) {if(!(pre[e[i].v]^pre[x])) v[x].push_back(e[i].v),v[e[i].v].push_back(x),flag=1;continue;}v[x].push_back(e[i].v);v[e[i].v].push_back(x);pre[e[i].v]=pre[x]^1;dfs(e[i].v,x);}
}
inline int find(int x) {return fa[x]==x?x:fa[x]=find(fa[x]);}
inline void merge(int x,int y) {if(sz[x]<sz[y]) fa[y]=x,sz[y]+=sz[x];else fa[x]=y,sz[x]+=sz[y];
}
int main() {n=read(),m=read(),Q=read();scanf("%s",S+1);for(re int x,y,i=1;i<=m;i++) {x=read(),y=read(),add(x,y),add(y,x);if(S[x]!=S[y]) a[++tot]=x,b[tot]=y;}for(re int i=1;i<=n;i++) {if(vis[i]) continue;flag=0;dfs(i,0);}for(re int i=1;i<=n;i++) fa[i]=i,sz[i]=1;for(re int i=1;i<=tot;i++) {int xx=find(a[i]),yy=find(b[i]);if(xx==yy) continue;merge(xx,yy);v[a[i]].push_back(b[i]);v[b[i]].push_back(a[i]);}for(re int i=1;i<=n;i++) f[i][i]=dp[i][i]=1,q[0].push(i),q[1].push(i);for(re int i=1;i<=n;i++) for(re int j=0;j<v[i].size();j++) {int x=i,y=v[i][j];if(x>y) std::swap(x,y);if(S[x]!=S[y]||f[x][y]) continue;f[x][y]=dp[x][y]=1,q[0].push(x),q[1].push(y);}while(!q[0].empty()) {int x=q[0].front(),y=q[1].front();q[0].pop(),q[1].pop();for(re int i=0;i<v[x].size();i++)for(re int j=0;j<v[y].size();j++) {int xx=v[x][i],yy=v[y][j];if(xx>yy) std::swap(xx,yy);if(S[xx]!=S[yy]||f[xx][yy]) continue;f[xx][yy]=dp[xx][yy]=1;q[0].push(xx),q[1].push(yy);}}for(re int x,y,i=1;i<=Q;i++) {x=read(),y=read();if(x>y) std::swap(x,y);puts(dp[x][y]?"YES":"NO");}return 0;
}