bzoj2001: [Hnoi2010]City 城市建设

传送门
这定是极其蛋疼的题目。
考虑分治询问。
我们可以先求出当询问在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]);
}