2019牛客暑期多校训练营(第八场)E:Explorer(LCT裸题 也可用线段树模拟并查集维护连通性)

题意:给定N,M,然后给出M组信息(u,v,l,r),表示u到v有[l,r]范围的通行证有效。问有多少种通行证可以使得1和N连通。

思路:和bzoj魔法森林有点像,LCT维护最小生成树。  开始和队友在想维护连通性,而不是维护树,这样好像会很麻烦。

队友yy了一个算法:用线段树模拟并查集维护连通性。(发现和标程有点像?

我的代码:LCT维护最小生成树。

...先给代码,后面补一下题解。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
vector<int>e[];
int n,m,t[],q[],siz[],ans,efree,a[],cnt,l[],r[],u[],v[],f[];
template<class T>
inline void read(T&a){
char c=getchar();
for(a=;(c<''||c>'')&&c!='-';c=getchar());
bool f=;if(c=='-')f=,c=getchar();
for(;c>=''&&c<='';c=getchar())a=a*+c-'';
if(f)a=-a;
}
inline int gf(int x){return x==f[x]?f[x]:gf(f[x]);}
void ins(int k,int nl,int nr,int l,int r,int x){
if(nl==l&&nr==r)return (void)e[k].push_back(x);
int mid=(nl+nr)>>;
if(r<=mid)ins(k<<,nl,mid,l,r,x);
else if(l>mid)ins(k<<|,mid+,nr,l,r,x);
else ins(k<<,nl,mid,l,mid,x),ins(k<<|,mid+,nr,mid+,r,x);
}
void dfs(int k,int l,int r){
vector<int>lastf;
int m=e[k].size(),mid=(l+r)>>;
for(register int i=;i<m;i++){
int x=u[e[k][i]],y=v[e[k][i]];
x=gf(x),y=gf(y);
if(siz[x]>siz[y])swap(x,y);
lastf.push_back(x);f[x]=y;siz[y]+=siz[x];
}
if(gf()==gf(n))ans+=a[r+]-a[l];
else if(l<r)dfs(k<<,l,mid),dfs(k<<|,mid+,r);
m=lastf.size();
for(register int i=m-;i>=;i--)f[lastf[i]]=lastf[i];
lastf.clear();
}
int main(){
read(n),read(m);
for(int i=;i<=n;i++)f[i]=i,siz[i]=;
for(int i=;i<=m;i++){
read(u[i]),read(v[i]),read(l[i]),read(r[i]);
a[++cnt]=l[i];a[++cnt]=r[i]+;
}
sort(a+,a++cnt);cnt=unique(a+,a++cnt)-a-;
for(int i=;i<=m;i++)
ins(,,cnt-,lower_bound(a+,a++cnt,l[i])-a,lower_bound(a+,a++cnt,r[i]+)-a-,i);
dfs(,,cnt-);printf("%d\n",ans);
return ;
}

LCT代码:

#include<bits/stdc++.h>
using namespace std;
const int maxn=;
const int inf=;
struct edge{
int x,y,a,b;
}e[maxn];
bool cmp (edge w,edge v){
return w.b>v.b;
}
void read(int &x){
char c=getchar(); x=;
for(;c>''||c<'';c=getchar());
for(;c<=''&&c>='';c=getchar()) x=(x<<)+(x<<)+c-'';
}
struct LCT
{
int Max[maxn],rev[maxn],ch[maxn][],fa[maxn],stc[maxn],top;
int isroot(int x){
return ch[fa[x]][]!=x&&ch[fa[x]][]!=x;
}
int get(int x){
return ch[fa[x]][]==x;
}
void pushdown(int x)
{
if(!rev[x]||!x) return ;
swap(ch[x][],ch[x][]);
if(ch[x][]) rev[ch[x][]]^=;
if(ch[x][]) rev[ch[x][]]^=;
rev[x]=;
}
void pushup(int x)
{
Max[x]=x;
if(ch[x][]&&e[Max[ch[x][]]].a>e[Max[x]].a) Max[x]=Max[ch[x][]];
if(ch[x][]&&e[Max[ch[x][]]].a>e[Max[x]].a) Max[x]=Max[ch[x][]];
}
void rotate(int x)
{
int old=fa[x],fold=fa[old],opt=get(x);
if(!isroot(old)) ch[fold][get(old)]=x;
fa[x]=fold;
ch[old][opt]=ch[x][opt^]; fa[ch[old][opt]]=old;
ch[x][opt^]=old; fa[old]=x;
pushup(old); pushup(x);
}
void splay(int x)
{
int top=; stc[++top]=x;
for(int i=x;!isroot(i);i=fa[i]) stc[++top]=fa[i];
for(int i=top;i;i--) pushdown(stc[i]);
for(int f;!isroot(x);rotate(x)){
if(!isroot(f=fa[x]))
rotate(get(x)==get(f)?f:x);
}
}
void access(int x)
{
int rson=;
for(;x;rson=x,x=fa[x]){
splay(x);
ch[x][]=rson;
pushup(x);
}
}
int find(int x){ access(x); splay(x); while(ch[x][]) x=ch[x][]; return x;}
int query(int x,int y) { make_root(y); access(x); splay(x); return Max[x]; }
void make_root(int x) { access(x); splay(x); rev[x]^=; }
void link(int x,int y) { make_root(x); fa[x]=y; splay(x); }
void cut(int x,int y) { make_root(x); access(y); splay(y); fa[x]=ch[y][]=; }
}S;
pair<int,int>fcy[maxn]; int tot;
int main()
{
int N,M,i,ans=;
read(N); read(M);
for(i=;i<=M;i++){
read(e[i].x); read(e[i].y);
read(e[i].a); read(e[i].b);
}
sort(e+,e+M+,cmp);
for(i=;i<=M;i++){
if(S.find(M+e[i].x)!=S.find(M+e[i].y)) {
S.link(i,M+e[i].x);
S.link(i,M+e[i].y);
}
else{
int tmp=S.query(M+e[i].x,M+e[i].y);
if(e[tmp].a>e[i].a) {
S.cut(tmp,M+e[tmp].x); S.cut(tmp,M+e[tmp].y);
S.link(i,M+e[i].x); S.link(i,M+e[i].y);
}
}
if(S.find(M+)==S.find(M+N)){
int tmp=S.query(M+,M+N);
if(e[tmp].a<=e[i].b){
tot++;
fcy[tot].first=e[tmp].a; fcy[tot].second=e[i].b;
}
}
}
sort(fcy+,fcy+tot+);
for(i=;i<=tot;i++){
int j=i,Mx=fcy[i].second;
while(j+<=tot&&fcy[j+].first<=Mx) {
j++; Mx=max(Mx,fcy[j].second);
}
ans+=Mx-fcy[i].first+;
i=j;
}
printf("%d\n",ans);
return ;
}
上一篇:关于ubuntu16.4 中安装最新的eclipse或者是STS出现页面特卡,且新建项目没有提示,preference选项中点击左侧标签右侧没反应的解决办法,参照google, 排版不太好,希望对一些小伙伴有所帮助


下一篇:java操作Hbase