第一问的做法好像不太一样

首先第二问非常简单,直接在主席树上二分就好了,单次查询的复杂度\(O(logn)\)

第一问并没有想到有二分这种神仙操作,依旧用的是主席树

我们可以对矩阵建出主席树,也就是像二维前缀和那样的主席树

但是众所周知我们写二维前缀和的时候是这么写的

	pre[x][y]+=pre[x-1][y]+pre[x][y-1]-pre[x-1][y-1]

我们发现这个样子我们根本没有办法优秀的建出主席树,因为这个样子还需要容斥

我们可以按照高维前缀和的思路来处理

for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++) c[i][j]+=c[i][j-1];
for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++) c[i][j]+=c[i-1][j];

也就是我们可以直接按照上面的方式来处理,所一我们每次需要的只有像一个主席树里添加一个点和把两个主席树合并的操作了

添加一个点是常规操作,合并两个主席树自然是需要线段树合并了

不过这里启发式合并好像更能保证复杂度的样子

单次查询像第二问一样就可以了,就是拎出四棵主席树来差分就好了

尽管现在预处理变得有些慢,但是现在单次查询的复杂度还是非常优秀的\(O(log(nm))=O(logn+logm)\)

这个题值域范围非常小,于是可以直接用二维前缀和进行二分,好像也是非常优秀的\(O(log(\text{值域范围}))\)

于是这个方法被吊打了

还是放上代码吧

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define re register
#define maxn 500005
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
inline int read()
{
char c=getchar();
int x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();
return x;
}
int n,m,Q;
namespace solve1
{
int rt[201][201];
int a[201][201],b[201*201],c[201][201];
int l[201*201*801],r[201*201*801],d[201*201*801],sum[201*201*801];
int sz,cnt;
inline int find(int x)
{
int l=1,r=sz;
while(l<=r)
{
int mid=l+r>>1;
if(b[mid]==x) return mid;
if(b[mid]<x) l=mid+1;
else r=mid-1;
}
return 0;
}
int change(int pre,int pos,int val,int x,int y)
{
int root=++cnt;
d[root]=d[pre]+1,sum[root]=sum[pre]+val;
if(x==y) return root;
l[root]=l[pre],r[root]=r[pre];
int mid=x+y>>1;
if(pos<=mid) l[root]=change(l[pre],pos,val,x,mid);
else r[root]=change(r[pre],pos,val,mid+1,y);
return root;
}
int merge(int a,int b,int x,int y)
{
if(!a) return b;if(!b) return a;
int root=++cnt;
if(x==y)
{
sum[root]=sum[a]+sum[b];
d[root]=d[a]+d[b];
return root;
}
int mid=x+y>>1;
l[root]=merge(l[a],l[b],x,mid),r[root]=merge(r[a],r[b],mid+1,y);
d[root]=d[l[root]]+d[r[root]],sum[root]=sum[l[root]]+sum[r[root]];
return root;
}
int query(int p1,int p2,int p3,int p4,int x,int y,int k)
{
if(x==y)
{
int t=-b[x];
if(k%t==0) return k/t;
return k/t+1;
}
int now=sum[l[p1]]+sum[l[p2]]-sum[l[p3]]-sum[l[p4]];
int mid=x+y>>1;
if(now<k) return query(r[p1],r[p2],r[p3],r[p4],mid+1,y,k-now)+d[l[p1]]+d[l[p2]]-d[l[p3]]-d[l[p4]];
return query(l[p1],l[p2],l[p3],l[p4],x,mid,k);
}
inline void solve()
{
for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++) a[i][j]=c[i][j]=read(),a[i][j]=-1*a[i][j],b[++sz]=a[i][j];
std::sort(b+1,b+sz+1);
sz=std::unique(b+1,b+sz+1)-b-1;
for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++) a[i][j]=find(a[i][j]);
for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++) rt[i][j]=change(rt[i][j-1],a[i][j],c[i][j],1,sz);
for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++) rt[i][j]=merge(rt[i-1][j],rt[i][j],1,sz);
for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++) c[i][j]+=c[i][j-1];
for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++) c[i][j]+=c[i-1][j];
int sx,sy,dx,dy,K;
while(Q--)
{
sx=read(),sy=read(),dx=read(),dy=read(),K=read();
if(c[dx][dy]+c[sx-1][sy-1]-c[dx][sy-1]-c[sx-1][dy]<K) puts("Poor QLW");
else printf("%d\n",query(rt[dx][dy],rt[sx-1][sy-1],rt[dx][sy-1],rt[sx-1][dy],1,sz,K));
}
}
}
namespace solve2
{
int rt[maxn];
int l[maxn*40],r[maxn*40],d[maxn*40],sum[maxn*40];
int a[maxn],b[maxn],c[maxn];
int sz,cnt;
inline int find(int x)
{
int l=1,r=sz;
while(l<=r)
{
int mid=l+r>>1;
if(b[mid]==x) return mid;
if(b[mid]<x) l=mid+1;
else r=mid-1;
}
return 0;
}
int change(int pre,int pos,int val,int x,int y)
{
int root=++cnt;
d[root]=d[pre]+1,sum[root]=sum[pre]+val;
if(x==y) return root;
l[root]=l[pre],r[root]=r[pre];
int mid=x+y>>1;
if(pos<=mid) l[root]=change(l[pre],pos,val,x,mid);
else r[root]=change(r[pre],pos,val,mid+1,y);
return root;
}
int query(int p1,int p2,int x,int y,int k)
{
if(x==y)
{
int t=-b[x];
if(k%t==0) return k/t;
return k/t+1;
}
int now=sum[l[p1]]-sum[l[p2]];
int mid=x+y>>1;
if(now<k) return query(r[p1],r[p2],mid+1,y,k-now)+d[l[p1]]-d[l[p2]];
return query(l[p1],l[p2],x,mid,k);
}
inline void solve()
{
std::swap(n,m);
for(re int i=1;i<=n;i++) a[i]=-1*read(),b[++sz]=a[i],c[i]=-1*a[i];
std::sort(b+1,b+sz+1);
sz=std::unique(b+1,b+sz+1)-b-1;
for(re int i=1;i<=n;i++) a[i]=find(a[i]);
for(re int i=1;i<=n;i++) rt[i]=change(rt[i-1],a[i],c[i],1,sz);
for(re int i=1;i<=n;i++) c[i]+=c[i-1];
int o,x,y,K;
while(Q--)
{
o=read(),x=read(),o=read(),y=read(),K=read();
if(c[y]-c[x-1]<K) puts("Poor QLW");
else printf("%d\n",query(rt[y],rt[x-1],1,sz,K));
}
}
}
int main()
{
n=read(),m=read(),Q=read();
if(n<=200&m<=200) solve1::solve();
else solve2::solve();
return 0;
}

【[SDOI2010]粟粟的书架】的更多相关文章

  1. bzoj1926[Sdoi2010]粟粟的书架 二分 主席树

    1926: [Sdoi2010]粟粟的书架 Time Limit: 30 Sec  Memory Limit: 552 MBSubmit: 1064  Solved: 421[Submit][Stat ...

  2. bzoj 1926: [Sdoi2010]粟粟的书架 (主席树+二分)

    链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1926 题面; 1926: [Sdoi2010]粟粟的书架 Time Limit: 30 Se ...

  3. [bzoj1926][Sdoi2010]粟粟的书架_二分_主席树

    粟粟的书架 bzoj-1926 Sdoi-2010 题目大意:题目链接 注释:略 想法:分成两个题 前面的我们可以二分,直接二分出来检验即可. 对于R=1的,相当一个数列,我们在上面建立主席树. 然后 ...

  4. Bzoj 1926: [Sdoi2010]粟粟的书架(二分答案+乱搞+主席树)

    1926: [Sdoi2010]粟粟的书架 Time Limit: 30 Sec Memory Limit: 552 MB Description 幸福幼儿园 B29 班的粟粟是一个聪明机灵.乖巧可爱 ...

  5. [SDOI2010]粟粟的书架 [主席树]

    [SDOI2010]粟粟的书架 考虑暴力怎么做 显然是提取出来 (x2-x1+1)*(y2-y1+1) 个数字拿出来 然后从大到小排序 然后就可以按次取数了- 然而接下来看数据范围 \(50\%\ r ...

  6. [BZOJ1926][SDOI2010]粟粟的书架

    BZOJ Luogu Description 幸福幼儿园 B29 班的粟粟是一个聪明机灵.乖巧可爱的小朋友,她的爱好是画画和读书,尤其喜欢 Thomas H. Cormen 的文章.粟粟家中有一个 R ...

  7. BZOJ1926[Sdoi2010]粟粟的书架——二分答案+主席树

    题目描述 幸福幼儿园 B29 班的粟粟是一个聪明机灵.乖巧可爱的小朋友,她的爱好是画画和读书,尤其喜欢 Thomas H. Co rmen 的文章.粟粟家中有一个 R行C 列的巨型书架,书架的每一个位 ...

  8. 【刷题】BZOJ 1926 [Sdoi2010]粟粟的书架

    Description 幸福幼儿园 B29 班的粟粟是一个聪明机灵.乖巧可爱的小朋友,她的爱好是画画和读书,尤其喜欢 Thomas H. Cormen 的文章.粟粟家中有一个 R行C 列的巨型书架,书 ...

  9. BZOJ1926:[SDOI2010]粟粟的书架——题解

    http://www.lydsy.com/JudgeOnline/problem.php?id=1926 https://www.luogu.org/problemnew/show/P2468 幸福幼 ...

  10. 【BZOJ1926】【SDOI2010】粟粟的书架 [主席树]

    粟粟的书架 Time Limit: 30 Sec  Memory Limit: 552 MB[Submit][Status][Discuss] Description 幸福幼儿园 B29 班的粟粟是一 ...

随机推荐

  1. nginx 服务器配置文件指令

    localtion 配置        语法结构: location [ =  ~  ~* ^~ ] uri{ ... }        uri 变量是带匹配的请求字符, 可以是不含正则表达的字符串, ...

  2. DP Intro - OBST

    http://radford.edu/~nokie/classes/360/dp-opt-bst.html Overview Optimal Binary Search Trees - Problem ...

  3. Thinking in java源码下载链接

    Thinking in java书上显示的下载源码到www.mindview.net站点,但是这个站点打不开了,后来找到真正的下载地址,贴于此. http://www.mindviewinc.com/ ...

  4. TOJ 3651 确定比赛名次

    描述 有N个比赛队(1<=N<=500),编号依次为1,2,3,....,N进行比赛,比赛结束后,裁判委员会要将所有参赛队伍从前往后依次排 名,但现在裁判委员会不能直接获得每个队的比赛成绩 ...

  5. Unity3D之OnGUI知识总结

    相对位置参考   http://blog.csdn.net/sunny__chen/article/details/51323265 自适应屏幕收缩  http://www.360doc.com/co ...

  6. windows下openssl config failed

    老大让我们去学学webpack, 我在一开始就遇到了坎, 调侃自己"webpack从入门到放弃", 在windows下, 报错是"openssl config faile ...

  7. Bootstrap table使用知识-转

    http://www.cnblogs.com/landeanfen/p/5005367.html 官方文档:http://bootstrap-table.wenzhixin.net.cn/zh-cn/ ...

  8. ASP.NET 使用 AjaxPro 实现前端跟后台交互

    使用 AjaxPro 进行交互,很多人都写过文章了,为什么还要继续老生常谈呢.因为有一些细节上的东西我们需要注意,因为这些细节如果不注意的话,那么程序会报错,而且维护性较差. 引言 一.首先,还是那句 ...

  9. 【学习笔记】JDBC数据库连接技术(Java Database Connectivity)

    一.JDBC简介 Java是通过JDBC技术实现对各种数据库的访问的,JDBC是Java数据库连接技术的简称.它可以把数据持久保存,是一种持久化机制. 1.持久化 持久化就是将程序中的数据在瞬时状态和 ...

  10. Scala 知识点掌握1

    Scala知识点巩固 1.Scala简介 Scala是一门面向对象和面向函数的编程语言,是一门静态编程语言,如 Java Scala(变量类型在编译阶段确定):源码文件需要基于 JVM 运行的. 动态 ...