zoukankan      html  css  js  c++  java
  • P4747 [CERC2017]Intrinsic Interval

    题目

    分析

    析合树模板题。

    先建出析合树。

    就是要统计这样的点对的个数,那么在析合树上其实就是对于每一个点这样统计:

    (1.) 对于合点:任意两个儿子之间可以任意选择,那么方案加上 (num imes (num-1)/2)

    (2.) 对于析点:任何两个儿子都不能选,于是只有它本身这个节点是一个合法的,答案加一。

    然后就完了。

    代码

    #include<bits/stdc++.h>
    using namespace std;
    template <typename T>
    inline void read(T &x){
    	x=0;char ch=getchar();bool f=false;
    	while(!isdigit(ch)){if(ch=='-'){f=true;}ch=getchar();}
    	while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
    	x=f?-x:x;
    	return ;
    }
    template <typename T>
    inline void write(T x){
    	if(x<0) putchar('-'),x=-x;
    	if(x>9) write(x/10);
    	putchar(x%10^48);
    	return ;
    }
    #define PII pair<int,int>
    #define mp make_pair
    const int N=2e5+5;
    int n,m,a[N];
    namespace RMQ1{
    	int Max[21][N],lg[N];
    	inline void Init(int n){
    		for(int i=2;i<=n;i++) lg[i]=lg[i>>1]+1;
    		for(int i=1;i<=n;i++) Max[0][i]=a[i];
    		for(int i=1;i<=20;i++) for(int j=1;j+(1<<i)-1<=n;j++) Max[i][j]=max(Max[i-1][j],Max[i-1][j+(1<<(i-1))]);
    		return ;
    	}
    	inline int Query(int l,int r){
    		int tmp=lg[r-l+1];
    		return max(Max[tmp][l],Max[tmp][r-(1<<tmp)+1]);
    	}
    }
    namespace RMQ2{
    	int Min[21][N],lg[N];
    	inline void Init(int n){
    		for(int i=2;i<=n;i++) lg[i]=lg[i>>1]+1;
    		for(int i=1;i<=n;i++) Min[0][i]=a[i];
    		for(int i=1;i<=20;i++) for(int j=1;j+(1<<i)-1<=n;j++) Min[i][j]=min(Min[i-1][j],Min[i-1][j+(1<<(i-1))]);
    		return ;
    	}
    	inline int Query(int l,int r){
    		int tmp=lg[r-l+1];
    		return min(Min[tmp][l],Min[tmp][r-(1<<tmp)+1]);
    	}
    }
    int Min[N<<2],ad[N<<2];
    inline void Pushup(int x){Min[x]=min(Min[x<<1],Min[x<<1|1]);return ;}
    inline void PushDown(int x){
    	if(!ad[x]) return ;
    	ad[x<<1]+=ad[x],ad[x<<1|1]+=ad[x];
    	Min[x<<1]+=ad[x],Min[x<<1|1]+=ad[x];ad[x]=0;
    	return ;
    }
    void Modify(int x,int l,int r,int ql,int qr,int v){
    	if(ql<=l&&r<=qr) return Min[x]+=v,ad[x]+=v,void();
    	PushDown(x);int mid=l+r>>1;
    	if(ql<=mid) Modify(x<<1,l,mid,ql,qr,v);
    	if(qr>mid) Modify(x<<1|1,mid+1,r,ql,qr,v);
    	Pushup(x);
    	return;
    }
    int Query(int x,int l,int r){
    	if(l==r) return l;
    	PushDown(x);int mid=l+r>>1;
    	if(!Min[x<<1]) return Query(x<<1,l,mid);
    	return Query(x<<1|1,mid+1,r);
    }
    inline bool Check(int l,int r){return RMQ1::Query(l,r)-RMQ2::Query(l,r)==r-l;}
    int head[N],nex[N<<1],to[N<<1],idx;
    inline void add(int u,int v){
    	nex[++idx]=head[u];
    	to[idx]=v;
    	head[u]=idx;
    	return ;
    }
    int fa[N][21],dep[N];
    void dfs(int x,int f){
    	fa[x][0]=f;dep[x]=dep[f]+1;
    	for(int i=1;i<=20;i++) fa[x][i]=fa[fa[x][i-1]][i-1];
    	for(int i=head[x];i;i=nex[i]){
    		int y=to[i];
    		if(y==f) continue;
    		dfs(y,x);
    	}
    	return ;
    }
    inline int QueryKth(int u,int k){
    	for(int i=0;i<=20&&k;i++) if(k&(1<<i)) u=fa[u][i],k-=(1<<i);
    	return u;
    }
    inline int QueryLca(int u,int v){
    	if(dep[u]<dep[v]) swap(u,v);
    	u=QueryKth(u,dep[u]-dep[v]);
    	if(u==v) return u;
    	for(int i=20;i>=0;i--) if(fa[u][i]!=fa[v][i]) u=fa[u][i],v=fa[v][i];
    	return fa[u][0];
    }
    int rt,top1,top2,sta1[N],sta2[N],id[N],L[N],R[N],top,sta[N],cnt,typ[N],lc[N];
    void Build(){
    	for(int i=1;i<=n;i++){
    		while(top1&&a[i]<=a[sta1[top1]]) Modify(1,1,n,sta1[top1-1]+1,sta1[top1],a[sta1[top1]]),top1--;
    		while(top2&&a[i]>=a[sta2[top2]]) Modify(1,1,n,sta2[top2-1]+1,sta2[top2],-a[sta2[top2]]),top2--;
    		Modify(1,1,n,sta1[top1]+1,i,-a[i]),sta1[++top1]=i;
    		Modify(1,1,n,sta2[top2]+1,i,a[i]),sta2[++top2]=i;
    		id[i]=++cnt;
    		L[cnt]=R[cnt]=i;
    		int tmp=Query(1,1,n),now=cnt;
    		while(top&&L[sta[top]]>=tmp){
    			if(typ[sta[top]]&&Check(lc[sta[top]],i)) R[sta[top]]=i,lc[sta[top]]=L[now],add(sta[top],now),now=sta[top--];
    			else if(Check(L[sta[top]],i)) typ[++cnt]=1,L[cnt]=L[sta[top]],R[cnt]=i,lc[cnt]=L[now],add(cnt,sta[top--]),add(cnt,now),now=cnt;
    			else{
    				add(++cnt,now);
    				do{
    					add(cnt,sta[top--]);
    				}while(top&&!Check(L[sta[top]],i));
    				L[cnt]=L[sta[top]],R[cnt]=i,add(cnt,sta[top--]),now=cnt;
    			}
    		}
    		sta[++top]=now;
    		Modify(1,1,n,1,i,-1);
    	}
    	rt=sta[1];
    	return ;
    }
    inline PII Query(int l,int r){
    	int x=id[l],y=id[r],Lca=QueryLca(x,y),posl=1,posr=1;
    	if(typ[Lca]&1) posl=L[QueryKth(x,dep[x]-dep[Lca]-1)],posr=R[QueryKth(y,dep[y]-dep[Lca]-1)];
    	else posl=L[Lca],posr=R[Lca];
    	return mp(posl,posr);
    }
    signed main(){
    	read(n);
    	for(int i=1;i<=n;i++) read(a[i]);
    	RMQ1::Init(n),RMQ2::Init(n);
    	Build();
    	dfs(rt,0);
    	read(m);
    	for(int i=1;i<=m;i++){
    		int x,y;
    		read(x),read(y);
    		PII t=Query(x,y);
    		write(t.first),putchar(' '),write(t.second),putchar('
    ');
    	}
    	return 0;
    }
    

    总结

    暂无。

  • 相关阅读:
    python的三个函数map、filter、reduce详解
    python 对文件进行切分行
    sed 命令匹配行操作
    停止学习框架
    JS中双等号“==”使用情况
    JavaScript中数组迭代方法
    JavaScript中replace()方法的第二个参数解析
    Vue生命周期
    JS中浮点数精度误差解决
    大量DOM操作的解决方案
  • 原文地址:https://www.cnblogs.com/Akmaey/p/14932402.html
Copyright © 2011-2022 走看看