[BZOJ 1135][POI2009]Lyz

题意

初始时滑冰俱乐部有 \(1\) 到 \(n\) 号的溜冰鞋各 \(k\) 双。已知 \(x\) 号脚的人可以穿 \(x\) 到 \(x+d\) 的溜冰鞋。

有 \(m\) 次操作,每次包含两个数 \(r_i,x_i\) 代表来了 \(x_i\) 个 \(r_i\) 号脚的人。\(x_i\) 为负,则代表走了这么多人。 对于每次操作,输出溜冰鞋是否足够。足够输出 TAK, 否则输出 NIE.

\(n\le 2\times 10^5,m\le5\times 10^5,k\le1\times 10^9,d\in[0,n],1\le r_i\le n-d,|x_i|\le1\times 10^9\)

题解

被某Robbery少许加强后丢到了胡策里...

由霍尔定理, 二分图存在完美匹配当且仅当一侧的任意一个子集中的点连接的另一侧的点的数量都不小于这个子集的大小. 虽然看上去要枚举子集但是实际上我们发现不连续的几段更有可能满足霍尔定理的要求, 所以只要连续区间都满足霍尔定理的要求那么所有子集就都满足了. 感性证明如下:

假设选中了一个不连续的子集, 那么显然可以在不改变另外一侧的邻接情况下在当前子集中添加新的点. 如果不能添加新的点的话可以把子集拆成若干部分, 每个部分是一个连续段. 显然拆开后或者添加新点后更可能会破坏霍尔定理的要求.

也就是说如果设 \(i\) 号脚的人共有 \(s_i\) 个, 那么溜冰鞋不足当且仅当存在任意一个区间 \([l,r]\) 满足下式:

\[\sum_{i=l}^rs_i>k(r-l+d+1)
\]

那么我们拆开移项就可以得到:

\[\begin{aligned}
\sum_{i=l}^rs_i-k(r-l+1)&>kd \\
\sum_{i=1}^r(s_i-k)&>kd
\end{aligned}
\]

于是就变成了一个支持单点加法的动态区间最大子段和问题. 线段树动态DP经典操作.

加强版里不保证 \(r\le n-d\), 需要继续考虑 \(r>n-d\) 的情况. 此时溜冰鞋不足的充要条件相当于:

\[\sum_{i=l}^rs_i>k(n-l+1)
\]

显然当 \(r=n\) 的时候左侧取到最大值, 我们只计算 \(r=n\) 时是否满足条件即可. 此时相当于:

\[\sum_{i=l}^ns_i>k(n-l+1)
\]

设 \(S\) 是 \(\langle s_i\rangle\) 的前缀和, 那么我们可以发现上式等价于:

\[\begin{aligned}
S_n-S_{l-1}&>kn-k(l-1) \\
k(l-1)-S_{l-1}&>kn-S_n
\end{aligned}
\]

左侧的最大值显然也可以用线段树维护出来.

参考代码

  1. #include <bits/stdc++.h>
  2. const int MAXN=1e5+10;
  3. typedef long long intEx;
  4. struct Node{
  5. struct Data{
  6. intEx sum;
  7. intEx lmax;
  8. intEx maxs;
  9. intEx rmax;
  10. Data(){}
  11. Data(intEx val){
  12. this->sum=val;
  13. this->lmax=this->rmax=this->maxs=std::max(this->sum,0ll);
  14. }
  15. Data friend operator+(const Data& a,const Data& b){
  16. Data ans;
  17. ans.sum=a.sum+b.sum;
  18. ans.lmax=std::max(a.lmax,a.sum+b.lmax);
  19. ans.rmax=std::max(a.rmax+b.sum,b.rmax);
  20. ans.maxs=std::max(a.rmax+b.lmax,std::max(a.maxs,b.maxs));
  21. return ans;
  22. }
  23. };
  24. int l;
  25. int r;
  26. Data val;
  27. Node* lch;
  28. Node* rch;
  29. Node(int,int);
  30. void Maintain();
  31. void Add(int,int);
  32. };
  33. int n;
  34. int q;
  35. int k;
  36. int d;
  37. int main(){
  38. scanf("%d%d%d%d",&n,&q,&k,&d);
  39. Node* N=new Node(1,n);
  40. for(int i=0;i<q;i++){
  41. int p,x;
  42. scanf("%d%d",&p,&x);
  43. N->Add(p,x);
  44. if(N->val.maxs>1ll*k*d)
  45. puts("NIE");
  46. else
  47. puts("TAK");
  48. }
  49. return 0;
  50. }
  51. Node::Node(int l,int r):l(l),r(r){
  52. if(l==r)
  53. this->val=Data(-k);
  54. else{
  55. int mid=(l+r)>>1;
  56. this->lch=new Node(l,mid);
  57. this->rch=new Node(mid+1,r);
  58. this->Maintain();
  59. }
  60. }
  61. void Node::Add(int x,int d){
  62. if(this->l==this->r)
  63. this->val=Data(this->val.sum+d);
  64. else{
  65. if(x<=this->lch->r)
  66. this->lch->Add(x,d);
  67. else
  68. this->rch->Add(x,d);
  69. this->Maintain();
  70. }
  71. }
  72. inline void Node::Maintain(){
  73. this->val=this->lch->val+this->rch->val;
  74. }

加强版写的比较蠢...写了一个维护最大子段和的和一个区间加法区间最值的线段树...

  1. #include <bits/stdc++.h>
  2. const int MAXN=1e5+10;
  3. typedef long long intEx;
  4. struct Node{
  5. struct Data{
  6. intEx sum;
  7. intEx lmax;
  8. intEx maxs;
  9. intEx rmax;
  10. Data(){}
  11. Data(intEx val){
  12. this->sum=val;
  13. this->lmax=this->rmax=this->maxs=std::max(this->sum,0ll);
  14. }
  15. Data friend operator+(const Data& a,const Data& b){
  16. Data ans;
  17. ans.sum=a.sum+b.sum;
  18. ans.lmax=std::max(a.lmax,a.sum+b.lmax);
  19. ans.rmax=std::max(a.rmax+b.sum,b.rmax);
  20. ans.maxs=std::max(a.rmax+b.lmax,std::max(a.maxs,b.maxs));
  21. return ans;
  22. }
  23. };
  24. int l;
  25. int r;
  26. Data val;
  27. Node* lch;
  28. Node* rch;
  29. Node(int,int);
  30. void Maintain();
  31. void Add(int,int);
  32. };
  33. struct NodeX{
  34. int l;
  35. int r;
  36. intEx add;
  37. intEx max;
  38. NodeX* lch;
  39. NodeX* rch;
  40. NodeX(int,int);
  41. void PushDown();
  42. void Maintain();
  43. void Add(int,int,int);
  44. void Add(const intEx&);
  45. };
  46. int n;
  47. int q;
  48. int k;
  49. int d;
  50. int main(){
  51. scanf("%d%d%d%d",&n,&q,&k,&d);
  52. Node* N=new Node(1,n);
  53. NodeX* K=new NodeX(0,n-1);
  54. for(int i=0;i<q;i++){
  55. int p,x;
  56. scanf("%d%d",&p,&x);
  57. N->Add(p,x);
  58. if(p!=n)
  59. K->Add(p,n-1,-x);
  60. // printf("%lld %lld\n",N->val.maxs,K->max);
  61. if(N->val.maxs>1ll*k*d||K->max>-N->val.sum)
  62. puts("No");
  63. else
  64. puts("Yes");
  65. }
  66. return 0;
  67. }
  68. NodeX::NodeX(int l,int r):l(l),r(r),add(0){
  69. if(l==r)
  70. this->max=1ll*l*k;
  71. else{
  72. int mid=(l+r)>>1;
  73. this->lch=new NodeX(l,mid);
  74. this->rch=new NodeX(mid+1,r);
  75. this->Maintain();
  76. }
  77. }
  78. void NodeX::Add(const intEx& d){
  79. this->add+=d;
  80. this->max+=d;
  81. }
  82. void NodeX::Add(int l,int r,int d){
  83. if(l<=this->l&&this->r<=r)
  84. this->Add(d);
  85. else{
  86. this->PushDown();
  87. if(l<=this->lch->r)
  88. this->lch->Add(l,r,d);
  89. if(this->rch->l<=r)
  90. this->rch->Add(l,r,d);
  91. this->Maintain();
  92. }
  93. }
  94. void NodeX::PushDown(){
  95. if(this->add){
  96. this->lch->Add(this->add);
  97. this->rch->Add(this->add);
  98. this->add=0;
  99. }
  100. }
  101. void NodeX::Maintain(){
  102. this->max=std::max(this->lch->max,this->rch->max);
  103. }
  104. Node::Node(int l,int r):l(l),r(r){
  105. if(l==r)
  106. this->val=Data(-k);
  107. else{
  108. int mid=(l+r)>>1;
  109. this->lch=new Node(l,mid);
  110. this->rch=new Node(mid+1,r);
  111. this->Maintain();
  112. }
  113. }
  114. void Node::Add(int x,int d){
  115. if(this->l==this->r)
  116. this->val=Data(this->val.sum+d);
  117. else{
  118. if(x<=this->lch->r)
  119. this->lch->Add(x,d);
  120. else
  121. this->rch->Add(x,d);
  122. this->Maintain();
  123. }
  124. }
  125. inline void Node::Maintain(){
  126. this->val=this->lch->val+this->rch->val;
  127. }

[BZOJ 1135][POI2009]Lyz的更多相关文章

  1. bzoj 1135 [POI2009]Lyz 线段树+hall定理

    1135: [POI2009]Lyz Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 573  Solved: 280[Submit][Status][ ...

  2. 1135: [POI2009]Lyz

    1135: [POI2009]Lyz https://lydsy.com/JudgeOnline/problem.php?id=1135 分析: hall定理+线段树连续区间的最大的和. 首先转化为二 ...

  3. 【BZOJ】1135: [POI2009]Lyz

    题意 有\(1\)到\(n(1 \le n \le 200000)\)号的溜冰鞋各\(k(1 \le k \le 10^9)\)双.已知\(x\)号脚的人可以穿\(x\)到\(x+d\)的溜冰鞋. 有 ...

  4. BZOJ1135: [POI2009]Lyz

    1135: [POI2009]Lyz Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 264  Solved: 106[Submit][Status] ...

  5. 【BZOJ1135】[POI2009]Lyz 线段树

    [BZOJ1135][POI2009]Lyz Description 初始时滑冰俱乐部有1到n号的溜冰鞋各k双.已知x号脚的人可以穿x到x+d的溜冰鞋. 有m次操作,每次包含两个数ri,xi代表来了x ...

  6. BZOJ 1115: [POI2009]石子游戏Kam

    1115: [POI2009]石子游戏Kam Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 883  Solved: 545[Submit][Stat ...

  7. BZOJ 1142: [POI2009]Tab

    1142: [POI2009]Tab Time Limit: 40 Sec  Memory Limit: 162 MBSubmit: 213  Solved: 80[Submit][Status][D ...

  8. bzoj 1133: [POI2009]Kon dp

    1133: [POI2009]Kon Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 242  Solved: 81[Submit][Status][D ...

  9. bzoj 1138: [POI2009]Baj 最短回文路 dp优化

    1138: [POI2009]Baj 最短回文路 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 161  Solved: 48[Submit][Sta ...

随机推荐

  1. Java subList、toArray、asList 注意点

    1. ArrayList的subList 结果不可以强转成ArrayList,否则抛出ClassCastException异常,原因是subList返回的是ArrayList的内部类SubList,并 ...

  2. 基于SpringBoot+SSM实现的Dota2资料库智能管理平台

    Dota2资料库智能管理平台的设计与实现 摘    要 当今社会,游戏产业蓬勃发展,如PC端的绝地求生.坦克世界.英雄联盟,再到移动端的王者荣耀.荒野行动的火爆.都离不开科学的游戏管理系统,游戏管理系 ...

  3. 【转载】Java8 HashMap之tableSizeFor

    Java8对许多内置的容器进行了优化与拓展,其中对HashMap的改变尤其大.之后将进行总结. 最近在看HashMap的源码时,发现了里面好多很不错的算法,相比Java7从性能上提高了许多.其中tab ...

  4. Docker基础-搭建本地私有仓库

    1.使用registry镜像创建私有仓库 安装Docker后,可以通过官方提供的registry镜像来简单搭建一套本地私有仓库环境: docker run -d -p 5000:5000 regist ...

  5. Spring源码分析:非懒加载的单例Bean初始化过程(下)

    上文Spring源码分析:非懒加载的单例Bean初始化过程(上),分析了单例的Bean初始化流程,并跟踪代码进入了主流程,看到了Bean是如何被实例化出来的.先贴一下AbstractAutowireC ...

  6. 【学习笔记】--- 老男孩学Python,day3 while 循环 运算符 逻辑、赋值运算

    1. 循环. while循环 while 条件: 代码块(循环体) 执行流程: 1. 判断条件是否为真. 如果真. 执行代码块 2. 再次判断条件是否为真...... 3. 当条件为假.执行else ...

  7. js-ES6学习笔记-正则的扩展

    1.在ES5中,RegExp构造函数的参数有两种情况.第一种情况是,参数是字符串,这时第二个参数表示正则表达式的修饰符(flag).第二种情况是,参数是一个正则表示式,这时会返回一个原有正则表达式的拷 ...

  8. Windows10设置

    按下Win+R键,输入gpedit.msc,打开组策略窗口

  9. 百度Ocr文字识别

    简述 最近开发一个项目需要用到Ocr文字识别技术来识别手写文字,在评估过程中体验了百度的文字识别和腾讯的文字识别.查找官方开发文档,发现它们都有印刷体和手写体两种符合项目需求的识别模式,但是腾讯的手写 ...

  10. KeyPress 和KeyDown 、KeUp之间的区别

    前几天,在写完一个功能模块上线测试的时候,出现了一个诡异的问题.input 框在输入查询内容之后,按回车键居然有两种不同的表现形式(input 框没有绑定键盘事件),谷歌和火狐功能正常,但IE在按了回 ...