ll@(++N,M,K),@(X,Y)[N],@(A,B)[K],L=1d3,x,y;double c;auto D=[&](ll&x,ll&y,ll i){return(x-X[i])**2+(y-Y[i])**2;};VLL o(K,1e18);VVI l(N);Timer t;t.set();Rand r;while(t.get()<5){c=r.get(0.0,3.2);x=cos(c)*L;y=sin(c)*L;rep(j,N)l[j]={X[j]*x+Y[j]*y,j};sortV(l);rep(i,K){M=A[i]*x+B[i]*y;M=bsearch_min[ll,m,0,N](M<l[m][0]);REP(j,max(0,M-L),min(N,M+L))o[i]<?=D(A[i],B[i],l[j][1]);}}c=sqrt(min[j,1,N](D(*X,*Y,j)));rep(i,K)A[i]=D(A[i],B[i],0);Wt(sum[i,0,K](min(sqrt(A[i]),c+sqrt(o[i]))*2));