码迷,mamicode.com
首页 > 其他好文 > 详细

HDU 1540 Tunnel Warfare 平衡树 / 线段树:单点更新,区间合并

时间:2016-07-16 22:49:57      阅读:400      评论:0      收藏:0      [点我收藏+]

标签:

                                      Tunnel Warfare

                                 Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)


Problem Description
During the War of Resistance Against Japan, tunnel warfare was carried out extensively in the vast areas of north China Plain. Generally speaking, villages connected by tunnels lay in a line. Except the two at the ends, every village was directly connected with two neighboring ones.

Frequently the invaders launched attack on some of the villages and destroyed the parts of tunnels in them. The Eighth Route Army commanders requested the latest connection state of the tunnels and villages. If some villages are severely isolated, restoration of connection must be done immediately!
 

 

Input
The first line of the input contains two positive integers n and m (n, m ≤ 50,000) indicating the number of villages and events. Each of the next m lines describes an event.

There are three different events described in different format shown below:

D x: The x-th village was destroyed.

Q x: The Army commands requested the number of villages that x-th village was directly or indirectly connected with including itself.

R: The village destroyed last was rebuilt.
 

 

Output
Output the answer to each of the Army commanders’ request in order on a separate line.
 

 

Sample Input
7 9 D 3 D 6 D 5 Q 4 Q 5 R Q 4 R Q 4
 

 

Sample Output
1 0 2 4
 

 

Source
 
题解:区间合并基础题
技术分享
//1085422276
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<cmath>
#include<map>
#include<bitset>
#include<set>
#include<vector>
#include<stack>
///#include<bits/stdc++.h>
using namespace std ;
typedef long long ll;
#define mem(a) memset(a,0,sizeof(a))
#define meminf(a) memset(a,127,sizeof(a));
#define memfy(a) memset(a,-1,sizeof(a))
#define TS printf("111111\n");
#define FOR(i,a,b) for( int i=a;i<=b;i++)
#define FORJ(i,a,b) for(int i=a;i>=b;i--)
#define READ(a,b) scanf("%d%d",&a,&b)
#define mod 530600414
#define eps 0.0000000001
#define inf 1000000000.0
inline ll read()
{
    ll x=0,f=1;
    char ch=getchar();
    while(ch<0||ch>9)
    {
        if(ch==-)f=-1;
        ch=getchar();
    }
    while(ch>=0&&ch<=9)
    {
        x=x*10+ch-0;
        ch=getchar();
    }
    return x*f;
}
//****************************************

#define maxn 500000+5
struct ss
{
    int l,r,L,R,v;
}tr[maxn];
int n,m,lflag,rflag;
void pushup(int k)
{
    if(tr[k<<1].L+tr[k<<1|1].R==(tr[k].r-tr[k].l+1))
    {
        tr[k].L=tr[k].R=(tr[k].r-tr[k].l+1);
    }
    else {
        if(tr[k<<1].L>=(tr[k<<1].r-tr[k<<1].l+1))
            tr[k].L=tr[k<<1].L+tr[k<<1|1].L;
        else tr[k].L=tr[k<<1].L;
        if(tr[k<<1|1].R>=tr[k<<1|1].r-tr[k<<1|1].l+1)
            tr[k].R=tr[k<<1|1].R+tr[k<<1].R;
        else tr[k].R=tr[k<<1|1].R;
    }
}
void build(int k,int s,int t)
{
    tr[k].l=s;
    tr[k].r=t;
    tr[k].L=tr[k].R=t-s+1;
    if(s==t)
    {
        return ;
    }
    int mid=(s+t)>>1;
    build(k<<1,s,mid);
    build(k<<1|1,mid+1,t);
}
void update(int k,int x,int c)
{
    if(tr[k].l==x&&tr[k].r==x)
    {
        tr[k].L=tr[k].R=c;return ;
    }
    int mid=(tr[k].l+tr[k].r)>>1;
    if(x<=mid)update(k<<1,x,c);
    else update(k<<1|1,x,c);
    pushup(k);
}
int ask(int k,int x)
{
    if(tr[k].l==tr[k].r&&tr[k].l==x)
    {
        if(tr[k].L)rflag=lflag=1;
        return tr[k].L;
    }
    int mid=(tr[k].l+tr[k].r)>>1;
    int A=0,B=0;
    if(x<=mid)
    {
        A=ask(k<<1,x);
    }
    else {
         B=ask(k<<1|1,x);
    }
    int ans=A+B;
     if(A && rflag)
     {
         ans+=tr[k<<1|1].L;
         if(tr[k<<1|1].L<tr[k<<1|1].r-tr[k<<1|1].l+1)
            rflag=0;
     }
     if(B && lflag)
     {
         ans+=tr[k<<1].R;
         if(tr[k<<1].R<tr[k<<1].r-tr[k<<1].l+1)
            lflag=0;
     }
     return ans;
}
int main()
{
  char ch[321];
  int x;

    while(scanf("%d%d",&n,&m)!=EOF)
    {
        stack<int >q;
        int last=0;
        build(1,1,n);
        FOR(i,1,m)
        {
            scanf("%s",ch);
            if(ch[0]==D)
            {cin>>x;
            q.push(x);
                update(1,x,0);
            }
            else if(ch[0]==Q)
            {
                cin>>x;
                rflag=lflag=0;
               cout<<ask(1,x)<<endl;
            }
            else {
                last=q.top();
                q.pop();
                update(1,last,1);
            }
        }
    }
    return 0;
}
代码

 

补个 treap写法

技术分享
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<map>
using namespace std;
struct data{
    int l,r,v,size,rnd,w;
}tr[5200005];
int n,size,root,ans;
int last[5000000];
void update(int k)//更新结点信息
{
    tr[k].size=tr[tr[k].l].size+tr[tr[k].r].size+tr[k].w;
}
void rturn(int &k)
{
    int t=tr[k].l;tr[k].l=tr[t].r;tr[t].r=k;
    tr[t].size=tr[k].size;update(k);k=t;
}
void lturn(int &k)
{
    int t=tr[k].r;tr[k].r=tr[t].l;tr[t].l=k;
    tr[t].size=tr[k].size;update(k);k=t;
}
void insert(int &k,int x)
{
    if(k==0)
    {
        size++;k=size;
        tr[k].size=tr[k].w=1;tr[k].v=x;tr[k].rnd=rand();
        return;
    }
    tr[k].size++;
    if(tr[k].v==x)tr[k].w++;
    else if(x>tr[k].v)
    {
        insert(tr[k].r,x);
        if(tr[tr[k].r].rnd<tr[k].rnd)lturn(k);
    }
    else
    {
        insert(tr[k].l,x);
        if(tr[tr[k].l].rnd<tr[k].rnd)rturn(k);
    }
}
void del(int &k,int x)
{
    if(k==0)return;
    if(tr[k].v==x)
    {
        if(tr[k].w>1)
        {
            tr[k].w--;tr[k].size--;return;
        }
        if(tr[k].l*tr[k].r==0)k=tr[k].l+tr[k].r;
        else if(tr[tr[k].l].rnd<tr[tr[k].r].rnd)
            rturn(k),del(k,x);
        else lturn(k),del(k,x);
    }
    else if(x>tr[k].v)
        tr[k].size--,del(tr[k].r,x);
    else tr[k].size--,del(tr[k].l,x);
}

int query_rank(int k,int x)
{
    if(k==0)return 0;
    if(tr[k].v==x)return tr[tr[k].l].size+1;
    else if(x>tr[k].v)
        return tr[tr[k].l].size+tr[k].w+query_rank(tr[k].r,x);
    else return query_rank(tr[k].l,x);
}

int query_num(int k,int x)
{
    if(k==0)return 0;
    if(x<=tr[tr[k].l].size)
        return query_num(tr[k].l,x);
    else if(x>tr[tr[k].l].size+tr[k].w)
        return query_num(tr[k].r,x-tr[tr[k].l].size-tr[k].w);
    else return tr[k].v;
}


void query_pre(int k,int x)
{
    if(k==0)return;
    if(tr[k].v == x) {
        ans = k;
        return ;
    }
    if(tr[k].v<x)
    {
        ans=k;query_pre(tr[k].r,x);
    }
    else query_pre(tr[k].l,x);
}
void query_sub(int k,int x)
{
    if(k==0)return;
     if(tr[k].v == x) {
        ans = k;
        return ;
    }
    if(tr[k].v>x)
    {
        ans=k;query_sub(tr[k].l,x);
    }
    else query_sub(tr[k].r,x);
}
int m;
map<int,int > mp;
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF) {
        int cnt = 0 ,t1, t2 , x;
        root = 0;
        mp.clear();
        tr[root].size = 0;
        char ch[10];
        insert(root,0);
        insert(root,n+1);
        for(int i=1;i<=m;i++) {
            scanf("%s",ch);
            if(ch[0]==D) {
                scanf("%d",&x) ;
                if(mp[x]==0) {
                    insert(root,x);
                    mp[x] = 1;
                }last[++cnt] = x;

            }
            if(ch[0]==R) {
                if(!cnt) continue;
                if(mp[last[cnt]] == 1) {
                    del(root,last[cnt]);
                    mp[last[cnt]] = 0;
                }
                cnt--;
            }
            if(ch[0]==Q) {
                scanf("%d",&x);
                ans = 0;
                query_sub(root,x);
                t1= tr[ans].v ;
                ans = 0 ;
                query_pre(root,x) ;
                t2 = tr[ans].v;
                printf("%d\n", max(t1 -  t2 - 1,0));
            }
        }
    }
}
View Code

 

HDU 1540 Tunnel Warfare 平衡树 / 线段树:单点更新,区间合并

标签:

原文地址:http://www.cnblogs.com/zxhl/p/4843726.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!