此题切记不要从1开始记录每组数据,因为在取模那里很容易出错!!!还是回归到0吧,因为这个WA了好几次了!
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <set>
#include <algorithm>
using namespace std;
double x[40005],y[40005],r[40005];
vector<int>res;
int N;
set<pair<double,int> >BaoLiu;
set<pair<double,int> >::iterator it;
pair<double,int> All[80010];
bool Judge(int i,int j) //此函数为包含判断
{
if((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j])>r[j]*r[j])
{
return false;
}
else return true;
}
void solve()
{
int i,j,id;
j=0;
for(i=0;i<N;i++) //保留左右边界
{
All[j++]=make_pair(x[i]-r[i],i);
All[j++]=make_pair(x[i]+r[i],i+N);
}
sort(All,All+2*N); //边界排序
for(i=0;i<2*N;i++)
{
id=All[i].second%N;
if(All[i].second<N) //扫描到左边界
{
it=BaoLiu.lower_bound(make_pair(y[id],id)); //实际上只用判断离该圆最近的两个圆是否包含该圆
if((it!=BaoLiu.end()&&Judge(id,it->second))||(it!=BaoLiu.begin()&&Judge(id,(--it)->second)))//判断是否被包含
continue;
res.push_back(id); //未被包含则保存下来
BaoLiu.insert(make_pair(y[id],id)); //保留到与扫描线相交并且不被其他圆包含的圆的集合中
}
else
BaoLiu.erase(make_pair(y[id],id)); //扫描到右边界直接把圆从集合中删除就可以,提高效率
}
}
int main()
{
int i;
scanf("%d",&N);
for(i=0;i<N;i++)
{
scanf("%lf%lf%lf",&r[i],&x[i],&y[i]);
}
solve();
sort(res.begin(),res.end()); //序号的排序
printf("%d\n",res.size());
printf("%d",res[0]+1); //序号记得加1,因为题目中是从1开始计序号的
for(i=1;i<res.size();i++)
{
printf(" %d",res[i]+1);
}
return 0;
}