#include <bits/stdc++.h>

#define maxn 200005
#define db long double
#define int long long

using namespace std;

int n, m;

struct node{
	int l, r;
	db prb;
	node () {}
	node (int x, int y, db z = 0) {
		l = x, r = y, prb = z;
	}
	int len() {
		return r - l;
	}
	db mid() {
		return (l + r) / 2.0;
	}
	bool operator ==(const node x) {
		return l == x.l && r == x.r;
	}
};

db ans;

void solve(vector<node> &a, vector<node> &b, int l, int r) {
	if (a.empty() || b.empty()) return;
	if (a.size() == 1 && b.size() == 1) {
		if (a[0] == b[0]) {
//			cout<<"AD:"<<a[0].len() * a[0].prb * b[0].prb / 3.0<<endl;
			ans += a[0].len() * a[0].prb * b[0].prb / 3.0;
			return;
		}
		if (a[0].len() == 0 && b[0].len() == 0) {
			ans += abs(a[0].l - b[0].l) * a[0].prb * b[0].prb;
			return;
		}
		if (a[0].len() == 0) {
			ans += fabs(a[0].l - b[0].mid()) * a[0].prb * b[0].prb;
			return;
		}
		if (b[0].len() == 0) {
			ans += fabs(b[0].l - a[0].mid()) * a[0].prb * b[0].prb;
			return;
		}
	}
	int mid = (l + r) >> 1;
//	cout<<mid<<" "<<a[0].l<<" "<<a[0].r<<" "<<b[0].l<<" "<<b[0].r<<endl;
	vector<node> la, ra, lb, rb;
	for (auto x : a) {
		if (x.r <= mid) la.push_back(x);
		else if (x.l >= mid) ra.push_back(x);
		else {
			auto s = node(x.l, mid, x.prb * (mid - x.l) / (x.r - x.l));
			auto t = node(mid, x.r, x.prb * (x.r - mid) / (x.r - x.l));
			la.push_back(s);
			ra.push_back(t);
		}
	}
	for (auto x : b) {
		if (x.r <= mid) lb.push_back(x);
		else if (x.l >= mid) rb.push_back(x);
		else {
			auto s = node(x.l, mid, x.prb * (mid - x.l) / (x.r - x.l));
			auto t = node(mid, x.r, x.prb * (x.r - mid) / (x.r - x.l));
			lb.push_back(s);
			rb.push_back(t);
		}
	}
	db s = 0, t = 0;
	db sp = 0, tp = 0;
	for (auto x : la) s += x.mid() * x.prb, sp += x.prb; 
	for (auto x : rb) t += x.mid() * x.prb, tp += x.prb;
//	cout<<"ED:"<<s<<" "<<t<<" "<<sp<<" "<<tp<<endl;
	ans += t * sp - s * tp;
	
	s = t = sp = tp = 0;
	for (auto x : lb) s += x.mid() * x.prb, sp += x.prb; 
	for (auto x : ra) t += x.mid() * x.prb, tp += x.prb;
//	cout<<"ED:"<<s<<" "<<t<<" "<<sp<<" "<<tp<<endl;
	ans += t * sp - s * tp;
	
	solve(la, lb, l, mid);
	solve(ra, rb, mid, r);
}

vector<node> a, b;

signed main() {
	scanf("%lld%lld", &n, &m);
	int sa = 0, sb = 0;
	for (int i = 1; i <= n; i++) {
		int x, y;
		scanf("%lld%lld", &x, &y);
		a.push_back(node(x, y));
		sa += y - x;
	}
	for (int i = 1; i <= m; i++) {
		int x, y;
		scanf("%lld%lld", &x, &y);
		b.push_back(node(x, y));
		sb += y - x;
	}
	if (sa == 0) for (auto &x : a) x.prb = 1000.0 / n;
	else for (auto &x : a) x.prb = x.len() * 1000.0 / sa;
	
	if (sb == 0) for (auto &x : b) x.prb = 1000.0 / m;
	else for (auto &x : b) x.prb = x.len() * 1000.0 / sb;
	
	solve(a, b, (int)-1e9, (int)1e9);
	
	printf("%.12lf\n", (double)ans / 1000000.0);
	return 0;
}