【二维莫队】【二维分块】bzoj2639 矩形计算
<法一>二维莫队，对n和m分别分块后，对块从上到下从左到右依次编号，询问以左上角所在块编号为第一关键字，以右下角标号为第二关键字排序，转移时非常厉害。
O(q*n*sqrt(n))。
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define N 201
#define M 100001
struct LiSan{int p,v;}t[N*N];
bool operator < (const LiSan &a,const LiSan &b){return a.v<b.v;}
int n,m,nm,zy,a[N][N],b[N*N],num[N][N],id[N][N];
struct ASK{int x1,y1,x2,y2,p;}Q[M];
bool operator < (const ASK &a,const ASK &b)
{return num[a.x1][a.y1]!=num[b.x1][b.y1] ?
num[a.x1][a.y1]<num[b.x1][b.y1] :
id[a.x2][a.y2]<id[b.x2][b.y2];}
void lisan()
{
scanf("%d%d",&n,&m); nm=n*m;
for(int i=1;i<=nm;++i)
{
scanf("%d",&t[i].v);
t[i].p=i;
}
sort(t+1,t+nm+1);
b[t[1].p]=++zy;
for(int i=2;i<=nm;++i)
{
if(t[i].v!=t[i-1].v) ++zy;
b[t[i].p]=zy;
}
int pen=0;
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j)
a[i][j]=b[++pen];
}
void makeblock()
{
int sun=sqrt(n),sum=sqrt(m);
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j)
num[i][j]=(i-1)/sun*sum+(j-1)/sum+1;
}
int ans,T[N*N];
bool vis[N][N];
inline void Insert(const int &x,const int &y)
{
if(!vis[x][y])
{
vis[x][y]=1;
ans+=((T[a[x][y]]++)<<1|1);
}
}
inline void Delete(const int &x,const int &y)
{
if(vis[x][y])
{
vis[x][y]=0;
ans-=((--T[a[x][y]])<<1|1);
}
}
inline void trans(const ASK &a,const ASK &b)
{
int t=min(a.x1-1,b.x2);
for(int i=b.x1;i<=t;++i)
for(int j=b.y1;j<=b.y2;++j)
Insert(i,j);
for(int i=max(b.x1,a.x2+1);i<=b.x2;++i)
for(int j=b.y1;j<=b.y2;++j)
Insert(i,j);
t=min(a.y1-1,b.y2);
for(int j=b.y1;j<=t;++j)
for(int i=b.x1;i<=b.x2;++i)
Insert(i,j);
for(int j=max(b.y1,a.y2+1);j<=b.y2;++j)
for(int i=b.x1;i<=b.x2;++i)
Insert(i,j);
t=min(b.x1-1,a.x2);
for(int i=a.x1;i<=t;++i)
for(int j=a.y1;j<=a.y2;++j)
Delete(i,j);
for(int i=max(a.x1,b.x2+1);i<=a.x2;++i)
for(int j=a.y1;j<=a.y2;++j)
Delete(i,j);
t=min(b.y1-1,a.y2);
for(int j=a.y1;j<=t;++j)
for(int i=a.x1;i<=a.x2;++i)
Delete(i,j);
for(int j=max(a.y1,b.y2+1);j<=a.y2;++j)
for(int i=a.x1;i<=a.x2;++i)
Delete(i,j);
}
int q,anss[M];
int main()
{
// freopen("bzoj2639.in","r",stdin);
// freopen("bzoj2639.out","w",stdout);
lisan();
makeblock();
int pen=0;
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j)
id[i][j]=++pen;
scanf("%d", &q);
for(int i=1;i<=q;++i)
{
scanf("%d%d%d%d",&Q[i].x1,&Q[i].y1,&Q[i].x2,&Q[i].y2);
if(Q[i].x1>Q[i].x2) swap(Q[i].x1,Q[i].x2);
if(Q[i].y1>Q[i].y2) swap(Q[i].y1,Q[i].y2);
Q[i].p=i;
}
sort(Q+1,Q+1+q);
for(int i=Q[1].x1;i<=Q[1].x2;++i)
for(int j=Q[1].y1;j<=Q[1].y2;++j)
Insert(i,j);
anss[Q[1].p]=ans;
for(int i=2;i<=q;++i)
{
trans(Q[i-1],Q[i]);
anss[Q[i].p]=ans;
}
for(int i=1;i<=q;++i)
printf("%d\n",anss[i]);
return 0;
}
<法二>二维分块，分块在矩阵中的扩展，预处理前缀矩阵每个数出现的次数，以及所有“子矩块”的答案。询问的时候整块的部分直接获取，零散的部分暴力转移。TLE。
复杂度O(n^2+q)*n*sqrt(n)，常数较大。
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
int f,C;
inline void R(int &x){
C=0;f=1;
for(;C<'0'||C>'9';C=getchar())if(C=='-')f=-1;
for(x=0;C>='0'&&C<='9';C=getchar())(x*=10)+=(C-'0');
x*=f;
}
inline void P(int x){
if(x<10)putchar(x+'0');
else{P(x/10);putchar(x%10+'0');}
}
#define N 202
#define BN 17
int n,m,nm;
struct Point{int x,y;}num[N][N];
inline bool operator == (const Point &a,const Point &b){return a.x==b.x && a.y==b.y;}
struct LiSan{int p,v;}t[N*N];
inline bool operator < (const LiSan &a,const LiSan &b){return a.v<b.v;}
int b[N*N],a[N][N],zy,szn,szm,ln[BN],rn[BN],lm[BN],rm[BN],sun=1,sum=1;
void makeblock()
{
szn=sqrt(n);
for(;sun*szn<n;++sun)
{
ln[sun]=rn[sun-1]+1;
rn[sun]=sun*szn;
}
ln[sun]=rn[sun-1]+1;
rn[sun]=n;
//
szm=sqrt(m);
for(;sum*szm<m;++sum)
{
lm[sum]=rm[sum-1]+1;
rm[sum]=sum*szm;
}
lm[sum]=rm[sum-1]+1;
rm[sum]=m;
//
for(int i=1;i<=sun;++i)
for(int j=1;j<=sum;++j)
for(int k=ln[i];k<=rn[i];++k)
for(int l=lm[j];l<=rm[j];++l)
num[k][l]=(Point){i,j};
}
int anss[BN][BN][BN][BN];
int T[N*N];
int sumv[BN][BN][N*N];
int haha[BN],hehe[BN];
void Init_Sumv()
{
for(int i=1;i<=sun;++i)
for(int j=1;j<=sum;++j)
{
memcpy(sumv[i][j]+1,sumv[i][j-1]+1,sizeof(int)*zy);
for(int k=ln[i];k<=rn[i];++k)
for(int l=lm[j];l<=rm[j];++l)
++sumv[i][j][a[k][l]];
}
for(int j=1;j<=sum;++j)
for(int i=1;i<=sun;++i)
for(int k=1;k<=zy;++k)
sumv[i][j][k]+=sumv[i-1][j][k];
}
inline int Calc(int X1,int Y1,int X2,int Y2,int v)
{
if(!(X2>=X1 && Y2>=Y1)) return 0;
// int a[20][20];
static int* b=&sumv[0][0][0];
int *x=b+haha[X2]+hehe[Y2]+v;
int *y=b+haha[X1-1]+hehe[Y2]+v;
int *z=b+haha[X2]+hehe[Y1-1]+v;
int *l=b+haha[X1-1]+hehe[Y1-1]+v;
// if(*z!=sumv[X2][Y1-1][v])puts("jhsdhe");
// if(*(sumv+x))!=sumv[X2][Y2][v])puts("BaoJingLa");
return *x-*y-*z+*l;
}
void Init_Ans()
{
for(int i=1;i<=sun;++i)
for(int j=1;j<=sum;++j)
for(int k=i;k<=sun;++k)
for(int l=j;l<=sum;++l)
{
anss[i][j][k][l]=anss[i][j][k-1][l];
for(int o=j;o<=l;++o)
for(int p=ln[k];p<=rn[k];++p)
for(int q=lm[o];q<=rm[o];++q)
{
anss[i][j][k][l]+=((Calc(i,j,k-1,l,a[p][q])+T[a[p][q]])<<1|1);
++T[a[p][q]];
}
for(int o=j;o<=l;++o)
for(int p=ln[k];p<=rn[k];++p)
for(int q=lm[o];q<=rm[o];++q)
--T[a[p][q]];
}
}
int q;
int main()
{
// freopen("bzoj2639.in","r",stdin);
// freopen("bzoj2639.out","w",stdout);
R(n); R(m); nm=n*m;
for(int i(1);i<BN;i++) haha[i]=haha[i-1]+BN*N*N;
for(int i(1);i<BN;i++) hehe[i]=hehe[i-1]+N*N;
for(int i=1;i<=nm;++i)
{
R(t[i].v);
t[i].p=i;
}
sort(t+1,t+nm+1);
b[t[1].p]=++zy;
for(int i=2;i<=nm;++i)
{
if(t[i].v!=t[i-1].v) ++zy;
b[t[i].p]=zy;
}
int pen=0;
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j)
a[i][j]=b[++pen];
makeblock();
Init_Sumv();
Init_Ans();
int X1,Y1,X2,Y2;
R(q);
for(;q;--q)
{
R(X1); R(Y1); R(X2); R(Y2);
if(X1>X2)
swap(X1,X2);
if(Y1>Y2)
swap(Y1,Y2);
Point zs=(Point){num[X1][Y1].x+(num[X1][Y1]==num[X1-1][Y1]),
num[X1][Y1].y+(num[X1][Y1]==num[X1][Y1-1])};
Point yx=(Point){num[X2][Y2].x-(num[X2][Y2]==num[X2+1][Y2]),
num[X2][Y2].y-(num[X2][Y2]==num[X2][Y2+1])};
int ans=0;
if(!(yx.x>=zs.x && yx.y>=zs.y))
{
for(int i=X1;i<=X2;++i)
for(int j=Y1;j<=Y2;++j)
{
ans+=((T[a[i][j]])<<1|1);
++T[a[i][j]];
}
for(int i=X1;i<=X2;++i)
for(int j=Y1;j<=Y2;++j)
--T[a[i][j]];
}
else
{
ans=anss[zs.x][zs.y][yx.x][yx.y];
int TX=ln[zs.x]-1;
int TY=rm[yx.y];
for(int i=X1;i<=TX;++i)
for(int j=Y1;j<=TY;++j)
{
ans+=((T[a[i][j]]+Calc(zs.x,zs.y,yx.x,yx.y,a[i][j]))<<1|1);
++T[a[i][j]];
}
int TX2=rn[yx.x];
int TY2=rm[yx.y]+1;
for(int i=X1;i<=TX2;++i)
for(int j=TY2;j<=Y2;++j)
{
ans+=((T[a[i][j]]+Calc(zs.x,zs.y,yx.x,yx.y,a[i][j]))<<1|1);
++T[a[i][j]];
}
int TX3=rn[yx.x]+1;
int TY3=lm[zs.y];
for(int i=TX3;i<=X2;++i)
for(int j=TY3;j<=Y2;++j)
{
ans+=((T[a[i][j]]+Calc(zs.x,zs.y,yx.x,yx.y,a[i][j]))<<1|1);
++T[a[i][j]];
}
int TX4=ln[zs.x];
int TY4=lm[zs.y]-1;
for(int i=TX4;i<=X2;++i)
for(int j=Y1;j<=TY4;++j)
{
ans+=((T[a[i][j]]+Calc(zs.x,zs.y,yx.x,yx.y,a[i][j]))<<1|1);
++T[a[i][j]];
}
for(int i=X1;i<=TX;++i)
for(int j=Y1;j<=TY;++j)
--T[a[i][j]];
for(int i=X1;i<=TX2;++i)
for(int j=TY2;j<=Y2;++j)
--T[a[i][j]];
for(int i=TX3;i<=X2;++i)
for(int j=TY3;j<=Y2;++j)
--T[a[i][j]];
for(int i=TX4;i<=X2;++i)
for(int j=Y1;j<=TY4;++j)
--T[a[i][j]];
}
P(ans),puts("");
}
return 0;
}
【二维莫队】【二维分块】bzoj2639 矩形计算的更多相关文章
- BZOJ.2639.矩形计算(二维莫队)
题目链接 二维莫队,按x,y坐标一起分块.(x,y)的所属的块为 x/sq(n)*sq(m) + y/sq(m) 排序时按照(左下点所在块,右上点的标号)排序 排序后 先得出一个询问的答案,然后利用上 ...
- [CSP-S模拟测试]:蔬菜（二维莫队）
题目描述 小$C$在家中开垦了一块菜地,可以抽象成一个$r\times c$大小的矩形区域,菜地的每个位置都种着一种蔬菜.秋天到了,小$C$家的菜地丰收了. 小$C$拟定了$q$种采摘蔬菜的计划,计划 ...
- csp-s模拟测试50(9.22)「施工（单调栈优化DP）」·「蔬菜（二维莫队？？？）」·「联盟（树上直径）」
改了两天,终于将T1,T3毒瘤题改完了... T1 施工(单调栈优化DP) 考场上只想到了n*hmaxn*hmaxn的DP,用线段树优化一下变成n*hmaxn*log但显然不是正解 正解是很**的单调 ...
- 洛谷 P4887 -【模板】莫队二次离线（第十四分块(前体)）（莫队二次离线）
题面传送门 莫队二次离线 mol ban tea,大概是这道题让我第一次听说有这东西? 首先看到这类数数对的问题可以考虑莫队,记 \(S\) 为二进制下有 \(k\) 个 \(1\) 的数集,我们实时 ...
- luogu P4887 莫队二次离线
珂朵莉给了你一个序列$a$,每次查询给一个区间$[l,r]$ 查询$l≤i<j≤r$,且$ai⊕aj$的二进制表示下有$k$个$1$的二元组$(i,j)$的个数.$⊕$是指按位异或. 直接暴力莫 ...
- luogu P4887 模板 莫队二次离线 莫队 离线
LINK:模板莫队二次离线 很早以前学的知识点 不过 很久了忘了. 考虑暴力 :每次莫队更新的时候 尝试更新一个点到一个区间的答案 可以枚举二进制下位数为k的数字 看一下区间内的这种数字有多少个. 不 ...
- [BZOJ2639]矩形计算
[BZOJ2639]矩形计算 题目大意: 给定一个\(n\times m(n,m\le200)\)的矩阵.\(q(q\le10^5)\)次询问,每次询问一个子矩阵中所有数字出现次数的平方和. 思路: ...
- 2019.01.08 bzoj3809: Gty的二逼妹子序列（莫队+权值分块）
传送门 题意:多组询问,问区间[l,r]中权值在[a,b]间的数的种类数. 看了一眼大家应该都知道要莫队了吧. 然后很容易想到用树状数组优化修改和查询做到O(mnlogamax)O(m\sqrt nl ...
- BZOJ 4241: 历史研究——莫队 二叉堆
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4241 题意:N个int范围内的数,M次询问一个区间最大的(数字*出现次数)(加权众数),可以 ...
随机推荐
- Codeforces 937.C Save Energy!
C. Save Energy! time limit per test 1 second memory limit per test 256 megabytes input standard inpu ...
- SICAU-OJ： 三角关系
三角关系 题意: 给出两个数n和k,统计(a,b,c)三元组满足(a+b)%k=0,(b+c)%k=0,(a+c)%k=0且1<=a,b,c<=n的数量. 题解: 由(a+b)%k=0,( ...
- centos 安装mysql 笔记
1.查询已安装软件的目录 rpm -ql mysql 2.mysql的安装卸载 a. 查找已安装的myslq 版本: #rpm -qa | grep mysql (注意大小写,如果mysql 不行 ...
- MySQL 8.0.11(zip)安装及配置
(1)下载MySQL8.0.11: (2)解压zip文件: 我解压到了D:/MySQL/mysql-8.0.11-winx64 (3)配置环境变量: 右键此电脑->属性 高级系统设置 环境变 ...
- xml数据格式
<?xml version="1.0" encoding="utf-8"?> <country> <name>中国</ ...
- 关于fragment点击能穿透问题
本人在做项目的过程中遇到的这个问题,然后就在网上百度了一下,之后也是在csdn上看到博友发过此类问题的解决办法,所以特此重新总结一下,顺便也给自己提个醒,避免出现此类问题.好!下面我们说一下问题: 举 ...
- bzoj 3208 暴力
对于每个操作,直接暴力做就行了,询问的话搜一遍,然后 就这么水过去了. /************************************************************** ...
- 基于x64的处理器意思
基于x64的处理器意思是CPU的架构是X64的,也是64位的CPU. 基本简介: "x86-64",有时会简称为"x64",是64位微处理器架构及其相应指令集的 ...
- tmux下vim颜色不正常问题
在解决了tmux下,make menuconfig颜色不正常问题https://www.cnblogs.com/zqb-all/p/9702582.html后,引入了新的问题,vim颜色错乱. 尝试了 ...
- 【转】使用者角度看bionic pthread_mutex和linux futex实现
使用者角度看bionic pthread_mutex和linux futex实现 本文所大篇幅引用的参考文章主要描述针对glibc和pthread实现:而本文的考察代码主要是android的bioni ...