Mad scientist Mike does not use slow hard disks. His modification of a hard drive has not one, but n different heads that can read data in parallel.

When viewed from the side, Mike's hard drive is an endless array of tracks. The tracks of the array are numbered from left to right with integers, starting with 1. In the initial state the i-th reading head is above the track number h_{i}. For each of the reading heads, the hard drive's firmware can move the head exactly one track to the right or to the left, or leave it on the current track. During the operation each head's movement does not affect the movement of the other heads: the heads can change their relative order; there can be multiple reading heads above any of the tracks. A track is considered read if at least one head has visited this track. In particular, all of the tracks numbered h_1, h_2, ..., h_{n} have been read at the beginning of the operation. [Image] 

Mike needs to read the data on m distinct tracks with numbers p_1, p_2, ..., p_{m}. Determine the minimum time the hard drive firmware needs to move the heads and read all the given tracks. Note that an arbitrary number of other tracks can also be read.


-----Input-----

The first line of the input contains two space-separated integers n, m (1 ≤ n, m ≤ 10^5) — the number of disk heads and the number of tracks to read, accordingly. The second line contains n distinct integers h_{i} in ascending order (1 ≤ h_{i} ≤ 10^10, h_{i} < h_{i} + 1) — the initial positions of the heads. The third line contains m distinct integers p_{i} in ascending order (1 ≤ p_{i} ≤ 10^10, p_{i} < p_{i} + 1) - the numbers of tracks to read.

Please, do not use the %lld specifier to read or write 64-bit integers in С++. It is recommended to use the cin, cout streams or the %I64d specifier.


-----Output-----

Print a single number — the minimum time required, in seconds, to read all the needed tracks.


-----Examples-----
Input
3 4
2 5 6
1 3 6 8

Output
2

Input
3 3
1 2 3
1 2 3

Output
0

Input
1 2
165
142 200

Output
81



-----Note-----

The first test coincides with the figure. In this case the given tracks can be read in 2 seconds in the following way:   during the first second move the 1-st head to the left and let it stay there;  move the second head to the left twice;  move the third head to the right twice (note that the 6-th track has already been read at the beginning). 

One cannot read the tracks in 1 second as the 3-rd head is at distance 2 from the 8-th track.

Difficulty: competition

Starter code:

Solutions:
import sys
from itertools import *
from math import *
def solve():
    n, m = map(int, input().split())
    h = list(map(int, input().split()))
    p = list(map(int, input().split()))
    ss, ll = 0, int(2.2e10)
    while ss < ll:
        avg = (ss + ll) // 2
        works = True
        hidx = 0
        pidx = 0
        while hidx < len(h) and pidx < len(p):
            leftget = p[pidx]
            curpos = h[hidx]
            if curpos - leftget > avg:
                works = False
                break
            getbacktime = max(0, 2*(curpos - leftget))
            alsotoright = max(0, avg - getbacktime)
            leftime = max(0, curpos - leftget)
            remtime = max(0, (avg - leftime) // 2)
            furthestright = curpos + max(alsotoright, remtime)
            while pidx < len(p) and p[pidx] <= furthestright: pidx += 1
            hidx += 1
        if pidx != len(p): works = False
        if works: ll = avg
        else: ss = avg + 1
    print(ss)


if sys.hexversion == 50594544 : sys.stdin = open("test.txt")
solve()

import sys
from itertools import *
from math import *
def solve():
    n, m = map(int, input().split())
    h = list(map(int, input().split()))
    p = list(map(int, input().split()))
    ss, ll = 0, int(2.1e10)
    while ss < ll:
        avg = (ss + ll) // 2
        works = True
        hidx = 0
        pidx = 0
        while hidx < len(h) and pidx < len(p):
            leftget = p[pidx]
            curpos = h[hidx]
            if curpos - leftget > avg:
                works = False
                break
            getbacktime = max(0, 2*(curpos - leftget))
            alsotoright = max(0, avg - getbacktime)
            leftime = max(0, curpos - leftget)
            remtime = max(0, (avg - leftime) // 2)
            furthestright = curpos + max(alsotoright, remtime)
            while pidx < len(p) and p[pidx] <= furthestright: pidx += 1
            hidx += 1
        if pidx != len(p): works = False
        if works: ll = avg
        else: ss = avg + 1
    print(ss)


if sys.hexversion == 50594544 : sys.stdin = open("test.txt")
solve()

import sys
from itertools import *
from math import *
def solve():
    n, m = map(int, input().split())
    h = list(map(int, input().split()))
    p = list(map(int, input().split()))
    ss, ll = 0, int(2.1e10)
    while ss < ll:
        avg = (ss + ll) // 2
        works = True
        hidx = 0
        pidx = 0
        while hidx < len(h) and pidx < len(p):
            leftget = p[pidx]
            curpos = h[hidx]
            if curpos - leftget > avg:
                works = False
                break
            getbacktime = max(0, 2*(curpos - leftget))
            alsotoright = max(0, avg - getbacktime)
            leftime = max(0, curpos - leftget)
            remtime = max(0, (avg - leftime) // 2)
            furthestright = curpos + max(alsotoright, remtime)
            while pidx < len(p) and p[pidx] <= furthestright: pidx += 1
            hidx += 1
        if pidx != len(p): works = False
        if works: ll = avg
        else: ss = avg + 1
    print(ss)


if sys.hexversion == 50594544 : sys.stdin = open("test.txt")
solve()

def can(d,a,b):
    d1=d
    mi=a[-1]
    ma=a[-1]
    while len(a)>0 and len(b)>0:
        if b[-1]<=mi:
            if abs(b[-1]-ma)<=d1:
                a.pop()
                if len(a)==0:
                    break
                ma=a[-1]
            else:
                b.pop()
                mi=a[-1]
                ma=a[-1]
        elif b[-1]>=ma:
            if abs(b[-1]-mi)<=d1:
                a.pop()
                if len(a)==0:
                    break
                ma=a[-1]
            else:
                b.pop()
                mi=a[-1]
                ma=a[-1]
        else:
            if abs(ma-mi)+min(abs(b[-1]-mi),abs(b[-1]-ma))<=d1:
                a.pop()
                if len(a)==0:
                    break
                ma=a[-1]
            else:
                b.pop()
                mi=a[-1]
                ma=a[-1]
    return len(a)==0  
n,m=list(map(int,input().split()))
s=list(map(int,input().split()))[::-1]
s1=list(map(int,input().split()))[::-1]
high=(10**12)+1
low=0
while high-low>1:
    mid=(high+low)//2
    if can(mid,s1.copy(),s.copy()):
        high=mid
    else:
        low=mid
if can(low,s1,s):
    print(low)
else:
    print(high)
    


def can(d,a,b):
    d1=d
    mi=a[-1]
    ma=a[-1]
    while len(a)>0 and len(b)>0:
        if b[-1]<=mi:
            if abs(b[-1]-ma)<=d1:
                a.pop()
                if len(a)==0:
                    break
                ma=a[-1]
            else:
                b.pop()
                mi=a[-1]
                ma=a[-1]
        elif b[-1]>=ma:
            if abs(b[-1]-mi)<=d1:
                a.pop()
                if len(a)==0:
                    break
                ma=a[-1]
            else:
                b.pop()
                mi=a[-1]
                ma=a[-1]
        else:
            if abs(ma-mi)+min(abs(b[-1]-mi),abs(b[-1]-ma))<=d1:
                a.pop()
                if len(a)==0:
                    break
                ma=a[-1]
            else:
                b.pop()
                mi=a[-1]
                ma=a[-1]
    return len(a)==0  
n,m=list(map(int,input().split()))
s=list(map(int,input().split()))[::-1]
s1=list(map(int,input().split()))[::-1]
high=(10**11)+1
low=0
while high-low>1:
    mid=(high+low)//2
    if can(mid,s1.copy(),s.copy()):
        high=mid
    else:
        low=mid
if can(low,s1,s):
    print(low)
else:
    print(high)
    


def can(d,a,b):
    d1=d
    mi=a[-1]
    ma=a[-1]
    x=len(a)-1
    y=len(b)-1
    while x>=0 and y>=0:
        if b[y]<=mi:
            if abs(b[y]-ma)<=d1:
                x-=1
                if x==-1:
                    break
                ma=a[x]
            else:
                y-=1
                mi=a[x]
                ma=a[x]
        elif b[y]>=ma:
            if abs(b[y]-mi)<=d1:
                x-=1
                if x==-1:
                    break
                ma=a[x]
            else:
                y-=1
                mi=a[x]
                ma=a[x]
        else:
            if abs(ma-mi)+min(abs(b[y]-mi),abs(b[y]-ma))<=d1:
                x-=1
                if x==-1:
                    break
                ma=a[x]
            else:
                y-=1
                mi=a[x]
                ma=a[x]
    return x==-1 
n,m=list(map(int,input().split()))
s=list(map(int,input().split()))[::-1]
s1=list(map(int,input().split()))[::-1]
high=(10**10)*3
low=0
while high-low>1:
    mid=(high+low)//2
    if can(mid,s1,s):
        high=mid
    else:
        low=mid
if can(low,s1,s):
    print(low)
else:
    print(high)
    


import sys
from itertools import *
from math import *
def solve():
    n, m = list(map(int, input().split()))
    h = list(map(int, input().split()))
    p = list(map(int, input().split()))
    ss, ll = 0, int(2.1e10)
    while ss < ll:
        avg = (ss + ll) // 2
        works = True
        hidx = 0
        pidx = 0
        while hidx < len(h) and pidx < len(p):
            leftget = p[pidx]
            curpos = h[hidx]
            if curpos - leftget > avg:
                works = False
                break
            getbacktime = max(0, 2*(curpos - leftget))
            alsotoright = max(0, avg - getbacktime)
            leftime = max(0, curpos - leftget)
            remtime = max(0, (avg - leftime) // 2)
            furthestright = curpos + max(alsotoright, remtime)
            while pidx < len(p) and p[pidx] <= furthestright:
                pidx += 1
            hidx += 1
        if pidx != len(p):
            works = False
        if works:
            ll = avg
        else:
            ss = avg + 1
    print(ss)
 
 

solve()