bzoj 2827: 千山鸟飞绝

该博客介绍了如何运用离散化和Splay数据结构解决一道关于鸟的威武值与士气值计算的问题。文章阐述了题意,即根据鸟的位置变化计算每只鸟士气值和团结值的最大乘积,并提供了题解思路和相应的代码实现。

题意:有n只鸟,每只鸟有个威武值。假如几只鸟在同一个位置,就可以互相鼓励,获得士气值(同一位置的其他鸟的威武值的最大值)和团结值(同一位置其他鸟的数量)。鸟的位置会改变,问每只鸟士气值的最大值和团结值的最大值的乘积。
题解:离散化+splay。每个点建一个splay,键值为威武值,同时记录士气值和团结值。移动鸟时,这只鸟的士气值就变成目标点splay威武值的最大值,士气值就是目标点splay的大小;目标splay打上标记就好了。
代码:

#include<bits/stdc++.h>
#define pa pair<int,int>
using namespace std;

int n,t,root[330010],w[30010];
struct order
{
    int v;
    pa pos;
}b[300010];
struct pnt
{
    int f,s[2],sz,c,c1,c2;//威武 士气 团结 (最大值) 
    int pos,cx,cy;//士气 团结 (当前值) 
}p[30010];
pa pos[30010];
map<pa,int>lsh;

void print(int x)
{
    printf("x=%d f=%d s={%d %d} sz=%d c=%d c1=%d c2=%d pos=%d cx=%d cy=%d\n",x,p[x].f,p[x].s[0],p[x].s[1],p[x].sz,p[x].c,p[x].c1,p[x].c2,p[x].pos,p[x].cx,p[x].cy);
}
void pushdown(int x)
{
    int lc=p[x].s[0],rc=p[x].s[1];
    if(lc)
    {
        p[lc].c1=max(p[lc].c1,p[x].cx);
        p[lc].c2=max(p[lc].c2,p[x].cy);
        if(p[x].cx!=-1)
        {
            p[lc].cx=max(p[lc].cx,p[x].cx);
            p[lc].cy=max(p[lc].cy,p[x].cy);
        }
    }
    if(rc)
    {
        p[rc].c1=max(p[rc].c1,p[x].cx);
        p[rc].c2=max(p[rc].c2,p[x].cy);
        if(p[x].cx!=-1)
        {
            p[rc].cx=max(p[rc].cx,p[x].cx);
            p[rc].cy=max(p[rc].cy,p[x].cy);
        }
    }
    p[x].cx=p[x].cy=-1;
}
void pushup(int x)
{
    int lc=p[x].s[0],rc=p[x].s[1];
    p[x].sz=1;
    if(lc)
    p[x].sz+=p[lc].sz;
    if(rc)
    p[x].sz+=p[rc].sz;
}
void rot(int x)
{
    int f=p[x].f,ff=p[f].f,R,r,w=x==p[f].s[0];
//  printf("%d %d\n",f,ff);
    R=f,r=p[x].s[w];
    p[R].s[1-w]=r;
    if(r)
    p[r].f=R;
    R=ff,r=x;
    p[R].s[f==p[ff].s[1]]=r;
    p[r].f=R;
    R=x,r=f;
    p[R].s[w]=r;
    p[r].f=R;
    pushup(f);
}
void preserve(int x)
{
    static stack<int>sta;
    while(x)
    {
        sta.push(x);
        x=p[x].f;
    }
    while(!sta.empty())
    {
        pushdown(sta.top());
        sta.pop();
    }
}
void splay(int x,int rt)
{
    preserve(x);/*
        for(int i=1;i<=n;i++)
        print(i);*/
    while(p[x].f!=rt)
    {
        int f=p[x].f,ff=p[f].f;
        if(ff==rt)
        {
            rot(x);
            continue;
        }
        if((p[ff].s[0]==f)==(p[f].s[0]==x))
        {
            rot(f);
            rot(x);
        }
        else
        {
            rot(x);
            rot(x);
        }/*
        for(int i=1;i<=n;i++)
        print(i);
        system("pause");*/
    }
    pushup(x);
    if(rt==0)
    root[p[x].pos]=x;
}
int getmax(int pos)
{
    int x=root[pos];
    while(p[x].s[1])
    {
        pushdown(x);
        x=p[x].s[1];
    }
    return x;
}
void add(int x,int f)
{
    p[x].f=f;
    if(p[x].c<p[f].c)
    p[f].s[0]=x;
    else
    p[f].s[1]=x;
    pushup(f);
    splay(x,0);
}
void insert(int x,int pos)
{
//  printf("%d %d\n",x,pos);
    p[x].pos=pos;
    if(!root[pos])
    {
        root[pos]=x;
        add(x,0);
        p[x].pos=pos;
        return;
    }
    int now=root[pos],f;
    while(now)
    {
        pushdown(now);
        f=now;
        if(p[x].c<p[now].c)
        now=p[now].s[0];
        else
        now=p[now].s[1];
    }
    int cx=p[getmax(pos)].c,cy=p[root[pos]].sz;
    p[x].c1=max(p[x].c1,cx);
    p[x].c2=max(p[x].c2,cy);
    p[x].cx=p[x].cy=-1;
    add(x,f);
    p[x].cx=p[x].c;
    p[x].cy=cy;
}
void del(int x)
{
    splay(x,0);
    if(!p[x].s[0]&&!p[x].s[1])
    {
        root[p[x].pos]=0;
    }
    else if(p[x].s[0]&&!p[x].s[1])
    {
        root[p[x].pos]=p[x].s[0];
        p[p[x].s[0]].f=0;
    }
    else if(!p[x].s[0]&&p[x].s[1])
    {
        root[p[x].pos]=p[x].s[1];
        p[p[x].s[1]].f=0;
    }
    else
    {
        int l,r;
        pushdown(x);
        l=p[x].s[0];
        while(p[l].s[1])
        {
            pushdown(l);
            l=p[l].s[1];
        }
        splay(l,x);
        p[l].s[1]=p[x].s[1];
        if(p[x].s[1])
        p[p[x].s[1]].f=l;
        p[l].f=0;
        root[p[x].pos]=l;
    }
    pushup(root[p[x].pos]);
    p[x].pos=0;
    p[x].s[0]=p[x].s[1]=0;
    p[x].sz=1;
    p[x].cx=p[x].cy=-1;
}
void wk(int x)
{
    if(!x)
    return;
    pushdown(x);
    wk(p[x].s[0]);
    wk(p[x].s[1]);
}
int main()
{
    memset(root,0,sizeof(root));
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d%d",&w[i],&pos[i].first,&pos[i].second);
        lsh[pos[i]]=1;
    }
    scanf("%d",&t);
    for(int i=1;i<=t;i++)
    {
        scanf("%d%d%d",&b[i].v,&b[i].pos.first,&b[i].pos.second);
        lsh[b[i].pos]=1;
    }
    int nn=0;
    for(map<pa,int>::iterator it=lsh.begin();it!=lsh.end();it++)
    it->second=++nn;/*
    for(map<pa,int>::iterator it=lsh.begin();it!=lsh.end();it++)
    printf("(%d,%d)=%d\n",it->first.first,it->first.second,it->second);*/
    for(int i=1;i<=n;i++)
    {
        p[i].c=w[i];
        p[i].cx=p[i].cy=-1;
        p[i].sz=1;
        p[i].pos=lsh[pos[i]];
        insert(i,p[i].pos);/*
        puts("root:");
        for(int i=1;i<=nn;i++)
        printf("%d ",root[i]);
        puts("\np");
        for(int i=1;i<=n;i++)
        print(i);*/
    }
//  puts("=======");
    for(int i=1;i<=t;i++)
    {
        del(b[i].v);
//      puts("r");
        insert(b[i].v,lsh[b[i].pos]);/*
        puts("root:");
        for(int i=1;i<=nn;i++)
        printf("%d ",root[i]);
        puts("\np");
        for(int i=1;i<=n;i++)
        print(i);*/
    }
    for(int i=1;i<=n;i++)
//  wk(i);
    splay(i,0);
    for(int i=1;i<=n;i++)
    printf("%lld\n",(long long)p[i].c1*p[i].c2);
}
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值