bzoj3514Codechef MARCH14 GERALD07加强版
- #include <iostream>
- #include <cstdio>
- #include <cstring>
- #include <cmath>
- #include <algorithm>
- #define maxn 400005
- #define maxm 200005
- #define maxk 4000005
- using namespace std;
- int n,m,k,type,ans,Ans,size,fa[maxn],son[maxn][],sm[maxn],Val[maxm],val[maxn],sum[maxk],root[maxm],lc[maxk],rc[maxk];
- bool rev[maxn];
- struct date{
- int u,v;
- }wi[maxm];
- struct note{
- bool which(int x){
- return son[fa[x]][]==x;
- }
- bool isroot(int x){
- return son[fa[x]][]!=x&&son[fa[x]][]!=x;
- }
- void update(int x){
- sm[x]=val[x];
- if (son[x][]) sm[x]=min(sm[x],sm[son[x][]]);
- if (son[x][]) sm[x]=min(sm[x],sm[son[x][]]);
- }
- void pushdown(int x){
- if (!rev[x]) return;
- rev[x]^=,swap(son[x][],son[x][]);
- if (son[x][]) rev[son[x][]]^=;
- if (son[x][]) rev[son[x][]]^=;
- }
- void relax(int x){
- if (!isroot(x)) relax(fa[x]);
- pushdown(x);
- }
- void rotata(int x){
- int y=fa[x],d=which(x),dd=which(y);
- if (!isroot(y)) son[fa[y]][dd]=x; fa[x]=fa[y];
- fa[son[x][d^]]=y,son[y][d]=son[x][d^];
- fa[y]=x,son[x][d^]=y;
- update(y);
- }
- void splay(int x){
- relax(x);
- while (!isroot(x)){
- if (isroot(fa[x])) rotata(x);
- else if (which(x)==which(fa[x])) rotata(fa[x]),rotata(x);
- else rotata(x),rotata(x);
- }
- update(x);
- }
- void access(int x){
- for (int p=;x;x=fa[x]){
- splay(x);
- son[x][]=p;
- p=x;
- update(x);
- }
- }
- void make_root(int x){
- access(x);
- splay(x);
- rev[x]^=;
- }
- void link(int x,int y){
- make_root(x);
- fa[x]=y;
- }
- void cut(int x,int y){
- make_root(x);
- access(y);
- splay(y);
- son[y][]=fa[x]=;
- update(y);
- }
- void split(int x,int y){
- make_root(x);
- access(y);
- splay(y);
- }
- int find_root(int x){
- access(x);
- splay(x);
- while (son[x][]) x=son[x][];
- return x;
- }
- int find(int x,int y){
- split(x,y);
- return sm[y];
- }
- }lct;
- void insert(int &k,int p,int l,int r,int x){
- k=++size,sum[k]=sum[p]+;
- int mid=(l+r)/;
- if (l==r) return;
- if (x<=mid) rc[k]=rc[p],insert(lc[k],lc[p],l,mid,x);
- else lc[k]=lc[p],insert(rc[k],rc[p],mid+,r,x);
- }
- void insert(int id,int x){
- insert(root[id],root[id-],,m+,x);
- }
- void Query(int k1,int k2,int l,int r,int x,int y){
- if (!k1&&!k2) return;
- if (l>=x&&r<=y){
- Ans+=(sum[k2]-sum[k1]);
- return;
- }
- int mid=(l+r)/;
- if (x<=mid) Query(lc[k1],lc[k2],l,mid,x,y);
- if (y>mid) Query(rc[k1],rc[k2],mid+,r,x,y);
- }
- int query(int x,int y){
- Ans=;
- int u=root[x-],v=root[y],l=,r=m+,mid;
- Query(u,v,l,r,,x-);
- return Ans;
- }
- int main(){
- int u,v,t1,t2,temp;
- scanf("%d%d%d%d",&n,&m,&k,&type);
- memset(sum,,sizeof(sum));
- memset(Val,,sizeof(Val));
- memset(fa,,sizeof(fa));
- memset(son,,sizeof(son));
- memset(rev,,sizeof(rev));
- for (int i=;i<=m;i++) scanf("%d%d",&wi[i].u,&wi[i].v);
- for (int i=;i<=n;i++) val[i]=sm[i]=m+;
- for (int i=;i<=m;i++) val[n+i]=sm[n+i]=i;
- for (int i=;i<=m;i++){
- u=wi[i].u,v=wi[i].v;
- if (u==v) Val[i]=m+;
- if (u==v) continue;
- if (lct.find_root(u)!=lct.find_root(v)){
- Val[i]=;
- lct.link(n+i,u),lct.link(n+i,v);
- }else{
- temp=lct.find(u,v);
- Val[i]=temp,temp+=n;
- lct.cut(temp,wi[temp-n].u),lct.cut(temp,wi[temp-n].v);
- lct.link(n+i,u),lct.link(n+i,v);
- }
- }
- memset(root,,sizeof(root));
- for (int i=;i<=m;i++){
- insert(i,Val[i]);
- }
- ans=;
- for (int i=;i<=k;i++){
- scanf("%d%d",&u,&v);
- if (type==) u^=ans,v^=ans;
- if (u>v) swap(u,v);
- ans=n-query(u,v);
- printf("%d\n",ans);
- }
- return ;
- }
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3514
题目大意:N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数。
做法:如果要求整个图中联通块的个数,记加入该边后不形成环的边数记为x,则答案为n-x,这题问的是加入编号为[l~r]中的边后联通块的个数,答案便稍微变化了一下,记能使联通块数目-1的边集为{S},答案为n-|S|,问题便成为了:有多少条边属于{S},仔细想想,可以先预处理一个数组val[i],怎么预处理呢?可以按编号依次加入边,有两种情况:*1.如果加入的这条边不形成环,则这条边的权值赋值为0,并加入这条边;*2.加入这条边形成环,记环上编号最小的边的编号为y,这加入的这条边权值为y,并删去编号最小的那条边,加入该边。这个过程用lct模拟即可,比较基础的操作。
得出每条边的权值val[i]后,对于一个询问[L~R],如果一条边的val<l,则该边属于{S},问题便简化为一个数组val,问在区间L~R中权值在0~L-1的个数,显而易见,可持久化线段树轻松搞定。问题得以圆满解决。(Lct+可持久化线段树)
bzoj3514Codechef MARCH14 GERALD07加强版的更多相关文章
- 【BZOJ-3514】Codechef MARCH14 GERALD07加强版 LinkCutTree + 主席树
3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 1288 Solved: 490 ...
- BZOJ 3514: Codechef MARCH14 GERALD07加强版( LCT + 主席树 )
从左到右加边, 假如+的边e形成环, 那么记下这个环上最早加入的边_e, 当且仅当询问区间的左端点> _e加入的时间, e对答案有贡献(脑补一下). 然后一开始是N个连通块, 假如有x条边有贡献 ...
- 【LCT+主席树】BZOJ3514 Codechef MARCH14 GERALD07加强版
3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 2023 Solved: 778 ...
- BZOJ 3514: Codechef MARCH14 GERALD07加强版 [LCT 主席树 kruskal]
3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 1312 Solved: 501 ...
- BZOJ_3514_Codechef MARCH14 GERALD07加强版_主席树+LCT
BZOJ_3514_Codechef MARCH14 GERALD07加强版_主席树+LCT Description N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数. I ...
- [BZOJ 3514]Codechef MARCH14 GERALD07加强版 (CHEF AND GRAPH QUERIES)
[BZOJ3514] Codechef MARCH14 GERALD07加强版 (CHEF AND GRAPH QUERIES) 题意 \(N\) 个点 \(M\) 条边的无向图,\(K\) 次询问保 ...
- [BZOJ3514]CodeChef MARCH14 GERALD07加强版(LCT+主席树)
3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 2177 Solved: 834 ...
- bzoj3514 Codechef MARCH14 GERALD07加强版 lct预处理+主席树
Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 1951 Solved: 746[Submi ...
- BZOJ3514 : Codechef MARCH14 GERALD07加强版
以边编号为权值 用Link-cut Tree维护最大生成树 对于新加的第i条边(u,v) a[i]表示当a[i]这条边加入后连通块个数会减少 若u==v则a[i]=m 若u与v不连通则连上,a[i]= ...
随机推荐
- 准备NOIP2017 最长公共子序列(模版)
一些概念: (1)子序列: 一个序列A = a1,a2,--an,中任意删除若干项,剩余的序列叫做A的一个子序列.也可以认为是从序列A按原顺序保留任意若干项得到的序列.例如: 对序列 1,3,5, ...
- BZOJ 1030 【JSOI2007】 文本生成器
Description JSOI交给队员ZYX一个任务,编制一个称之为"文本生成器"的电脑软件:该软件的使用者是一些低幼人群,他们现在使用的是GW文本生成器v6版.该软件可以随机生 ...
- LeetCode 01 Two Sum swift
class TwoSum { func sumTow(nums: [Int], target: Int)->[Int]{ ,]; ;x<nums.count;x++){ ;y<num ...
- 【转】$_POST 与 php://input的区别分析
$data = file_get_contents("php://input"); php://input 是个可以访问请求的原始数据的只读流. POST 请求的情况下,最好使用 ...
- noi题库(noi.openjudge.cn) 1.5编程基础之循环控制T36——T45
T36 计算多项式的值 描述 假定多项式的形式为xn+xn-1+-+x2+x+1,请计算给定单精度浮点数x和正整数n值的情况下这个多项式的值. 输入 输入仅一行,包括x和n,用单个空格隔开.x在flo ...
- string to char* and char* to string 玩转 String 和 Char*
char 类型是c语言中常见的一个数据类型,string是c++中的一个,它的定义为 Strings are objects that represent sequences of character ...
- WPF 3D模型 3D场景
1.首先得说明的是这并不是真正的3D,模型被导出为一系列的单个图片,例如一个3D户型图,以某个视角旋转360°,渲染出一系列连续的单个图片文件. 2.在Image.MouseMove事件中添加相应代码 ...
- 理解AngularJS生命周期:利用ng-repeat动态解析自定义directive
ng-repeat是AngularJS中一个非常重要和有意思的directive,常见的用法之一是将某种自定义directive和ng-repeat一起使用,循环地来渲染开发者所需要的组件.比如现在有 ...
- CI(CodeIgniter)框架入门教程——第二课 初始MVC
本文转载自:http://www.softeng.cn/?p=53 今天的主要内容是,使用CodeIgniter框架完整的MVC内容来做一个简单的计算器,通过这个计算器,让大家能够体会到我在第一节课中 ...
- Realm Java的学习、应用、总结
从React Native珠三角沙龙会议了解到Realm这个开源库,然后开始学习.理解和使用Realm.Realm是跨平台.支持多种主流语言,这里主要是对Realm Java结合实际项目的一些情况进行 ...