hdu4562 守护雅典娜

先筛选出来两种圆，一种是包含雅典娜，不包含怪兽的，一种是不包含雅典娜，包含怪兽的，然后圆的半径大小排序，分别对两种圆dp，最后合并两种DP，求出最大的就是结果

#include <cstdio>
#include <set>
#include <iostream>
#include <string>
#include <queue>
#include <cstring>
#include <algorithm>
#define LL long long
#define Inf (1<<30)
#define max(x,y)  x>y?x:y
using namespace std;
int n,X,Y;
int cntp,cntq;
struct node
{
	int x,y,r;
	int num;
}p[1100],q[1100];
bool cmp(node a,node b)
{
	return a.r<b.r;
}
int judge(node a,node b)//判断两个圆是否内含或者相离，是的话返回对
{
	int d=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
	if((a.r-b.r)*(a.r-b.r)<=d&&d<=(a.r+b.r)*(a.r+b.r))return 0;
	return 1;
}
int main()
{	
	int i,j,t,o=1;
	//freopen("D:\\oo.txt","r",stdin);
	scanf("%d",&t);
	while(t--)
	{
		cntp=cntq=0;
		scanf("%d%d%d",&n,&X,&Y);		
		for(i=0;i<n;i++)
		{
			int x,y,r;
			scanf("%d%d%d",&x,&y,&r);
			int d1=x*x+y*y;
			int d2=(x-X)*(x-X)+(y-Y)*(y-Y);
			if(d1<r*r&&d2>r*r)//包含雅典娜，不包含怪兽
			{
				p[cntp].x=x;p[cntp].y=y;p[cntp].r=r;p[cntp].num=1;
				cntp++;
			}
			if(d1>r*r&&d2<r*r)//不包含雅典娜，包含怪兽
			{
				q[cntq].x=x;q[cntq].y=y;q[cntq].r=r;q[cntq].num=1;
				cntq++;
			}
		}
		int ans=0;
		sort(p,p+cntp,cmp);
		sort(q,q+cntq,cmp);
		for(i=1;i<cntp;i++)
		{
			for(j=0;j<i;j++)
				if(judge(p[i],p[j]))p[i].num=max(p[i].num,p[j].num+1);
			ans=max(ans,p[i].num);
		}
		if(cntp==1)ans=max(p[0].num,ans);//比较坑，因为这里错了很多次
		for(i=1;i<cntq;i++)
		{
			for(j=0;j<i;j++)
				if(judge(q[i],q[j]))q[i].num=max(q[i].num,q[j].num+1);
			ans=max(ans,q[i].num);
		}		
		if(cntq==1)ans=max(q[0].num,ans);//比较坑，因为这里错了很多次
		for(i=0;i<cntp;i++)
		{
			for(j=0;j<cntq;j++)
				if(judge(p[i],q[j]))ans=max(ans,p[i].num+q[j].num);
		}
		printf("Case %d: %d\n",o++,ans);
	}
    return 0;
}