hdu 3511

  题目意思实际上就是找嵌套数最大的圆并输出其嵌套数。

  一开始见时间有5000ms，就果断暴力，然后果断超时— — 之后用二分，终于不tle啦，但wa。。。TAT。上网找题解，发现有两种方法，一个是用随机数，另一个是计算几何中的扫描法，其中用了set，没见过，于是好好学习了一下，发现思想还是很好懂的，该解法的所在地址是：

http://hi.baidu.com/bobobry/blog/item/e270bf097b3bd9376b60fbbe.html

 

以下是计算几何的解法：

 

#include<cstdio>
#include<cstring>
#include<stdlib.h>
#include<set>
#include<cmath>
using namespace std;
const int M=50100;
int max(int a,int b) {return a>b?a:b;}

int n,t;
struct cirs
{
 int x,y;
 int r,w;
}cir[M];
struct event
{
 int x;
 int id,mark;   //mark为1，则圆为运算对象，反之，圆要删出set
}eve[2*M];
int nowx;

double cacul(int id,int mark)
{
 double d;
 d=sqrt((double)cir[id].r*cir[id].r-(double)(nowx-cir[id].x)*(nowx-cir[id].x));
 return mark?(double)cir[id].y+d:(double)cir[id].y-d;
}

bool equal(double a,double b)
{
 if(fabs(a-b)<1e-8) return 1;
 else return 0;
}

struct nodes
{
 int id;
 int mark;
 bool operator<(const nodes a)const   //set的<要重载
 {
  double y1=cacul(id,mark);
  double y2=cacul(a.id,a.mark);
  return y1>y2 || equal(y1,y2) && mark>a.mark;
 }
};
set<nodes> line;
set<nodes>::iterator it,pre,succ;

int cmp(const void *a,const void *b)
{
 int t;
 struct event *x,*y;
 x=(struct event *)a;
 y=(struct event *)b;
 t=x->x-y->x;
 if(t!=0) return t;
 else t=cir[y->id].y-cir[x->id].y;
 return t;
}

void move_line()     //计算几何中的扫描法
{
 line.clear();
 int i,j;
 for(i=0;i<t;i++)
 {
  nowx=eve[i].x;
  if(eve[i].mark==1)    //这个圆没运算过
  {
   struct nodes node;
   node.id=eve[i].id;
   node.mark=1;
   it=line.insert(node).first;
   succ=pre=it;
   succ++;
   if(it==line.begin() || succ==line.end())
   {
    cir[it->id].w=1;           //此时是最外层的圆
   }
   else               //当前圆嵌套与前一个圆中
   {
    pre--;
    if(pre->id==succ->id)
    {
     cir[it->id].w=cir[pre->id].w+1;
    }
    else           //处于同一嵌套级别的圆的嵌套数
    {
     cir[it->id].w=max(cir[pre->id].w,cir[succ->id].w);
    }
   }
   struct nodes node2;
   node2.id=eve[i].id;
   node2.mark=0;
   line.insert(node2).second;
  }
  else              //该圆的右点以被扫描到，该离开set了
  {
   struct nodes node,node2;
   node.id=eve[i].id;
   node.mark=1;
   line.erase(node);
   node2.mark=0;
   node2.id=eve[i].id;
   line.erase(node2);
  }
 }
}

int main()
{
// freopen("in.txt","r",stdin);
 int i,j;
 while(scanf("%d",&n)==1)
 {
  t=0;
  for(i=0;i<n;i++)
  {
   scanf("%d%d%d",&cir[i].x,&cir[i].y,&cir[i].r);
   cir[i].w=0;
   eve[t].x=cir[i].x-cir[i].r;
   eve[t].id=i;
   eve[t].mark=1;
   t++;
   eve[t].x=cir[i].x+cir[i].r;
   eve[t].id=i;
   eve[t].mark=0;
   t++;
  }
  qsort(eve,t,sizeof(eve[0]),cmp);
  move_line();
  int ans=-1;
  for(i=0;i<n;i++)
  {
   if(cir[i].w>ans)
    ans=cir[i].w;
  }
  printf("%d\n",ans);
 }
 return 0;
}

 

接下的这份随机函数的代码就真的很考验rp的哈，我交了8次才过了3次。。。两次是直接拿那位大牛的来交的。。。。

 

#include<cstdio>
#include<cstring>
#include <cstdlib>
#include<ctime>
#include<iostream>
using namespace std;
const int M=50000;

struct cir
{
 double x,y,r;
}c[M];

bool vis[M];
int n;
int main()
{

 srand((unsigned)time(NULL));
 while(scanf("%d",&n)==1)
 {
  int i,j,ans=-1;
  for(i=0;i<n;i++)
   scanf("%lf%lf%lf",&c[i].x,&c[i].y,&c[i].r);
  memset(vis,0,sizeof(vis));
  int tmp=1;

  for(j=1;j<n;j++)
  {
   if((c[0].x-c[j].x)*(c[0].x-c[j].x)+(c[0].y-c[j].y)*(c[0].y-c[j].y)<c[j].r*c[j].r)
    tmp++;
  }

  if(tmp>ans) tmp=ans;

  int rd;
  for(i=0;i<n && i<5000 && ans!=n;)
  {
   rd=rand()%n;
   if(vis[rd]) continue;
   vis[rd]=1;
   i++;
   tmp=1;
   for(j=0;j<n;j++)
   {
    if(rd==j) continue;
    if((c[rd].x-c[j].x)*(c[rd].x-c[j].x)+(c[rd].y-c[j].y)*(c[rd].y-c[j].y)<c[j].r*c[j].r)
     tmp++;
   }
   if(tmp>ans) ans=tmp;
  }

  printf("%d\n",ans);
 }
 return 0;
}