http://acm.hunnu.edu.cn/online/?action=problem&type=show&id=11287&courseid=109
题意:
在地球上给一些一矩形,再给一些查询矩形,问哪些矩形与查询矩形相交。
思路:
考虑地球是圆的,变成平面坐标系后经度小于-180和大于180的都要扩展下,
题目说给的纬度、经度的矩形是上顶点和其对角线。换成坐标后是Y坐标为
经度。通过—+360来移到一个平面来。
我的代码:
/*
program:hunnu_11287
author:BlackAndWhile
*/
#include<stdio.h>
#include<math.h>
struct point
{
int x,y;
}a[100001][2],b[100001][2],tb[2],ansa[100001][2];
int T,n,m,i,j,ans;
int bl[100001];
int cross(point a[2],point b[2])//判断矩形交
{
int i,j;
if(a[1].y<=b[1].y&&a[1].y>=b[0].y&&((a[1].x<=b[1].x&&a[0].x>=b[0].x)||(a[1].x>=b[1].x&&a[1].x<=b[0].x)||(a[0].x<=b[0].x&&a[0].x>=b[1].x))) return 1;
if(a[0].y<=b[1].y&&a[0].y>=b[0].y&&((a[1].x<=b[1].x&&a[0].x>=b[0].x)||(a[1].x>=b[1].x&&a[1].x<=b[0].x)||(a[0].x<=b[0].x&&a[0].x>=b[1].x))) return 1;
if(a[1].x>=b[1].x&&a[1].x<=b[0].x&&((a[1].y>=b[1].y&&a[0].y<=b[0].y)||(a[1].y<=b[1].y&&a[1].y>=b[0].y)||(a[0].y>=b[0].y&&a[0].y<=b[1].y))) return 1;
if(a[0].x>=b[1].x&&a[0].x<=b[0].x&&((a[1].y>=b[1].y&&a[0].y<=b[0].y)||(a[1].y<=b[1].y&&a[1].y>=b[0].y)||(a[0].y>=b[0].y&&a[0].y<=b[1].y))) return 1;
return 0;
}
int main()
{
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
for(i=0;i<n;i++)
{
scanf("%d,%d %d,%d",&a[i][0].x,&a[i][0].y,&a[i][1].x,&a[i][1].y);
ansa[i][0]=a[i][0];
ansa[i][1]=a[i][1];
if(a[i][0].y>a[i][1].y) a[i][1].y+=360;//犯2了,我把原始数据改了,一直WA这里
}
scanf("%d",&m);
for(i=0;i<m;i++)
{
scanf("%d,%d %d,%d",&b[i][0].x,&b[i][0].y,&b[i][1].x,&b[i][1].y);
if(b[i][0].y>b[i][1].y) b[i][1].y+=360;
bl[i]=0;
}
ans=0;
for(i=0;i<n;i++)
{
for(j=0;j<m;j++)
{
if(cross(a[i],b[j])||cross(b[j],a[i]))
{
bl[i]=1;
ans=1;
break;
}
else
{
tb[0]=b[i][0];tb[1]=b[i][1];
tb[0].y-=360;tb[1].y-=360;//查询区间左移360
if(cross(a[i],tb)||cross(tb,a[i]))
{
bl[i]=1;
ans=1;
break;
}
tb[0].y+=720;tb[1].y+=720;//右移360
if(cross(a[i],tb)||cross(tb,a[i]))
{
bl[i]=1;
ans=1;
break;
}
}
}
}
if(ans)
{
for(i=0;i<n;i++)
if(bl[i]) printf("%d,%d %d,%d\n",ansa[i][0].x,ansa[i][0].y,ansa[i][1].x,ansa[i][1].y);
}
else
printf("No data found.\n");
printf("\n");
}
return 0;
}