首页 > 代码库 > BZOJ3514 Codechef MARCH14 GERALD07加强版
BZOJ3514 Codechef MARCH14 GERALD07加强版
Description
N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数。
Input
第一行四个整数N、M、K、type,代表点数、边数、询问数以及询问是否加密。
接下来M行,代表图中的每条边。
接下来K行,每行两个整数L、R代表一组询问。对于type=0的测试点,读入的L和R即为询问的L、R;对于type=1的测试点,每组询问的L、R应为L xor lastans和R xor lastans。
Output
K行每行一个整数代表该组询问的联通块个数。
Sample Input
3 5 4 0
1 3
1 2
2 1
3 2
2 2
2 3
1 5
5 5
1 2
1 3
1 2
2 1
3 2
2 2
2 3
1 5
5 5
1 2
Sample Output
2
1
3
1
1
3
1
HINT
对于100%的数据,1≤N、M、K≤200,000。
正解:LCT加主席树
解题报告:当一条边加进来的时候,如果产生了环,我们把环上最小的一条边删掉并记录,如果没有环,记录为0,到时候进行区间查询的时候,我们可以访问这个数组,如果被删掉的边小于l(0也包括),说明这条边联通了两个本来不连通的区间,然后答案-1即可,这样题目就变成了给定一个区间,求小于这个数的元素的数量,删加边显然用lct做。
1 #include <iostream> 2 #include <iomanip> 3 #include <cstdlib> 4 #include <cstdio> 5 #include <cmath> 6 #include <string> 7 #include <cstring> 8 #include <algorithm> 9 #define MIN(a,b) (a<b?a:b) 10 #define RG register 11 const int N = 500000; 12 const int M = 15000000; 13 const int inf = 2147483641; 14 15 using namespace std; 16 17 int gi(){ 18 RG char ch=getchar();RG int x=0; 19 while(ch<‘0‘ || ch>‘9‘) ch=getchar(); 20 while(ch>=‘0‘ && ch<=‘9‘) x=x*10+ch-‘0‘,ch=getchar(); 21 return x; 22 } 23 24 int fa[N],c[2][N],rev[N],st[N],lt[N],mn[N],a[N],b[N],rt[N]; 25 int n,m,k,type,cnt; 26 27 struct date{ 28 int l,r,s; 29 }tr[M]; 30 31 int isroot(int x){ 32 return c[0][fa[x]]!=x && c[1][fa[x]]!=x; 33 } 34 35 void pushdown(int x){ 36 if (rev[x]==0) return; 37 rev[x]^=1,rev[c[0][x]]^=1,rev[c[1][x]]^=1; 38 swap(c[0][x],c[1][x]); 39 return; 40 } 41 42 void update(int x){ 43 mn[x]=x<=n?inf:x; 44 mn[x]=MIN(mn[x],MIN(mn[c[0][x]],mn[c[1][x]])); 45 return; 46 } 47 48 void rotate(int x){ 49 int l,r,y=fa[x],z=fa[y]; 50 if (c[0][y]==x) l=0; 51 else l=1;r=l^1; 52 if (!isroot(y)) 53 if (c[0][z]==y) c[0][z]=x; 54 else c[1][z]=x; 55 fa[x]=z,fa[y]=x,fa[c[r][x]]=y; 56 c[l][y]=c[r][x],c[r][x]=y; 57 update(y),update(x); 58 return; 59 } 60 61 void splay(int x){ 62 int tot=0;st[++tot]=x; 63 for (RG int i=x; !isroot(i); i=fa[i]) st[++tot]=fa[i]; 64 for (RG int i=tot; i; --i) pushdown(st[i]); 65 while(!isroot(x)){ 66 int y=fa[x],z=fa[y]; 67 if (!isroot(y)) 68 if (c[0][y]==x ^ c[0][z]==y) rotate(x); 69 else rotate(y); 70 rotate(x); 71 } 72 return; 73 } 74 75 void access(int x){ 76 int t=0; 77 while(x){ 78 splay(x); 79 c[1][x]=t;update(x); 80 t=x,x=fa[x]; 81 } 82 return; 83 } 84 85 void rever(int x){ 86 access(x),splay(x),rev[x]^=1; 87 return; 88 } 89 90 void link(int l,int r){ 91 rever(l),fa[l]=r; 92 return; 93 } 94 95 void cut(int l,int r){ 96 rever(l),access(r),splay(r); 97 c[0][r]=fa[l]=0; 98 return; 99 } 100 101 int query(int x){ 102 access(x),splay(x); 103 while(c[0][x]) x=c[0][x]; 104 return x; 105 } 106 107 void build(int &t,int l,int r,int s){ 108 tr[++cnt]=tr[t];t=cnt; 109 ++tr[t].s; 110 if (l==r) return; 111 int mid=(l+r)>>1; 112 if (s<=mid) build(tr[t].l,l,mid,s); 113 else build(tr[t].r,mid+1,r,s); 114 return; 115 } 116 117 int query(int x,int l,int r,int ql,int qr){ 118 if (ql>qr) return 0; 119 if (ql<=l && r<=qr) return tr[x].s; 120 int ans=0,mid=(l+r)>>1; 121 if (ql<=mid) ans+=query(tr[x].l,l,mid,ql,qr); 122 if (qr>mid) ans+=query(tr[x].r,mid+1,r,ql,qr); 123 return ans; 124 } 125 126 int main(){ 127 freopen("splay.in","r",stdin); 128 freopen("splay.out","w",stdout); 129 n=gi(),m=gi(),k=gi(),type=gi(); 130 for (RG int i=0; i<=n; ++i) mn[i]=inf; 131 for (RG int i=1; i<=m; ++i){ 132 int l=gi(),r=gi(); 133 if (l==r) {lt[i]=m;continue;} 134 a[i]=l,b[i]=r; 135 if (query(l)==query(r)){ 136 rever(l),access(r),splay(r); 137 int z=a[mn[r]-n],y=b[mn[r]-n],k=mn[r]; 138 lt[i]=mn[r]-n; 139 cut(z,k),cut(k,y); 140 } 141 link(l,i+n),link(i+n,r); 142 } 143 for (RG int i=1; i<=m; ++i){ 144 int g=lt[i]; 145 rt[i]=rt[i-1]; 146 build(rt[i],0,m,g); 147 } 148 int jl=0; 149 for (RG int i=1; i<=k; ++i){ 150 int l=gi()^jl,r=gi()^jl,ans=n; 151 ans-=query(rt[r],0,m,0,l-1)-query(rt[l-1],0,m,0,l-1); 152 printf("%d\n",ans); 153 if (type) jl=ans; 154 } 155 return 0; 156 }
BZOJ3514 Codechef MARCH14 GERALD07加强版
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。