HDU4973 A simple simulation problem

最新推荐文章于 2020-11-20 16:07:49 发布
最新推荐文章于 2020-11-20 16:07:49 发布
阅读量375 收藏
点赞数
分类专栏： 2017暑假 线段树 线段树分治
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xzxxzx401/article/details/75668476
版权

HDU4973 A simple simulation problem

线段树分治

传送门：HustOJ

题意

设计数据结构，支持操作：
初始时给定n，表示序列123…n
D l r：将序列下标lr之间的数复制一份，比如12345 D 2 4，变为12233445，在D 1 4变为112222333445。
Q l r：查询下标lr之间的数中，出现次数最多的数的出现次数。

思路

是个线段树分治的题。

线段树开到n，每个叶子节点表示一个数，节点内存这个数出现了多少次。其余节点维护一个sum，表示他管控区间内共有多少个数，以及num，表示出现次数最多的数出现了多少次。

注意到区间可能不是整区间，比如122333445，我操作3,7。表示操作23334，2和4就不是一个完整的节点。这就是线段树分治的思路，将这个操作分开，分成前，中，后三部分。前后分别是一个单点操作，中是一个区间操作，可以lazy。

但是这又一个问题：如何将它给的实际序列下标信息转化成线段树的区间信息并分开前后？这需要写一个find函数：

int zhao(LL pos, LL &num, LL &tim, int l, int r, int rt)
/*
返回值是实际序列中pos位对应数
num是查询的pos位数字是实际序列中的第几个同样的数字
tim是该数字在实际序列中出现的总次数
*/
{
    if(l==r)
    {
        num=pos;
        tim=stree[rt].num;
        return l;
    }
    pushdown(rt);
    int m=(l+r)>>1;
    if(stree[rt<<1].sum>=pos) return zhao(pos, num, tim, lson);
    else return zhao(pos-stree[rt<<1].sum, num, tim, rson);
}

比如对于序列：12233334444567

我find10，返回值是4，num是3，tim是4。num作用在于在于单点加倍时，该点有一部分数字落在带加倍区间内，另一部分不落在。对于区间左端点，我肯定加倍tim-num+1个数字，右端点加倍num个数字。

返回值得作用在于扣去左右端点对应的数字。这样在区间加倍和区间查询时，直接查询[左侧返回值+1，右侧返回值-1]内。

恩。。。那个find那块比较不好理解，而且变量特别多，不要搞混。

另外。。这种多组数据的。。一定注意线段树要清干净。我刚开始疯狂wa，因为lazy没清干净，下一组数据查询时直接pushdown了>_<

另外据正雄大佬说，这题不用分治，完全单点加一个小优化就能过。不过这样有坑。

void db(LL L, LL R, int l, int r, int rt)
{
    if(l==r)
    {
        stree[rt].sum+=(R-L+1);
        stree[rt].num=stree[rt].sum;
        return;
    }
    int m=(l+r)>>1;
    if(R<=stree[rt<<1].sum) db(L, R, lson);
    else if(L>stree[rt<<1].sum) db(L-stree[rt<<1].sum, R-stree[rt<<1].sum, rson);
    else
    {
        LL cnm=stree[rt<<1].sum;
        db(L, stree[rt<<1].sum, lson);
        db(1, R-cnm, rson);

        //坑！！！
        //db(L, stree[rt<<1].sum, lson);
        //db(1, R-stree[rt<<1].sum+1, rson);
    }
    pushup(rt);
}

如果按注释里面那么写，在第一个db完成后，左儿子的sum值就变了。R-他就不对了。如果没注意到这里，样例会输出2 2。

代码

线段树分治：注意线段树清干净！！，直接在外面置零就行

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <iomanip>
#define _ ios_base::sync_with_stdio(0),cin.tie(0)
#define M(a,b) memset(a,b,sizeof(a))
#define N n
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1

using namespace std;

const int MAXN=100007;
const int oo=0x3f3f3f3f;
typedef long long LL;
const LL loo=4223372036854775807ll;
typedef long double LB;
const LL mod=1e9+7;
const double eps=1e-6;
struct SegTree
{
    LL sum; LL num;
    int lazy;
}stree[4*MAXN];
void pushup(int rt)
{
    stree[rt].sum=stree[rt<<1].sum+stree[rt<<1|1].sum;
    stree[rt].num=max(stree[rt<<1].num, stree[rt<<1|1].num);
}
void pushdown(int rt)
{
    if(stree[rt].lazy)
    {
        stree[rt<<1].lazy+=stree[rt].lazy;
        stree[rt<<1|1].lazy+=stree[rt].lazy;
        for(int i=0;i<stree[rt].lazy;i++)
        {
            stree[rt<<1].sum*=2;
            stree[rt<<1].num*=2;
            stree[rt<<1|1].sum*=2;
            stree[rt<<1|1].num*=2;
        }
        stree[rt].lazy=0;
    }

}
void build(int l, int r, int rt)
{
    stree[rt].lazy=stree[rt].num=stree[rt].sum=0;
    //注意区间外清空！！！否则非叶子节点就有了lazy
    if(l==r)
    {
        stree[rt].num=1;

        stree[rt].sum=1;
        stree[rt].lazy=0;
        return;
    }
    int m=(l+r)>>1;
    build(lson);
    build(rson);
    pushup(rt);
}
void update1(int pos, LL num, int l, int r, int rt)//单点更新
{
    if(l==r)
    {
        stree[rt].sum+=num;
        stree[rt].num+=num;
        return;
    }
    pushdown(rt);
    int m=(l+r)>>1;
    if(pos<=m) update1(pos, num, lson);
    else update1(pos, num, rson);
    pushup(rt);
}
void update2(LL L, LL R, int l, int r, int rt)//区间更新
{
    if(L<=l&&r<=R)
    {
        stree[rt].lazy++;
        stree[rt].num*=2;
        stree[rt].sum*=2;
        return;
    }
    pushdown(rt);
    int m=(l+r)>>1;
    if(L<=m) update2(L, R, lson);
    if(R> m) update2(L, R, rson);
    pushup(rt);
}
LL query(int L, int R, int l, int r, int rt)
{
    if(L<=l&&r<=R)
    {
        return stree[rt].num;
    }
    pushdown(rt);
    int m=(l+r)>>1;
    LL ma=0;
    if(L<=m)  ma=max(ma, query(L, R, lson));
    if(R> m)  ma=max(ma, query(L, R, rson));
    return ma;
}
int zhao(LL pos, LL &num, LL &tim, int l, int r, int rt)
/*
返回值是实际序列中pos位对应数
tim是该数字在实际序列中出现的总次数
num是查询的pos位数字是实际序列中的第几个同样的数字
*/
{
    if(l==r)
    {
        num=pos;
        tim=stree[rt].num;
        return l;
    }
    pushdown(rt);
    int m=(l+r)>>1;
    if(stree[rt<<1].sum>=pos) return zhao(pos, num, tim, lson);
    else return zhao(pos-stree[rt<<1].sum, num, tim, rson);
}
int main()
{
    int T;scanf("%d", &T);
    for(int t=0;t<T;t++)
    {
        printf("Case #%d:\n", t+1);
        int n, q;scanf("%d%d", &n, &q);
        build(1, n, 1);
        for(int i=0;i<q;i++)
        {
            char c[2];LL l, r;
            scanf("%s%lld%lld", c, &l, &r);
            LL ln, rn;int lp, rp;LL tmp;
            rp=zhao(r, rn, tmp, 1, n, 1);
            lp=zhao(l, ln, tmp, 1, n, 1);
            if(c[0]=='D')
            {
                if(lp==rp) update1(lp, rn-ln+1, 1, n, 1);
                else
                {
                    update1(rp, rn, 1, n, 1);
                    update1(lp, tmp-ln+1, 1, n, 1);
                    if(1<rp-lp)
                        update2(lp+1, rp-1, 1, n, 1);
                }
            }
            else
            {
                if(lp==rp)
                {
                    printf("%lld\n", rn-ln+1);
                }
                else
                {
                    LL ans=0;
                    if(rn>ans) ans=rn;
                    if(tmp-ln+1>ans) ans=tmp-ln+1;
                    if(rp-lp>1)
                    {
                        ans=max(ans, query(lp+1, rp-1, 1, n, 1));
                    }
                    printf("%lld\n", ans);
                }
            }
        }
    }

    return 0;
}

单点更新+蜜汁优化：注意两次更新之间的耦合关系！！

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <iomanip>
#define _ ios_base::sync_with_stdio(0),cin.tie(0)
#define M(a,b) memset(a,b,sizeof(a))
#define N n
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1

using namespace std;

const int MAXN=100007;
const int oo=0x3f3f3f3f;
typedef long long LL;
const LL loo=4223372036854775807ll;
typedef long double LB;
const LL mod=1e9+7;
const double eps=1e-6;
struct SegTree
{
    LL sum; LL num;
}stree[4*MAXN];
void pushup(int rt)
{
    stree[rt].sum=stree[rt<<1].sum+stree[rt<<1|1].sum;
    stree[rt].num=max(stree[rt<<1].num, stree[rt<<1|1].num);
}
void build(int l, int r, int rt)
{
    if(l==r)
    {
        stree[rt].num=1;
        stree[rt].sum=1;
        return;
    }
    int m=(l+r)>>1;
    build(lson);
    build(rson);
    pushup(rt);
}
void db(LL L, LL R, int l, int r, int rt)
{
    if(l==r)
    {
        stree[rt].sum+=(R-L+1);
        stree[rt].num=stree[rt].sum;
        return;
    }
    int m=(l+r)>>1;
    if(R<=stree[rt<<1].sum) db(L, R, lson);
    else if(L>stree[rt<<1].sum) db(L-stree[rt<<1].sum, R-stree[rt<<1].sum, rson);
    else
    {
        LL cnm=stree[rt<<1].sum;
        db(L, stree[rt<<1].sum, lson);
        db(1, R-cnm, rson);
    }
    pushup(rt);
}
LL query(LL L, LL R, int l, int r, int rt)
{
    if(l==r)
    {
        return R-L+1;
    }
    if(L==1&&R==stree[rt].sum) return stree[rt].num;

    int m=(l+r)>>1;
    LL res=-1;
    if(R<=stree[rt<<1].sum) return query(L, R, lson);
    else if(L>stree[rt<<1].sum) return query(L-stree[rt<<1].sum, R-stree[rt<<1].sum, rson);
    else
    {
        return max(query(L, stree[rt<<1].sum, lson),
            query(1, R-stree[rt<<1].sum, rson));
    }
}
int main()
{
    int T;scanf("%d", &T);
    for(int t=0;t<T;t++)
    {
        printf("Case #%d:\n", t+1);
        int n, q;scanf("%d%d", &n, &q);
        build(1, n, 1);
        for(int i=0;i<q;i++)
        {
            char c;LL l, r;
            scanf(" %c %lld %lld", &c, &l, &r);
            if(c=='D')
            {
                db(l, r, 1, n, 1);
            }
            else
            {
                printf("%lld\n", query(l, r, 1, n, 1));
            }
        }
    }

    return 0;
}
xzxxzx401
关注
专栏目录
HDU 5974 A Simple Math Problem 2016大连区域赛
LuRiCheng的博客
11-06 904
A Simple Math Problem Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 14 Accepted Submission(s): 12 Problem Description Given two positive in
HDU1757 - A Simple Math Problem - 矩阵快速幂
终有绿洲摇曳在沙漠
03-05 391
1.题目描述 A Simple Math Problem Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 4507    Accepted Submission(s): 2716 Problem Descripti
A simple problem
七刀的专栏
07-12 798
A simple problem 题目描述 Dreamone has a lovely cat. The cat can comfort you if you are very dejected, she can also play with you if  you are very bored. So Dreamone loves her very much. Of cours
hdu4973
xcszbdnl的专栏
08-21 624
线段树，每个点记录该点出现的次数，对于l,r的更新， 找到l所对应的shu
HDU4973 【几何。】
weixin_30469895的博客
10-27 88
题意： 给你一个以原点为圆心的两个圆，一个大圆，一个小圆，然后给你一个硬币和他的速度，问你经过大圆的时间； 思路： 直接杠。。 然后wa的怀疑人生，后面wa在了速度的方向，如果我说一个点在两个圆的左上方，那么速度也是左上方，那么就是0； 其他都还简单； 贴一发比赛时的挫代码。。。 #include <iostream> #include <stdio.h&g...
H - A simple problem（挖）
Lzzz
11-20 107
#include <stdio.h> int a[10000]; int b[10000]; int main() { int t; scanf("%d", &t); while (t--) { int n; scanf("%d", &n); if (n < 0) { printf("-"); n = n * (-1); } if (n == 1) printf("1\n"); else { int flag = 0.
H - A simple problem （循环小数的判断）
weixin_30257433的博客
12-03 226
H - A simple problem HDU - 2522 Zty很痴迷数学问题.。一天，yifenfei出了个数学题想难倒他，让他回答1 / n。但Zty却回答不了^_^. 请大家编程帮助他. Input第一行整数T，表示测试组数。后面T行，每行一个整数 n (1<=|n|<=10^5).Output输出1/n. ...
hdu-api:A simple SDK for HDU. 一个提供一卡通服务、考试、课表、选课和公共信息等 API 的 SDK
05-04
A simple SDK for HDU. hdu-api 是一个集结 HDU 所有教务管理服务的 SDK，提供了一卡通服务、考试、课表、选课和一些公共信息如空闲教室、上课时间等信息的 API。 hdu-api 主要基于 Requests 库和 Beautiful Soup 库...
hdu1757 - A Simple Math Problem（矩阵快速幂模板题）
red_red_red的博客
05-14 994
Lele now is thinking about a simple function f(x). If x < 10 f(x) = x. If x >= 10 f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + …… + a9 * f(x-10); And ai(0<=i<=9) can only be 0 or 1 ...
A Simple Problem
just_gsj的博客
09-28 485
题面描述： 求 0~n 中有多少个整数 k 满足(2 * k) XOR (3 * k) == k, 其中 XOR 为异或操作。 输入数据： 输入数据的第一行为一个正整数 T(T&amp;lt;=30)，表示测试数据的组数。 接下来的 T 行中，每行为一组测试数据，包含一个正整数 n(n&amp;lt;=1000000000)。 输出数据： 对每一组输入数据，输出一行结果“Case #id: M”，表示第 id 组...
hdu4973 A simple simulation problem.（成段更新+下标转化）
Iguodala的博客
05-19 448
http://acm.hdu.edu.cn/showproblem.php?pid=4973 题意：给你一串连续升序数字和两种操作，D表示将[l, r]区间内的元素个数加倍，Q表示查询[l,r]区间内的最大元素个数。 ps:看的别人的博客。。草。。 思路：既然元素再怎么加倍也总是连续的，那么每个不同元素当做线段树子节点，开一个sum域表示同种元素的个数。于是难点就转化为了
hdu4973(线段树)
u013509299的专栏
08-22 777
题意：给1-n的数字，m次操作，每次操作将区间[l,r]中的数字翻倍，比如1,2,3,4,5,D 2,3,得到的结果是1,2,2,3,3,4,5， 求区间[l,r]中同一个数字出现的最多次数 解题思路:自己AC的一道线段树，感觉很不错。以[1,n]建立线段树，线段树中维护三个信息，max即当前区间内该区间内同一个数字出现的最多次数，sum即当前区间内数字的总个数，flag延迟标记。对于每次操
Just a Simple Problem
无限迭代中......
12-09 249
http://whu.contest.codeforces.com/group/LqmZzHyOQN/contest/234192/problem/E 题解：虽然不知道题目在干嘛 但是跑了几分钟的暴力，似乎是素数，就AC了 /* *@Author: STZG *@Language: C++ */ #include &lt;bits/stdc++.h&gt; #include&lt;io...
hdu 4973 A simple simulation problem.（2014 Multi-University Training Contest 10）
ACway的专栏
08-23 6141
A simple dynamic programming problem Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 455    Accepted Submission(s): 187 Problem Descri
hdu4973 A simple simulation problem.（多校第十场1003）
huantwofat
08-22 482
可以直接用树状数组求和加二分判断l
2018 ACM 国际大学生程序设计竞赛上海大都会赛 H [A Simple Problem with Integers]
ACM2017
08-05 487
题目链接：https://www.nowcoder.com/acm/contest/163/H 四糸智乃的题解: https://www.nowcoder.com/discuss/90721?type=101&amp;order=0&amp;pos=1&amp;page=0 来源：牛客网 时间限制：C/C++ 2秒，其他语言4秒 空间限制：C/C++ 262144K，其他语言524288...
2014多校10（1003）hdu4973（简单线段树区间操作）
cq_phqg
09-04 704
A simple simulation problem. Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 679    Accepted Submission(s): 269 Problem Description
hdu4973 线段树
江上白衣凌云
08-23 476
题意：给定的序列 有两个操作 2将其中一段加倍 3询问
A - LOOPS HDU - 3853
最新发布
07-28
所以，问题A - LOOPS HDU - 3853是关于解决迷宫问题的一个题目。 #### 引用[.reference_title] - *1* *3* [期望DP——C - LOOPS HDU - 3853](https://blog.csdn.net/GreyBtfly/article/details/80534173)[target="_...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值