题目





像素有点低啊~

算了凑合一下就好啦~

题目大意

给你一个首尾相接的数列,每次对一个区间进行操作:

顺时针操作,如果当前值比vvv大,就交换。输出最后的vvv。


比赛思路

首先这题的时限这么仁慈,一定有天大玄机。

并且一看这题,就感觉像是一个数据结构。

首先在想链表,但显然这题链表不可做。

然后一直在想带修主席树。

没想出来……

最后弃疗,直接打了个暴力。

WTF?15分?说好的25分呢?

然而实际上这题很没良心地捆绑数据,将10分和另外15分绑在一起了。

出题人,你怎么能这样子啊?你忍心吗?


正解

这题WMY大佬说可以用带修主席树做。

刚了一个下午,最终,他弃疗了……

原因是标记不好下传。

实际上正解是分块。

首先看到时间复杂度,我们就应该想到这题可以随意给你搞事情。

然而我就是没有想到分块!!!

首先,对于一个区间,如果有一个操作经过了这个区间,设区间中的最大值为mxmxmx。若mx>vmx>vmx>v,则交换,否则继续。

这个结论是很显然的,依靠这个结论可以再拿15分。

我们可以将其分块,每个块的大小为KKK。对于每个块,我们维护一个大根堆,存下这个块里面的所有值。

如果处理整块,就直接和最大值比较,然后像之前一样操作。并且,在这个块上打一个标记。

如果处理散块,就要将这个散块还原,然后暴力搞一遍。

如何还原呢?

首先,对于每个块,我们将标记存在一个小根堆里面。

在还原的时候,我们从左到右扫,对于每个值,用小根堆的堆顶操作。如果ai>va_i>vai​>v,就交换(也就是将vvv弹出,将aia_iai​加入,并且改变aia_iai​的值)

最后将标记清空。

这就还原了整个块了,然后暴力搞一遍,重构大根堆。

那么问题来了,为什么每次用小根堆的堆顶操作?

可以感性地理解一下:

对于第一个,这些标记的操作都会对它有操作。如果当前的这个值大于vvv,那么就要被交换。而交换那么多遍,最后真正能对它做出影响的是最小的vvv,其它的东西都会传到后面去。

然后对于后面的,也是一样的道理。

和氧化还原反应好像啊!——ZJQ

所以整块的时间复杂度是O(qnKlg⁡K)O(q\frac{n}{K}\lg K)O(qKn​lgK),散块的时间复杂度是O(qKlg⁡q)O(qK\lg q)O(qKlgq)

然后平衡规划一下,得出KKK大概为n\sqrt nn​。

然而分块的常数是有差异的,所以KKK的取值据实际而定。

我取了800800800。当我取600600600时,程序就崩了,或许是堆太多了吧。(我用了STL的堆)


代码

  1. using namespace std;
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #include <queue>
  6. #define N 400000
  7. #define K 800
  8. int n,q;
  9. int a[N];
  10. int m;
  11. #define bel(x) ((x)/K)
  12. #define nxt(x) (((x)+1==m)?(0):((x)+1))
  13. priority_queue<int> h[N/K];
  14. priority_queue<int,vector<int>,greater<int> > bz[N/K];
  15. inline void pushdown(int b,int l,int r,int &v){//处理散块
  16. //还原
  17. if (!bz[b].empty()){
  18. for (int i=b*K;i<b*K+K && i<n;++i){
  19. int t=bz[b].top();
  20. if (a[i]>t){
  21. bz[b].pop();
  22. bz[b].push(a[i]);
  23. a[i]=t;
  24. }
  25. }
  26. while (!bz[b].empty())
  27. bz[b].pop();
  28. }
  29. //暴力处理
  30. for (int i=l;i<=r;++i)
  31. if (a[i]>v)
  32. swap(a[i],v);
  33. //重构
  34. while (!h[b].empty())
  35. h[b].pop();
  36. for (int i=b*K;i<b*K+K && i<n;++i)
  37. h[b].push(a[i]);
  38. }
  39. inline void getinto(int b,int &v){//表示处理整块,v进入b中,再出来
  40. int t=h[b].top();
  41. if (t>v){
  42. h[b].pop();
  43. h[b].push(v);
  44. bz[b].push(v);
  45. v=t;
  46. }
  47. }
  48. int main(){
  49. scanf("%d%d",&n,&q);
  50. for (int i=0;i<n;++i)
  51. scanf("%d",&a[i]);
  52. m=(n-1)/K+1;
  53. for (int i=0;i<m;++i)
  54. for (int j=0;j<K && i*K+j<n;++j)
  55. h[i].push(a[i*K+j]);
  56. while (q--){
  57. int l,r,v;
  58. scanf("%d%d%d",&l,&r,&v);
  59. l--,r--;
  60. int bl=bel(l),br=bel(r);
  61. if (bl==br){
  62. if (l<=r)
  63. pushdown(bl,l,r,v);
  64. else{
  65. pushdown(bl,l,min(bl*K+K-1,n-1),v);
  66. for (int i=nxt(bl);i!=br;i=nxt(i))
  67. getinto(i,v);
  68. pushdown(br,br*K,r,v);
  69. }
  70. }
  71. else{
  72. if (l==bl*K)
  73. getinto(bl,v);
  74. else
  75. pushdown(bl,l,min(bl*K+K-1,n-1),v);
  76. for (int i=nxt(bl);i!=br;i=nxt(i))
  77. getinto(i,v);
  78. if (r==min(br*K+K-1,n-1))
  79. getinto(br,v);
  80. else
  81. pushdown(br,br*K,r,v);
  82. }
  83. printf("%d\n",v);
  84. }
  85. return 0;
  86. }

我才发现原来要打个cpp才能有颜色,我之前打的都是C++,天哪,博客更新之后就是不一样!


总结

看见时限大的题,往分块方面想一想,或许就能很好解决了。

分块的优点,在于它只需要对整块和散块分块处理,也就是说,不像线段树那样下传时这么复杂。

还有,这题有没有其他的方法。比如,分块套分块(手动滑稽)

JZOJ5967 常数国的更多相关文章

  1. Before NOIP 2018

    目录 总结 刷题 2018 - 9 - 24 2018 - 9 - 25 2018 - 9 - 26 2018 - 9 - 27 2018 - 9 - 28 2018 - 9 - 29 2018 - ...

  2. [jzoj NOIP2018模拟10.23]

    丢分主要是下面几个方面: 1.T2代码交错了,有个特判没写丢了10分 2.T1线段树加等差数列写错了(其实二维差分就可以,但我当时不会) 3.T3思考再三还是为了10分写上了主席树,还是写错了 总体评 ...

  3. 解析大型.NET ERP系统 多国语言实现

    实现多国语言有许多种实现方案,无外乎是一种字符串替换技术,将界面控件的文本标签替换成相应语言的文字..NET Windows Forms实现多国语言的方法有以下几种: 1 .NET的方案,使用资源文件 ...

  4. [LeetCode] Insert Delete GetRandom O(1) - Duplicates allowed 常数时间内插入删除和获得随机数 - 允许重复

    Design a data structure that supports all following operations in average O(1) time. Note: Duplicate ...

  5. [LeetCode] Insert Delete GetRandom O(1) 常数时间内插入删除和获得随机数

    Design a data structure that supports all following operations in average O(1) time. insert(val): In ...

  6. Flex 1046: 找不到类型,或者它不是编译时常数;1180: 调用的方法 CompPropInfo 可能未定义

    导入项目之后一直报这个错误, 1046: 找不到类型,或者它不是编译时常数: 1180: 调用的方法 CompPropInfo 可能未定义 想这应该是没有把当前这个类编译进项目当中,找了半天也没有找到 ...

  7. 高质量,高效率的多国语言软件开发(Web/PC/Mobile),使用接口约束/调用不同语言资源

    偶然间翻出了几年前写的一个小程序,把当时的资料整理整理分享一下. 当时为了给自己的软件实现多国语言功能,而开发的辅助工具:SE String Resource. 这是当时基于自己另一款 IDE 软件抽 ...

  8. ARM的常数表达式

    ARM的常数表达式   如果说Intel指令中的立即数,相信大家都很熟悉.类似的,Arm指令中的“立即数”就是常数表达式.之所以称为常数表达式,而不称为立即数是有原因的. Intel指令属于CISC指 ...

  9. Atitti 存储引擎支持的国内点与特性attilax总结

    Atitti 存储引擎支持的国内点与特性attilax总结 存储引擎处理的事情: · 并发性:某些应用程序比其他应用程序具有很多的颗粒级锁定要求(如行级锁定). · 事务支持:并非所有的应用程序都需要 ...

随机推荐

  1. CSS3——过渡

    过渡(transition)是CSS3中具有颠覆性的特征之一,我们可以在不使用 Flash 动画或 JavaScript 的情况下,当元素从一种样式变换为另一种样式时为元素添加效果. 帧动画:通过一帧 ...

  2. LUOGU P2290 [HNOI2004]树的计数(组合数,prufer序)

    传送门 解题思路 \(prufer\)序,就是所有的不同的无根树,都可以转化为唯一的序列.做法就是每次从度数为\(1\)的点中选出一个字典序最小的,把这个点删掉,并把这个点相连的节点加入序列,直到只剩 ...

  3. H5 回到顶部按钮

    //这里我用的是svg,所以代码老长了,可以把svg部分换成自己的图片什么的$(function(){ var obj = $('<div style="width:1.3rem;po ...

  4. [NOI 2018]冒泡排序

    题意:求所有字典序大于给定序列且满足条件的排列个数之和. 思路: 考虑dp即可,打表出卡特兰数优化,直接dp可以44... #include <bits/stdc++.h> using n ...

  5. Android中visibility属性

    Android开发中,大部分控件都有visibility这个属性,其属性有3个分别为“visible ”.“invisible”.“gone”.主要用来设置控制控件的显示和隐藏. 1) 可见(visi ...

  6. windows 映射samba Linux服务器,输入正确的账号密码却提示“ 指定的网络密码不正确

    重启Linux samba服务也没用,重启Linux和windows系统也没用,急!!! 最佳答案 linux中要添加对应的系统用户和samba用户useradd titiansmbpasswd -a ...

  7. timestamp的自动更新 ON UPDATE CURRENT_TIMESTAMP

    最近有一个关于MySQL版本升级的事,涉及到一些关于时间类型的细节问题需要查明,因此到官网找到相关文章,翻出来比较方便自己理解,博客这里也贴一下. 参考官网网址: https://dev.mysql. ...

  8. Mac 解决硬盘插入不能写的问题

    软件解决 链接:https://pan.baidu.com/s/1H_zvPPpW0dp7aRUvjDnkQA  密码:8fit 有个NTFS的移动硬盘,默认Mac系统是不能写入,只能读取的,我们可以 ...

  9. 自己写一个依赖注入容器Container

    前言:在平时的写代码中为了解耦.方便扩展,经常使用一些DI容器(如:Autofac.Unity),那是特别的好用. 关于它的底层实现代码 大概是这样. 一.使用依赖注入的好处 关于使用依赖注入创建对象 ...

  10. Ubuntu 18.04/18.10快速开启Google BBR的方法

    说明:Ubuntu 18.04改变挺大的,内核直接升到了正式版4.15,而BBR内核要求为4.9,也就是说满足了,所以我们不需要换内核就可以很快的开启BBR,这里简单说下方法. 提示:Ubuntu 1 ...