链接
题解
不妨把非树边的端点称作”特殊点“,预处理200个特殊点到其它所有点的距离
那么一次查询 的答案就是枚举特殊点 , 和树上 长度取较小值就是答案
しかし,T了!
经过我的实测,每次 大概 ,那么 次就是
我好难啊…
そして、他の人のコードを見た
然后我就发现一个大佬是这样么写的:
他不枚举特殊点,而是枚举非树边,对于一条非树边 , 的时候把 的初始距离都设为 ,然后开始 ,再拿 更新每次查询,就过了!
我抄了他的代码,然后我也过了!
でも、どうして?
这样做的时候,有一种不存在的方案被考虑了,这种方案就是
+ 1(不存在的边) +
但是就全局来看,我们其实考虑过
+
所以那种不合法答案实际上会被更新掉
代码
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#define iinf 0x3f3f3f3f
#define linf (1ll<<60)
#define eps 1e-8
#define maxn 100010
#define maxe 210010
#define maxt 210
#define maxk 17
#define cl(x) memset(x,0,sizeof(x))
#define rep(i,a,b) for(i=a;i<=b;i++)
#define drep(i,a,b) for(i=a;i>=b;i--)
#define em(x) emplace(x)
#define emb(x) emplace_back(x)
#define emf(x) emplace_front(x)
#define fi first
#define se second
#define de(x) cerr<<#x<<" = "<<x<<endl
using namespace std;
using namespace __gnu_pbds;
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
ll read(ll x=0)
{
ll c, f(1);
for(c=getchar();!isdigit(c);c=getchar())if(c=='-')f=-f;
for(;isdigit(c);c=getchar())x=x*10+c-0x30;
return f*x;
}
struct Graph
{
int etot, head[maxn], to[maxe], next[maxe], w[maxe];
void clear(int N)
{
for(int i=1;i<=N;i++)head[i]=0;
etot=0;
}
void adde(int a, int b, int c=0){to[++etot]=b;w[etot]=c;next[etot]=head[a];head[a]=etot;}
#define forp(_,__) for(auto p=__.head[_];p;p=__.next[p])
}G;
struct UnionFind
{
ll f[maxn], size[maxn];
void init(ll n)
{
for(auto i=1;i<=n;i++)f[i]=i, size[i]=1;
}
ll find(ll x){return f[x]==x?x:f[x]=find(f[x]);}
void merge(ll x, ll y)
{
auto fx=find(x), fy=find(y);
f[fx]=fy;
if(fx!=fy)size[fy]+=size[fx];
}
}uf;
struct Doubling_LCA
{
int f[maxn][maxk+1], depth[maxn];
void clear(int n){for(int i=1;i<=n;i++)depth[i]=0, cl(f[i]);}
void dfs(Graph &G, int pos, int pre)
{
for(auto k=1;(1<<k)<=depth[pos];k++)f[pos][k]=f[f[pos][k-1]][k-1];
for(auto p(G.head[pos]);p;p=G.next[p])
if(G.to[p]!=pre)
{
f[G.to[p]][0]=pos;
depth[G.to[p]]=depth[pos]+1;
dfs(G,G.to[p],pos);
}
}
void run(Graph &G, int root)
{
depth[root]=1;
dfs(G,root,0);
}
int q(int x, int y)
{
if(depth[x]<depth[y])swap(x,y);
for(auto k(maxk);~k;k--)
if(depth[f[x][k]]>=depth[y])
x=f[x][k];
if(x==y)return x;
for(auto k(maxk);~k;k--)
if(f[x][k]!=f[y][k])
x=f[x][k], y=f[y][k];
return f[x][0];
}
int jp(int x, int b)
{
for(auto k=0;k<=maxk;k++)
if(b&(1<<k))x=f[x][k];
return x;
}
}db;
vector<pii> e;
int d[maxn], n, m, num[maxn], t, ans[maxn], a[maxn], b[maxn];
int q[maxn];
void bfs(int S1, int S2)
{
for(int i=1;i<=n;i++)d[i]=iinf;
int l=0, r=0;
auto to = G.to, nex = G.next, head = G.head;
d[S1]=d[S2]=0;
q[r++]=S1, q[r++]=S2;
while(l<r)
{
auto x=q[l++];
for(int p=head[x];p;p=nex[p])
{
if(d[to[p]]==iinf)
{
d[to[p]]=d[x]+1;
q[r++]=to[p];
}
}
}
}
int dis(int x, int y)
{
int lca = db.q(x,y);
return db.depth[x] + db.depth[y] - 2*db.depth[lca];
}
int main()
{
int u, v, q, i;
n=read(), m=read();
uf.init(n);
rep(i,1,m)
{
u=read(), v=read();
if(uf.find(u)==uf.find(v))
{
e.emb( pii(u,v) );
}
else
{
G.adde(u,v);
G.adde(v,u);
uf.merge(u,v);
}
}
db.run(G,1);
for(auto pr:e)
{
G.adde(pr.fi,pr.se);
G.adde(pr.se,pr.fi);
}
q=read();
rep(i,1,q)
{
a[i]=read(), b[i]=read();
ans[i] = dis(a[i],b[i]);
}
for(auto pr:e)
{
bfs(pr.fi,pr.se);
rep(i,1,q)
{
ans[i] = min( ans[i], d[a[i]]+1+d[b[i]]);
}
}
rep(i,1,q)
{
printf("%d\n",ans[i]);
}
return 0;
}