AGC 002 D Stamp Rally - 分块

显然可以整体二分加可持久化并查集但是注意到撤销不需要可持久化可以直接按秩合并;其次注意到二分之后可以直接kruskal重构树;然而我都不想写,所以就随手敲了一个分块,然后精细的实现可以做到空间线性。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cmath>
#define gc getchar()
#define N 100010
#define cpy(a,b,n) memcpy(a,b,sizeof(int)*(n+1))
#define debug(x) cerr<<#x<<"="<<x
#define sp <<" "
#define ln <<endl
using namespace std;
inline int inn()
{
    int x,ch;while((ch=gc)<'0'||ch>'9');
    x=ch^'0';while((ch=gc)>='0'&&ch<='9')
        x=(x<<1)+(x<<3)+(ch^'0');return x;
}
int fa[2][N],sz[2][N],u[N],v[N],x[N],y[N],z[N],pre[N],nxt[N],ans[N];
vector<int> ql;
inline int findf(int *fa,int x)
{
    int fx=x,y;while(fx^fa[fx]) fx=fa[fx];
    while(x^fx) y=fa[x],fa[x]=fx,x=y;return fx;
}
inline int merge_union(int *fa,int *sz,int x,int y)
{
    int fx=findf(fa,x),fy=findf(fa,y);
    if(fx!=fy) fa[fx]=fy,sz[fy]+=sz[fx];
    return 0;
}
inline int ok(int *fa,int *sz,int x,int y,int z)
{
    int fx=findf(fa,x),fy=findf(fa,y);
    if(fx==fy) return sz[fx]>=z;
    return sz[fx]+sz[fy]>=z;
}
int main()
{
    int n=inn(),m=inn(),s=sqrt(m+0.5),bc=(m-1)/s+1;
    for(int i=1;i<=m;i++) u[i]=inn(),v[i]=inn();
    int q=inn(),las=0,now=1;
    for(int i=1;i<=n;i++) fa[las][i]=i,sz[las][i]=1;
    for(int i=1;i<=q;i++) x[i]=inn(),y[i]=inn(),z[i]=inn();
    for(int i=1;i<=q;i++) pre[i]=i-1,nxt[i]=i+1;
    nxt[0]=1,pre[q+1]=q;
    for(int i=1,L,R;i<=bc;i++)
    {
        L=(i-1)*s+1,R=min(m,i*s),cpy(fa[now],fa[las],n),cpy(sz[now],sz[las],n);
        for(int j=L;j<=R;j++) merge_union(fa[now],sz[now],u[j],v[j]);ql.clear();
        for(int j=nxt[0];j<=q;j=nxt[j])
            if(ok(fa[now],sz[now],x[j],y[j],z[j]))
                ql.push_back(j),j[nxt][pre]=j[pre],j[pre][nxt]=j[nxt];
        for(int j=L;j<=R;j++)
        {
            merge_union(fa[las],sz[las],u[j],v[j]);
            for(int k=0,p;k<ql.size();k++)
            {
                if(ans[p=ql[k]]) continue;
                if(ok(fa[las],sz[las],x[p],y[p],z[p]))
                    ans[p]=j;
            }
        }
    }
    for(int i=1;i<=q;i++) printf("%d\n",ans[i]);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值