2015多校第一场1001 hdu 5298 Solid Geometry Homework

题意：很多圆球和平面把空间分成了多个区域，要求你给分布每个区域的点染色，使相邻的不同，有些点是已经给好颜色的（有时候没有），要求仅用两种颜色，然后输出被染色的点染了什么颜色。

做法：实际上，暴力的按他在面或者球的哪一边，把点分成两类，然后二分图暴力染色，把冲突的情况判掉就好了

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <iostream>
using namespace std;
struct plane{
	__int64  a,b,c,d;
}pla[4111];
struct sphere{
	__int64 x,y,z,r;
}sph[4111];
struct point{
	__int64 x,y,z;
	__int64 color;
	__int64 tag;
}poi[4111];
struct pointtag{
	point po;
	int Id;
}p1[2111],p2[2111];
int m,n,p,q;
int platest(plane pl,point po)
{
	int a=pl.a,b=pl.b,c=pl.c,d=pl.d;
	if(a*po.x+b*po.y+c*po.z+d<=0)return 0;
	else return 1;
}
int sphtest(sphere sp,point po)
{
	__int64 x1=sp.x,y1=sp.y,z1=sp.z,r=sp.r;
	if((x1-po.x)*(x1-po.x)+(y1-po.y)*(y1-po.y)+(z1-po.z)*(z1-po.z)-r*r<=0)return 0;
	else return 1;
}
int main()
{
	int T;
//	freopen("1011.in","r",stdin);
//	freopen("1.out","w",stdout);
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d%d%d",&m,&n,&p,&q);
		for(int i=1;i<=m;i++)
		{
			scanf("%I64d%I64d%I64d%I64d",&pla[i].a,&pla[i].b,&pla[i].c,&pla[i].d);
		}
		for(int i=1;i<=n;i++)
		{
			scanf("%I64d%I64d%I64d%I64d",&sph[i].x,&sph[i].y,&sph[i].z,&sph[i].r);
		}
		for(int i=1;i<=p;i++)
		{
			scanf("%I64d%I64d%I64d",&poi[i].x,&poi[i].y,&poi[i].z);
			poi[i].color=1;
			poi[i].tag=0;
		}
		for(int i=p+1;i<=p+q;i++)
		{
			scanf("%I64d%I64d%I64d",&poi[i].x,&poi[i].y,&poi[i].z);
			poi[i].color=0;
			poi[i].tag=0;
		}
		if(p==0)
		{
			poi[1].color=1;
		}
		for(int i=1;i<=p+q;i++)
		{
			for(int j=1;j<=m;j++)
			{
				poi[i].tag^=platest(pla[j],poi[i]);
			}
		}
		for(int i=1;i<=p+q;i++)
		{
			for(int j=1;j<=n;j++)
			{
				poi[i].tag^=sphtest(sph[j],poi[i]);
			}
		}
		int col[3];
		memset(col,0,sizeof(col));
		for(int i=1;i<=p+q;i++)
		{
			col[poi[i].tag]|=poi[i].color;
		}
		if((col[1]&col[0])!=0)
		{
			printf("Impossible\n");
			if(T>0)printf("\n");
			continue;
		}
		else {
			if(col[1]==0)col[1]=3^col[0];
			else col[0]=col[1]^3;
		}
		int Y;
		bool key=true;
		for(int i=1;i<=p+q;i++)
		{
			if(poi[i].color==1)
			{
				Y=poi[i].tag;
				break;
			}	
		}
		if(!key)printf("Impossible\n");
		else{
			if(p==0)
			for(int i=p+1;i<=p+q;i++)
			{
				printf("Both\n");
			}
			else{
				for(int i=p+1;i<=p+q;i++)
				{
					if(poi[i].tag==Y)printf("Y\n");
					else printf("R\n");
				}
			}
		}
		if(T>0)printf("\n");	
	}
}
/*
1
1 1 0 2
0 0 1 0
0 0 0 2
0 0 4
0 0 -1 

1
0 10 3 10
0 0 0 20000
0 0 0 40000
0 0 0 60000
0 0 0 80000
0 0 0 100000
0 0 0 120000
0 0 0 140000
0 0 0 160000
0 0 0 180000
0 0 0 200000
0 0 0
999999 999999 999999
0 0 91111
0 0 10000
0 0 30000
0 0 50000
0 0 70000
0 0 90000
0 0 110000
0 0 130000
0 0 150000
0 0 170000
0 0 190000
*/