Problem
Solution
- 做模板使用,不证明(
不会证 。 - 原图中两点间的最小割为最小割树上两点间最小割,可通过树上倍增询问。
- 最小割树的建立:任意两点间的最小割将当前点集分成两部分,以此为基础分治,每次从当前点集中任取两点求最小割,将最小割作为Gomory-Hu Tree上两点间的边权,并根据最小割将当前点集分成两部分,s那边的与t那边的,边界就是当前点集中只有一个点。
- 相当于跑了n-1次dinic。
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn = 505,maxm = 1505,inf = 1e9;
int cnt = 1,last[maxn],ceng[maxn],cnt0=1,belong[maxn];
int cur[maxn],n,m,cnt1,f[maxn][20],dep[maxn],front[maxn],Q,gm[maxn][20];
struct edge{
int v,next,w;
}e[maxm<<2];
edge te[maxn<<1];
inline void add(int u,int v,int w)
{
e[++cnt] = {v,last[u],w};
last[u] = cnt;
e[++cnt] = {u,last[v],0};
last[v] = cnt;
}
inline void add_tree(int u,int v,int w)
{
te[++cnt1] = {v,front[u],w};
front[u] = cnt1;
te[++cnt1] = {u,front[v],w};
front[v] = cnt1;
}
bool bfs(int s,int t)
{
queue<int>q;
memset(ceng,0,sizeof(ceng));
ceng[s]=1;
q.push(s);
while(!q.empty())
{
int u=q.front();q.pop();
for(int i=last[u];i;i=e[i].next)
{
int v=e[i].v,w=e[i].w;
if(ceng[v]==0&&w>0)
{
ceng[v]=ceng[u]+1;
q.push(v);
}
}
}
return ceng[t];
}
int dfs(int u,int dis,int t)
{
if(u==t)
return dis;
int minf = dis;
for(int i=cur[u];i;i=e[i].next)
{
cur[u]=e[i].next;
int v=e[i].v,w=e[i].w;
if(ceng[v]==ceng[u]+1&&w!=0)
{
int di=dfs(v,min(dis,w),t);
if(di>0)
{
e[i].w-=di;
e[i^1].w+=di;
dis-=di;
}
}
}
return minf-dis;
}
int dinic(int s,int t)
{
int ans=0;
while(bfs(s,t))
{
for(int i=1;i<=n;i++)
cur[i]=last[i];
ans+=dfs(s,inf,t);
}
return ans;
}
void init()
{
for(int i=2;i<=cnt;i+=2)
{
e[i].w+=e[i^1].w;
e[i^1].w = 0;
}
}
void dfs0(int u,int fa,int from)
{
if(belong[u]==from)
belong[u] = cnt0;
for(int i=last[u];i;i=e[i].next)
{
int v=e[i].v,w=e[i].w;
if(v!=fa&&w>0&&e[i^1].w>0)
dfs0(v,u,from);
}
}
int num[maxn],tmp1[maxn],tmp2[maxn];
void build(int l,int r)
{
if(l==r) return ;
init();
int s = num[l],t = num[l+1];
int flow = dinic(s,t),num1=0,num2=0;
add_tree(s,t,flow);
for(int i=l;i<=r;i++)
{
if(!ceng[num[i]]) tmp1[++num1]=num[i];
else tmp2[++num2]=num[i];
}
for(int i=l;i<=l+num1-1;i++) num[i] = tmp1[i-l+1];
for(int i=l+num1;i<=r;i++) num[i] = tmp2[i-l-num1+1];
build(l,l+num1-1);
build(l+num1,r);
}
void dfs_tree(int u,int fa)
{
for(int i=1;i<=15;i++)
{
f[u][i] = f[f[u][i-1]][i-1];
gm[u][i] = min(gm[u][i-1],gm[f[u][i-1]][i-1]);
}
for(int i=front[u];i;i=te[i].next)
{
int v=te[i].v,w=te[i].w;
if(v!=fa)
{
dep[v] = dep[u]+1;
f[v][0] = u;
gm[v][0] = w;
dfs_tree(v,u);
}
}
}
int query(int u,int v)
{
int ans = inf;
if(dep[u]<dep[v]) swap(u,v);
for(int i=15;i>=0;i--)
{
if(dep[f[u][i]]>=dep[v])
{
ans = min(ans,gm[u][i]);
u = f[u][i];
}
}
for(int i=15;i>=0;i--)
{
if(f[u][i]!=f[v][i])
{
ans=min(ans,min(gm[u][i],gm[v][i]));
u = f[u][i];v = f[v][i];
}
}
if(u!=v)
ans=min(ans,min(gm[u][0],gm[v][0]));
return ans;
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1,u,v,w;i<=m;i++)
{
scanf("%d%d%d",&u,&v,&w);
add(u,v,w);
add(v,u,w);
}
for(int i=1;i<=n;i++)
num[i] = i;
build(1,n);
dep[1] = 1;
dfs_tree(1,0);
scanf("%d",&Q);
for(int i=1,u,v;i<=Q;i++)
{
scanf("%d%d",&u,&v);
cout<<query(u,v)<<endl;
}
return 0;
}