律师行业协会网站建设,上海什么公司最有名,单位网站建设典型材料,网站建设与制作培训通知传送门 文章目录题意#xff1a;题意#xff1a;题意#xff1a;
给你一个n∗mn*mn∗m的矩阵#xff0c;每次询问一个矩形#xff0c;左上角是(x1,y1)(x_1,y_1)(x1,y1)#xff0c;右上角是(x2,y2)(x_2,y_2)(x2,y2)#xff0c;和一个hhh#xff0c;问是否能选出…传送门
文章目录题意题意题意
给你一个n∗mn*mn∗m的矩阵每次询问一个矩形左上角是(x1,y1)(x_1,y_1)(x1,y1)右上角是(x2,y2)(x_2,y_2)(x2,y2)和一个hhh问是否能选出矩阵中的若干个数使其和≥h\ge h≥h能的话求出选的个数最少是多少个。 数据范围
题意
初看这题没看范围还以为是什么高级的数据结构比如二维主席树什么的听说二维主席树随便就能卡成nnn\sqrt nnn的所以一直卡着没什么想法。最后看到了数据范围直呼上当了。 可以将这个题分成两个题来做。 先看第二种n1,m≤5e5n1,m \le 5e5n1,m≤5e5。 这个显然就变成一个序列我们只需要在序列上建立主席树让后查询即可变成了主席树裸题。 看第一种n,m≤200n,m\le200n,m≤200。 有一种贪心的策略那就是尽可能选p(i,j)p(i,j)p(i,j)大的数加上这样选出来的一定个数最少所以我们二分选的最小的p(i,j)p(i,j)p(i,j)让后现在我们需要快速求出来矩阵中≥p(i,j)\ge p(i,j)≥p(i,j)的数的和以及数量这两个显然可以预处理出来具体细节可以看代码。 注意求数量的时候由于可能有多个都是p(i,j)p(i,j)p(i,j)所以要减去多余的。
// Problem: P2468 [SDOI2010]粟粟的书架
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P2468
// Memory Limit: 500 MB
// Time Limit: 3000 ms
//
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize(Ofast,no-stack-protector,unroll-loops,fast-math)
//#pragma GCC target(sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tunenative)
//#pragma GCC optimize(2)
#includecstdio
#includeiostream
#includestring
#includecstring
#includemap
#includecmath
#includecctype
#includevector
#includeset
#includequeue
#includealgorithm
#includesstream
#includectime
#includecstdlib
#define X first
#define Y second
#define L (u1)
#define R (u1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].ltr[u].r1)
#define Len(u) (tr[u].r-tr[u].l1)
#define random(a,b) ((a)rand()%((b)-(a)1))
#define db puts(---)
using namespace std;//void rd_cre() { freopen(d://dp//data.txt,w,stdout); srand(time(NULL)); }
//void rd_ac() { freopen(d://dp//data.txt,r,stdin); freopen(d://dp//AC.txt,w,stdout); }
//void rd_wa() { freopen(d://dp//data.txt,r,stdin); freopen(d://dp//WA.txt,w,stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pairint,int PII;const int N1000010,M2000,mod1e97,INF0x3f3f3f3f;
const double eps1e-6;int n,m,q;
int a[201][201],g[201][201][1001],f[201][201][1001];
int root[N],tot;
struct Node {int l,r;int cnt;int sum;
}tr[N*40];void insert(int p,int q,int l,int r,int pos) {qtot; tr[q]tr[p];tr[q].cnt; tr[q].sumpos;if(lr) return;int mid(lr)1;if(posmid) insert(tr[p].l,tr[q].l,l,mid,pos);else insert(tr[p].r,tr[q].r,mid1,r,pos);
}int query(int p,int q,int l,int r,int sum) {if(lr) {int cnttr[q].cnt-tr[p].cnt;int one(tr[q].sum-tr[p].sum)/cnt;return (sumone-1)/one;}LL alltr[q].sum-tr[p].sum;//coutall sumendl;if(allsum) return -1;int mid(lr)1;LL rstr[tr[q].r].sum-tr[tr[p].r].sum;if(rssum) return query(tr[p].r,tr[q].r,mid1,r,sum);else return query(tr[p].l,tr[q].l,l,mid,sum-rs)tr[tr[q].r].cnt-tr[tr[p].r].cnt;
}void solve1() {nm;for(int i1;in;i) {int x; scanf(%d,x);insert(root[i-1],root[i],1,M,x);}while(q--) {int x,xx,l,r,h; scanf(%d%d%d%d%d,x,l,xx,r,h);int ansquery(root[l-1],root[r],1,M,h);if(ans-1) puts(Poor QLW);else printf(%d\n,ans);}
}bool check(int mid,int x1,int y1,int x2,int y2,int h) {LL sumf[x2][y2][mid]-f[x1-1][y2][mid]-f[x2][y1-1][mid]f[x1-1][y1-1][mid];if(sumh) return false;int ansg[x2][y2][mid1]-g[x1-1][y2][mid1]-g[x2][y1-1][mid1]g[x1-1][y1-1][mid1];h-f[x2][y2][mid1]-f[x1-1][y2][mid1]-f[x2][y1-1][mid1]f[x1-1][y1-1][mid1];totans(hmid-1)/mid;return true;
}void solve2() {for(int i1;in;i) for(int j1;jm;j) scanf(%d,a[i][j]);for(int i1;in;i) {for(int j1;jm;j) {LL sum1,sum2; sum1sum20;for(int k1000;k1;k--) {f[i][j][k]f[i-1][j][k]f[i][j-1][k]-f[i-1][j-1][k];g[i][j][k]g[i-1][j][k]g[i][j-1][k]-g[i-1][j-1][k];if(a[i][j]k) f[i][j][k]k,g[i][j][k]1;}}}for(int i1;in;i) {for(int j1;jm;j) {LL sum1,sum2; sum1sum20;for(int k1000;k1;k--) {f[i][j][k]f[i][j][k1];g[i][j][k]g[i][j][k1];}}}while(q--) {int x1,y1,x2,y2,h; scanf(%d%d%d%d%d,x1,y1,x2,y2,h);int l1,r1000,ans-1; tot-1;while(lr) {int mid(lr)1;if(check(mid,x1,y1,x2,y2,h)) ansmid,lmid1;else rmid-1;}if(ans-1) puts(Poor QLW);else printf(%d\n,tot);}
}int main()
{
// ios::sync_with_stdio(false);
// cin.tie(0); scanf(%d%d%d,n,m,q);if(n1) { solve1(); }else { solve2(); }return 0;
}
/**/
