bzoj 2827: 千山鸟飞绝

2827: 千山鸟飞绝

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 802  Solved: 228
[Submit][Status][Discuss]

Description

话说有一天doyouloveme和vfleaking到山里玩。谁知doyouloveme刚刚进山,所有的鸟儿竟被他的神犇气场给惊得全部飞走了。vfleaking顿时膜拜不已。

这时鸟王用鸟语说道:“[email protected]#$%……?”安抚了一下众鸟的情绪。鸟王生性好斗,作出了一个决定——要排鸟布阵把刚才吓到它们的人类赶出山去。

每只鸟都有一个编号,都有一个威武值。每秒钟鸟王都会发一个命令,编号为v的鸟飞到(x,y)去(坐标系原点是山顶,坐标单位为鸟爪)。鸟飞得很快,一秒之内就飞到了,可以看作是瞬间移动。如果编号为v的鸟和编号为u的鸟某一时刻处在同一位置,它们就会互相鼓励,增加各自的士气值和团结值。一只鸟的士气值等于此刻与它处在同一位置的鸟中的威武值的最大值,团结值等于此刻与它处在同一位置的鸟的只数。如果每一时刻都没有鸟与它处在同一位置,则士气值和团结值都为0。要注意自己不能鼓励自己,计算士气值和团结值时不能算上自己。

t秒钟后,doyouloveme目测出了现在每只鸟的战斗力,于是感叹了一句:“不妙,我们得走了。”

正所谓团结的鸟儿一个顶俩,所以doyouloveme这样描述战斗力:一只鸟战斗力值等于它在0到t秒中士气值的最大值与团结值的最大值的乘积。注意不是乘积的最大值,而是最大值的乘积。

vfleaking很想知道现在每只鸟的战斗力,但是他当然不会啦,于是他把这个任务交给了你来完成。

Input

第一行一个数n,代表鸟的只数。(鸟王那家伙你可以完全忽视掉)

接下来n行,每行三个整数w,x,y描述每只鸟的威武值和初始坐标。第i+1行描述编号为i的鸟。

接下来一行有一个数t,代表经过时间ts。

接下来t行,每行三个整数v,x,y描述鸟王每秒的命令。

Output

一共n行,每行一个数,代表每只鸟的战斗力。

Sample Input

5
1 1 1
3 1 2
4 4 4
2 0 1
2 2 3
5
1 1 2
2 4 4
2 4 3
3 0 1
5 0 1

Sample Output

3
4
6
8
8

HINT

对于样例的解释:

首先5只鸟的位置为(1,1),(1,2),(4,4),(0,1),(2,3),士气和团结值都是0。

鸟1飞到了(1,2),于是鸟1和鸟2互相鼓励,鸟1士气变为3,鸟2士气变为1。鸟1鸟2的团结值变为1。

然后鸟2飞到(4,4),与鸟3互相鼓励,鸟2士气变为4,鸟3士气变为3。鸟2与鸟3的团结值变为1。

鸟2然后飞到了(4,3),一个没有鸟的地方。于是士气和团结值都变为了0。

接下来鸟3和鸟5都飞到了鸟4的位置,于是三只鸟互相鼓励,鸟4、鸟5士气变为4,鸟3士气仍为3。鸟3、鸟4、鸟5的团结值都变为2。

于是大家的战斗力:

鸟1:3 * 1 = 3

鸟2:4 * 1 = 4

鸟3:3 * 2 = 6

鸟4:4 * 2 = 8

鸟5:4 * 2 = 8

1≤n≤30000   0≤t≤300000   坐标范围为整数,且不超过INT_MIN~INT_MAX

威武值为不超过INT_MAX的非负整数。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100010
#define M 500010
using namespace std;
int ch[M][2];
int rnd[M],key[M],sz[M],cnt[M],d1[M],d2[M],tot;
int z[N],rt[M],sum;
int mx[N],tog[N],w[N];
struct query{
    int x,y,id;
    void input(){scanf("%d%d%d",&id,&x,&y);}
}q[M];
struct point{
    int x,y;
    point(int x=0,int y=0):x(x),y(y){}
    bool operator < (const point &b)const{
        return x<b.x||x==b.x&&y<b.y;
    }
    bool operator == (const point &b)const{
        return x==b.x&&y==b.y;
    }
}p[N],san[M];
int cmp(int x,int val){
    if(w[key[x]]==w[val])return key[x]<val;
    return w[key[x]]<w[val];
}
int cmp2(int x,int val){
    if(key[x]==val)return -1;
    return w[key[x]]<=w[val];
}
void mark_down(int x,int u,int v){
    d1[x]=max(d1[x],u);
    d2[x]=max(d2[x],v);
    mx[key[x]]=max(mx[key[x]],u);
    tog[key[x]]=max(tog[key[x]],v);
}
void push_up(int x){
    sz[x]=cnt[x]+sz[ch[x][1]]+sz[ch[x][0]];
}
void push_down(int x){
    if(!d1[x]&&!d2[x]||!x)return;
    if(ch[x][0])mark_down(ch[x][0],d1[x],d2[x]);
    if(ch[x][1])mark_down(ch[x][1],d1[x],d2[x]);
    d1[x]=d2[x]=0;
}
void rot(int &x,int d){
    int k=ch[x][d^1];
    ch[x][d^1]=ch[k][d];
    ch[k][d]=x;
    push_up(x);push_up(k);
    x=k;
}
void insert(int &x,int val){
    if(x==0){
        tot++;
        ch[tot][0]=ch[tot][1]=0;
        key[tot]=val;
        rnd[tot]=rand();
        sz[tot]=cnt[tot]=1;
        x=tot;
        return;
    }
    push_down(x);
    int d=cmp(x,val);
    insert(ch[x][d],val);
    if(rnd[ch[x][d]]>rnd[x])rot(x,d^1);
    push_up(x);
}
int get_max(int x){
    if(!x)return -1;
    push_down(x);
    if(ch[x][1])return get_max(ch[x][1]);
    return w[key[x]];
}
void gao(int id,int i){//位置,编号
    if(rt[id]){
        mx[i]=max(get_max(rt[id]),mx[i]);
        mark_down(rt[id],w[i],0);
    }
    insert(rt[id],i);
    mark_down(rt[id],0,sz[rt[id]]-1);
}
int haxi(point v){
    return lower_bound(san+1,san+1+sum,v)-san;
}
void del(int &x,int val){
    if(x==0)return;
    push_down(x);
    int d=cmp(x,val);
    if(key[x]==val)d=-1;
    if(d==-1){
        if(ch[x][0]*ch[x][1]==0)x=ch[x][0]+ch[x][1];
        else{
            int d2=rnd[ch[x][0]]>rnd[ch[x][1]];
            push_down(ch[x][d2^1]);
            rot(x,d2);
            del(ch[x][d2],val);
        }
    }
    else del(ch[x][d],val);
    if(x)push_up(x);
}
int main(){
    int n;scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d%d%d",&w[i],&p[i].x,&p[i].y);
        san[++sum]=p[i];
    }
    int t;scanf("%d",&t);
    for(int i=0;i<t;i++){
        q[i].input();
        san[++sum]=point(q[i].x,q[i].y);
    }
    sort(san+1,san+sum+1);
    sum=unique(san+1,san+sum+1)-san-1;
    for(int i=1;i<=n;i++){//枚举每只鸟
        z[i]=haxi(p[i]);//z[i]是初始时的位置
        gao(z[i],i);
    }
    for(int i=0;i<t;i++){
        int u=q[i].id,xx=q[i].x,yy=q[i].y;
        del(rt[z[u]],u);
        z[u]=haxi(point(xx,yy));
        gao(z[u],u);
    }
    for(int i=1;i<=n;i++)del(rt[z[i]],i);
    for(int i=1;i<=n;i++)printf("%lld\n",1LL*mx[i]*tog[i]);
    return 0;
}

原文地址:https://www.cnblogs.com/thmyl/p/8278008.html

时间: 2024-10-11 07:15:04

bzoj 2827: 千山鸟飞绝的相关文章

2827: 千山鸟飞绝 非旋treap

国际惯例的题面:看起来很不可做的样子,我们先来整理一下题意吧.就是,维护每个点曾经拥有过的最大的两个属性值,支持把点的位置移动.我们用map对每个位置进行离散化,对每个位置建立一个平衡树.为了方便分离,这个平衡树的关键字是节点编号.然后把每个点当做一个节点,放进其所在位置的平衡树里.剩下要做的就是平衡树分离出一个点,合并一个点,还有打标记了.对于士气值的标记,我们维护平衡树中的max,每次合并的时候,用这个新点的威武值去给整棵树打标记,再用树中的max给这个新点打标记.团结值的标记,合并后一起打

BZOJ 1013: [JSOI2008]球形空间产生器sphere

二次联通门 : BZOJ 1013: [JSOI2008]球形空间产生器sphere /* BZOJ 1013: [JSOI2008]球形空间产生器sphere 高斯消元 QAQ SB的我也能终于能秒题了啊 设球心的坐标为(x,y,z...) 那么就可以列n+1个方程,化化式子高斯消元即可 */ #include <cstdio> #include <iostream> #include <cstring> #define rg register #define Max

bzoj 3309 DZY Loves Math - 莫比乌斯反演 - 线性筛

对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0. 给定正整数a,b,求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b). Input 第一行一个数T,表示询问数. 接下来T行,每行两个数a,b,表示一个询问. Output 对于每一个询问,输出一行一个非负整数作为回答. Sample Input 4 7558588 9653114 6514903 445

【BZOJ】[HNOI2009]有趣的数列

[算法]Catalan数 [题解] 学了卡特兰数就会啦>_<! 因为奇偶各自递增,所以确定了奇偶各自的数字后排列唯一. 那么就是给2n个数分奇偶了,是不是有点像入栈出栈序呢. 将做偶数标为-1,做奇数标为+1,显然当偶数多于奇数时不合法,因为它压不住后面的奇数. 然后其实这种题目,打表就可知啦--QAQ 然后问题就是求1/(n+1)*C(2n,n)%p了,p不一定是素数. 参考bzoj礼物的解法. 看到网上清一色的素数筛+分解质因数解法,不解了好久,感觉写了假的礼物-- 后来觉得礼物的做法才比

洛谷 P2709 BZOJ 3781 小B的询问

题目描述 小B有一个序列,包含N个1~K之间的整数.他一共有M个询问,每个询问给定一个区间[L..R],求Sigma(c(i)^2)的值,其中i的值从1到K,其中c(i)表示数字i在[L..R]中的重复次数.小B请你帮助他回答询问. 输入输出格式 输入格式: 第一行,三个整数N.M.K. 第二行,N个整数,表示小B的序列. 接下来的M行,每行两个整数L.R. 输出格式: M行,每行一个整数,其中第i行的整数表示第i个询问的答案. 输入输出样例 输入样例#1: 6 4 3 1 3 2 1 1 3

BZOJ 1012: [JSOI2008]最大数maxnumber(线段树)

012: [JSOI2008]最大数maxnumber Time Limit: 3 Sec  Memory Limit: 162 MB Description 现在请求你维护一个数列,要求提供以下两种操作:1. 查询操作.语法:Q L 功能:查询当前数列中末尾L个数中的最大的数,并输出这个数的值.限制:L不超过当前数列的长度.2. 插入操作.语法:A n 功能:将n加上t,其中t是最近一次查询操作的答案(如果还未执行过查询操作,则t=0),并将所得结果对一个固定的常数D取模,将所得答案插入到数列

【BZOJ】【1016】【JSOI2008】最小生成树计数

Kruskal/并查集+枚举 唉我还是too naive,orz Hzwer 一开始我是想:最小生成树删掉一条边,再加上一条边仍是最小生成树,那么这两条边权值必须相等,但我也可以去掉两条权值为1和3的,再加上权值为2和2的,不也满足题意吗?事实上,如果这样的话……最小生成树应该是1和2,而不是1和3或2和2!!! 所以呢?所以对于一个图来说,最小生成树有几条边权为多少的边,都是固定的!所以我们可以做一遍Kruskal找出这些边权,以及每种边权出现的次数.然后,对于每种边权,比方说出现了$v_i$

【BZOJ】【2844】albus就是要第一个出场

高斯消元解XOR方程组 srO  ZYF  Orz 膜拜ZYF…… http://www.cnblogs.com/zyfzyf/p/4232100.html 1 /************************************************************** 2 Problem: 2844 3 User: Tunix 4 Language: C++ 5 Result: Accepted 6 Time:252 ms 7 Memory:2052 kb 8 *******

[BZOJ 1066] [SCOI2007] 蜥蜴 【最大流】

题目链接:BZOJ - 1066 题目分析 题目限制了高度为 x 的石柱最多可以有 x 只蜥蜴从上面跳起,那么就可以用网络流中的边的容量来限制.我们把每个石柱看作一个点,每个点拆成 i1, i2,从 i1 到 i2 连一条边,容量为这个石柱 i 的高度,即跳跃次数限制.来到这个石柱就是向 i1 连边,从这个石柱跳起就是从 i2 向外连边,这样只要从石柱 i 跳起就一定会消耗 i1 到 i2 的边的容量.如果 i 有蜥蜴,就从 S 到 i1 连一条容量为 1 的边,如果从石柱 i 能跳出边界,就从