#include<iostream>
#include<vector>
#include<cmath>
#include<algorithm>
#include<iomanip>
#define int long long
using namespace std;
const int N = 1e3+5;
int t=1, n, m, p;
int a[N], mp[N][N];
double ans=0;
struct note {
    double x, y;
}jud[N], fea[N];
struct edge {
    int a, b;
    double l;
};
bool cmp (edge x, edge y) {
    if(x.l==y.l) {
        if(x.a==y.a) {
            return x.b < y.b;
        }
        return x.a < y.a;
    }
    return x.l < y.l;
}
bool vis[N], vis1[N];
void work(int x)
{
    vector<edge> vec;
    int cnt=0;
    for(int i=1; i<=x; i++)
    {
        vis[i]=0;
        vis1[i]=0;
        cin>>fea[i].x>>fea[i].y;
    }
    
    for(int i=1; i<=n; i++)
    {
        for(int j=1; j<=x; j++)
        {
            double dis = sqrt((jud[i].x-fea[j].x)*(jud[i].x-fea[j].x)+(jud[i].y-fea[j].y)*(jud[i].y-fea[j].y));
            vec.push_back({i, j, dis});
        }
    }
    sort(vec.begin(), vec.end(), cmp);
    for(edge i: vec)
    {
//        cout<<i.a<<' '<<i.b<<' '<<i.l<<'\n';
        if(cnt==n) break;
        if(vis[i.a]==0 && vis1[i.b]==0)
        {
            cnt++;
            ans += i.l;
            vis[i.a]=1;
            vis1[i.b]=1;
        }
    }
}
void solve()
{
    cin>>n>>m>>p;
    for(int i=1; i<=n; i++)
    {
        cin>>jud[i].x>>jud[i].y;
    }
    work(m);
    work(p);
    cout<<fixed<<setprecision(6)<<ans;
}
signed main()
{
    ios::sync_with_stdio(0);
    cin.tie(0), cout.tie(0);
//    cin>>t;
    while(t--) {
        solve();
    }
}
/*
for(int i=1; i<=n; i++)
3 3 3
1 0
2 0
1 3
0 0
3 0
2 2
1 1
2 1
3 2
7.650282
*/