好题+数据结构神题+感人肺腑pascal被卡系列,我下面的代码几乎写到最优
可怎耐bzoj上pascal开的是O1,c++开的是O2,这怎么可能跑得过!!!
还是说说方法吧,这是一道算贡献的好题,因为我们不可能把边加进去依次算连通块个数
我们考虑在区间[l,r]中每条边i是否能使两个连通块变成1个,什么样的边可以呢?
显然,这条边端点所在的连通块在i加入之前不连通或者是使这两个连通块相连的边号<l
这样算法就呼之欲出了,我们首先依次把边i加入,如果构成环,则弹出环上编号最小的边并记录,否则记为0
这里我们可以用link cut tree解决
下面每次询问,我们只要求[l,r]有多少个数小于l即可
这题有是否强制在线两种情况,在线可以直接rush,但是我第一次TLE之后又写了一遍离线,发现比原来快但还是远远慢于只在线的c++,实在没办法只好弃疗,说一下两种情况
在线:显然这是主席树的经典问题
离线:注意这个个数是满足区间减法的,于是我们可以拆成两个求前缀询问并排序(注意这里用计数排序才有显著的提升),然后按边编号的顺序依次加入然后用树状数组统计,具体见程序

 const inf=;
type node=record
l,r,s:longint;
end;
way=record
s,e:longint;
end;
point=record
be,op,y:longint;
end;
xxx=record
po,next:longint;
end; var son:array[..,..] of longint;
q,w,v,fa:array[..] of longint;
qq:array[..] of point;
ee:array[..] of xxx;
rev:array[..] of boolean;
tree:array[..*] of node;
e:array[..] of way;
p,f,a,h,anss:array[..] of longint;
len,j,k1,k2,n,m,k,ty,t,i,ans,l,r,mm:longint; function min(a,b:longint):longint;
begin
if w[a]>w[b] then exit(b) else exit(a);
end; function getf(x:longint):longint;
begin
if f[x]<>x then f[x]:=getf(f[x]);
exit(f[x]);
end; procedure swap(var a,b:longint);
var c:longint;
begin
c:=a;
a:=b;
b:=c;
end; function root(x:longint):boolean;
begin
exit((son[fa[x],]<>x) and (son[fa[x],]<>x));
end; procedure update(x:longint);
begin
v[x]:=min(x,min(v[son[x,]],v[son[x,]]));
end; procedure push(x:longint);
begin
if rev[x] then
begin
swap(son[x,],son[x,]);
rev[son[x,]]:=not rev[son[x,]];
rev[son[x,]]:=not rev[son[x,]];
rev[x]:=false;
end;
end; procedure rotate(x,w:longint);
var y:longint;
begin
y:=fa[x];
if not root(y) then
begin
if son[fa[y],]=y then son[fa[y],]:=x
else son[fa[y],]:=x;
end;
fa[x]:=fa[y];
son[y,-w]:=son[x,w];
if son[x,w]<> then fa[son[x,w]]:=y;
son[x,w]:=y;
fa[y]:=x;
update(y);
end; procedure splay(x:longint);
var i,s,y:longint;
fl:boolean;
begin
i:=x;
s:=;
while not root(i) do
begin
inc(s);
q[s]:=i;
i:=fa[i];
end;
inc(s);
q[s]:=i;
for i:=s downto do
push(q[i]);
if s= then exit;
fl:=true;
while fl do
begin
y:=fa[x];
if y=q[s] then
begin
if son[y,]=x then rotate(x,)
else rotate(x,);
fl:=false;
end
else begin
if fa[y]=q[s] then fl:=false;
if son[fa[y],]=y then
begin
if son[y,]=x then rotate(y,)
else rotate(x,);
rotate(x,);
end
else begin
if son[y,]=x then rotate(x,)
else rotate(y,);
rotate(x,);
end;
end;
end;
update(x);
end; procedure access(x:longint);
var y:longint;
begin
y:=;
repeat
splay(x);
son[x,]:=y;
update(x);
y:=x;
x:=fa[x];
until x=;
end; procedure makeroot(x:longint);
begin
access(x);
splay(x);
rev[x]:=not rev[x];
end; procedure link(x,y:longint);
begin
makeroot(x);
fa[x]:=y;
end; procedure path(x,y:longint);
begin
makeroot(x);
access(y);
splay(y);
end; procedure cut(x,y:longint);
begin
path(x,y);
son[y,]:=;
fa[x]:=;
end; function build(l,r:longint):longint;
var m,q:longint;
begin
inc(t);
if l=r then exit(t)
else begin
m:=(l+r) shr ;
q:=t;
tree[q].l:=build(l,m);
tree[q].r:=build(m+,r);
exit(q);
end;
end; function add(l,r,last,x:longint):longint;
var m,q:longint;
begin
inc(t);
if l=r then
begin
tree[t].s:=tree[last].s+;
exit(t);
end
else begin
m:=(l+r) shr ;
q:=t;
if x<=m then
begin
tree[q].r:=tree[last].r;
tree[q].l:=add(l,m,tree[last].l,x);
end
else begin
tree[q].l:=tree[last].l;
tree[q].r:=add(m+,r,tree[last].r,x);
end;
tree[q].s:=tree[tree[q].l].s+tree[tree[q].r].s;
exit(q);
end;
end; function ask(l,r,x,y:longint):longint;
var m,s:longint;
begin
if l=r then exit(tree[y].s-tree[x].s)
else begin
m:=(l+r) shr ;
if k1<=m then exit(ask(l,m,tree[x].l,tree[y].l))
else begin
s:=tree[tree[y].l].s-tree[tree[x].l].s;
exit(s+ask(m+,r,tree[x].r,tree[y].r));
end;
end;
end; procedure cha(var a,b:point);
var c:point;
begin
c:=a;
a:=b;
b:=c;
end; function lowbit(x:longint):longint;
begin
exit(x and (-x));
end; procedure jia(x,y:longint);
begin
inc(len);
ee[len].po:=y;
ee[len].next:=p[x];
p[x]:=len;
end; procedure ins(x:longint);
begin
while x<=m+ do
begin
inc(h[x]);
x:=x+lowbit(x);
end;
end; function an(x:longint):longint;
begin
an:=;
while x> do
begin
an:=an+h[x];
x:=x-lowbit(x);
end;
end; begin
readln(n,m,k,ty);
for i:= to n do
begin
f[i]:=i;
w[i]:=inf;
end;
t:=n;
w[]:=inf;
for i:= to m do
begin
readln(e[i].s,e[i].e);
if e[i].s=e[i].e then
begin
a[i]:=i;
continue;
end;
k1:=getf(e[i].s);
k2:=getf(e[i].e);
if k1=k2 then //并查集判连通
begin
path(e[i].s,e[i].e);
l:=v[e[i].e];
a[i]:=w[l];
cut(e[w[l]].s,l);
cut(e[w[l]].e,l);
end
else f[k1]:=k2;
inc(t); w[t]:=i;
link(e[i].s,t);
link(e[i].e,t);
end;
if ty= then //在线做法
begin
t:=;
h[]:=build(,m);
for i:= to m do
begin
h[i]:=add(,m,h[i-],a[i]);
end;
ans:=;
for i:= to k do
begin
readln(l,r);
l:=l xor ans;
r:=r xor ans;
if l>r then swap(l,r);
k1:=l-;
ans:=n-ask(,m,h[l-],h[r]);
writeln(ans);
end;
end
else begin //离线做法
t:=;
for i:= to k do
begin
readln(l,r);
inc(t);
qq[t].op:=-;
qq[t].y:=l;
qq[t].be:=i;
jia(r,t);
inc(t);
qq[t].op:=;
jia(l-,t);
qq[t].y:=l;
qq[t].be:=i;
end;
j:=;
a[m+]:=inf;
for i:= to m do
begin
while j<=i do
begin
ins(a[j]+); //树状数组编号从1开始,所以+
inc(j);
end;
l:=p[i];
while l<> do
begin
r:=ee[l].po;
inc(anss[qq[r].be],qq[r].op*an(qq[r].y));
l:=ee[l].next;
end;
end;
for i:= to k do
writeln(anss[i]+n);
end;
end.

bzoj3514的更多相关文章

  1. BZOJ3514 Codechef MARCH14 GERALD07加强版 LCT

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ3514 题意概括 N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数. N ...

  2. [BZOJ3514]CodeChef MARCH14 GERALD07加强版(LCT+主席树)

    3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 2177  Solved: 834 ...

  3. 【BZOJ3514】Codechef MARCH14 GERALD07加强版(LCT)

    题意:N个点M条边的无向图,q次询问保留图中编号在[l,r]的边的时候图中的联通块个数. 询问加密,强制在线 n,m,q<=200000 题意:RYZ作业 以下转载自hzwer http://h ...

  4. 【BZOJ3514】Codechef MARCH14 GERALD07加强版(LCT_主席树)

    题目: BZOJ3514 分析: 看到这题真的是一脸懵逼无从下手,只好膜题解.看到「森林的联通块数 = 点数 - 边数」这一句话就立刻什么都会了 QAQ . 这题最重要的就是意识到上面那个式子(正确性 ...

  5. 【BZOJ-3514】Codechef MARCH14 GERALD07加强版 LinkCutTree + 主席树

    3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 1288  Solved: 490 ...

  6. BZOJ3514 : Codechef MARCH14 GERALD07加强版

    以边编号为权值 用Link-cut Tree维护最大生成树 对于新加的第i条边(u,v) a[i]表示当a[i]这条边加入后连通块个数会减少 若u==v则a[i]=m 若u与v不连通则连上,a[i]= ...

  7. 【bzoj3514】Codechef MARCH14 GERALD07加强版

    hzwer上少有的几道需要建一下模的 要不是有这么几道题 我都觉得lct只会考裸题了呢 题解看hzwer吧 http://hzwer.com/4358.html 唯一蛋疼的就是为了处理0这个呵呵的位置 ...

  8. 【LCT+主席树】BZOJ3514 Codechef MARCH14 GERALD07加强版

    3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 2023  Solved: 778 ...

  9. bzoj3514(LCT+主席树)

    题目描述 N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数. 题解 对于一个截止时间来说,越晚的变越好. 所以我们可以维护一颗以边的序号为关键字的最大生成树,然后用主席树维 ...

随机推荐

  1. Update msi using vbscript

    参考: http://stackoverflow.com/questions/1609250/how-do-i-add-update-a-property-inside-an-msi-from-the ...

  2. Generate GUID using vbscript

    在 .msi 中 的 Component table,查看 ComponentId 列,是一个16进制数的字符串, 用 InstallShield IDE 添加一个 component ,Compon ...

  3. 存储占用:Memory Map 汉化去广告版

    转载说明 本篇文章可能已经更新,最新文章请转:http://www.sollyu.com/storage-occupancy-memory-map-localization-to-billboards ...

  4. PHP连接sqlserver的两种方法,向sqlserver2000中写入数据,中文乱码

    项目环境是php5.3.28 项目用的ThinkPHP3.2.3  已经mysql5.5数据库,要和另一个项目对接,需要连接sqlsever2000数据库进行一些操作. 第一种用php自带扩展连接数据 ...

  5. 七牛云覆盖上传 php

    使用七牛云过程中遇到了需要上传覆盖的情况,最终解决,分享给大家. 七牛云sdk上传示例中是这样写的 <?php require_once 'path_to_sdk/vendor/autoload ...

  6. Android开发系列之搭建开发环境

    接触Android好久了,记得09年刚在中国大陆有点苗头的时候,我就知道了google有个Android,它是智能机操作系统.后来在Android出1.5版本之后,我第一时间下载了eclipse开发工 ...

  7. Cassandra1.2文档学习(5)—— Snitch

    参考资料:http://www.datastax.com/documentation/cassandra/1.2/webhelp/index.html#cassandra/architecture/a ...

  8. hdu 3234 Exclusive-OR

    Exclusive-OR 题意:输入n个点和Q次操作(1 <= n <= 20,000, 2 <= Q <= 40,000).操作和叙述的点标号k(0 < k < ...

  9. ios字符串操作

    string的操作应用 NSRange range = [self.general rangeOfString:@"."]; NSString *str = [self.gener ...

  10. [BZOJ 2212] [Poi2011] Tree Rotations 【线段树合并】

    题目链接:BZOJ - 2212 题目分析 子树 x 内的逆序对个数为 :x 左子树内的逆序对个数 + x 右子树内的逆序对个数 + 跨越 x 左子树与右子树的逆序对. 左右子树内部的逆序对与是否交换 ...