//#pragma GCC optimize(2)
#include <bits/stdc++.h>
#define F first
#define S second
#define endl '\n'
#define int ll
#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<__int128,__int128> pii;
//typedef pair<long double,long double> pii;
// Custom input and output operators for __int128
istream& operator>>(istream& in, __int128& x) {
	string s;
	in >> s;
	x = 0;
	bool neg = false;
	int i = 0;
	if (s[0] == '-') {
		neg = true;
		i = 1;
	}
	for (; i < s.size(); i++) {
		x = x * 10 + (s[i] - '0');
	}
	if (neg) x = -x;
	return in;
}

ostream& operator<<(ostream& out, __int128 x) {
	if (x == 0) {
		out << "0";
		return out;
	}
	if (x < 0) {
		out << "-";
		x = -x;
	}
	string s;
	while (x > 0) {
		s += '0' + x % 10;
		x /= 10;
	}
	reverse(s.begin(), s.end());
	out << s;
	return out;
}
const int N = 4e5+10;
int n,m;
vector<pii> work(vector<pii> a, vector<pii> b)
{
	vector<pii> now;
	vector<int> pos;
	for(auto [l,r]:b)
	{
		pos.push_back(l);
		pos.push_back(r);
	}
	int j = 0;
	for(auto [l,r] : a)
	{
		while(j<pos.size() && pos[j] < l)
		{
			j++;
		}
		if(j<pos.size() && pos[j]>=l && pos[j]<=r)
		{
			while(j<pos.size() && pos[j]<=r)
			{
				now.push_back({l,pos[j]});
				l = pos[j];
				j++;
			}
			if(l<=r) now.push_back({l,r});
		}
		else if(l!=r) now.push_back({l,r});
	}
	return now;
}
long double sufl[N];
long double suflm[N];
long double prel[N], prelm[N];
void solve()
{
	cin>>n>>m;
	vector<pii> a(n),b(m);
	for(int i=0;i<n;i++)  cin>>a[i].F>>a[i].S;
	for(int i=0;i<m;i++) cin>>b[i].F>>b[i].S;
	sort(a.begin(), a.end()); sort(b.begin(), b.end());
	int flaga = 1;
	for(auto [l,r]:a)
	{
		if(l!=r) flaga = 0;
	}
	int flagb = 1;
	for(auto [l,r]:b)
	{
		if(l!=r) flagb = 0;
	}
	vector<pii> na = a;
	vector<pii> nb = b;
	
	if(flaga == 0) na = work(a,b);
	if(flagb == 0) nb = work(b,a);
	
	
	
	na.erase(unique(na.begin(), na.end()), na.end());
	nb.erase(unique(nb.begin(), nb.end()), nb.end());


	n = na.size();
	m = nb.size();
	
	for(int i=m-1;i>=0;i--)
	{
		sufl[i] = nb[i].S - nb[i].F;
		sufl[i]+=flagb;
		
		suflm[i] = sufl[i]*(nb[i].S + nb[i].F)/2 ;
		
		if(i!=m-1) 
		{
			sufl[i] += sufl[i+1];
			suflm[i] += suflm[i+1];
		}
	}

	long double ans = 0;
	for(int i=0,j=0;i<n;i++)
	{
		while(j<m && nb[j].F < na[i].F)
		{
			j++;
		}
		
		if(j>=m) break;
		auto [l,r] = na[i];
		if(j<m && nb[j].F == na[i].F )
		{
			if(nb[j].S == na[i].S) ans += ((r*r*r - l*l*l)/3.0 + (long double)r*l*l - (long double)l*r*r);
			j++;
//			cout<<ans<<endl;
		}

		if(j<m) ans+= (r-l+flaga)*(suflm[j] - sufl[j]*(r+l)/2.0);
//		cout<<i<<' '<<j<<endl;
//		cout<<ans<<endl;
	}
	
	for(int i=0;i<m;i++)
	{
		prel[i] = nb[i].S - nb[i].F;
		prel[i] += flagb;
		prelm[i] = prel[i]*(nb[i].S + nb[i].F)/2 ;
		if(i!=0) 
		{
			prel[i] += prel[i-1];
			prelm[i] += prelm[i-1];
		}
	}
//	ans = 0;
//	cout<<endl;
	for(int i=0,j=0;i<n;i++)
	{
		while(j<m && nb[j].S < na[i].F)
		{
			j++;
		}
		j--;
		auto [l,r] = na[i];
//		cout<<i<<' '<<j<<endl;
		if(j>=0) ans+= -(r-l+flaga)*(prelm[j] - prel[j]*(r+l)/2.0);
//		cout<<ans<<endl;
		if(j==-1) j=0;
	}
	long double suma = 0;
	long double sumb = 0;
	for(int i=0;i<n;i++)
	{
		suma += na[i].S - na[i].F + flaga;
	}
	for(int i=0;i<m;i++)
	{
		sumb += nb[i].S - nb[i].F + flagb;
	}
//		cout<<endl;
//	for(auto [l,r]:na) cout<<l<<' '<<r<<endl;
//	cout<<endl;	
//	cout<<endl;
//	for(auto [l,r]:nb) cout<<l<<' '<<r<<endl;
//	cout<<endl;	
	ans/=(suma*sumb);
	cout<<fixed<<setprecision(11)<<ans;
}
signed main()
{
 	IOS;
	int T = 1;
//	cin>>T;
	while(T--)
	{
		solve();
	}
}