【BZOJ 2738】 矩阵乘法|分块|离线

好像还有其他的高大上的算法

不过分块最简单了

直接分成n块 暴力搞就行

详见popoQQQ http://blog.csdn.net/popoqqq/article/details/41356899

卡过去的 跑的好慢

错过几次

1.范围500不是200

2.最后输出q个 不是n个

3.中间有变量写错

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

const int MAXN=60000+10;
bool b[MAXN];
int x1[MAXN],x2[MAXN],y1[MAXN],y2[MAXN],qk[MAXN];
struct H
{
	int x,y;
	int v;
}num[500*500+10];
int sum[510][510],cnt[510][510];
int n,q;
bool cmp(H a,H b)
{
	return a.v<b.v;
}
int ans[MAXN];
int main()
{
	cin>>n>>q;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			int x,tmp=i*n-n+j;
			scanf("%d",&x);
			num[tmp].x=i;
			num[tmp].y=j;
			num[tmp].v=x;
		}
	sort(num+1,num+n*n+1,cmp);
	for(int i=1;i<=q;i++) scanf("%d %d %d %d %d",&x1[i],&y1[i],&x2[i],&y2[i],&qk[i]);
	for(int i=1;i<=n;i++)
	{
		for(int j=i*n-n+1;j<=i*n;j++)
		{
			cnt[num[j].x][num[j].y]=1;//cout<<num[j].x<<" "<<num
		}
		for(int j=1;j<=n;j++)
			for(int k=1;k<=n;k++)
				sum[j][k]=sum[j][k-1]+sum[j-1][k]-sum[j-1][k-1]+cnt[j][k];//cout<<sum[n][n]<<endl;
		for(int j=1;j<=q;j++)
		{
			if(b[j]) continue;
			int su=sum[x2[j]][y2[j]]-sum[x2[j]][y1[j]-1]-sum[x1[j]-1][y2[j]]+sum[x1[j]-1][y1[j]-1];//cout<<i<<' '<<j<<' '<<sum[n][n]<<endl;
			if(su<qk[j]) continue;
			int k=su-qk[j];//cout<<'*'<<k<<endl;
			int w=i*n;
			for(;k;w--)
				if(num[w].x<=x2[j]&&num[w].x>=x1[j]&&num[w].y<=y2[j]&&num[w].y>=y1[j])
					k--;
			for(;true;w--)
				if(num[w].x<=x2[j]&&num[w].x>=x1[j]&&num[w].y<=y2[j]&&num[w].y>=y1[j])
				{
					ans[j]=num[w].v;
					break;
				}
			b[j]=true;	
		}
	}
	for(int i=1;i<=q;i++) printf("%d\n",ans[i]);
	return 0;
}