P1527 [国家集训队]矩阵乘法

题意：

给你一个 n×n 的矩阵，每次询问一个子矩形的第 k 小数。

题解：

整体二分稍微加强化
模板题是一个序列，现在升级成一个矩阵求，做法和原理都是一样的
使用整体二分解决的题目，首先要将所有的操作，包括插入、修改、删除和询问按照时间顺序打一个时间戳，存放在操作序列中，不妨称这个序列为 q 。
结果函数solve(l,r,L,R)表示操作序列上第[l,r]个询问的答案已经被确定在[L，R]的值域范围内
这题和模板的不同在于是二维的，所以树状数组也要是二维的
二维前缀和公式:
左上角是(x,y),右下角是(u,v)
ans=query(u,v)-query(x-1,v)-query(u,y-1)+query(x,y)
剩下过程和一维写法差不多，详细看代码
复杂度：O(Qlog_Nlog_Nlog_1e9)

代码：

#include<bits/stdc++.h>
#define debug(a,b) printf("%s = %d\n",a,b)
typedef long long ll;
using namespace std;inline int read(){int s=0,w=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();//s=(s<<3)+(s<<1)+(ch^48);return s*w;
}
const int M=1010000;
const int N=1e3+9;
struct node{int x,y,u,v,k,id;//第k小的数，第id个查询操作 node(int x=0,int y=0,int u=0,int v=0,int k=0,int id=0):x(x),y(y),u(u),v(v),k(k),id(id){}
}q[M],q1[M],q2[M];
int tot=0,n,m,tree[N][N],ans[M];
inline int lowbit(int x)
{return x&(-x);
}
inline void update(int x,int y,int v){for(int i=x;i<=n;i+=lowbit(i)){for(int j=y;j<=n;j+=lowbit(j)){tree[i][j]+=v;}}return ;
}
inline int query(int x,int y){int res=0;for(int i=x;i;i-=lowbit(i)){for(int j=y;j;j-=lowbit(j)){res+=tree[i][j];}}return res;
}
int Matquery(int x)
{/*     (u,y-1)(x,y)(x-1,v)     (u,v)*/return query(q[x].u,q[x].v)-query(q[x].x-1,q[x].v)-query(q[x].u,q[x].y-1)+query(q[x].x-1,q[x].y-1);
}
inline void solve(int ql,int qr,int L,int R)
{if(ql>qr)return ;if(L==R){for(int i=ql;i<=qr;i++)if(q[i].id!=0)ans[q[i].id]=L;return ;}int len1=0,len2=0;int mid=L+R>>1;for(int i=ql;i<=qr;i++){if(q[i].id==0)//插入操作 {if(q[i].k<=mid){update(q[i].x,q[i].y,1);q1[++len1]=q[i];}else q2[++len2]=q[i];}else //查询操作 {int tmp=Matquery(i);if(tmp>=q[i].k){q1[++len1]=q[i];}else {q[i].k-=tmp;q2[++len2]=q[i];}}}for(int i=1;i<=len1;i++)q[ql+i-1]=q1[i];for(int i=1;i<=len2;i++)q[ql+len1+i-1]=q2[i];for(int i=ql;i<=ql+len1-1;i++){if(q[i].id==0&&q[i].k<=mid)update(q[i].x,q[i].y,-1);}solve(ql,ql+len1-1,L,mid);solve(ql+len1,qr,mid+1,R);
}
int main()
{cin>>n>>m;for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){int x;cin>>x;q[++tot]=(node){i,j,0,0,x,0};}}for(int i=1;i<=m;i++){int x,y,a,b,k;cin>>x>>y>>a>>b>>k;q[++tot]=(node){x,y,a,b,k,i};}solve(1,tot,INT_MIN,INT_MAX);for(int i=1;i<=m;i++)printf("%d\n",ans[i]);return 0;
}

P1527 [国家集训队]矩阵乘法相关推荐

洛谷P1527 [国家集训队] 矩阵乘法 [整体二分，二维树状数组]
题目传送门矩阵乘法题目描述给你一个N*N的矩阵,不用算矩阵乘法,但是每次询问一个子矩形的第K小数. 输入输出格式输入格式: 第一行两个数N,Q,表示矩阵大小和询问组数: 接下来N行N列一共N* ...
P1527 [国家集训队]矩阵乘法整体二分 + 二维树状数组
传送门题意: 思路: 算是个整体二分的板子啦,不过这个是二维的矩阵,我们只需要把一位树状数组改成二维的,让后动态维护单点加,区间查询前缀和即可. //#pragma GCC optimize(2) ...
【洛谷1527】 [国家集训队]矩阵乘法（整体二分）
传送门洛谷 Solution 考虑看到什么k小就整体二分套上去试一下. 矩形k小整体二分+二维树状数组就好了. 代码实现 // luogu-judger-enable-o2 /*mail: mlea ...
Luogu4451 [国家集训队]整数的lqp拆分
题目链接:洛谷题目大意:求对于所有$n$的拆分$a_i$,使得$\sum_{i=1}^ma_i=n$,$\prod_{i=1}^mf_{a_i}$之和.其中$f_i$为斐波那契数列的第$i$项. 数 ...
zcmu-1184（矩阵乘法）
1184: 帮我求算一下斐波那契数吧 Time Limit: 1 Sec Memory Limit: 128 MB Submit: 221 Solved: 42 [Submit][Status][ ...
国家集训队论文分类整理[转]
国家集训队论文分类整理转自这里 dalao写的东西,是非常有学习价值的.反正noip不一定会用,但是对以后肯定有用的. 组合数学计数与统计 2001 - 符文杰:<Pólya原理及其应用&g ...
国家集训队论文分类整理
国家集训队论文分类整理组合数学计数与统计 2001 - 符文杰:<Pólya原理及其应用> 2003 - 许智磊:<浅谈补集转化思想在统计问题中的应用> 2007 - 周冬 ...
国家集训队论文分类整理（转）
国家集训队论文分类整理 ----------转自https://www.cnblogs.com/AbandonZHANG/archive/2012/07/21/2601889.html 距离ACM/I ...
国家集训队论文集题目
这些论文还是很不错的,可以当成专题到谷歌上搜着看看,权当拓展下知识. 国家集训队1999论文集陈宏:<数据结构的选择与算法效率--从IOI98试题PICTURE谈起> 来煜坤:<把 ...

P1527 [国家集训队]矩阵乘法

题意：

题解：

代码：

P1527 [国家集训队]矩阵乘法相关推荐

最新文章

热门文章