传送门
这定是极其蛋疼的题目。
考虑分治询问。
我们可以先求出当询问在l到r之间是,必定在MST上的边。
(未确定权值的边即使取到-inf也在MST上的边)
同时,我们可以缩点。
然后删除必定不在MST上的边。
(未确定权值的边即使取到-inf也在MST上的边)
分析时间复杂度:
点数最多为r-l+1(其他都被缩掉了)
边数最多为2*(r-l+1)(MST上r-l+1条和未确定的r-l+1条)
每一次分治时间复杂度O(mlogm)(m表示边数)
然后就是T(n)=2*T(n/2)+O(nlogn)
显然是O(nlog^2n)的。
#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define M 50505
#define N 20202
#define inf 0x3f3f3f3f
#define ll long long
using namespace std;
struct edge{int x,y,w,pos;}e[25][M],d[M],t[M];
struct query{int x,y;}q[M];
int sum[N],f[N],sz[N],c[M],a[M];
int n,m,Q,x,y;
ll ans[M];
inline bool cmp(const edge &a,const edge &b){
return a.w<b.w;
}
inline int get(int x){
return f[x]==x?x:f[x]=get(f[x]);
}
inline void merge(int x,int y){
if (sz[x]>sz[y]) swap(x,y);
sz[y]+=sz[x]; f[x]=y;
}
inline void clear(int tot){
for (int i=1;i<=tot;i++){
f[d[i].x]=d[i].x;
f[d[i].y]=d[i].y;
sz[d[i].x]=sz[d[i].y]=1;
}
}
inline void cont(int &tot,ll &cnt){
int tmp=0;
clear(tot);
sort(d+1,d+tot+1,cmp);
for (int i=1;i<=tot;i++)
if (get(d[i].x)!=get(d[i].y))
merge(f[d[i].x],f[d[i].y]),t[++tmp]=d[i];
for (int i=1;i<=tmp;i++){
f[t[i].x]=t[i].x;
f[t[i].y]=t[i].y;
sz[t[i].x]=sz[t[i].y]=1;
}
for (int i=1;i<=tmp;i++)
if (t[i].w!=-inf&&get(t[i].x)!=get(t[i].y))
merge(f[t[i].x],f[t[i].y]),cnt+=t[i].w;
tmp=0;
for (int i=1;i<=tot;i++)
if (get(d[i].x)!=get(d[i].y)){
t[++tmp]=d[i];
c[d[i].pos]=tmp;
t[tmp].x=f[d[i].x];
t[tmp].y=f[d[i].y];
}
for (int i=1;i<=tmp;i++) d[i]=t[i];
tot=tmp;
}
inline void reduct(int &tot){
int tmp=0;
clear(tot);
sort(d+1,d+tot+1,cmp);
for (int i=1;i<=tot;i++)
if (get(d[i].x)!=get(d[i].y)){
merge(f[d[i].x],f[d[i].y]);
t[++tmp]=d[i];
c[d[i].pos]=tmp;
}
else if (d[i].w==inf){
t[++tmp]=d[i];
c[d[i].pos]=tmp;
}
for (int i=1;i<=tmp;i++) d[i]=t[i];
tot=tmp;
}
void solve(int l,int r,int now,ll cnt){
int tot=sum[now];
if (l==r) a[q[l].x]=q[l].y;
for (int i=1;i<=tot;i++)
e[now][i].w=a[e[now][i].pos];
for (int i=1;i<=tot;i++)
d[i]=e[now][i],c[d[i].pos]=i;
if (l==r){
ans[l]=cnt;
clear(tot);
sort(d+1,d+tot+1,cmp);
for (int i=1;i<=tot;i++)
if (get(d[i].x)!=get(d[i].y))
merge(f[d[i].x],f[d[i].y]),ans[l]+=d[i].w;
return;
}
for (int i=l;i<=r;i++)
d[c[q[i].x]].w=-inf;
cont(tot,cnt);
for (int i=l;i<=r;i++)
d[c[q[i].x]].w=inf;
reduct(tot);
for (int i=1;i<=tot;i++)
e[now+1][i]=d[i];
sum[now+1]=tot;
int mid=(l+r)/2;
solve(l,mid,now+1,cnt);
solve(mid+1,r,now+1,cnt);
}
int main(){
scanf("%d%d%d",&n,&m,&Q);
for (int i=1;i<=m;i++){
scanf("%d%d%d",&x,&y,&a[i]);
e[0][i].x=x; e[0][i].y=y;
e[0][i].w=a[i]; e[0][i].pos=i;
}
for (int i=1;i<=Q;i++)
scanf("%d%d",&q[i].x,&q[i].y);
sum[0]=m;
solve(1,Q,0,0);
for (int i=1;i<=Q;i++)
printf("%lld\n",ans[i]);
}