ZOJ 3993Safest Buildings （数学题）

来源：互联网发布：淘宝一口价设置技巧编辑：程序博客网时间：2024/06/06 02:35

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

转化为线性的来做，-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;}/*1004 10 50 03 43 53 64 10 40 0-7 -64 55 45 10 40 0-5 6-6 5-7 8-8 75 10000 50000 010000 100009999 11     -10000100  250*/

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

阅读全文

0 0