bzoj 2648 SJY摆棋子 cdq分治+树状数组

题面

题目传送门

解法

bzoj2716

自己cdq写的还是丑啊,别人A掉了我T飞了

代码

#include <bits/stdc++.h>
#define inf 1 << 30
#define N 1000010
using namespace std;
template <typename node> void chkmax(node &x, node y) {x = max(x, y);}
template <typename node> void chkmin(node &x, node y) {x = min(x, y);}
template <typename node> void read(node &x) {
    x = 0; int f = 1; char c = getchar();
    while (!isdigit(c)) {if (c == ‘-‘) f = -1; c = getchar();}
    while (isdigit(c)) x = x * 10 + c - ‘0‘, c = getchar(); x *= f;
}
struct Node {
    int op, x, y, tim, id;
} a[N], t[N];
int mx, f[N], ans[N];
bool cmp1(Node a, Node b) {return a.x < b.x;}
bool cmp2(Node a, Node b) {return a.y > b.y;}
bool cmp3(Node a, Node b) {return a.x > b.x;}
int lowbit(int x) {return x & -x;}
void modify(int x, int v) {
    for (int i = x; i <= mx; i += lowbit(i))
        if (v == -inf) f[i] = v;
            else chkmax(f[i], v);
}
int query(int x) {
    int ret = -inf;
    for (int i = x; i; i -= lowbit(i))
        ret = max(ret, f[i]);
    return ret;
}
void cdq1(int l, int r) {
    if (l >= r) return;
    int mid = (l + r) >> 1, tx = l, ty = mid + 1;
    for (int i = l; i <= r; i++)
        if (a[i].tim <= mid) t[tx++] = a[i];
            else t[ty++] = a[i];
    for (int i = l; i <= r; i++) a[i] = t[i];
    tx = l, ty = mid + 1;
    while (ty <= r) {
        while (tx <= mid && a[tx].x <= a[ty].x) {
            if (a[tx].op == 1) modify(a[tx].y, a[tx].x + a[tx].y);
            tx++;
        }
        if (a[ty].op == 2) chkmin(ans[a[ty].id], a[ty].x + a[ty].y - query(a[ty].y));
        ty++;
    }
    for (int i = l; i < tx; i++) if (a[i].op == 1) modify(a[i].y, -inf);
    cdq1(l, mid); cdq1(mid + 1, r);
}
void cdq2(int l, int r) {
    if (l >= r) return;
    int mid = (l + r) >> 1, tx = l, ty = mid + 1;
    for (int i = l; i <= r; i++)
        if (a[i].tim <= mid) t[tx++] = a[i];
            else t[ty++] = a[i];
    for (int i = l; i <= r; i++) a[i] = t[i];
    tx = l, ty = mid + 1;
    while (ty <= r) {
        while (tx <= mid && a[tx].y >= a[ty].y) {
            if (a[tx].op == 1) modify(a[tx].x, a[tx].x - a[tx].y);
            tx++;
        }
        if (a[ty].op == 2) chkmin(ans[a[ty].id], a[ty].x - a[ty].y - query(a[ty].x));
        ty++;
    }
    for (int i = l; i < tx; i++) if (a[i].op == 1) modify(a[i].x, -inf);
    cdq2(l, mid); cdq2(mid + 1, r);
}
void cdq3(int l, int r) {
    if (l >= r) return;
    int mid = (l + r) >> 1, tx = l, ty = mid + 1;
    for (int i = l; i <= r; i++)
        if (a[i].tim <= mid) t[tx++] = a[i];
            else t[ty++] = a[i];
    for (int i = l; i <= r; i++) a[i] = t[i];
    tx = l, ty = mid + 1;
    while (ty <= r) {
        while (tx <= mid && a[tx].x >= a[ty].x) {
            if (a[tx].op == 1) modify(a[tx].y, a[tx].y - a[tx].x);
            tx++;
        }
        if (a[ty].op == 2) chkmin(ans[a[ty].id], a[ty].y - a[ty].x - query(a[ty].y));
        ty++;
    }
    for (int i = l; i < tx; i++) if (a[i].op == 1) modify(a[i].y, -inf);
    cdq3(l, mid); cdq3(mid + 1, r);
}
void cdq4(int l, int r) {
    if (l >= r) return;
    int mid = (l + r) >> 1, tx = l, ty = mid + 1;
    for (int i = l; i <= r; i++)
        if (a[i].tim <= mid) t[tx++] = a[i];
            else t[ty++] = a[i];
    for (int i = l; i <= r; i++) a[i] = t[i];
    tx = l, ty = mid + 1;
    while (ty <= r) {
        while (tx <= mid && a[tx].x >= a[ty].x) {
            if (a[tx].op == 1) modify(mx - a[tx].y + 1, -a[tx].x - a[tx].y);
            tx++;
        }
        if (a[ty].op == 2) chkmin(ans[a[ty].id], -a[ty].x - a[ty].y - query(mx - a[ty].y + 1));
        ty++;
    }
    for (int i = l; i < tx; i++) if (a[i].op == 1) modify(mx - a[i].y + 1, -inf);
    cdq4(l, mid); cdq4(mid + 1, r);
}
int main() {
    int n, m, tot = 0; read(n), read(m);
    for (int i = 1; i <= n; i++) {
        int x, y; read(x), read(y); chkmax(mx, max(x, y));
        a[++tot] = (Node) {1, x, y, tot, 0};
    }
    int cntq = 0;
    for (int i = 1; i <= m; i++) {
        read(a[++tot].op), read(a[tot].x), read(a[tot].y);
        chkmax(mx, max(a[tot].x, a[tot].y));
        if (a[tot].op == 2) a[tot].id = ++cntq; a[tot].tim = tot;
    }
    for (int i = 1; i <= cntq; i++) ans[i] = inf;
    for (int i = 1; i <= mx; i++) f[i] = -inf;
    sort(a + 1, a + tot + 1, cmp1); cdq1(1, tot);
    for (int i = 1; i <= mx; i++) f[i] = -inf;
    sort(a + 1, a + tot + 1, cmp2); cdq2(1, tot);
    for (int i = 1; i <= mx; i++) f[i] = -inf;
    sort(a + 1, a + tot + 1, cmp3); cdq3(1, tot);
    for (int i = 1; i <= mx; i++) f[i] = -inf;
    sort(a + 1, a + tot + 1, cmp3); cdq4(1, tot);
    for (int i = 1; i <= cntq; i++) cout << ans[i] << "\n";
    return 0;
}

原文地址:https://www.cnblogs.com/copperoxide/p/9478329.html

时间: 2024-10-01 02:20:41

bzoj 2648 SJY摆棋子 cdq分治+树状数组的相关文章

BZOJ 2244: [SDOI2011]拦截导弹 [CDQ分治 树状数组]

传送门 题意:三维最长不上升子序列以及每个元素出现在最长不上升子序列的概率 $1A$了好开心 首先需要从左右各求一遍,长度就是$F[0][i]+F[1][i]-1$,次数就是$G[0][i]*G[1][i]$ 我们可以用一些转换来简化代码 反转之后变成$LIS$,然后再反转并且$x,y$取反还是$LIS$,写一遍就可以啦 然后本题的树状数组需要维护最大值以及最大值的数量,还有一个时间戳 #include <iostream> #include <cstdio> #include &

BZOJ 2648 SJY摆棋子(KD树)

[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2716 [题目大意] 给出一些点,同时不断插入点和询问某点离插入点最近距离 [题解] 我们对于给定的点直接建树,之后动态插入查询即可,重建会超时, 直接插入就可以过了 [代码] #include <cstdio> #include <algorithm> using namespace std; const int N=1500000,INF=1e9; inline in

BZOJ 2683 简单题 cdq分治+树状数组

题意:链接 **方法:**cdq分治+树状数组 解析: 首先对于这道题,看了范围之后,二维的数据结构是显然不能过的,于是我们可能会考虑把一维排序之后另一位上数据结构什么的,然而cdq分治却能够很好的体现它的作用. 首先,对于每一个询问求和,显然是x在它左边的并且出现时间在它之前的所有的change对他可能会有影响. 我们按照x第一关键字,y第二关键字,操作第三关键字来排序所有的询问,然后在cdq的时候,每次递归处理左半区间,按照x动态的将y这一列的值加到树状数组里,来更新右半边的所有询问,注意这

HDU 5618:Jam&#39;s problem again(CDQ分治+树状数组处理三维偏序)

http://acm.hdu.edu.cn/showproblem.php?pid=5618 题意:-- 思路:和NEUOJ那题一样的.重新写了遍理解了一下,算作处理三维偏序的模板了. 1 #include <cstdio> 2 #include <algorithm> 3 #include <iostream> 4 #include <cstring> 5 using namespace std; 6 #define INF 0x3f3f3f3f 7 #d

XJOI NOIP2015模拟赛Day1 T2 ctps bitset优化 或 排序+cdq分治+树状数组+平衡树

题意: 4维空间中有1个点集A,|A|=n,用(a,b,c,d)表示每个点. 共有m个询问,每次询问输入一个点(a,b,c,d),求最大的S,其中S={p|p∈A且ap<=a,bp<=b,cp<=c,dp<=d},输出|S| 输入格式: 第一行n 接下来n行有n个4维点对 第n+2行有一个数m 再接下来m行每行有一个四维点对,表示每个询问 输出格式: 对于每个询问输出一个数 **方法:**bitset优化 或 排序+cdq分治+树状数组+平衡树 解析: 神题,考场不会,暴力骗40,

【BZOJ4553】[Tjoi2016&amp;Heoi2016]序列 cdq分治+树状数组

[BZOJ4553][Tjoi2016&Heoi2016]序列 Description 佳媛姐姐过生日的时候,她的小伙伴从某宝上买了一个有趣的玩具送给他.玩具上有一个数列,数列中某些项的值可能会变化,但同一个时刻最多只有一个值发生变化.现在佳媛姐姐已经研究出了所有变化的可能性,她想请教你,能否选出一个子序列,使得在任意一种变化中,这个子序列都是不降的?请你告诉她这个子序列的最长长度即可.注意:每种变化最多只有一个值发生变化.在样例输入1中,所有的变化是: 1 2 3 2 2 3 1 3 3 1

BZOJ2683: 简单题(CDQ分治 + 树状数组)

BZOJ2683: 简单题(CDQ分治 + 树状数组) 题意: 你有一个\(N*N\)的棋盘,每个格子内有一个整数,初始时的时候全部为\(0\),现在需要维护两种操作: 命令 参数限制 内容 \(1\ x\ y\ A\) \(1\le x,y \le N\),A是正整数 将格子\(x,y\)里的数字加上\(A\) \(2\ x1\ y1\ x2\ y2\) \(1\le x1\le x2\le N,1\le y1\le y2\le N\) 输出\(x1\ y1\ x2\ y2\)这个矩形内的数字

BZOJ 2648(SJY摆棋子-KD_Tree)

2648: SJY摆棋子 Time Limit: 20 Sec  Memory Limit: 128 MB Submit: 1180  Solved: 391 [Submit][Status][Discuss] Description 这天,SJY显得无聊.在家自己玩.在一个棋盘上,有N个黑色棋子.他每次要么放到棋盘上一个黑色棋子,要么放上一个白色棋子,如果是白色棋子,他会找出距离这个白色棋子最近的黑色棋子.此处的距离是 曼哈顿距离 即(|x1-x2|+|y1-y2|) .现在给出N<=5000

[cdq分治][树状数组] Bzoj P3262 陌上花开

Description 有n朵花,每朵花有三个属性:花形(s).颜色(c).气味(m),用三个整数表示. 现在要对每朵花评级,一朵花的级别是它拥有的美丽能超过的花的数量. 定义一朵花A比另一朵花B要美丽,当且仅Sa>=Sb,Ca>=Cb,Ma>=Mb. 显然,两朵花可能有同样的属性.需要统计出评出每个等级的花的数量. Input 第一行为N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分别表示花的数量和最大属性值. 以下N