【BZOJ】2738 矩阵乘法 整体二分+树状数组

题目传送门

伏地膜LYF大佬，连整体二分都会，不愧是一代神犇。

向大佬低头，这位大佬讲的很详细，并且很好懂。

整体二分，是一种离线的算法，就是对所有状态和询问一起二分。

递归，每次对所有状态进行统计，然后根据当前所有状态的一半将当前所有询问进行划分。

当所有状态只剩下一种时，就可计算当前所有询问的答案。

当然，现在我还是初学整体二分，需要做更多的题目来熟练这种算法。

附上AC代码：

#include <cstdio>
#include <cctype>
#include <algorithm>
#define N 510
#define M 60010
#define mid ((l+r)>>1)
#define lowbit(x) (x&-x)
using namespace std;

struct dian{
	int x,y,w;
	bool operator < (const dian lyf) const {return w<lyf.w;}
}a[N*N];
struct note{
	int sx,sy,ex,ey,k,id;
}q[M],q1[M],q2[M];
int n,m,num,ans[M],t[N][N];

inline char nc(){
	static char ch[100010],*p1=ch,*p2=ch;
	return p1==p2&&(p2=(p1=ch)+fread(ch,1,100010,stdin),p1==p2)?EOF:*p1++;
}

inline void read(int& a){
	static char c=nc(),f=1;
	for (;!isdigit(c);c=nc()) if (c=='-') f=-1;
	for (a=0;isdigit(c);a=a*10+c-'0',c=nc());
	a*=f;return;
}

inline void change(int x,int y,int c){
	for (int i=x; i<=n; i+=lowbit(i))
		for (int j=y; j<=n; j+=lowbit(j))
			t[i][j]+=c;
}

inline int query(int x,int y){
	int sum=0;
	for (int i=x; i; i-=lowbit(i))
		for (int j=y; j; j-=lowbit(j))
			sum+=t[i][j];
	return sum;
}

inline int ask(note Q){
	return query(Q.ex,Q.ey)-query(Q.sx-1,Q.ey)-query(Q.ex,Q.sy-1)+query(Q.sx-1,Q.sy-1);
}

inline void solve(int l,int r,int ql,int qr){
	if (ql>qr) return;
	if (l==r){
		for (int i=ql; i<=qr; ++i) ans[q[i].id]=a[l].w;
		return;
	}
	for (int i=l; i<=mid; ++i) change(a[i].x,a[i].y,1);
	int ll=0,rr=0,tmp;
	for (int i=ql; i<=qr; ++i)
		if ((tmp=ask(q[i]))>=q[i].k) q1[++ll]=q[i];
			else q[i].k-=tmp,q2[++rr]=q[i];
	for (int i=l; i<=mid; ++i) change(a[i].x,a[i].y,-1);
	int len=ql;
	for (int i=1; i<=ll; ++i) q[len++]=q1[i];
	for (int i=1; i<=rr; ++i) q[len++]=q2[i];
	solve(l,mid,ql,ql+ll-1);
	solve(mid+1,r,ql+ll,qr);
	return;
}

int main(void){
	read(n),read(m);
	for (int i=1; i<=n; ++i)
		for (int j=1; j<=n; ++j)
			a[++num].x=i,a[num].y=j,read(a[num].w);
	sort(a+1,a+1+num);
	for (int i=1; i<=m; ++i)
		read(q[i].sx),read(q[i].sy),read(q[i].ex),read(q[i].ey),read(q[i].k),q[i].id=i;
	solve(1,num,1,m);
	for (int i=1; i<=m; ++i) printf("%d\n",ans[i]);
	return 0;
}