[Ahoi2005]LANE 航线规划
题目描述

输入
输出
样例输入
样例输出
首先,题意是求动态图的两点桥数。
将边删除会很麻烦,所以先将要删去的边去掉,形成最后的图。
求边双联通分量,缩点,形成一棵树(一定是一棵树)树边为关键边。
再根据数据加边,每加入一边,则形成一个环,则环上路不再为关键边,最后将答案逆序输出。
加入边后,两点之间的路径中的关键边可用树剖和线段树维护。
线段树中,1表是关键边,0表不是,维护和。每加入一边,求两点LCA,并将线段树相应区间归零(要做延迟表记)。查询就输出线段树相应区间的和
注意细节,代码有点长
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<map>
#include<set>
using namespace std;
struct Messi1
{
int next,to;
}edge1[],edge2[];
struct Messi2
{
int p,q,c;
}a[],b[];
int isbridge[],son[],size[],fa[],dep[],top[],pos[],tot,m;
int head1[],head2[],num1,num2,dfn[],low[],dfscnt,sccno[],n,scccnt,ans[];
int lazy[],c[],cx;
typedef pair<int,int> dou;//定义数对
set<dou> vis,vis2;//此处用set可节约内存,且方便。
void add1(int u,int v)//原图
{
num1++;
edge1[num1].next=head1[u];
edge1[num1].to=v;
head1[u]=num1;
}
void add2(int u,int v)//缩点后的图
{
num2++;
edge2[num2].next=head2[u];
edge2[num2].to=v;
head2[u]=num2;
}
void tarjan(int x,int fa)
{int i;
dfn[x]=low[x]=++dfscnt;
int child=;
for (i=head1[x]; i; i=edge1[i].next)
{
int v=edge1[i].to;
if (dfn[v]==)
{
child++;
tarjan(v,x);
low[x]=min(low[x],low[v]);
if (dfn[x]<low[v]) isbridge[i]=;
}
else if (v!=fa) low[x]=min(low[x],dfn[v]);
}
}
void dfs(int x,int v)
{int i;
//cout<<x<<' '<<v<<endl;
dfn[x]=;
sccno[x]=v;
for (i=head1[x]; i; i=edge1[i].next)
if (isbridge[i]==&&dfn[edge1[i].to]==)
{
dfs(edge1[i].to,v);
}
}
void make_tree()
{int i,j,x1,y1;
tarjan(,-);
memset(dfn,,sizeof(dfn));
for (i=; i<=n; i++)
if (dfn[i]==)
{
scccnt++;
dfs(i,scccnt);
}
for (i=; i<=n; i++)
{
for (j=head1[i]; j; j=edge1[j].next)
{
int v=edge1[j].to;
if (sccno[i]!=sccno[v])
{if (sccno[i]>sccno[v]) x1=sccno[v],y1=sccno[i];
else y1=sccno[v],x1=sccno[i];
if (vis2.find(make_pair(x1,y1))==vis2.end())
//注:set的find函数返回指针,不存在时反回最后节点子节点,
所以不能写!find()。用end函数反回最后节点子节点
{
add2(x1,y1);
add2(y1,x1);
vis2.insert(make_pair(x1,y1));
}
}
}
}
}
void dfs1(int u,int pa,int depth)
{
son[u]=;
size[u]=;
fa[u]=pa;
dep[u]=depth;
for (int j=head2[u]; j; j=edge2[j].next)
{
int v=edge2[j].to;
if (v!=pa)
{
dfs1(v,u,depth+);
size[u]+=size[v];
if (size[v]>size[son[u]]) son[u]=v;
}
}
}
void dfs2(int u,int tp)
{
top[u]=tp;
pos[u]=++tot;
if (son[u]) dfs2(son[u],tp);
for (int j=head2[u]; j; j=edge2[j].next)
if (edge2[j].to!=son[u]&&edge2[j].to!=fa[u])
{
dfs2(edge2[j].to,edge2[j].to);
}
}
void pushdown(int rt)
{
if (lazy[rt])
{
c[rt*]=;
lazy[rt*]=;
c[rt*+]=;
lazy[rt*+]=;
lazy[rt]=;
}
}
int ask(int rt,int l,int r,int L,int R)
{
if (l!=r) pushdown(rt);
if (l>=L&&r<=R)
{
return c[rt];
}
pushdown(rt);
int mid=(l+r)/,s=;
if (L<=mid) s+=ask(rt*,l,mid,L,R);
if (R>mid) s+=ask(rt*+,mid+,r,L,R);
return s;
}
int query(int x,int y)
{
int s=;
while (top[x]!=top[y])
{
if (dep[top[x]]<dep[top[y]]) swap(x,y);
s+=ask(,,n,pos[top[x]],pos[x]);
x=fa[top[x]];
}
if (dep[x]>dep[y]) swap(x,y);
if (x!=y) s+=ask(,,n,pos[x]+,pos[y]);
return s;
}
void update(int rt,int l,int r,int L,int R)
{
if (l!=r)pushdown(rt);
if (l>=L&&r<=R)
{
c[rt]=;
lazy[rt]=;
return;
}
pushdown(rt);
int mid=(l+r)/;
if (L<=mid) update(rt*,l,mid,L,R);
if (R>mid) update(rt*+,mid+,r,L,R);
c[rt]=c[rt*]+c[rt*+];
}
void change(int x,int y)
{
while (top[x]!=top[y])
{
if (dep[top[x]]<dep[top[y]]) swap(x,y);
update(,,n,pos[top[x]],pos[x]);
x=fa[top[x]];
}
if (dep[x]>dep[y]) swap(x,y);
if (x!=y) update(,,n,pos[x]+,pos[y]);
}
void build(int rt,int l,int r)
{
if (l==r)
{
c[rt]=;
return;
}
int mid=(l+r)/;
build(rt*,l,mid);
build(rt*+,mid+,r);
c[rt]=c[rt*]+c[rt*+];
}
int main()
{int i,j,x,y;
//freopen("lane.in","r",stdin);
//freopen("lane.out","w",stdout);
cin>>n>>m;
for (i=; i<=m; i++)
{
scanf("%d%d",&a[i].p,&a[i].q);
vis.insert(make_pair(a[i].p,a[i].q));
}
int tot1=;
scanf("%d",&cx);
while (cx!=-)
{
scanf("%d%d",&x,&y);
if (cx==)
{
tot1++;
b[tot1].p=x;
b[tot1].q=y;
b[tot1].c=;
vis.erase(make_pair(x,y));
}
else
{
tot1++;
b[tot1].c=;
b[tot1].p=x;
b[tot1].q=y;
}
scanf("%d",&cx);
}
for (i=; i<=m; i++)
{
if (vis.find(make_pair(a[i].p,a[i].q))!=vis.end())
{
add1(a[i].p,a[i].q);
add1(a[i].q,a[i].p);
}
}
make_tree();
dfs1(,,);
dfs2(,);
n=scccnt;
//cout<<scccnt<<endl;
build(,,n);
memset(ans,-,sizeof(ans));
for (i=tot1; i>=; i--)
{
if (b[i].c==)
{
int u=sccno[b[i].p],v=sccno[b[i].q];
change(u,v);
}
else
{
int u=sccno[b[i].p],v=sccno[b[i].q];
ans[i]=query(u,v);
}
}
for (i=; i<=tot1; i++)
if (ans[i]!=-) printf("%d\n",ans[i]);
}
[Ahoi2005]LANE 航线规划的更多相关文章
- BZOJ 1969: [Ahoi2005]LANE 航线规划( 树链剖分 )
首先我们要时光倒流, 倒着做, 变成加边操作维护关键边. 先随意搞出一颗树, 树上每条边都是关键边(因为是树, 去掉就不连通了)....然后加边(u, v)时, 路径(u, v)上的所有边都变成非关键 ...
- 【BZOJ 1969】 1969: [Ahoi2005]LANE 航线规划 (树链剖分+线段树)
1969: [Ahoi2005]LANE 航线规划 Description 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系——一个巨大的由千百万星 ...
- 【BZOJ1969】[Ahoi2005]LANE 航线规划 离线+树链剖分+线段树
[BZOJ1969][Ahoi2005]LANE 航线规划 Description 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系——一个巨大的由 ...
- BZOJ1969: [Ahoi2005]LANE 航线规划(LCT)
Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 587 Solved: 259[Submit][Status][Discuss] Description ...
- 【刷题】BZOJ 1969 [Ahoi2005]LANE 航线规划
Description 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系--一个巨大的由千百万星球构成的Samuel星系. 星际空间站的Samuel ...
- 【bzoj1959】[Ahoi2005]LANE 航线规划 树链剖分+线段树
题目描述 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系——一个巨大的由千百万星球构成的Samuel星系. 星际空间站的Samuel II巨型计算 ...
- [bzoj1969] [Ahoi2005]LANE 航线规划
tarjan.并查集.树状数组.树链剖分. 时间倒流,变删边为加边. 先求一波边双联通分量,缩点. 题目保证最后还是整张图联通的..所以就是一棵树. 现在的操作就是,将路径上的边权置0(加边时),查询 ...
- BZOJ 1969: [Ahoi2005]LANE 航线规划 [树链剖分 时间倒流]
题意: 一张图,删除边,求两点之间的割边数量.保证任意时刻图连通 任求一棵生成树,只有树边可能是割边 时间倒流,加入一条边,就是两点路径上的边都不可能是割边,区间覆盖... 然后本题需要把边哈希一下, ...
- ●BZOJ 1969 [Ahoi2005]LANE 航线规划
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=1969 题解: 线段树,树链剖分,反向考虑思路是很巧妙,但是感觉代码真的恶心.. 反着考虑,先 ...
随机推荐
- 201621123040《Java程序设计》第八周学习总结
1.本周学习总结 2.书面作业 2.1ArrayList代码分析 2.1.1解释ArrayList的contains源代码 ArrayList.contain()方法通过调用indexOf()来判断元 ...
- C程序设计-----第1次作业
一. PTA作业. 在完成PTA作业的时候我没有认真读题.每次都是提交完整代码 6-1(1) #include <stdio.h> //P++等价于(p)++还是等价于*(p++)? ...
- C++类型萃取
stl中的迭代器和C++中的类型萃取: http://www.itnose.net/detail/6487058.html 赐教!
- JAVA中GridBagLayout布局管理器应用详解
很多情况下,我们已经不需要通过编写代码来实现一个应用程序的图形界面,而是通过强大的IDE工具通过拖拽辅以简单的事件处理代码即可很轻松的完成.但是我们不得不面对这样操作存在的一些问题,有时候我们希望能够 ...
- android 广播安装指定下载的apk
// 广播出去,由广播接收器来处理下载完成的文件 Intent sendIntent = new Intent("com.test.downloadComplete"); ...
- Python split()方法
Python split()方法 描述 Python split()通过指定分隔符对字符串进行切片,如果参数num 有指定值,则仅分隔 num 个子字符串 语法 split()方法语法: str.sp ...
- Hibernate实体类注解解释
Hibernate注解1.@Entity(name="EntityName")必须,name为可选,对应数据库中一的个表2.@Table(name="",cat ...
- thinkphp中定义自己的函数
可以在前台和后台的公共文件夹中common.php中定义自己的函数,这样就可以在控制器中调用,而不需要调用对象了 /** * @name addvtorandp * @author 黄峰1664253 ...
- C++ 实现一个信号量
C++ 实现一个信号量 信号量有很多应用场景,事实上只要是生产者-消费者模型,一般都需要一个信号量来控制. POSIX接口是有PV信号量API的.但C++标准没有.下面是一个PV信号量的简单实现.有些 ...
- LXC学习实践(3)快速体验第一个容器
1.搭建第一个 LXC 虚拟计算机 #yum install lxc* 2.安装软件包后要检查 Linux 发行版的内核对 LXC 的支持情况,可以使用下面命令 #lxc-checkconfig #l ...