[BZOJ 1120]POI2009 str

最新推荐文章于 2025-03-10 03:51:21 发布
最新推荐文章于 2025-03-10 03:51:21 发布
阅读量1k 收藏
点赞数
分类专栏： POI BZOJ 文章标签： 2009 poi
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/Zeyu_King/article/details/41443203
版权

吐槽一下，好啃的题目。。。。。

我们可以看到在两个点之间的矩形内分界线是直线x+y+C=0或x-y+C=0。矩形之外都是矩形的分界线。

然后思路很明显了，我们将要询问的矩形全记下来，最后用扫描线加树状数组求答案即可。

程序。。。。还是贴上来吧微笑

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int Maxn=100005;
#define lowbit(x) ((x)&(-(x)))
int n,m,z,p,i,j,x0,x1,y0,y1,tot[3];
int px,py,qx,qy,l,r,N,mid,tr[Maxn],a[Maxn],b[Maxn],ANS[Maxn][2];
struct arr
{
  int x,y0,y1,num;
  bool flag;
  void make(int X,int Y0,int Y1,bool Flag,int NUM)
   { x=X; y0=Y0; y1=Y1; flag=Flag; num=NUM; }
  bool operator <(const arr &a)const
    { return x<a.x; }
}Qy[3][Maxn*6];

struct Point
{
  int x,y;
  bool operator <(const Point &a)const
    { return x<a.x||(x==a.x&&y<a.y); }
} dot[Maxn];

void work0(int K){   //x1>x0 y1>y0
  if (x1-x0==y1-y0){
  	Qy[0][tot[0]++].make(x0,  1,y1-1,1,i*100+K);
  	Qy[0][tot[0]++].make(n,   y0+1,m,1,i*100+(K^1));
  	Qy[0][tot[0]++].make(x1-1,y0+1,m,0,i*100+(K^1));
  	
  	Qy[1][tot[1]++].make(x1-1,1,x1+y0-1,1,i*100+K);
  	Qy[1][tot[1]++].make(x0,  1,x1+y0-1,0,i*100+K);
  	Qy[1][tot[1]++].make(x1-1,x1+y0+1,n+m,1,i*100+(K^1));
  	Qy[1][tot[1]++].make(x0,  x1+y0+1,n+m,0,i*100+(K^1));
  	
  } else
  if (x1-x0>y1-y0){
  	px=(x0+x1+y0-y1)/2;
  	qx=(x0+x1+y1-y0)/2;
  	if ( (x0+x1+y0-y1)&1 ){
 	  Qy[0][tot[0]++].make(px,1,m,1,i*100+K);
  	  Qy[0][tot[0]++].make(n, 1,m,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(qx,1,m,0,i*100+(K^1));
  	  
  	  Qy[1][tot[1]++].make(qx,1,px+y1,1,i*100+K);
  	  Qy[1][tot[1]++].make(px,1,px+y1,0,i*100+K);
  	  Qy[1][tot[1]++].make(qx,px+y1+1,n+m,1,i*100+(K^1));
  	  Qy[1][tot[1]++].make(px,px+y1+1,n+m,0,i*100+(K^1));
  	  
  	} else
  	{
  	  Qy[0][tot[0]++].make(px-1,y1,m,1,i*100+K);
  	  Qy[0][tot[0]++].make(px,1,y1-1,1,i*100+K);
  	  Qy[0][tot[0]++].make(n,   1,y0,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(qx,  1,y0,0,i*100+(K^1));
  	  Qy[0][tot[0]++].make(n,   y0+1,m,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(qx-1,y0+1,m,0,i*100+(K^1));
  	  
  	  Qy[1][tot[1]++].make(qx-1,1,px+y1-1,1,i*100+K);
  	  Qy[1][tot[1]++].make(px,  1,px+y1-1,0,i*100+K);
  	  Qy[1][tot[1]++].make(qx-1,px+y1+1,n+m,1,i*100+(K^1));
  	  Qy[1][tot[1]++].make(px,  px+y1+1,n+m,0,i*100+(K^1));
  	  
  	}
  } else
  if (x1-x0<y1-y0){
  	py=(y0+y1+x1-x0)/2;
  	qy=(y0+y1+x0-x1)/2;
  	if ( (y0+y1+x1-x0)&1 ){
  	  Qy[0][tot[0]++].make(x0,1,py,1,i*100+K);
  	  Qy[0][tot[0]++].make(n, 1,qy,1,i*100+K);
  	  Qy[0][tot[0]++].make(x1,1,qy,0,i*100+K);
  	  Qy[0][tot[0]++].make(x0,py+1,m,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(n, qy+1,m,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(x1,qy+1,m,0,i*100+(K^1));
  	  
  	  Qy[1][tot[1]++].make(x1,1,py+x0,1,i*100+K);
  	  Qy[1][tot[1]++].make(x0,1,py+x0,0,i*100+K);
  	  Qy[1][tot[1]++].make(x1,py+x0+1,n+m,1,i*100+(K^1));
  	  Qy[1][tot[1]++].make(x0,py+x0+1,n+m,0,i*100+(K^1));
  	  
    } else
  	{
  	  Qy[0][tot[0]++].make(x0,1,py-1,1,i*100+K);
  	  Qy[0][tot[0]++].make(n, 1,qy-1,1,i*100+K);
  	  Qy[0][tot[0]++].make(x1,1,qy-1,0,i*100+K);
  	  Qy[0][tot[0]++].make(x0,py+1,m,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(n, qy+1,m,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(x1,qy+1,m,0,i*100+(K^1));
  	  
  	  Qy[1][tot[1]++].make(x1,1,py+x0-1,1,i*100+K);
  	  Qy[1][tot[1]++].make(x0,1,py+x0-1,0,i*100+K);
  	  Qy[1][tot[1]++].make(x1,py+x0+1,n+m,1,i*100+(K^1));
  	  Qy[1][tot[1]++].make(x0,py+x0+1,n+m,0,i*100+(K^1));
  	}
  	
  }
}

void work1(int K){   //x1>x0 y0>y1
  if (x1-x0==y0-y1){
  	Qy[0][tot[0]++].make(x0,  y1+1,m,1,i*100+K);
  	Qy[0][tot[0]++].make(n,   1,y0-1,1,i*100+(K^1));
  	Qy[0][tot[0]++].make(x1-1,1,y0-1,0,i*100+(K^1));
  	
  	Qy[2][tot[2]++].make(x1-1,-m,x0-y1-1,1,i*100+K);
  	Qy[2][tot[2]++].make(x0,  -m,x0-y1-1,0,i*100+K);
  	Qy[2][tot[2]++].make(x1-1, x0-y1+1,n,1,i*100+(K^1));
  	Qy[2][tot[2]++].make(x0,   x0-y1+1,n,0,i*100+(K^1));
  } else
  if (x1-x0>y0-y1){
  	px=(x0+x1+y1-y0)/2;
  	qx=(x0+x1+y0-y1)/2;
  	if ( (x0+x1+y1-y0)&1 ){
  	  Qy[0][tot[0]++].make(px,1,m,1,i*100+K);
  	  Qy[0][tot[0]++].make(n, 1,m,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(qx,1,m,0,i*100+(K^1));
  	  
  	  Qy[2][tot[2]++].make(qx,-m,qx-y0,1,i*100+K);
  	  Qy[2][tot[2]++].make(px,-m,qx-y0,0,i*100+K);
  	  Qy[2][tot[2]++].make(qx,qx-y0+1,n,1,i*100+(K^1));
  	  Qy[2][tot[2]++].make(px,qx-y0+1,n,0,i*100+(K^1));
    } else
    {
      Qy[0][tot[0]++].make(px-1,1,y1,1,i*100+K);
      Qy[0][tot[0]++].make(px,y1+1,m,1,i*100+K);
  	  Qy[0][tot[0]++].make(n,   y0,m,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(qx,  y0,m,0,i*100+(K^1));
  	  Qy[0][tot[0]++].make(n,   1,y0-1,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(qx-1,1,y0-1,0,i*100+(K^1));
  	  
  	  Qy[2][tot[2]++].make(qx-1,-m,qx-y0-1,1,i*100+K);
  	  Qy[2][tot[2]++].make(px,  -m,qx-y0-1,0,i*100+K);
  	  Qy[2][tot[2]++].make(qx-1,qx-y0+1,n,1,i*100+(K^1));
  	  Qy[2][tot[2]++].make(px,  qx-y0+1,n,0,i*100+(K^1));
    }
  } else
  {
  	py=(y0+y1+x0-x1)/2;
  	qy=(y0+y1+x1-x0)/2;
  	if ( (y0+y1+x1-x0)&1 ){
  	  Qy[0][tot[0]++].make(x0,py+1,m,1,i*100+K);
  	  Qy[0][tot[0]++].make(n, qy+1,m,1,i*100+K);
  	  Qy[0][tot[0]++].make(x1,qy+1,m,0,i*100+K);
  	  Qy[0][tot[0]++].make(x0,1,py,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(n, 1,qy,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(x1,1,qy,0,i*100+(K^1));
  	  
  	  Qy[2][tot[2]++].make(x1,-m,x0-py-1,1,i*100+K);
  	  Qy[2][tot[2]++].make(x0,-m,x0-py-1,0,i*100+K);
  	  Qy[2][tot[2]++].make(x1,x0-py,n,1,i*100+(K^1));
  	  Qy[2][tot[2]++].make(x0,x0-py,n,0,i*100+(K^1));
  	} else
  	{
  	  Qy[0][tot[0]++].make(x0,py+1,m,1,i*100+K);
  	  Qy[0][tot[0]++].make(n, qy+1,m,1,i*100+K);
  	  Qy[0][tot[0]++].make(x1,qy+1,m,0,i*100+K);
  	  Qy[0][tot[0]++].make(x0,1,py-1,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(n, 1,qy-1,1,i*100+(K^1));
  	  Qy[0][tot[0]++].make(x1,1,qy-1,0,i*100+(K^1));
  	  
  	  Qy[2][tot[2]++].make(x1,-m,x0-py-1,1,i*100+K);
  	  Qy[2][tot[2]++].make(x0,-m,x0-py-1,0,i*100+K);
  	  Qy[2][tot[2]++].make(x1,x0-py+1,n,1,i*100+(K^1));
  	  Qy[2][tot[2]++].make(x0,x0-py+1,n,0,i*100+(K^1));
  	}
  }
}

void work2(int K){
  if ( (y0+y1)&1 ){
  	Qy[0][tot[0]++].make(n,1,(y0+y1)/2,1,i*100+K);
  	Qy[0][tot[0]++].make(n,(y0+y1)/2+1,m,1,i*100+(K^1));
  } else
  {
  	Qy[0][tot[0]++].make(n,1,(y0+y1)/2-1,1,i*100+K);
  	Qy[0][tot[0]++].make(n,(y0+y1)/2+1,m,1,i*100+(K^1));
  }
}

void work3(int K){
  if ( (x0+x1)&1 ){
  	Qy[0][tot[0]++].make((x0+x1)/2,1,m,1,i*100+K);
  	Qy[0][tot[0]++].make(n,1,m,1,i*100+(K^1));
  	Qy[0][tot[0]++].make((x0+x1)/2,1,m,0,i*100+(K^1));
  } else
  {
  	Qy[0][tot[0]++].make((x0+x1)/2-1,1,m,1,i*100+K);
  	Qy[0][tot[0]++].make(n,1,m,1,i*100+(K^1));
  	Qy[0][tot[0]++].make((x0+x1)/2,1,m,0,i*100+(K^1));
  }
}

int low(int x){
  int ret=N+1;
  l=1; r=N;
  while (l<=r){
  	mid=(l+r)>>1;
  	if (a[mid]>=x) ret=mid,r=mid-1;
  	  else l=mid+1;
  }
  return ret;
}

int up(int x){
  int ret=0;
  l=1; r=N;
  while (l<=r){
  	mid=(l+r)>>1;
  	if (a[mid]<=x) ret=mid,l=mid+1;
  	  else r=mid-1;
  }
  return ret;
}

void add(int x){
  for (int i=x;i<=N;i+=lowbit(i))
    tr[i]++;
}

int query(int x){
  int ret=0;
  for (int i=x;i>0;i-=lowbit(i))
    ret+=tr[i];
  return ret;
}

int query(int l,int r){
  if (l>r) return 0;
  return query(r)-query(l-1);
}

void calc(int K){
  memset(tr,0,sizeof(tr));
  sort(Qy[K],Qy[K]+tot[K]);
  for (i=0;i<tot[K];i++){
  	Qy[K][i].y0=low(Qy[K][i].y0);
  	Qy[K][i].y1=up(Qy[K][i].y1);
  }
  for (i=1,j=0;i<=z;i++){
  	for (;j<tot[K]&&Qy[K][j].x<dot[i].x;j++){
  	  int tmp=query(Qy[K][j].y0,Qy[K][j].y1);
  	  if (Qy[K][j].flag) ANS[Qy[K][j].num/100][Qy[K][j].num%100]+=tmp;
  	    else ANS[Qy[K][j].num/100][Qy[K][j].num%100]-=tmp;
  	}
  	add(b[i]);
  	while (i<z && dot[i+1].x==dot[i].x)
  	  add(b[++i]);
  }
  for (;j<tot[K];j++){
  	int tmp=query(Qy[K][j].y0,Qy[K][j].y1);
  	if (Qy[K][j].flag) ANS[Qy[K][j].num/100][Qy[K][j].num%100]+=tmp;
  	  else ANS[Qy[K][j].num/100][Qy[K][j].num%100]-=tmp;
  }

}

void calc(){
  for (i=1;i<=z;i++) a[i]=dot[i].y;
  sort(a+1,a+z+1);
  N=unique(a+1,a+z+1)-a-1;
  for (i=1;i<=z;i++) b[i]=lower_bound(a+1,a+N+1,dot[i].y)-a;
  calc(0);

  for (i=1;i<=z;i++) a[i]=dot[i].x+dot[i].y;
  sort(a+1,a+z+1);
  N=unique(a+1,a+z+1)-a-1;
  for (i=1;i<=z;i++) b[i]=lower_bound(a+1,a+N+1,dot[i].x+dot[i].y)-a;
  calc(1);

  for (i=1,j=0;i<=z;i++) a[i]=dot[i].x-dot[i].y;
  sort(a+1,a+z+1);
  N=unique(a+1,a+z+1)-a-1;
  for (i=1;i<=z;i++) b[i]=lower_bound(a+1,a+N+1,dot[i].x-dot[i].y)-a;
  calc(2);
}

int main(){
  scanf("%d%d%d%d",&n,&m,&z,&p);
  for (i=1;i<=z;i++)
  	scanf("%d%d",&dot[i].x,&dot[i].y);
  sort(dot+1,dot+z+1);
  
  for (i=1;i<=p;i++){
  	scanf("%d%d",&x0,&y0);
  	scanf("%d%d",&x1,&y1);
  	
  	if (x0==x1){
  	  if (y0<y1) work2(0);
	  else
  	  {
  	  	swap(y1,y0);
		work2(1);
  	  }
  	  continue;
  	}
  	if (y0==y1){
  	  if (x0<x1) work3(0);
	  else
  	  {
  	  	swap(x1,x0);
		work3(1);
  	  }
  	  continue;
  	}
  	
  	if (x0<x1&&y0<y1) work0(0);
  	else if (x1<x0&&y1<y0){
  	  swap(x0,x1);
  	  swap(y0,y1);
	  work0(1);
    }
  	else if (x0<x1&&y1<y0) work1(0);
  	else if (x1<x0&&y0<y1){
  	  swap(x0,x1);
  	  swap(y0,y1);
	  work1(1);
    }
    
  }
  calc();
  for (i=1;i<=p;i++)
    printf("%d %d %d\n",ANS[i][0],ANS[i][1],z-ANS[i][0]-ANS[i][1]);
  return 0;
}


google hack database（最新版本）
03-07
本文档是最新版本的google hack 语法查阅文档，包含了从最开始到2019.3.7 所搜集的所有google hack 语法示例
google hacking
自由无限
09-19 1079
google hacking其实并算不上什么新东西,在早几年我在一些国外站点上就看见过相关的介绍,但是由于当时并没有重视这种技术,认为最多就只是用来找找未改名的mdb或者别人留下的webshell什么的,并无太大实际用途.但是前段时间仔细啃了些资料才猛然发觉google hacking其实并非如此简单...google hacking的简单实现    记得以前看见过一篇文章写的就是简单的通
谷歌黑客数据库指南
最新发布
weixin_44893706的博客
03-10 385
Exploits不仅限于 Web 服务，它可以涉及各种类型的漏洞，包括操作系统、应用程序等。Shellcodes确实是用于攻击系统的一种技术，通常与漏洞利用结合使用。GHDB更多是通过查询暴露的信息来发现潜在的攻击目标，主要用于信息收集。
Google Hacking的实现以及应用
游侠之魂的专栏
10-12 787
Google Hacking其实并算不上什么新东西,在早几年我在一些国外站点上就看见过相关的介绍,但是由于当时并没有重视这种技术,认为最多就只是用来找找未改名的mdb或者别人留下的webshell什么的,并无太大实际用途.但是前段时间仔细啃了些资料才猛然发觉Google Hack
Google-Hacking-Database (GHDB) 项目使用教程
gitblog_00057的博客
05-27 423
Google-Hacking-Database (GHDB) 项目使用教程 Google-Hacking-Database-GHDB The GHDB is an index of search queries (we call them dorks) used to find publicly available in...
Google Hacking Database (GHDB) 开源项目教程
gitblog_00426的博客
08-31 1328
Google Hacking Database (GHDB) 开源项目教程 Google-Hacking-Database-GHDBThe GHDB is an index of search queries (we call them dorks) used to find publicly available information, intended for pentesters and ...
google hacking 技术手册
10-31
### Google Hacking技术手册知识点详解 #### 一、Google Hacking概述 Google Hacking是一种利用Google搜索引擎的强大功能来发现安全漏洞或敏感信息的技术。通过掌握特定的搜索语法和技巧，安全专家可以找到可能被...
Google Hacking从0到1
weixin_45626840的博客
04-02 1550
什么是Google Hack搜索，它和普通的搜索由什么区别？深入理解Google Hacking，学会如何利用Google Hacking进行信息收集。就像百度，在搜索框输入我们要检索的东西/概念/关键字，百度匹配关键字，按照一定的逻辑/规则/算法，返回与之相关的信息。这是现代软件程序提供给我们的一个服务。Google搜索是一个过程，它的目标是找到关于某一主题的信息。这个过程从一个基本搜索开始，通过各种方式对其进行修正，直到只返回相关信息的页面。
Google Hacking语法总结
carmi的博客
01-17 2262
Google Hacking是利用谷歌搜索的强大，来在浩瀚的互联网中搜索到我们需要的信息。轻量级的搜索可以搜素出一些遗留后门，不想被发现的后台入口，中量级的搜索出一些用户信息泄露，源代码泄露，未授权访问等等，重量级的则可能是mdb文件下载，CMS未被锁定install页面，网站配置密码，php远程文件包含漏洞等重要信息。基本搜索:逻辑与：and逻辑或：or逻辑非：–完整匹配：“关键词”通配符：*?
Google hacking
qq_43371350的博客
04-11 3638
如何使用Google hacking
黑客的数据库：Google进阶搜索技巧
Cairo_A的博客
04-18 633
这些黑客的搜索技巧，非常值得借鉴，可以解决一些棘手的搜索问题，但切记不要用于非法的目的。
web渗透--6--Google Hacking
武天旭的博客
09-11 2000
3.3 搜索引擎信息收集 Google是一款功能强大的搜索引擎，能够根据robots协议抓取互联网上几乎所有页面，其中包括大量账号密码等敏感信息。GHDB(Google Hacking Database)存放了大量敏感信息，Google Hacking是指利用Google等搜索引擎搜索信息来进行入侵的技术和行为，攻击者利用搜索引擎高级指令能够轻松获取指定的敏感信息。以下是利用GHDB结合搜索引...
【安全】【信息搜集】Google Hacking
qqchaozai的专栏
05-25 1229
山永远都在那，躲不开，只有持续的登高，然后，踩下它！ ~~ 怕山太高？怕自己能力不够？怕技术太深？一辈子怕下去？ ~~ 再复杂的技术，归根结底就是0 和 1，何惧哉！
Google Hacking使用详解
热门推荐
lucashu的博客
04-10 2万+
GoogleHacking：利用搜索引擎有争对性的搜索信息来对网络入侵的技术和行为。搜索引擎对于搜索的关键字提供了很多种语法，构造出特殊的关键字，能够快速全面的让攻击者挖掘到有价值的信息。 轻量级的搜索可搜索出一些遗留后门，不想被发现的后台入口，中量级的搜索出一些用户信息泄露，源代码泄露，未授权访问等等，重量级的则可能是mdb文件下载，CMS 未被锁定install页面，网站配置密码，SQL注入、远程文件包含漏洞等重要信息。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值