poj 1819 Disks

http://poj.org/problem?id=1819

 #include <cstdio>
 #include <cstring>
 #include <cmath>
 #include <cstdlib>
 #include <algorithm>
 using namespace std;

 const double eps=1e-;
 struct point
 {
     double x;
     double y,r;
 }pp[];
 bool flag[];
 int stack1[],n;

 int main()
 {
     while(scanf("%d",&n)!=EOF)
     {
         double x;
         int p=,i,last;
         for(last=i=; i<=n; i++){
             flag[i]=true;
             scanf("%lf",&pp[i].r);
             int k=;
             pp[i].x=pp[i].y=pp[i].r;
             for(int j=; j<=p; j++)
             {
                 x=pp[stack1[j]].x+*sqrt(pp[i].r*pp[stack1[j]].r);
                 if(x-pp[i].x>eps)
                 {
                     pp[i].x=x;
                     k=j;
                 }
             }
             p=k;
             stack1[++p]=i;
             if(pp[i].x+pp[i].r-pp[last].x-pp[last].r>eps)
             {
                 last=i;
             }
         }
         for(int i=; i<=p; i++) flag[stack1[i]]=;
         for(int j=last+; j<=n; j++)
         {
             flag[j]=true;
         }
         int num=;
         for(int i=; i<=n; i++)
         {
             if(flag[i])
             {
                 num++;
             }
         }
         printf("%d\n",num);
         for(int i=; i<=n; i++)
         {
             if(flag[i])
             {
                 printf("%d\n",i);
             }
         }
     }
     return ;
 }