#include<bits/stdc++.h>
using namespace std;
using db=long double;
const db eps=1e-22;
const int B=1200;
int sgn(db x)
{
	if(x<-eps) return -1;
	else if(x>eps) return 1;
	else return 0;
}
struct Vec
{
	db x,y;
	Vec(db a=0,db b=0) {x=a,y=b;}
	db norm() {return x*x+y*y;}
	db abs() {return sqrt(norm());}
};
Vec operator + (Vec x,Vec y) {return Vec(x.x+y.x,x.y+y.y);}
Vec operator - (Vec x,Vec y) {return Vec(x.x-y.x,x.y-y.y);}
Vec operator * (Vec x,db y) {return Vec(x.x*y,x.y*y);}
db cdot(Vec x,Vec y) {return x.x*y.x+x.y*y.y;}
db cross(Vec x,Vec y) {return x.x*y.y-x.y*y.x;}
int ccw(Vec x,Vec y,Vec z) {return sgn(cross(y-x,z-x));}
vector<Vec> convex(vector<Vec> a)
{
	int n=a.size();
	int mnx=0;
	for(int i=0;i<n;i++) if(sgn(a[i].x-a[mnx].x)==-1||(sgn(a[i].x-a[mnx].x)==0&&sgn(a[i].y-a[mnx].y)==-1)) mnx=i;
	swap(a[mnx],a[0]);
	sort(a.begin()+1,a.end(),[&](Vec x,Vec y){
		return ccw(a[0],x,y)==-1||(ccw(a[0],x,y)==0&&(x-a[0]).norm()<(y-a[0]).norm());
	});
	vector<int> st; int tp=0;
	st.push_back(0),tp++;
	for(int i=1;i<n;i++)
	{
		if(sgn(a[i].x-a[i-1].x)==0&&sgn(a[i].y-a[i-1].y)==0) continue;
		while(tp>=2&&ccw(a[st[tp-2]],a[st[tp-1]],a[i])>=0) tp--,st.pop_back();
		st.push_back(i),tp++;
	}
	vector<Vec> ans;
	for(int i=tp-1;i>=0;i--) ans.push_back(a[st[i]]);
	return ans;
}
struct Seg
{
	Vec s,t;
	Seg() {}
	Seg(Vec a,Vec b) {s=a,t=b;}
	Vec point(db p) {return s+(t-s)*p;}
};
Vec proj(Seg l,Vec x)
{
	db p=cdot(l.t-l.s,x-l.s)/(l.t-l.s).norm();
	return l.point(p);
}
#define N 100005
int n,m;
Vec h[N],d[N];
vector<Vec> bads;
vector<Vec> v;
int findmin(const vector<Vec> &a,auto cmp)
{
	int l=0,r=a.size()-1,d=1;
	if(cmp(a.back(),a[0])) swap(l,r),d=-1;
	for(;(r-l)*d>1;)
	{
		int mid=(l+r)/2;
		if(cmp(a[mid],a[mid-d])&&cmp(a[mid],a[l])) l=mid;
		else r=mid;
	}
	return l;
}
void rebuild()
{
	for(auto i:bads) v.push_back(i);
	v=convex(v);
	// for(auto [x,y]:v) printf("%Lf %Lf\n",x,y);
	bads.clear();
}
signed main()
{
#ifdef xay5421
	freopen("b.in","r",stdin);
#endif
std::fesetround(FE_TOWARDZERO);
	cin>>n>>m;
	for(int i=1;i<=n;i++) scanf("%Lf %Lf",&h[i].x,&h[i].y),scanf("%Lf %Lf",&d[i].x,&d[i].y);
	v.resize(m);
	for(int i=0;i<m;i++) scanf("%Lf %Lf",&v[i].x,&v[i].y);
	v=convex(v);
	db ans=0;
	for(int i=n;i>=1;i--)
	{
		if(i%B==0) rebuild();
		auto cmp1=[&](Vec x,Vec y){
			return cdot(d[i],x-h[i])<cdot(d[i],y-h[i]);
		};
		Vec L=h[i],R=h[i];
		int id1=findmin(v,cmp1);
		if(cmp1(v[id1],L)) L=v[id1];
		for(auto i:bads) if(cmp1(i,L)) L=i;
		auto cmp2=[&](Vec x,Vec y){
			return cdot(d[i],x-h[i])>cdot(d[i],y-h[i]);
		};
		int id2=findmin(v,cmp2);
		if(cmp2(v[id2],R)) R=v[id2];
		for(auto i:bads) if(cmp2(i,R)) R=i;
		Vec A=proj(Seg(h[i],h[i]+d[i]),L);
		Vec B=proj(Seg(h[i],h[i]+d[i]),R);
		ans+=(B-A).abs();
		// printf("%d : %Lf %Lf %Lf %Lf\n",i,A.x,A.y,B.x,B.y);
		bads.push_back(A);
		bads.push_back(B);
	}
	printf("%.11Lf\n",ans);
	return 0;
}