import sys

# sys.setrecursionlimit(200005)
# sys.set_int_max_str_digits(1000005)
int1 = lambda x: int(x)-1
pDB = lambda *x: print(*x, end="\n", file=sys.stderr)
p2D = lambda x: print(*x, sep="\n", end="\n\n", file=sys.stderr)
def II(): return int(sys.stdin.readline())
def LI(): return list(map(int, sys.stdin.readline().split()))
def LLI(rows_number): return [LI() for _ in range(rows_number)]
def LI1(): return list(map(int1, sys.stdin.readline().split()))
def LLI1(rows_number): return [LI1() for _ in range(rows_number)]
def SI(): return sys.stdin.readline().rstrip()

dij = [(0, 1), (-1, 0), (0, -1), (1, 0)]
# dij = [(0, 1), (-1, 0), (0, -1), (1, 0), (1, 1), (1, -1), (-1, 1), (-1, -1)]
inf = (1 << 63)-1
# md = 10**9+7
md = 998244353

def ccw(s, t, p):
    x1 = p[0]-s[0]
    y1 = p[1]-s[1]
    x2 = t[0]-s[0]
    y2 = t[1]-s[1]
    if x1*y2-x2*y1 > 0: return 1
    if x1*y2-x2*y1 < 0: return -1
    if x1*x2+y1*y2 < 0: return -2
    if x1**2+y1**2 > x2**2+y2**2: return 2
    return False

# 共有点の有無のみ
def common_point(a, b, s, t):
    cc = []
    cc.append(ccw(a, b, s))
    cc.append(ccw(a, b, t))
    cc.append(ccw(s, t, a))
    cc.append(ccw(s, t, b))
    if any(c == 0 for c in cc): return True
    if any(abs(c) == 2 for c in cc): return False
    if cc[0] == cc[1] or cc[2] == cc[3]: return False
    return True

co = 10**10

n, m = LI()
xy, st = [], []
for _ in range(n):
    x, y, s, t = LI()
    xy.append((x*co, y*co))
    st.append((s*co, t*co))

dist = [[inf]*2*n for _ in range(2*n)]

def can(a, b, i, j):
    for k in range(n):
        if k == i or k == j: continue
        if common_point(a, b, xy[k], st[k]): return False
    return True

for i in range(n):
    dist[i*2][i*2] = 0
    dist[i*2+1][i*2+1] = 0
    if can(xy[i], st[i], i, i):
        x, y = xy[i]
        s, t = st[i]
        dist[i*2][i*2+1] = dist[i*2+1][i*2] = round(((s-x)**2+(t-y)**2)**0.5)
    for j in range(i):
        if can(xy[i], xy[j], i, j):
            x, y = xy[i]
            s, t = xy[j]
            dist[i*2][j*2] = dist[j*2][i*2] = round(((s-x)**2+(t-y)**2)**0.5)
        if can(xy[i], st[j], i, j):
            x, y = xy[i]
            s, t = st[j]
            dist[i*2][j*2+1] = dist[j*2+1][i*2] = round(((s-x)**2+(t-y)**2)**0.5)
        if can(st[i], xy[j], i, j):
            x, y = st[i]
            s, t = xy[j]
            dist[i*2+1][j*2] = dist[j*2][i*2+1] = round(((s-x)**2+(t-y)**2)**0.5)
        if can(st[i], st[j], i, j):
            x, y = st[i]
            s, t = st[j]
            dist[i*2+1][j*2+1] = dist[j*2+1][i*2+1] = round(((s-x)**2+(t-y)**2)**0.5)
# p2D(dist)

for j in range(n*2):
    for i in range(n*2):
        ij = dist[i][j]
        if ij == inf: continue
        for k in range(n*2):
            jk = dist[j][k]
            if jk == inf: continue
            dist[i][k] = min(dist[i][k], ij+jk)

for _ in range(m):
    i, b, j, d = LI1()
    ans = dist[i*2+b][j*2+d]
    print(ans/co)