nowcoder19814 最短路

链接

点击跳转

题解

不妨把非树边的端点称作”特殊点“，预处理200个特殊点到其它所有点的距离

那么一次查询$(a,b)$的答案就是枚举特殊点$i$，$min(dist[i][u]+dist[i][v])$和树上$(u,v)$长度取较小值就是答案

しかし，T了！

经过我的实测，每次$bfs$大概$10ms$，那么$200$次就是$2000ms=2s>1s$

我好难啊…

そして、他の人のコードを見た

然后我就发现一个大佬是这样么写的：

他不枚举特殊点，而是枚举非树边，对于一条非树边$(u,v)$，$bfs$的时候把$u,v$的初始距离都设为$0$，然后开始$bfs$，再拿$dist[u]+1+dist[v]$更新每次查询，就过了！

我抄了他的代码，然后我也过了！

でも、どうして？

这样做的时候，有一种不存在的方案被考虑了，这种方案就是

$a\to u$ + 1(不存在的边) + $a\to u$

但是就全局来看，我们其实考虑过

$a\to u$ + $a\to u$

所以那种不合法答案实际上会被更新掉

代码

#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;
}