float X[300],Y[300],x,y,d,s,t;N,i,j,m;C(c,k){c=0;for(k=N;k--;c>m?m=c:0)c+=k==i|k==j|hypot(X[k]-s,Y[k]-t)<1;}main(){for(;scanf("%d",&N)*N;printf("%d\n",m)){for(i=N;i--;)scanf("%f%f",X+i,Y+i);for(m=1;++i<N;)for(j=i;j--;d<4?s=(X[i]+X[j])/2-(y*=d=sqrt(1/d-.25)),C(t=(Y[i]+Y[j])/2+(x*=d)),s+=x+=y,C(t-=x):0)x=X[i]-X[j],y=Y[i]-Y[j],d=x*x+y*y;}}