BZOJ3083 遥远的国度

本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作。

 

 

本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!

 

题目链接:BZOJ3083

正解:树链剖分+线段树+分类讨论

解题报告:

  路径修改$+$查询就用链剖资瓷就好了,问题就在于换根操作怎么办。

  我们的做法是不真正的换根,而是选择分类讨论。

  假设每次的询问点为$x$,当前根为$rt$,如果$rt$在$x$子树中,那么实际的子树就是除了含$rt$的这棵$x$的子树之外的所有部分;

  如果$rt$就是$x$,那么子树是全部;

  否则的话$x$的子树不会发生变化。

  分类讨论一下就做完了==

  修改操作用线段树维护,上一个$set$标记就好了。

 

//It is made by ljh2000
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。
#include <algorithm>
#include <iostream>
#include <cstring>
#include <vector>
#include <cstdio>
#include <string>
#include <queue>
#include <cmath>
#include <ctime>
#define lc root<<1
#define rc root<<1|1
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define reg(i,x) for(int i=first[x];i;i=next[i])
using namespace std;
typedef long long LL;
const int MAXN = 200011;
const int MAXM = 200011; 
int n,m,rt,ecnt,first[MAXN],to[MAXM],nxt[MAXM],deep[MAXN],top[MAXN],size[MAXN],son[MAXN],father[MAXN],id[MAXN],end[MAXN],pre[MAXN],f[MAXN][18],bel[MAXN];
LL ans,tag[MAXN*3],a[MAXN*3],vv[MAXN];
bool hav[MAXN*3];//!!!
inline void link(int x,int y) { nxt[++ecnt]=first[x]; first[x]=ecnt; to[ecnt]=y; }

inline int getint(){
    int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
    if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
}

inline void pushdown(int root,int l,int r){ 
	if(!hav[root] || l==r)return ;
	tag[lc]=tag[rc]=a[lc]=a[rc]=tag[root]; 
	hav[lc]=hav[rc]=true;
	hav[root]=false; tag[root]=0; 
}

inline void build(int root,int l,int r){
	if(l==r) { a[root]=vv[ pre[l] ]; return ; }
	int mid=(l+r)>>1; build(lc,l,mid); build(rc,mid+1,r);
	a[root]=min(a[lc],a[rc]);
}

inline void modify(int root,int l,int r,int ql,int qr,LL val){
	pushdown(root,l,r); 
	if(ql<=l && r<=qr) { tag[root]=val; hav[root]=true; a[root]=val; return ; }
	int mid=(l+r)>>1;
	if(ql<=mid) modify(lc,l,mid,ql,qr,val);
	if(qr>mid) modify(rc,mid+1,r,ql,qr,val);
	a[root]=min(a[lc],a[rc]);
}

inline LL query(int root,int l,int r,int ql,int qr){
	pushdown(root,l,r);
	if(ql<=l && r<=qr) return a[root];
	int mid=(l+r)>>1;
	if(ql>mid) return query(rc,mid+1,r,ql,qr);
	else if(qr<=mid) return query(lc,l,mid,ql,qr);
	else return min( query(lc,l,mid,ql,qr) , query(rc,mid+1,r,ql,qr) );
}

inline void dfs(int x,int fa){
	size[x]=1;
	for(int i=first[x];i;i=nxt[i]) {
		int v=to[i]; if(v==fa) continue;
		deep[v]=deep[x]+1; f[v][0]=x; father[v]=x; dfs(v,x);
		size[x]+=size[v]; if(size[v]>size[son[x]]) son[x]=v;
	}
}

inline void dfs2(int x,int fa){
	id[x]=++ecnt; pre[ecnt]=x;
	if(son[x]) { top[son[x]]=top[x]; dfs2(son[x],x); }
	for(int i=first[x];i;i=nxt[i]) {
		int v=to[i]; if(v==fa || v==son[x]) continue;
		top[v]=v; dfs2(v,x);
	}
	end[x]=ecnt;
}

inline void lca(int x,int y,LL z){
	int f1=top[x],f2=top[y];
	while(f1!=f2) {
		if(deep[f1]<deep[f2]) swap(f1,f2),swap(x,y);
		modify(1,1,n,id[f1],id[x],z);
		x=father[f1]; f1=top[x];
	}
	if(deep[x]<deep[y]) swap(x,y);
	modify(1,1,n,id[y],id[x],z);
}

inline int getF(int x,int F){
	for(int i=bel[deep[x]];i>=0;i--) 
		if(deep[f[x][i]]>deep[F]) 
			x=f[x][i];
	return x;
}

inline void work(){
	n=getint(); m=getint(); int x,y,type,pos; LL z;
	for(int i=1;i<n;i++) { x=getint(); y=getint(); link(x,y); link(y,x); }
	for(int i=1;i<=n;i++) scanf("%lld",&vv[i]);
	rt=getint(); 
	deep[rt]=1; ecnt=0; dfs(rt,0);
	top[rt]=rt; dfs2(rt,0);
	for(int i=2;i<=n;i++) bel[i]=bel[i>>1]+1;
	for(int j=1;j<=17;j++) for(int i=1;i<=n;i++) f[i][j]=f[f[i][j-1]][j-1];
	build(1,1,n);
	while(m--) {
		type=getint(); x=getint();
		if(type==1) rt=x;
		else if(type==2) {
			y=getint(); scanf("%lld",&z);
			lca(x,y,z);
		}
		else {
			if(x==rt) ans=query(1,1,n,1,n);
			else if(id[x]<=id[rt] && end[rt]<=end[x]){
				ans=1LL<<60; 
				pos=getF(rt,x);//!!!
				if(id[pos]>1) ans=query(1,1,n,1,id[pos]-1);//!!!
				if(end[pos]<n) ans=min(ans,query(1,1,n,end[pos]+1,n));//!!!
			}
			else ans=query(1,1,n,id[x],end[x]);
			printf("%lld\n",ans);
		}
	}
}

int main()
{
#ifndef ONLINE_JUDGE
	freopen("3083.in","r",stdin);
	freopen("3083.out","w",stdout);
#endif
    work();
    return 0;
}
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。

  

 make:

//It is made by ljh2000
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <string>
#include <complex>
using namespace std;
typedef long long LL;
typedef long double LB;
typedef complex<double> C;
const double pi = acos(-1);

inline int getint(){
    int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
    if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
}

int main()
{
    freopen("3083.in","w",stdout);
    srand(time(NULL));
	int n,m,size,type;
	n=10; m=10;
	//n=10000; m=10000;
	size=10000000;
	//n=6; m=6;
	//size=100;
	printf("%d %d\n",n,m);
	for(int i=2;i<=n;i++) printf("%d %d\n",rand()%(i-1)+1,i);
	for(int i=1;i<=n;i++) printf("%d ",rand()%size+1);
	printf("\n%d\n",rand()%n+1);
	for(int i=1;i<=m;i++) {
		type=rand()%3+1; printf("%d ",type);
		if(type==1) {
			printf("%d\n",rand()%n+1);
		}
		else if(type==2) {
			printf("%d %d %d\n",rand()%n+1,rand()%n+1,rand()%size+1);
		}
		else {
			printf("%d\n",rand()%n+1);
		}
	}
    return 0;
}
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。

  

 

 pai:

//It is made by ljh2000
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <string>
#include <complex>
#include <bitset>
using namespace std;
typedef long long LL;
typedef long double LB;
typedef complex<double> C;
const double pi = acos(-1);

inline int getint(){
    int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
    if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
}

int main()
{
    int t=0;
	while(++t) {
		printf("%d : ",t);
		system("make.exe");
		system("3083.exe");
		system("2.exe");
		if(system("fc 3083.out 2.out")){ printf("WA\n"); break;  }
		else printf("AC\n");
	}
    return 0;
}
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。

  

 

 数据:

 1.in

10 10
1 2
1 3
1 4
2 5
1 6
4 7
2 8
6 9
7 10
22411 20657 24240 21525 12974 17337 3717 26377 23217 20558 
2
2 5 6 15012
3 5
3 9
3 7
1 3
3 10
3 2
2 8 1 29431
1 8
2 7 6 30507

1.out

15012
23217
3717
20558
15012

  

 

2.in

6 6
1 2
1 3
1 4
3 5
5 6
38 66 42 98 40 49 
6
1 4
3 5
3 3
3 3
3 5
3 5

  

2.out

40
40
40
40
40

  

posted @ 2017-05-05 09:44  ljh_2000  阅读(152)  评论(0编辑  收藏