[BZOJ2843] 极地旅行社

题目描述

Description

不久之前,Mirko建立了一个旅行社,名叫“极地之梦”。这家旅行社在北极附近购买了N座冰岛,并且提供观光服务。当地最受欢迎的当然是帝企鹅了,这些小家伙经常成群结队的游走在各个冰岛之间。Mirko的旅行社遭受一次重大打击,以至于观光游轮已经不划算了。旅行社将在冰岛之间建造大桥,并用观光巴士来运载游客。Mirko希望开发一个电脑程序来管理这些大桥的建造过程,以免有不可预料的错误发生。这些冰岛从1到N标号。一开始时这些岛屿没有大桥连接,并且所有岛上的帝企鹅数量都是知道的。每座岛上的企鹅数量虽然会有所改变,但是始终在[0, 1000]之间。你的程序需要处理以下三种命令:
1."bridge A B"——在A与B之间建立一座大桥(A与B是不同的岛屿)。由于经费限制,这项命令被接受,当且仅当A与B不联通。若这项命令被接受,你的程序需要输出"yes",之后会建造这座大桥。否则,你的程序需要输出"no"。
2."penguins A X"——根据可靠消息,岛屿A此时的帝企鹅数量变为X。这项命令只是用来提供信息的,你的程序不需要回应。
3."excursion A B"——一个旅行团希望从A出发到B。若A与B连通,你的程序需要输出这个旅行团一路上所能看到的帝企鹅数量(包括起点A与终点B),若不联通,你的程序需要输出"impossible"。

Input

第一行一个正整数N,表示冰岛的数量。
第二行N个范围[0, 1000]的整数,为每座岛屿初始的帝企鹅数量。
第三行一个正整数M,表示命令的数量。接下来M行即命令,为题目描述所示。
1<=N<=30000,1<=M<=100000

Output

对于每个bridge命令与excursion命令,输出一行,为题目描述所示。

Sample Input

5
4 2 4 5 6
10
excursion 1 1
excursion 1 2
bridge 1 2
excursion 1 2
bridge 3 4
bridge 3 5
excursion 4 5
bridge 1 3
excursion 2 4
excursion 2 5

Sample Output

4
impossible
yes
6
yes
yes
15
yes
15
16

题目分析

一群dalao用lct艹掉了......
但是我是蒟蒻怎么会lct? 树链剖分即视感
先离线所有询问 按照要求将给定的森林建出来 然后对森林进行树链剖分
操作更新查询神马的都是基础操作了 就不说了.
最后用并查集判断连通性

此题坑点 建树时不能盲目加边 因为不按要求来做最终树会建成图..... 也就是说加边要判连通性
代码丑哭QAQ

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
int n,m;
int val[40000],head[30010],net[30010*2],to[30010*2],value[30010];
int tot,cnt;
void add(int x,int y)
{
    net[++tot]=head[x];
    to[tot]=y;
    head[x]=tot;
}
int f[30010];
int find(int x) 
{
    return f[x]==x?x:f[x]=find(f[x]);
}
char s[20];
int deep[30010],size[30010],son[30010],fa[30010],top[30010],number[30010];
void dfs(int x,int temp)
{
    size[x]=1;
    deep[x]=deep[temp]+1;
    fa[x]=temp;
    for(int i=head[x];i;i=net[i])
    {
    	if(to[i]!=temp)
    	{
    		dfs(to[i],x);
    		size[x]+=size[to[i]];
    		if(size[son[x]]<size[to[i]]) son[x]=to[i];
    	}
    }
}
void dfs2(int x,int temp)
{
    number[x]=++cnt;
    value[cnt]=val[x];
    top[x]=temp;
    if(son[x]) dfs2(son[x],temp);
    for(int i=head[x];i;i=net[i])
        if(to[i]!=fa[x]&&to[i]!=son[x]) dfs2(to[i],to[i]);
}
struct your
{
    int x,y,sum,tmp,ans;
}a[5*30010],q[400000];
void build(int dx,int dy,int num)
{
    a[num].x=dx,a[num].y=dy;
    if(dx==dy)
    {
        a[num].sum=value[dx];
        return;
    }
    int mid=(dx+dy)>>1;
    build(dx,mid,num<<1),build(mid+1,dy,num<<1|1);
    a[num].sum=a[num<<1].sum+a[num<<1|1].sum;
}
void update(int dx,int c,int num)
{
    if(dx==a[num].x&&dx==a[num].y) 
    {
        a[num].sum=c;
        return ;
    }
    int mid=(a[num].x+a[num].y)>>1;
    if(dx<=mid) update(dx,c,num<<1);
    else update(dx,c,num<<1|1);
    a[num].sum=a[num<<1].sum+a[num<<1|1].sum;
}
int ask(int dx,int dy,int num)
{
    if(dx==a[num].x&&dy==a[num].y) return a[num].sum;
    int mid=(a[num].x+a[num].y)>>1;
    if(dx>mid) return ask(dx,dy,num<<1|1);
    if(dy<=mid) return ask(dx,dy,num<<1);
    return ask(dx,mid,num<<1)+ask(mid+1,dy,num<<1|1);
}
int getsum(int x,int y)
{
    int sum=0;
    while(top[x]!=top[y])
    {
        if(deep[top[x]]<deep[top[y]]) swap(x,y);
        sum+=ask(number[top[x]],number[x],1);
        x=fa[top[x]];
    }
    if(deep[x]>deep[y]) swap(x,y);
    sum+=ask(number[x],number[y],1);
    return sum;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&val[i]);    
    scanf("%d",&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%s",&s[0]);
        int x,y;
        scanf("%d%d",&x,&y);
        if(s[0]=='b') q[i].tmp=1;
     	else if(s[0]=='p') q[i].tmp=2;
        else q[i].tmp=3;
    	q[i].x=x,q[i].y=y;
    }
    for(int i=1;i<=n;i++) f[i]=i;
    for(int i=1;i<=m;i++)
    {
        if(q[i].tmp==1) 
        {
            if(find(q[i].x)!=find(q[i].y)) 
                add(q[i].x,q[i].y),add(q[i].y,q[i].x),f[find(q[i].x)]=find(q[i].y),q[i].x=-1;
        	else q[i].x=-2;
        }
        else if(q[i].tmp==3) 
        {
                if(find(q[i].x)!=find(q[i].y)) q[i].ans=-1; 
            else q[i].ans=-2;
        }
    }
    for(int i=1;i<=n;i++)
        if(!deep[i]) dfs(i,0),dfs2(i,i);
    build(1,n,1);
    for(int i=1;i<=m;i++)
    {
        if(q[i].tmp==1) 
        {
            if(q[i].x==-2) printf("no\n");
            else printf("yes\n");
        }
        else if(q[i].tmp==2) update(number[q[i].x],q[i].y,1);
        else 
        {
            if(q[i].ans==-1) printf("impossible\n");
            else printf("%d\n",getsum(q[i].x,q[i].y));
        }
    }
}

发表评论

邮箱地址不会被公开。 必填项已用*标注