●BZOJ 1969 [Ahoi2005]LANE 航线规划
题链：
http://www.lydsy.com/JudgeOnline/problem.php?id=1969
题解：
线段树，树链剖分，反向考虑
思路是很巧妙，但是感觉代码真的恶心。。
反着考虑，先按照给出的操作把所有该删的边都删掉，
那么剩下的也还是一个联通块。（题目保证了的）
然后在该联通块里随便选出某些边来形成一颗树。
那么显然，若只看树边的话，每一条树边都是一条关键边，
如果此时加一条其他的边 e(x,y)，使得形成了一个环，
那么显然 x ~ y 路径上的边都不再是关键边。
那么此时问题变为了树上的区间覆盖问题，树链剖分+线段树来解决就好了。
（另外，会用到 hash或者 map的）
代码：（三份，分别是用到map，hash和hash表的，hahaha）
map
#include<map>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 40500
#define MAXM 100500
using namespace std;
int fa[MAXN],bgs[MAXN],top[MAXN],dep[MAXN],pos[MAXN],ord[MAXN],ANS[MAXN];
int N,M,snt,ont,ant;
typedef pair<int,int>pii;
map<pii,bool>H;
struct Oper{
int type,x,y;
}stk[MAXN];
struct Edge{
int to[MAXM*2],nxt[MAXM*2],mark[MAXM*2],head[MAXN],ent;
Edge(){
ent=2; memset(head,0,sizeof(head));
}
void Adde(int u,int v){
to[ent]=v; nxt[ent]=head[u]; head[u]=ent++;
to[ent]=u; nxt[ent]=head[v]; head[v]=ent++;
}
}E;
struct SGT{
int ls[MAXN*2],rs[MAXN*2],rem[MAXN*2],sz,rt;
void pushup(int u){
rem[u]=rem[ls[u]]+rem[rs[u]];
}
void Build(int &u,int l,int r){
u=++sz;
if(l==r){rem[u]=1; return;}
int mid=(l+r)>>1;
Build(ls[u],l,mid);
Build(rs[u],mid+1,r);
pushup(u);
}
void Modify(int u,int l,int r,int al,int ar){
if(!rem[u]) return;
if(al<=l&&r<=ar){rem[u]=0; return;}
int mid=(l+r)>>1;
if(al<=mid) Modify(ls[u],l,mid,al,ar);
if(mid<ar) Modify(rs[u],mid+1,r,al,ar);
pushup(u);
}
int Query(int u,int l,int r,int al,int ar){
if(!rem[u]) return 0;
if(al<=l&&r<=ar) return rem[u];
int mid=(l+r)>>1,ret=0;
if(al<=mid) ret+=Query(ls[u],l,mid,al,ar);
if(mid<ar) ret+=Query(rs[u],mid+1,r,al,ar);
return ret;
}
}DT;
void read(int &x){
static int f; static char ch;
x=0; f=1; ch=getchar();
while(ch<'0'||'9'<ch){if(ch=='-')f=-1;ch=getchar();}
while('0'<=ch&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
x=x*f;
}
int dfs1(int u,int dad){
static bool vis[MAXN];
int bgn=0,num=1,son; vis[u]=1; fa[u]=dad;
for(int i=E.head[u],v;i;i=E.nxt[i]){
v=E.to[i]; if(vis[v]) continue;
if(H[make_pair(u,v)]) continue;
E.mark[i]=1; E.mark[i^1]=1;
son=dfs1(v,u);
if(bgn<son) bgn=son,bgs[u]=v;
num+=son;
}
return num;
}
void dfs2(int u,int tp){
ord[++ont]=u; pos[u]=ont;
top[u]=tp; dep[u]=dep[fa[u]]+1;
if(bgs[u]) dfs2(bgs[u],tp);
for(int i=E.head[u],v;i;i=E.nxt[i]) if(E.mark[i]){
v=E.to[i];
if(v==bgs[u]||v==fa[u]) continue;
dfs2(v,v);
}
}
void cover(int u,int v){
while(top[u]!=top[v]){
if(dep[top[u]]>dep[top[v]]) swap(u,v);
DT.Modify(DT.rt,1,N,pos[top[v]],pos[v]);
v=fa[top[v]];
}
if(u==v) return;
if(dep[u]>dep[v]) swap(u,v);
DT.Modify(DT.rt,1,N,pos[u]+1,pos[v]);
}
int answer(int u,int v){
int ret=0;
while(top[u]!=top[v]){
if(dep[top[u]]>dep[top[v]]) swap(u,v);
ret+=DT.Query(DT.rt,1,N,pos[top[v]],pos[v]);
v=fa[top[v]];
}
if(u==v) return ret;
if(dep[u]>dep[v]) swap(u,v);
ret+=DT.Query(DT.rt,1,N,pos[u]+1,pos[v]);
return ret;
}
int main(){
//freopen("1969.out","w",stdout);
read(N); read(M);
for(int i=1,u,v;i<=M;i++)
read(u),read(v),E.Adde(u,v);
while(1){
read(stk[snt+1].type); if(stk[snt+1].type==-1) break;
snt++; read(stk[snt].x); read(stk[snt].y);
if(stk[snt].type) continue;
H[make_pair(stk[snt].x,stk[snt].y)]=1;
H[make_pair(stk[snt].y,stk[snt].x)]=1;
}
dfs1(1,0);
dfs2(1,1);
DT.Build(DT.rt,1,N);
for(int u=1,v;u<=N;u++)
for(int i=E.head[u];i;i=E.nxt[i]){
v=E.to[i]; if(u>v||E.mark[i]) continue;
if(H[make_pair(u,v)]) continue;
cover(u,v);
}
while(snt){
if(stk[snt].type) ANS[++ant]=answer(stk[snt].x,stk[snt].y);
else cover(stk[snt].x,stk[snt].y);
snt--;
}
while(ant) printf("%d\n",ANS[ant--]);
return 0;
}
hash
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 40500
#define MAXM 100500
using namespace std;
int fa[MAXN],bgs[MAXN],top[MAXN],dep[MAXN],pos[MAXN],ord[MAXN],ANS[MAXN];
int N,M,snt,ont,ant;
struct Oper{
int type,x,y;
}stk[MAXN];
struct Edge{
int to[MAXM*2],nxt[MAXM*2],mark[MAXM*2],head[MAXN],ent;
Edge(){
ent=2; memset(head,0,sizeof(head));
}
void Adde(int u,int v){
to[ent]=v; nxt[ent]=head[u]; head[u]=ent++;
to[ent]=u; nxt[ent]=head[v]; head[v]=ent++;
}
}E;
struct Hash{
long long C[MAXN];int cnt;
void Push(int a,int b){
if(a>b) swap(a,b);
C[++cnt]=(1ll*a<<30)+b;
}
void Finish(){
sort(C+1,C+cnt+1);
}
bool Exist(int a,int b){
static int p; static long long w;
if(a>b) swap(a,b);
w=(1ll*a<<30)+b;
p=lower_bound(C+1,C+cnt+1,w)-C;
return C[p]==w;
}
}H;
struct SGT{
int ls[MAXN*2],rs[MAXN*2],rem[MAXN*2],sz,rt;
void pushup(int u){
rem[u]=rem[ls[u]]+rem[rs[u]];
}
void Build(int &u,int l,int r){
u=++sz;
if(l==r){rem[u]=1; return;}
int mid=(l+r)>>1;
Build(ls[u],l,mid);
Build(rs[u],mid+1,r);
pushup(u);
}
void Modify(int u,int l,int r,int al,int ar){
if(!rem[u]) return;
if(al<=l&&r<=ar){rem[u]=0; return;}
int mid=(l+r)>>1;
if(al<=mid) Modify(ls[u],l,mid,al,ar);
if(mid<ar) Modify(rs[u],mid+1,r,al,ar);
pushup(u);
}
int Query(int u,int l,int r,int al,int ar){
if(!rem[u]) return 0;
if(al<=l&&r<=ar) return rem[u];
int mid=(l+r)>>1,ret=0;
if(al<=mid) ret+=Query(ls[u],l,mid,al,ar);
if(mid<ar) ret+=Query(rs[u],mid+1,r,al,ar);
return ret;
}
}DT;
void read(int &x){
static int f; static char ch;
x=0; f=1; ch=getchar();
while(ch<'0'||'9'<ch){if(ch=='-')f=-1;ch=getchar();}
while('0'<=ch&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
x=x*f;
}
int dfs1(int u,int dad){
static bool vis[MAXN];
int bgn=0,num=1,son; vis[u]=1; fa[u]=dad;
for(int i=E.head[u],v;i;i=E.nxt[i]){
v=E.to[i]; if(vis[v]) continue;
if(H.Exist(u,v)) continue;
E.mark[i]=1; E.mark[i^1]=1;
son=dfs1(v,u);
if(bgn<son) bgn=son,bgs[u]=v;
num+=son;
}
return num;
}
void dfs2(int u,int tp){
ord[++ont]=u; pos[u]=ont;
top[u]=tp; dep[u]=dep[fa[u]]+1;
if(bgs[u]) dfs2(bgs[u],tp);
for(int i=E.head[u],v;i;i=E.nxt[i]) if(E.mark[i]){
v=E.to[i];
if(v==bgs[u]||v==fa[u]) continue;
dfs2(v,v);
}
}
void cover(int u,int v){
while(top[u]!=top[v]){
if(dep[top[u]]>dep[top[v]]) swap(u,v);
DT.Modify(DT.rt,1,N,pos[top[v]],pos[v]);
v=fa[top[v]];
}
if(u==v) return;
if(dep[u]>dep[v]) swap(u,v);
DT.Modify(DT.rt,1,N,pos[u]+1,pos[v]);
}
int answer(int u,int v){
int ret=0;
while(top[u]!=top[v]){
if(dep[top[u]]>dep[top[v]]) swap(u,v);
ret+=DT.Query(DT.rt,1,N,pos[top[v]],pos[v]);
v=fa[top[v]];
}
if(u==v) return ret;
if(dep[u]>dep[v]) swap(u,v);
ret+=DT.Query(DT.rt,1,N,pos[u]+1,pos[v]);
return ret;
}
int main(){
//freopen("1969.out","w",stdout);
read(N); read(M);
for(int i=1,u,v;i<=M;i++)
read(u),read(v),E.Adde(u,v);
while(1){
read(stk[snt+1].type); if(stk[snt+1].type==-1) break;
snt++; read(stk[snt].x); read(stk[snt].y);
if(stk[snt].type) continue;
H.Push(stk[snt].x,stk[snt].y);
}
H.Finish();
dfs1(1,0);
dfs2(1,1);
DT.Build(DT.rt,1,N);
for(int u=1,v;u<=N;u++)
for(int i=E.head[u];i;i=E.nxt[i]){
v=E.to[i]; if(u>v||E.mark[i]) continue;
if(H.Exist(u,v)) continue;
cover(u,v);
}
while(snt){
if(stk[snt].type) ANS[++ant]=answer(stk[snt].x,stk[snt].y);
else cover(stk[snt].x,stk[snt].y);
snt--;
}
while(ant) printf("%d\n",ANS[ant--]);
return 0;
}
hash表（第一次写哈希表）
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 40500
#define MAXM 100500
using namespace std;
int fa[MAXN],bgs[MAXN],top[MAXN],dep[MAXN],pos[MAXN],ord[MAXN],ANS[MAXN];
int N,M,snt,ont,ant;
struct Oper{
int type,x,y;
}stk[MAXN];
struct Edge{
int to[MAXM*2],nxt[MAXM*2],mark[MAXM*2],head[MAXN],ent;
Edge(){ent=2;}
void Adde(int u,int v){
to[ent]=v; nxt[ent]=head[u]; head[u]=ent++;
to[ent]=u; nxt[ent]=head[v]; head[v]=ent++;
}
}E;
struct Hash{
int to1[MAXN],to2[MAXN],nxt[MAXN],head[50000],ent,BIT,MOD;
Hash(){BIT=12356789; MOD=43517; ent=2;}
void Push(int a,int b){
static int u;
if(a>b) swap(a,b);
u=(1ll*a*BIT+b)%MOD;
to1[ent]=a; to2[ent]=b;
nxt[ent]=head[u]; head[u]=ent++;
}
bool Exist(int a,int b){
static int u;
if(a>b) swap(a,b);
u=(1ll*a*BIT+b)%MOD;
for(int i=head[u];i;i=nxt[i])
if(to1[i]==a&&to2[i]==b) return 1;
return 0;
}
}H;
struct SGT{
int ls[MAXN*2],rs[MAXN*2],rem[MAXN*2],sz,rt;
void pushup(int u){
rem[u]=rem[ls[u]]+rem[rs[u]];
}
void Build(int &u,int l,int r){
u=++sz;
if(l==r){rem[u]=1; return;}
int mid=(l+r)>>1;
Build(ls[u],l,mid);
Build(rs[u],mid+1,r);
pushup(u);
}
void Modify(int u,int l,int r,int al,int ar){
if(!rem[u]) return;
if(al<=l&&r<=ar){rem[u]=0; return;}
int mid=(l+r)>>1;
if(al<=mid) Modify(ls[u],l,mid,al,ar);
if(mid<ar) Modify(rs[u],mid+1,r,al,ar);
pushup(u);
}
int Query(int u,int l,int r,int al,int ar){
if(!rem[u]) return 0;
if(al<=l&&r<=ar) return rem[u];
int mid=(l+r)>>1,ret=0;
if(al<=mid) ret+=Query(ls[u],l,mid,al,ar);
if(mid<ar) ret+=Query(rs[u],mid+1,r,al,ar);
return ret;
}
}DT;
void read(int &x){
static int f; static char ch;
x=0; f=1; ch=getchar();
while(ch<'0'||'9'<ch){if(ch=='-')f=-1;ch=getchar();}
while('0'<=ch&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
x=x*f;
}
int dfs1(int u,int dad){
static bool vis[MAXN];
int bgn=0,num=1,son; vis[u]=1; fa[u]=dad;
for(int i=E.head[u],v;i;i=E.nxt[i]){
v=E.to[i]; if(vis[v]) continue;
if(H.Exist(u,v)) continue;
E.mark[i]=1; E.mark[i^1]=1;
son=dfs1(v,u);
if(bgn<son) bgn=son,bgs[u]=v;
num+=son;
}
return num;
}
void dfs2(int u,int tp){
ord[++ont]=u; pos[u]=ont;
top[u]=tp; dep[u]=dep[fa[u]]+1;
if(bgs[u]) dfs2(bgs[u],tp);
for(int i=E.head[u],v;i;i=E.nxt[i]) if(E.mark[i]){
v=E.to[i];
if(v==bgs[u]||v==fa[u]) continue;
dfs2(v,v);
}
}
void cover(int u,int v){
while(top[u]!=top[v]){
if(dep[top[u]]>dep[top[v]]) swap(u,v);
DT.Modify(DT.rt,1,N,pos[top[v]],pos[v]);
v=fa[top[v]];
}
if(u==v) return;
if(dep[u]>dep[v]) swap(u,v);
DT.Modify(DT.rt,1,N,pos[u]+1,pos[v]);
}
int answer(int u,int v){
int ret=0;
while(top[u]!=top[v]){
if(dep[top[u]]>dep[top[v]]) swap(u,v);
ret+=DT.Query(DT.rt,1,N,pos[top[v]],pos[v]);
v=fa[top[v]];
}
if(u==v) return ret;
if(dep[u]>dep[v]) swap(u,v);
ret+=DT.Query(DT.rt,1,N,pos[u]+1,pos[v]);
return ret;
}
int main(){
//freopen("1969.out","w",stdout);
read(N); read(M);
for(int i=1,u,v;i<=M;i++)
read(u),read(v),E.Adde(u,v);
while(1){
read(stk[snt+1].type); if(stk[snt+1].type==-1) break;
snt++; read(stk[snt].x); read(stk[snt].y);
if(stk[snt].type) continue;
H.Push(stk[snt].x,stk[snt].y);
}
dfs1(1,0);
dfs2(1,1);
DT.Build(DT.rt,1,N);
for(int u=1,v;u<=N;u++)
for(int i=E.head[u];i;i=E.nxt[i]){
v=E.to[i]; if(u>v||E.mark[i]) continue;
if(H.Exist(u,v)) continue;
cover(u,v);
}
while(snt){
if(stk[snt].type) ANS[++ant]=answer(stk[snt].x,stk[snt].y);
else cover(stk[snt].x,stk[snt].y);
snt--;
}
while(ant) printf("%d\n",ANS[ant--]);
return 0;
}
●BZOJ 1969 [Ahoi2005]LANE 航线规划的更多相关文章
- BZOJ 1969: [Ahoi2005]LANE 航线规划( 树链剖分 )
首先我们要时光倒流, 倒着做, 变成加边操作维护关键边. 先随意搞出一颗树, 树上每条边都是关键边(因为是树, 去掉就不连通了)....然后加边(u, v)时, 路径(u, v)上的所有边都变成非关键 ...
- 【刷题】BZOJ 1969 [Ahoi2005]LANE 航线规划
Description 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系--一个巨大的由千百万星球构成的Samuel星系. 星际空间站的Samuel ...
- BZOJ 1969: [Ahoi2005]LANE 航线规划 [树链剖分 时间倒流]
题意: 一张图,删除边,求两点之间的割边数量.保证任意时刻图连通 任求一棵生成树,只有树边可能是割边 时间倒流,加入一条边,就是两点路径上的边都不可能是割边,区间覆盖... 然后本题需要把边哈希一下, ...
- 【BZOJ 1969】 1969: [Ahoi2005]LANE 航线规划 （树链剖分+线段树）
1969: [Ahoi2005]LANE 航线规划 Description 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系——一个巨大的由千百万星 ...
- 【BZOJ】1969: [Ahoi2005]LANE 航线规划
题目链接: 传送~~ 题解: 老夫实在是码不动了…… 正着搞显然不好做,尝试倒着乱搞.先给被删除的边标记一个时间戳,先删除的时间戳大,同时维护询问时间戳,询问早的时间戳大.没被删除过的边时间戳都是0 ...
- 【BZOJ1969】[Ahoi2005]LANE 航线规划 离线+树链剖分+线段树
[BZOJ1969][Ahoi2005]LANE 航线规划 Description 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系——一个巨大的由 ...
- [Ahoi2005]LANE 航线规划
题目描述 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系——一个巨大的由千百万星球构成的Samuel星系. 星际空间站的Samuel II巨型计算 ...
- BZOJ1969: [Ahoi2005]LANE 航线规划(LCT)
Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 587 Solved: 259[Submit][Status][Discuss] Description ...
- 【bzoj1959】[Ahoi2005]LANE 航线规划 树链剖分+线段树
题目描述 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系——一个巨大的由千百万星球构成的Samuel星系. 星际空间站的Samuel II巨型计算 ...
随机推荐
- 团队项目7——团队冲刺（beta版本）
beta版本冲刺计划安排:http://www.cnblogs.com/ricardoCYF/p/8018413.html 12.06:http://www.cnblogs.com/ricardoCY ...
- node.js基础
//安装淘宝npm镜像 npm install -g cnpm --registry=https://registry.npm.taobao.org//require表示引包,引包就是引用自己的一个特 ...
- OO前三次作业总结
一.第一次作业 1.程序设计分析  图1 第一次作业类图  类图 程序的大致结构如下: 代码分析 ...
- LeetCode & Q219-Contains Duplicate II
Array Hash Table Description: Given an array of integers and an integer k, find out whether there ar ...
- api-gateway实践（15）3.6JL分支和3.7并行改造需求
一.名称改为"API网关" --哪个地方的名称?二.开发者视图中,API网关显示两个视图. 1. 服务分类视图:支持按照业务分为多个类别,分类方式参照应用服务化的分类:人像比对.自 ...
- Linux知识积累 （9） 创建用户、分配权限和更改所有者
一.useradd和adduser 1.useradd命令: 用于Linux中创建的新的系统用户. useradd可用来建立用户帐号.帐号建好之后,再用passwd设定帐号的密码． 而可用userde ...
- SpringCloud的EurekaClient ： 客户端应用访问注册的微服务（无断路器场景）
演示客户端应用如何访问注册在EurekaServer里的微服务 一.概念和定义 采用Ribbon或Feign方式访问注册到EurekaServer中的微服务.1.Ribbon实现了客户端负载均衡,2. ...
- ASP.NET CORE系列【三】使用Entity Framework Core进行增删改查
身份验证 以前我们熟悉的web.config中配置的form验证,现在没有了.我们来看看在Core里面如何配置: 首先需要NuGet安装一个包:Microsoft.AspNetCore.Authent ...
- Python的下载及安装
1.官网下载地址:https://www.python.org/downloads/ 2.python设置环境变量: 在系统变量里添加Python的安装位置 3.在cmd里输入python里即可