判断几何体是否相交 + 正方形的坐标确定+ 输入scanf格式控制

题目来源:http://poj.org/problem?id=3449

1:给定一个正方形 对角的坐标 （x0 ,y0）, (x2,y2)， 求出另2个点(x1,y1), (x3,y3)的坐标公式。(x1,y1)(x3,y3) 坐标可交换，不影响结果

x0+x2 = x1+x3

y0+ y2 = y1+y3

y2-y0 = x3 -x1

x2-x0 = y1-y3

x1 = (x0 + x2 -y2 + y0)/2

y1= (y0 + y2 + x2 -x0)/2

x3=(x0 + x2 + y2 -y0) /2

y3=(y0+y2 -x2 +x0) /2

 

2:判断几何体相交， 等价于 判断几何体的每条线段 与 另一个几何体的每条线段枚举

3：这题 没有考虑 误差 ，也AC, 可以写一个 考虑误差的版本

4：输入输出过分恶心

5： 排序， 对几何体名称进行排序，后，每次从小到大进行 匹配， 每次都是先匹配小的，然后可以实现从小到大的输出

代码如下：

#include <iostream>
#include <algorithm>
#include <stdlib.h>
#include <stack>
#include <iostream>
#include <stdio.h>
#include <string>
#include <string.h>
#include <algorithm>
#include <stdlib.h>
#include <vector>
#include <set>
#include <math.h>
#include <cmath>
#include <map>
#include <stack>
#include <queue>
#define max(a,b) ( (a) > (b)  ? (a) : (b))
#define min( a,b ) ( (a) < (b) ? (a) : (b)  )

using namespace std ;

typedef long long LL ;

struct Point
{
double x,y;
Point(){}
Point(double x,double y):x(x),y(y){}
Point operator -(Point p)
{
return Point(x-p.x , y- p.y);
}
Point operator +( Point p)
{
return Point( x + p.x , y  + p.y);
}
double operator^(Point p){
return x* p.y - y*p.x;
}
bool operator == (Point p)
{
return (x==p.x) && (y==p.y);
}
};
struct data  // 为每个几何体建一个数据结构
{
int n,m; // n 表示 几何体的边数， m，表示 相交几何体的数量。
char c, name[30] ; // c表示几何体的名称，  name 表示各个相交几何体的名称。
Point p[25]; // p为几何体的各个顶点
bool operator <( const data &d)const
{
return c<d.c;
}
void read()
{
printf("c=%c, n=%d, m=%d, name=%s\n", c,n,m,name);
}
}sha[30];
// 判断点p0是否在线段 p1p2内
int on_segment(Point p1, Point p2, Point p0)
{
if( (p1-p0).x  * ( p2-p0).x <= 0   && (p1-p0).y * (p2-p0).y <= 0 )
return 1;
return 0;
}
// 判断 线段p1p2与q1q2是否相交
int intersection( Point p1, Point p2, Point q1, Point q2)
{
double d1= (p2-p1)^(q1- p1);
double d2= (p2-p1)^(q2-p1);
double d3= (q2-q1)^(p1-q1);
double d4= (q2 - q1)^(p2 - q1);
if( (d1==0 && on_segment(p1,p2,q1) ) || ( d2==0 && on_segment( p1,p2,q2) ) ||

(d3==0 && on_segment( q1,q2,p1)) || ( d4==0 && on_segment(q1,q2,p2) )      )
return 1;
if(  d1 * d2 < 0 && d3* d4 <0)
return 1;
return 0;
}
// 判断几何体 a 与几何体 b 是否相交
bool index(int a,int b)
{
int i,j;
for(i=0; i<sha[a].n; i++)
{
for(j=0 ; j<sha[b].n ;j++)
{
if( intersection( sha[a].p[i]  , sha[a].p[(i+1) % sha[a].n], sha[b].p[j], sha[b].p[ (j+1) % sha[b].n ] )  )
return 1;
}
}
return 0;
}
// 已知正方形的(x0,y0)(x2,y2) 求(x1,y1)(x3,y3)的坐标
void square(int i)
{
sha[i].n=4;
sha[i].p[2].x=sha[i].p[1].x;
sha[i].p[2].y=sha[i].p[1].y;
double x0,y0,x2,y2;
x0= sha[i].p[0].x;
y0= sha[i].p[0].y;
x2= sha[i].p[2].x;
y2= sha[i].p[2].y;
sha[i].p[1].x= (x0+ x2 + y2 - y0) / 2.0;
sha[i].p[1].y=( y0 + y2 -x2 +x0 ) / 2.0;
sha[i].p[3].x=(x0+ x2 -y2 + y0) / 2.0;
sha[i].p[3].y = (y0 + y2 +x2 - x0) / 2.0;
}
//已知矩形 （x0,y0）(x1,y1),(x2,y2)求(x3,y3)
void rec(int i)
{
sha[i].n=4;
sha[i].p[3]= sha[i].p[0] + sha[i].p[2]-sha[i].p[1];
}

int  main(){
char str[15];
while(1)
{
int i=0;
while(   1  )
{
scanf("%s",str);
if (str[0] == '-' ||   str[0] == '.') break;

sha[i].c=str[0];
scanf("%s",str);
if(str[0] == 's' || str[0] == 'l')
sha[i].n=2;
else if(str[0] == 'r'  || str[0] == 't') sha[i].n=3;
else scanf("%d",&sha[i].n);
for(int j=0; j<sha[i].n; j++)
{
scanf(" (%lf,%lf)",&sha[i].p[j].x,&sha[i].p[j].y);  // 格式控制输入
}
if(str[0] == 's') square(i);
else if(str[0] == 'r') rec(i);
sha[i].m=0;
i++;
}
if(str[0] == '.') break;
int len =i;
sort(sha,sha+len);   // 排序
for(int j=0; j<len ;j++)  // 几何体匹配
{
for(int k = j+1 ; k< len ; k++)
{
if(index( j, k ))
{
sha[j].name[sha[j].m ++ ] = sha[k].c;
sha[k].name[sha[k].m++ ] = sha[j].c;
}
}
}

for(int j=0 ;j<len ; j++)   // 输出
{
printf("%c ",sha[j].c);
if(sha[j].m==0)
printf("has no intersections\n");
else
{
printf("intersects with ");
if(sha[j].m==1)
printf("%c\n",sha[j].name[0]);
else if( sha [j].m == 2)
printf("%c and %c\n",sha[j].name[0],sha[j].name[1]);
else
{
int k;
for( k=0 ; k< sha[j].m -1 ;k++)
printf("%c, ",sha[j].name[k]   );
printf("and %c\n",sha[j].name[k]);
}
}
}
printf("\n");
}
return 0 ;
}