BZOJ.1926.[SDOI2010]粟粟的书架(前缀和 主席树 二分)
题意: 在给定矩形区域内找出最少的数,满足和>=k。输出数的个数。两种数据范围。
0~50 注意到(真没注意到...)P[i,j]<=1000,我们可以利用前缀和预处理.
num[i][j][v]表示(1,1)(i,j)中值>=v的个数,val[i][j][v]表示(1,1)(i,j)值>=v的所有数的和。(不要被什么>=v坑,和二维前缀和一样,只是一个点的初始值为A[i,j]>=k)
依旧二分。(好像可以离线之类的优化空间?不管了)
60~100 在一个数列的区间中找 满足和>=k。用主席树+二分做。注意去重对sz[]有影响!P[i,j]<=1000,不用离散化。
//480556kb 16480ms
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
#define Failed puts("Poor QLW")
int r,c,Q;
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
namespace Problem1
{
const int N=202;
int A[N][N],num[N][N][1002],val[N][N][1002];
inline int Get_num(int x1,int y1,int x2,int y2,int v){
return num[x2][y2][v]-num[x1][y2][v]-num[x2][y1][v]+num[x1][y1][v];
}
inline int Get_val(int x1,int y1,int x2,int y2,int v){
return val[x2][y2][v]-val[x1][y2][v]-val[x2][y1][v]+val[x1][y1][v];
}
void Solve()
{
int Max=0;
for(int i=1; i<=r; ++i)
for(int j=1; j<=c; ++j) Max=std::max(Max,A[i][j]=read());
for(int k=1; k<=Max; ++k)
for(int i=1; i<=r; ++i)
for(int j=1; j<=c; ++j)//每个k直接递推就行啊,sb了。
num[i][j][k]=num[i-1][j][k]+num[i][j-1][k]-num[i-1][j-1][k]+(A[i][j]>=k),
val[i][j][k]=val[i-1][j][k]+val[i][j-1][k]-val[i-1][j-1][k]+(A[i][j]>=k?A[i][j]:0);//括号!要不成了val=(..+..+..)>=k?A:0;
int x1,y1,x2,y2,h;
while(Q--)
{
x1=read()-1,y1=read()-1,x2=read(),y2=read(),h=read();
int l=1,r=Max,mid,res=0;
while(l<=r)
{
if(Get_val(x1,y1,x2,y2,mid=l+r>>1)>=h) l=mid+1,res=mid;
else r=mid-1;
}
if(res) printf("%d\n",Get_num(x1,y1,x2,y2,res)-(Get_val(x1,y1,x2,y2,res)-h)/res/*重复多算*/);
else Failed;
}
}
}
namespace Problem2
{
#define lson son[rt][0]
#define rson son[rt][1]
const int N=5e5+5,S=500005*20,Max=1e3;
int A[N],root[N],tot,son[S][2],sz[S],sum[S];
inline void Update(int rt){
sum[rt]=sum[lson]+sum[rson];
}
void Insert(int rt,int &y,int l,int r,int p)
{
sz[y=++tot]=sz[rt]+1, sum[y]=sum[rt]/*no need?*/;
if(l==r) sum[y]+=p;//s_y=s_x+ref[p]
else{
int m=l+r>>1;
if(p<=m) son[y][1]=rson, Insert(lson,son[y][0],l,m,p);
else son[y][0]=lson, Insert(rson,son[y][1],m+1,r,p);
Update(y);
}
}
int Query(int rt,int y,int l,int r,int k)
{
if(sz[y]<=sz[rt]+k) return sum[y]-sum[rt];
if(l==r) return (sum[y]-sum[rt])/(sz[y]-sz[rt])*k;
if(sz[son[y][1]]>=sz[rson]+k) return Query(rson,son[y][1],(l+r>>1)+1,r,k);
return Query(lson,son[y][0],l,l+r>>1,k-sz[son[y][1]]+sz[rson])+sum[son[y][1]]-sum[rson];
}
void Solve()
{
for(int i=1; i<=c; ++i) A[i]=read();
int x,s,t,h;
for(int i=1; i<=c; ++i) Insert(root[i-1],root[i],1,Max,A[i]);
while(Q--)
{
x=read(),s=read()-1,x=read(),t=read(),h=read();
int l=1,r=t-s+1,mid;
int tst;
while(l<r)
{
if((tst=Query(root[s],root[t],1,Max,mid=l+r>>1))>=h) r=mid;
else l=mid+1;
}
if(r==t-s+1) Failed;
else printf("%d\n",l);
}
}
}
int main()
{
r=read(),c=read(),Q=read();
r==1?Problem2::Solve():Problem1::Solve();
return 0;
}
BZOJ.1926.[SDOI2010]粟粟的书架(前缀和 主席树 二分)的更多相关文章
- BZOJ 4556 [Tjoi2016&Heoi2016]字符串 ——后缀数组 ST表 主席树 二分答案
Solution 1: 后缀数组暴力大法好 #include <map> #include <cmath> #include <queue> #include &l ...
- BZOJ 3439: Kpm的MC密码 (trie+dfs序主席树)
题意 略 分析 把串倒过来插进trietrietrie上, 那么一个串的kpmkpmkpm串就是这个串在trietrietrie上对应的结点的子树下面的所有字符串. 那么像 BZOJ 3551/354 ...
- bzoj 1926: [Sdoi2010]粟粟的书架 (主席树+二分)
链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1926 题面; 1926: [Sdoi2010]粟粟的书架 Time Limit: 30 Se ...
- Bzoj 1926: [Sdoi2010]粟粟的书架(二分答案+乱搞+主席树)
1926: [Sdoi2010]粟粟的书架 Time Limit: 30 Sec Memory Limit: 552 MB Description 幸福幼儿园 B29 班的粟粟是一个聪明机灵.乖巧可爱 ...
- 【刷题】BZOJ 1926 [Sdoi2010]粟粟的书架
Description 幸福幼儿园 B29 班的粟粟是一个聪明机灵.乖巧可爱的小朋友,她的爱好是画画和读书,尤其喜欢 Thomas H. Cormen 的文章.粟粟家中有一个 R行C 列的巨型书架,书 ...
- BZOJ 1926: [Sdoi2010]粟粟的书架(主席树,二分答案)
BZOJ 1926: [Sdoi2010]粟粟的书架(主席树,二分答案) 题意 : 给你一个长为\(R\)宽为\(C\)的矩阵,第\(i\)行\(j\)列的数为\(P_{i,j}\). 有\(m\)次 ...
- bzoj 1926: [Sdoi2010]粟粟的书架
#include<cstdio> #include<iostream> #define N 201 #define M 500008 using namespace std; ...
- BZOJ1926 [Sdoi2010]粟粟的书架 【主席树 + 二分 + 前缀和】
题目 幸福幼儿园 B29 班的粟粟是一个聪明机灵.乖巧可爱的小朋友,她的爱好是画画和读书,尤其喜欢 Thomas H. Co rmen 的文章.粟粟家中有一个 R行C 列的巨型书架,书架的每一个位置都 ...
- 【BZOJ-1926】粟粟的书架 二分 + 前缀和 + 主席树
1926: [Sdoi2010]粟粟的书架 Time Limit: 30 Sec Memory Limit: 552 MBSubmit: 616 Solved: 238[Submit][Statu ...
随机推荐
- spring boot 事务配置
事务的作用这里不细说,相信很多人也在工作中使用过. 那么在spring-boot是如何配置事务的,事实上非常简便. 直接贴代码吧. 首先配置数据源 myqlDataSource,这个就不说了.之前的 ...
- Nginx ab压力测试
20-ab压力测试及nginx性能统计模块 优化的启发,打开的文件太多 Nginx 错误日志显示,打开文件数太多 系统层面 more /proc/sys/net/core/somaxconn 单个Ng ...
- 流媒体技术学习笔记之(十二)Linux(Ubuntu)环境运行EasyDarwin
Debug问题??? ./easydarwin -c ./easydarwin.xml & //这样的话是80端口 ./easydarwin -c ./easydarwin.xml -d // ...
- 经典设计模式-iOS的实现
最近看了<HeadFirst 设计模式>这本书,给组内伙伴准备一次分享,把这次分享记录下来,有需要的可以看看. 这本书主要介绍了四人帮23种经典设计模式中的的14种,也是常用的几种.看完这 ...
- centos7.2 安装 composer
安装Composer https://pkg.phpcomposer.com/#how-to-install-composer # 下载composer.phar curl -sS https:// ...
- 分享一个C#创建Barcode的DLL
用于工作需要产生Barcode,随手从网上找了一个DLL(原文地址忘了) http://files.cnblogs.com/panchunting/barcode_bin.zip 使用非常简单,只需添 ...
- apiCloud app调用浏览器打开网页的方法
在APP调用浏览器有两种方法:1.使用openApp2.使用openWin 两种方法调用浏览器后的效果有一点不同:1.使用openApp调用浏览器后,如果手机内有多个浏览器,会首先弹出选择浏览器的框2 ...
- Java 连接远程Linux 服务器执行 shell 脚本查看 CPU、内存、硬盘信息
pom.xml jar 包支持 <dependency> <groupId>com.jcraft</groupId> <artifactId>jsch& ...
- 在Scrapy项目【内外】使用scrapy shell命令抓取 某网站首页的初步情况
Windows 10家庭中文版,Python 3.6.3,Scrapy 1.5.0, 时隔一月,再次玩Scrapy项目,希望这次可以玩的更进一步. 本文展示使用在 Scrapy项目内.项目外scrap ...
- python实现广度优先搜索和深度优先搜索
图的概念 图表示的是多点之间的连接关系,由节点和边组成.类型分为有向图,无向图,加权图等,任何问题只要能抽象为图,那么就可以应用相应的图算法. 用字典来表示图 这里我们以有向图举例,有向图的邻居节点是 ...