给定 N 个格点 (Xi, Yi), i = 1, 2, ..., N 在二维坐标系上。您必须处理 Q 查询,每个查询将以“xy d”的形式给出。令 ABC 为顶点位于 A(x+d, y)、B(x, y) 和 C(x, y+d) 的三角形。对于每个查询,您必须找出有多少给定的格点位于三角形 ABC 的内部或边界上。输入
输入的第一行包含两个空格分隔的整数 N 和 Q。接下来的 N 行中的每一行都有两个空格分隔的整数 Xi、Yi,分别给出格点的 x 和 y 坐标。然后下面的 Q 行包含三个空格分隔的整数 x、y、d,给出一个查询。输出
对于每个查询,在一行上输出一个整数,表示位于三角形内部或边界上的给定格点的数量。约束
1 ≤ N ≤ 300000 (3 * 105)
1 ≤ Q ≤ 200000 (2 * 105)
1 ≤ Xi, Yi ≤ 300000 (3 * 105)
1 ≤ x, y, d ≤ 300000 (3 * 105)
样本
Input 1:
5 3
1 3
1 5
3 6
4 4
2 6
1 5 3
1 5 4
1 1 1
Output 1:
3
3
0
Input 2:
2 2
1 5
3 7
2 5 6
2 3 4
Output 2:
1
0
#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;
int sizer[300000]={0};//sizer[i] tells the number of points having i as the abscissa
int cmp(const void *a,const void *b)
{
return *(int *)a-*(int *)b;
}
int main()
{
int **X=(int **)malloc(300000*sizeof(int));//each pointer points to an array of ints
for(int i=0;i<300000;i++)
X[i]=NULL;
int N,Q;
int x,y,d;
scanf("%d %d",&N,&Q);
scanf("%d %d",&x,&y);
sizer[x-1]++;
for(int i=1;i<=N;i++)
{
scanf("%d %d",&x,&y);
sizer[x-1]++;
X[x-1]=(int *)realloc(X[x-1],sizer[x-1]*sizeof(int));
X[x-1][sizer[x-1]-1]=y;}
for(int i=0;i<300000;i++)
{
if(X[i]!=NULL)
qsort(X[i],sizer[i],sizeof(int),cmp);}//sorts all the ordinates
for(int i=1;i<=Q;i++)//query loop
{
scanf("%d %d %d",&x,&y,&d);
int count=0;
for(int j=x;j<=x+d;j++)//this loop works for each row from x to x+d
{
if(X[j-1]==NULL)//no points on X=j line
continue;
else if(X[j-1][sizer[j-1]-1]<y)//the largest ordinate point on X=j line is <y ie below the triangle
continue;
else if(X[j-1][0]+j>x+y+d)//the smallest ordinate point on X=j line is above the triangle
continue;
int k=0;
for(;X[j-1][k]<y&&k<sizer[j-1];k++);//k is the position from where counting of points begins.moving along the X=j line and counting the points
int v=k;
for(;X[j-1][v]+j<=x+y+d&&v<sizer[j-1];v++);
count=count+v-k;}
printf("%d\n",count);}
return 0;}
当输入案例非常大但在小型测试案例中可以正常工作时,我会在在线法官上获得 SIGSEGV。我哪里错了?