ZOJ 3993 Safest Buildings （数学题）

高中数学学得还凑合，但是开始就猜的心态必须要制止，开始就猜就过分了，以后还是必须有理论证明再去写。西安留下了阴影，是我功利心太强。想要数据的可以拿走下面的测试数据

转化为线性的来做，-r1到r1是圆心的范围。R到R是全集

#include <string.h>
#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <math.h>
using namespace std;
struct node
{
    int i;
    double x,y,l,val;
}a[105];
int b[105];
bool cmp(node x,node y)
{
    return x.val>y.val;
}
int main()
{
   int t;
   //ios::sync_with_stdio(false);
   cin>>t;
   while(t--)
   {
       int n;
       memset(b,0,sizeof(b));
       memset(a,0,sizeof(a));
       double R,r,r1;
       cin>>n>>R>>r;
       r1=R-r;
       for(int i=1;i<=n;i++)
       {
           cin>>a[i].x>>a[i].y;
           a[i].i=i;
           a[i].l=sqrt(a[i].x*a[i].x+a[i].y*a[i].y);
           if(a[i].l>R)
           {
               a[i].val=-1;
               continue;
           }
           if(a[i].l>=r1)
           {
               a[i].val=min(2*r1,r1-a[i].l+r);
           }
           else
           {
               //cout<<"R"<<R<<" "<<r<<" "<<r1<<endl;
               a[i].val+=min(r,r1-a[i].l);
               a[i].val+=min(r,r1+a[i].l);
             //  cout<<a[i].val<<endl;
           }

       }
       sort(a+1,a+1+n,cmp);
       int cnt=0;
       if(a[1].val>=0)b[cnt++]=a[1].i;
       for(int i=2;i<=n;i++)
       {
           if(a[i-1].val==a[i].val&&a[i].val>=0)
            b[cnt++]=a[i].i;
           else break;
       }
       sort(b,b+cnt);
       printf("%d\n",cnt);
       for(int i=0;i<cnt;i++)
       {
           if(i==cnt-1)
            printf("%d\n",b[i]);
            else
            printf("%d ",b[i]);
       }
   }
    return 0;
}
/*
100
4 10 5
0 0
3 4
3 5
3 6
4 10 4
0 0
-7 -6
4 5
5 4
5 10 4
0 0
-5 6
-6 5
-7 8
-8 7
5 10000 5000
0 0
10000 10000
9999 1
1     -10000
100  250
*/

-------------------------------希望让我坚持下来的是热爱，而不是仅仅功利心----------------------------