●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巨型计算 ...
随机推荐
- django搭建web (二) urls.py
URL模式: 在app下的urls.py中 urlpatterns=[ url(正则表达式,view函数,参数,别名,前缀)] urlpatterns=[ url(r'^hello/$',hello. ...
- 服务器数据恢复_服务器xfs数据丢失数据恢复
简介:太原一家公司的服务器出现故障,服务器是linux服务器,连接了一台某型号的存储,文件系统为xfs文件系统.管理员使用xfs_repair工具试图对文件系统进行修复但修复失败,linux服务器中所 ...
- asp.net web api 控制器
1控制器操作的参数 控制器操作的参数可以是内置类型也可以是自定义类型,无参也是允许的. 2控制器操作返回值 类型 说明 void 操作返回值为void时,Web API返回空HTTP响应,其状态码为2 ...
- python之路--day11---迭代器和生成器
迭代: 迭代是一个重复的过程,每次重复即一次迭代,并且每次迭代的结果都是下一次迭代的初始值 为什么要有迭代器: 数据类型的取值,字符串,列表,元组依靠索引可以取值,但是字典,集合,文件这些数据类型无法 ...
- $.ajax 中的contentType
$.ajax 中的contentType 在 cnodejs.org 论坛中有一个问题,让我也很奇怪,说是 $.ajax 设置数据类型 applicaiton/json之后,服务器端(express) ...
- MySQL Group Relication 部署环境入门篇
一:环境介绍 cenos 6.7 版本 数据库的版本5.7.19 二:部署规划单机多实例的部署 端口号 数据目录 group_repplicatoon 通信接口 3307 /data ...
- ajax中设置contentType: “application/json”的作用
最近在做项目交互的时候,刚开始向后台传递数据返回415,后来百度添加了 contentType:"application/json"之后返回400,然后把传输的数据格式改为json ...
- ELK学习总结（3-3）elk的组合查询
1.bool 查询: must: 必须 should: 可以满足,也可以不满足. must_not: minimum_should_match: 至少要x个匹配才算匹配成功 disable_coor ...
- Linux--慕课学习
刚开始接触Linux,很有幸的在慕课网上看到了Peter老师的Linux入门课程,老师讲课真的式行云流水,深入浅出,循循善诱,层层递进. 老师分享的都是自己多年来总结的经验.看完之后也学到了很多东西. ...
- Lua编写wireshark插件初探——解析Websocket上的MQTT协议
一.背景 最近在做物联网流量分析时发现, App在使用MQTT协议时往往通过SSL+WebSocket+MQTT这种方式与服务器通信,在使用SSL中间人截获数据后,Wireshark不能自动解析出MQ ...