2017暑期集训Day 25 树状数组

最新推荐文章于 2024-07-01 20:23:08 发布
最新推荐文章于 2024-07-01 20:23:08 发布
阅读量342 收藏
点赞数
分类专栏： SDU_ACM_Practice 文章标签： poj
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Yancy_/article/details/76934143
版权

poj 2352 && hdu 1541 Stars

[Solution]

对于每个点，询问在其左下角的点的数量，我们按照y方向排序，这样按照升序方向遍历可以确保后面的点的纵坐标比前面的大，我们把前面的点的横坐标存到树状数组里，这样每次取出横坐标比其小的数量即可

[Code]

#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
using namespace std;
#define N 52505
int n, m, l, r;
int a[N], c[N];
int lowbit(int x) {  return x & -x; }
int sum(int x)
{
    int ret = 0;
    while(x > 0){
        ret += c[x];
        x -= lowbit(x);
    }
    return ret;
}
void add(int x)
{
    while(x <= 32005){
        c[x]++;
        x += lowbit(x);
    }
}
int main()
{
    //freopen("b.in", "r", stdin);
    while(~scanf("%d", &n) && n){
        for(int i = 0; i <= 32005; i++)  {  c[i] = 0;  a[i] = 0;  }
        for(int i = 1; i <= n; i++){
            int x, y;
            scanf("%d%d", &x, &y);
            x++;
            y = sum(x);
            a[y]++;
            add(x);
        }
        for(int i = 0; i < n; i++)  printf("%d\n", a[i]);
    }
    return 0;
}

HDU 2838-Cow Sorting

[Solution]

树状数组求解逆序对

[Code]

#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
using namespace std;
typedef long long ll;
#define N 102505
int n, m, l, r, delta;
long long c[N], d[N];
int a[N];
int lowbit(int x) {  return x & -x; }
ll sum(int x)
{
    ll ret = 0;
    while(x > 0){
        ret += c[x];
        x -= lowbit(x);
    }
    return ret;
}
void add(int x)
{
    while(x <= n){
        c[x] += delta;
        x += lowbit(x);
    }
}
ll sumd(int x)
{
    ll ret = 0;
    while(x > 0){
        ret += d[x];
        x -= lowbit(x);
    }
    return ret;
}
void addd(int x)
{
    while(x <= n){
        d[x]++;
        x += lowbit(x);
    }
}
int main()
{
   // freopen("b.in", "r", stdin);
    scanf("%d", &n);
    for(int i = 1; i <= n; i++)  scanf("%d", a + i);
    long long ans = 0;
    for(int i = n; i >= 1; i--){
        ans += 1LL * a[i] * sumd(a[i] - 1) + sum(a[i] - 1);
        delta = a[i];
        add(a[i]);
        addd(a[i]);
    }
    cout << ans << endl;
    return 0;
}

HDU - 1556 Color the ball

[Solution]

区间修改，单点查询

[Code]

#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
using namespace std;
#define N 102505
int n, m, l, r, delta;
int a[N], c[N];
int lowbit(int x) {  return x & -x; }
int sum(int x)
{
    int ret = 0;
    while(x > 0){
        ret += c[x];
        x -= lowbit(x);
    }
    return ret;
}
void add(int x)
{
    while(x <= n){
        c[x] += delta;
        x += lowbit(x);
    }
}
int main()
{
   // freopen("b.in", "r", stdin);
    while(~scanf("%d", &n) && n){
        for(int i = 1; i <= n; i++)  c[i] = 0;
        for(int i = 1; i <= n; i++){
            scanf("%d%d", &l, &r);
            delta = 1;
            add(l);
            if (r + 1 <= n){
                delta = -1;
                add(r + 1);
            }
        }
        for(int i = 1; i < n; i++)
            printf("%d ", sum(i));
        printf("%d\n", sum(n));
    }
    return 0;
}

POJ 2299 Ultra-QuickSort (归并排序求逆序数)

[Solution]

需要离散化

[Code]

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
#define N 502000
typedef long long ll;
int n, m, l, r;
int a[N], c[N];
struct node{
  int w, num;
} b[N];
int lowbit(int x) {  return x & -x; }
int sum(int x)
{
    int ret = 0;
    while(x > 0){
        ret += c[x];
        x -= lowbit(x);
    }
    return ret;
}
void add(int x)
{
    while(x <= n){
        c[x] ++;
        x += lowbit(x);
    }
}
bool cmp(node a, node b){
   return a.w < b.w;
}
int main()
{
   // freopen("b.in", "r", stdin);
    while(~scanf("%d", &n) && n){
        for(int i = 1; i <= n ;i++)  c[i] = 0;
        for(int i = 1; i <= n; i++)  {
            scanf("%d", &b[i].w);
            b[i].num = i;
        }
        sort(b + 1, b + n + 1, cmp);
        int cnt = 0;
        a[b[1].num] = ++cnt;
        for(int i = 2; i <= n; i++)  {
            if (b[i].w == b[i - 1].w)
                a[b[i].num] = cnt;
            else
                a[b[i].num] = ++cnt;
        }
        ll ans = 0LL;
        for(int i = n; i >= 1; i--){
            int x = a[i];
            ans += sum(x - 1);
            add(x);
        }
        cout << ans << endl;
    }
    return 0;
}

HDU 5862 Counting Intersections

[Problem]

给定n(1e5)个平行于坐标轴的线段，询问交点个数。

[Solution]

扫描线+树状数组

我们按照x方向枚举平行与y轴的线段，我们需要得到的是在此时[y1, y2] 有多少个横向线段相交，因此我们对于每个横向线段，只保留两个点，在lx 的时候对应y 的位置数量+1， 在rx+1的位置上对应y数量-1，这样用树状数组维护y值对应的短点数量即可

[Code]

#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
#include<set>
#include<map>
using namespace std;
typedef long long ll;
#define N 402000
int n, m, k, cnt, T, num, delta;
int a[N], c[N];
struct node{
  int num, y, y1, y2, x;
} point[N];
int lowbit(int x) {  return x & -x; }
int bitsum(int x)
{
    int ret = 0;
    while(x > 0){
        ret += c[x];
        x -= lowbit(x);
    }
    return ret;
}
void add(int x)
{
    while(x <= cnt){
        c[x] += delta;
        x += lowbit(x);
    }
}
bool cmp(node a, node b)
{
    if (a.x == b.x)  return a.num > b.num;
    return a.x < b.x;
}
int main()
{
  //  freopen("b.in", "r", stdin);
    scanf("%d", &T);
    while(T--){
        scanf("%d", &n);
        int x1, x2, y1, y2;
        num = 0;
        set<int> dict;
        map<int, int> id;
        for(int i = 1; i <= n ;i++){
            scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
            if (!dict.count(y1))  dict.insert(y1);
            if (!dict.count(y2))  dict.insert(y2);
            if (x1 == x2){
                num++;
                point[num].num = 1;
                point[num].y1 = min(y1, y2);
                point[num].y2 = max(y1, y2);
                point[num].x = x1;
            }
            else  {
                num++;
                if (x1 > x2 ) { int t = x1; x1 = x2; x2 = t; }
                point[num].num = 3;
                point[num].y = y1;
                point[num].x = x1;
                num++;
                point[num].num = 2;
                point[num].y = y1;
                point[num].x = x2 + 1;
            }
        }
        ll ans = 0; cnt = 0;
        sort(point + 1, point + num + 1, cmp);
        for(set<int>::iterator it = dict.begin(); it != dict.end(); ++it)  id[*it] = ++cnt;
        for(int i = 1; i <= cnt; i++) c[i] = 0;
        for(int i = 1; i <= num; i++){
            if (point[i].num == 3){
                int d = id[point[i].y];
                delta = 1;
                add(d);
            }
            else if (point[i].num == 2){
                int d = id[point[i].y];
                delta = -1;
                add(d);
            }
            else{
                int y1 = point[i].y1, y2 = point[i].y2;
                y1 = id[y1];  y2 = id[y2];
                ll tot = bitsum(y2) - bitsum(y1 - 1);
                ans += 0LL + bitsum(y2) - bitsum(y1 - 1);
            }
        }
        cout << ans << endl;
    }
    return 0;
}

ZOJ3672:Gao The Sequence

[Problem]

给定n（2e5）个数，询问多少个子区间[l. r] 的算数平均值大于k

[Solution]

∑ai >= (r - l + 1) * k
∑(ai - k) >= 0
令b[i] = a[i] - k
令sum[i] = ∑b[i]
原式为sum[r] - sum[l - 1] >= 0  ->   sum[r] >= sum[l - 1]

这样我们从小到大枚举r, 每次取出小于sum[r]的数量，然后把sum[r] 加入到bit 中即可

令sum[i] =∑a[i]
sum[r] - sum[l - 1] >= (r - l + 1) * k
sum[r] - r * k >= sum[l - 1] - (l - 1) * k
b[i] = sum[i] - i * k;
b[r] >= b[l] (0 <= l < r)

[Code]

#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
#include<set>
#include<map>
using namespace std;
typedef long long ll;
#define N 202500
int n, m, k, cnt;
ll c[N];
int a[N];
ll sum[N];
int lowbit(int x) {  return x & -x; }
ll bitsum(int x)
{
    ll ret = 0;
    while(x > 0){
        ret += c[x];
        x -= lowbit(x);
    }
    return ret;
}
void add(int x)
{
    while(x <= cnt){
        c[x] ++;
        x += lowbit(x);
    }
}

int main()
{
  //  freopen("b.in", "r", stdin);
    scanf("%d%d", &n, &k);
    for(int i = 1; i <= n; i++)  scanf("%d", a + i);
    sum[0] = 0LL;
    set<ll> dict;
    dict.insert(0);
    for(int i = 1; i <= n; i++) {
        a[i] = a[i] - k;
        sum[i] = sum[i - 1] + a[i];
        if (!dict.count(sum[i])) dict.insert(sum[i]);
    }
    map<ll , int > id;
    cnt = 0;
    ll ans = 0LL;
    for(set<ll>::iterator it= dict.begin(); it != dict.end(); ++it)   id[*it] = ++cnt;
    add(id[0]);
    for(int i = 1; i <= n; i++){
        int x = id[sum[i]];
        ans += bitsum(x);
        add(x);
    }
    cout << ans;
    return 0;
}
YuHsin
关注
专栏目录
2017暑期集训Day 14 树形dp
Yancy_的博客
07-29 429
POJ 2486 Apple Tree POJ 3345 Bribing FIPA POJ 2378 Tree Cutting POJ 3107 Godfather POJ 3140 Contestants Division HDU 5326 Work CodeForces 696B Puzzles LightOJ 1382 The Queue
2017暑期集训 Day 3 搜索与并查集
Yancy_的博客
07-14 376
A - 食物链[solution] 并查集新操作！ 维护决策的矛盾信息，每次决策之间的选择是互相矛盾的 每个动物可能为三种情况：A、B、C，i代表第i种动物是A类动物，i+n代表第i种动物是B类，i+2*n是c类，这样 对于d=1时，x、y为同类动物，即如果x是a类，则y也为a类，即unite(x,y),同理，都为b类，unite(x+n,y +n),同理，unite(x+2*n,y+2*n
ZROI2019暑期集训 Day4
dfn8726的博客
08-08 128
由于我太菜，这篇文章咕了…… 转载于:https://www.cnblogs.com/p-b-p-b/p/11324438.html
2017暑假七林集训day14
整型的映山红
07-03 180
早上 从今天早上开始，为期一个星期的数据结构大作战就开始了。
2017暑假七林集训day15
整型的映山红
07-04 230
早上 因为zcy学长住的远，所以只要是他讲课考试都是在上午QAQ 拿到题，t1又读不懂题。发现t3昨晚刚好做过，但是估计会调试很久，为了避免出现以前把赌注压到一道题上解过翻车的现象，先把t2的暴力写了在说。t3昨晚 yy了一个单调栈的做法，先处理跨区间的点对，然后再递归搞 区间内部的点对。这个过程中可以用两个单调栈来O(1)做。t3一直调试到了11：00，搞了一会t2搞不出来，时间就到了。
2017暑假七林集训day11
整型的映山红
06-30 262
上午 今天因为一些特殊原因，考试的时间被调整到了上午。拿到t1，和昨天一样，并没有什么很好的思路，于是去看t2。t2是一个丢蛋的题，然后发现自己连样例都玩不动，跳过看t3。t3有一个挺容易拿的30分做法，果断写了，先避免爆零。倒回去看t2，准备拿n >= 100的部分分（n >= 100相当于无限蛋），结果测评的时候发现手写log蜜汁挂掉。又返回去看t3，写了个（可能能过一些点的）优化，结果交的
USST暑期集训day1
最新发布
2302_80108031的博客
07-01 388
找到2^j接近len max=max(st[i][j],st[][j]);2.向上枚举(二分，找到深度相同if d[x]>=d[u],v=x)3.（无根树）树的重心：mx最小的节点,mx=儿子的子树中最大的节点数。5.堆(优先队列)：二叉树，只能去堆顶拿，堆顶是最值。从上到下搜索，遇到包含，return 值，不包含，return 0。st[i][j]，i向后，pow(2,j)位[i,i+2^j];pre[i][j]：i向上跳2^j步到达的节点。--------以上都是线性--------
【ACM】—蓝桥杯大一暑期集训Day4
H1727548的博客
07-18 1767
ACM大一暑期集训Day4练习题
2021暑假集训 Day 2
SSL_hhw的博客
08-10 168
回顾：就不知道今天为什么那么想尝试富文本编辑器了…… A 今天的第一题非常的无语呵，想了90min？而且差一个longlong就A了，别吧，心态炸裂。 Decription: n个小球排成一列，第i个小球颜色为ai， 定义小球i与j的距离为|(i-j)*(ai-aj)|， 求所有小球间距离之和， i与j 和 j与i 不必重复计算。 Input: 第一行一个n 第二行给出a1,a2,...,an-1,an Output 按照题目要求作答 ...
【2020暑期海亮集训】Day 5：普及组模拟赛（7.21）
Admin-zx 的博客
07-21 481
目录引子题目T1-时间（time）【题目描述】【输入数据】【输出数据】【样例输入 1】【样例输出 1】【样例输入 2】【样例输出 2】【样例输入 3】【样例输出 3】【数据范围】【附加知识】【解题思路】【Code】T3-数学（math）【题目描述】【输入数据】【输出数据】【样例输入】【样例输出】【数据范围】【解题思路】【Code】T4-发际线（hair）【题目描述】【输入数据】【输出数据】【样例输入】【样例输出】【数据范围】【解题思路】【Code】 引子 怎么还有10天结束这恐怖的集训呀！整天只有考试，考试
17天正睿集训总结
qq_40681184的博客
07-28 772
【写在前面】   "正睿普及=IOI"是常识    这里开始亲身17天体验     ——来自蒟蒻的微笑(*￣︶￣) 持续更新中 - 已完成 1/17 7.28-Day0 集训拉开序幕。 首先签字（转：" 宁可做大佬也不可做CXK "，一群人纷纷在大佬区签了字）+ 拿CXK专属纪念T恤+开始点外卖。 8：00 开始为时半个小时的开会，老师洪亮的声音让我们认清现实： 准备好爆 0 吧。 老师继续...
暑假2019培训：Day2提高组测试赛
星暗宇的博客
07-07 333
概述 今天我们的Chty_syq放了点水，第一题来了个水题…… 但我依旧考的不好 题目顺序 1． 小刚传说 2． 中珂院的难题 3． 欧拉图 1． 小刚传说 (legendary.cpp) 【 时空限制】 时间限制： 2s 空间限制： 256MB 【 问题描述】 众所周知， 刘小刚是海亮中学的金牌教练， 然而世人并不知道他的传奇经历， 以及那个 曾经轰动世界的名字 sharpland。 20...
正睿17+2天集训总结TOT
qq_40681184的博客
08-16 326
Part1 Feeling 暑期开始踏入C班集训的时候，觉得C班的考试和讲课好难啊。但一对比正睿的B\C班，顿时无话可说。在海亮考试勉勉强强不会掉出前面的一半，到了正睿就是勉勉强强能不排上倒数（虽然最后一次就是倒数了）。 B班的速度太快，还是C班注重一些基本的模板和概念，更加具体。 题目质量很高，PPT上的题很少见但是很经典，学到的也有不少吧。 听得都差不多，像二分、贪心、线段树、树链剖分一些基础...
2021 7.5~7.20日集训总结
zero_orez6的博客
07-29 194
Day.1 集训开始的第一天,以数据结构为主体，在之前的基础上进行了不小的扩展。 并查集 启发式合并：一种根据子树大小为基准的合并集合的方式，能够使得“合并集合”操作和“询问”的复杂度降到O(logn)O(log n)O(logn) 路径压缩：在只考虑每个点与其根节点的情况下，我们可以将任意一个子树的父亲设置为根节点，使得查询操作降到O(1)O(1)O(1) ...
BZOJ 刷题总结
Yancy_的博客
01-29 1478
背景 辣鸡的人总要想法自救，便产生了寒假学些新算法、在Bzoj刷些题的想法。一来为明年省赛做准备…压力不小；二来寒假也可以有些事情做。 PS 1.不定时更新做题的思路和吐槽 2.按照hzw刷题顺序训练，具体依照BZOJ题表 3.希望寒假能够至少刷够100道题目吧，在此立个FLAG，希望别被青岛的妖风吹跑 Problem &amp; Solution BZOJ 1003 ...
SDU省赛选拔-ACM ICPC 2010–2011, NEERC, Northern Subregional Contest
Yancy_的博客
04-14 1252
A. Alien Communication Masterclass【题目】给定序列A、B，构造一个等式，使得对于任意的A[i]进制都满足，对于任意的B[i]都 不满足【分析i】同一个式子，不同进制的结果如果不相同，一般是在进位上，首先我们可以想到如果对于A序列的进制都满足呢，可以写成x1 * x2 * x3 *x4 = 0这种，每种进制可以构造一个等于0的式子，对于A[i]，可以构造（10 -
sduacm16级寒假训练 自测
Yancy_的博客
02-16 584
POJ 3579 Median【sol】显然，每个点与其他的点形成的差是有序的，这样相当于有n组升序的序列，询问中位数。二分中位数，然后判断比他小的数的个数，复杂度n* logn*logn#include<cstdio> #include<iostream> #include<cstring> #include<cmath> #include<algorithm> using namespace
Educational Codeforces Round 37 所有题目整理！
Yancy_的博客
02-05 548
A.Water The Garden int ct = a / p; if (a % p != 0) ++ct; 等价于 ct = (a + p - 1) / p; B. Tea Queue solution：按照题意模拟即可 C. Swap Adjacent Elements 题意：给定1-n的一个排列，存在一些i，其中1 思路：首先考虑数字1，假设其位置为pos1，必然
北京大学郭炜教授详解线段树与树状数组在ACM/ICPC竞赛中的应用
线段树和树状数组是数据结构中的两种重要工具，用于高效处理区间查询和更新问题，在ACM/ICPC等计算机编程竞赛中广泛应用。北京大学信息学院郭炜教授在暑期课程《ACM/ICPC竞赛训练》中详细讲解了这两种数据结构。 ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值