POJ 1819 Disks G++ 几何 没掌握

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
//英语    看博友分析    抄博友程序    几何      没掌握  
double r[1008];
double x[1008];//圆点所在位置的x坐标 
int touch[1008];//圆i与其后的圆j相切 
int flag[1008];
int main()
{
	while(1)
	{
		int n;
		cin>>n;
		if(cin.eof()==1)
		{
			break;
		}
		memset(touch,-1,sizeof(touch));//抄博友程序 
		memset(flag,0,sizeof(flag));
		double exp=1e-10;
		for(int i=0;i<n;i++)
		{
			cin>>r[i];
			x[i]=r[i];
		}
		double mx=0;
		int ti=0;
		for(int i=0;i<n;i++)
		{
			for(int j=0;j<i;j++)
			{
				if((x[j]+2*sqrt(r[i]*r[j]))-x[i]>exp)
				{
					x[i]=(x[j]+2*sqrt(r[i]*r[j]));
					touch[i]=j;//抄博友程序
				}
			}
			if(x[i]+r[i]-mx>exp)//抄博友程序 
			{
				mx=x[i]+r[i];
				ti=i;
			}
		} 
	
		int js=0;
		for(int i=0;i<n;i++)
		{
			//if(touch[i]!=-1) //加 wa 
			for(int j=touch[i]+1;j<i;j++)
			{
				if(flag[j]==0)
				{
					flag[j]=1;
					js++;				
				}
			}
		}
		for(int i=ti+1;i<n;i++)//抄博友程序 
		{
			if(flag[i]==0)
			{
				flag[i]=1;
				js++;
			}
		}
		cout<<js<<endl;
		for(int i=0;i<n;i++)
		{
			if(flag[i]==1)
			cout<<i+1<<endl;
		}		
	}
	return 0;
}