#include<bits/stdc++.h>
#define rep(i,a,b) for (register int i=(a);i<=(b);i++)
#define drep(i,a,b) for (register int i=(a);i>=(b);i--)
typedef long long ll;
using namespace std;

typedef long double db;
db eps=3e-11;

const int N=2000010;

int n,m,sa[N],sb[N],cnt;
db a[N][2],b[N][2],dis[N]; 

int main()
{
	ios::sync_with_stdio(false); cin.tie(0);
	
	cin>>n>>m;
	rep(i,1,n) 
	{
		rep(j,0,1) cin>>a[i][j];
		if (fabs(a[i][0]-a[i][1])<eps) a[i][1]+=eps;
		rep(j,0,1) dis[++cnt]=a[i][j];
	}
	rep(i,1,n) 
	{
		rep(j,0,1) cin>>b[i][j];
		if (fabs(b[i][0]-b[i][1])<eps) b[i][1]+=eps;
		rep(j,0,1) dis[++cnt]=b[i][j];
	}
	sort(dis+1,dis+cnt+1);
	int sz=unique(dis+1,dis+cnt+1)-dis-1;
	
	rep(i,1,n)
	{
		int l=lower_bound(dis+1,dis+sz+1,a[i][0])-dis;
		int r=lower_bound(dis+1,dis+sz+1,a[i][1])-dis;
		sa[l]++;
		sa[r]--;
	}
	rep(i,1,m)
	{
		int l=lower_bound(dis+1,dis+sz+1,b[i][0])-dis;
		int r=lower_bound(dis+1,dis+sz+1,b[i][1])-dis;
		sb[l]++;
		sb[r]--;
	}
	rep(i,1,sz) sa[i]+=sa[i-1],sb[i]+=sb[i-1];
	
	db lena=0,lenb=0;
	rep(i,1,sz-1)
	{
		db d=dis[i+1]-dis[i];
		if (sa[i]) lena+=d;
		if (sb[i]) lenb+=d;
	}
	
	db a0=0,a1=0,b0=0,b1=0,ans=0;
	rep(i,1,sz-1)
	{
		db d=dis[i+1]-dis[i];
		db pa=d/lena,pb=d/lenb;
		db mid=(dis[i]+dis[i+1])/2;
		if (sa[i])
		{
			ans+=pa*mid*b0;
			ans-=pa*b1;
		}
		if (sb[i])
		{
			ans+=pb*mid*a0;
			ans-=pb*a1;
		}
		if (sa[i]&&sb[i])
		{
			ans+=pa*pb*d*d*d/3;
		}
		if (sa[i])
		{
			a0+=pa;
			a1+=pa*mid;
		}
		if (sb[i])
		{
			b0+=pb;
			b1+=pb*mid;
		}
	}
	
	printf("%.15Lf\n",ans); 
	
	return 0;
}