Problem: There is a string of length n+1 of characters 'A' and 'B'. The first character and last character of the string are equal to 'A'.

You are given m indices p_1, p_2, …, p_m (0-indexation) denoting the other indices of characters 'A' in the string.

Let's denote the minimum distance between two neighboring 'A' as l, and maximum distance between neighboring 'A' as r.

For example, (l,r) of string "ABBAABBBA" is (1,4).

And let's denote the balance degree of a string as the value of r-l.

Now Dreamoon wants to change exactly k characters from 'B' to 'A', and he wants to make the balance degree of the string as small as possible.

Please calculate the required minimum possible value of balance degree.

Input

The first line contains one integer t denoting the number of test cases (1 ≤ t ≤ 400 000).

For each test case, the first line contains three integers n, m and k (1 ≤ n ≤ 10^{15}, 0 ≤ m ≤ 400 000, 0 ≤ k < n - m).

The second line contains m integers p_1, p_2, …, p_m, (0 < p_1 < p_2 < … < p_m < n).

The total sum of m is at most 400 000.

Output

For each test case, print one integer: the smallest possible value of balance degree after k changes of 'B' to 'A'.

Example

Input


5
80 3 5
11 24 50
81 7 12
4 10 17 26 37 48 61
25 10 14
3 4 7 12 13 15 17 19 21 23
1 0 0

10 2 0
2 4


Output


5
2
0
0
4

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const long long N = 4e5 + 7;
const long long INF = 1e15 + 7;
long long T, n, m, k;
long long a[N], b[N], u[N];
long long check1(long long v, long long t) {
  long long ans = 0;
  for (long long i = 1; i <= m; i++) {
    if (i >= t) v++;
    if (v % 2 == 0)
      ans += a[i] / (v / 2);
    else
      ans += (a[i] - 1) / ((v - 1) / 2);
    if (i >= t) v--;
  }
  return ans;
}
long long check2(long long v, long long t) {
  long long ans = 0;
  for (long long i = 1; i <= m; i++) {
    if (i <= t) v--;
    if (v % 2 == 0)
      ans += (a[i] - (a[i] + (v / 2) - 1) / (v / 2) + 1);
    else
      ans += (a[i] - (a[i] + ((v + 1) / 2)) / ((v + 1) / 2) + 1);
    if (i <= t) v++;
  }
  return ans;
}
void cb1(long long x) {
  long long sm = 0;
  for (long long i = 1; i <= m; i++) {
    b[i] = ((x + 1) & 1) ? ((a[i] - 1) / (x / 2)) : (a[i] / ((x + 1) / 2));
    u[i] = ((x & 1) ? ((a[i] - 1) / ((x - 1) / 2)) : (a[i] / (x / 2))) - b[i];
    sm += b[i];
  }
  for (long long i = 1; i <= m; i++) {
    if (sm + u[i] >= k + m) {
      b[i] += (k + m - sm);
      return;
    }
    b[i] += u[i];
    sm += u[i];
  }
}
void cb2(long long x) {
  long long sm = 0;
  for (long long i = 1; i <= m; i++) {
    long long x1 = x / 2, x2 = (x + 1) / 2;
    b[i] = ((x - 1) & 1) ? ((a[i] + x1) / x1) : ((a[i] + x1 - 1) / x1);
    u[i] = (x & 1) ? ((a[i] + x2) / x2) : ((a[i] + x2 - 1) / x2);
    b[i] = (a[i] - b[i] + 1);
    u[i] = (a[i] - u[i] + 1);
    u[i] -= b[i];
    sm += b[i];
  }
  for (long long i = 1; i <= m; i++) {
    if (sm + u[i] >= n - k) {
      b[i] += (n - k - sm);
      break;
    }
    b[i] += u[i];
    sm += u[i];
  }
  for (long long i = 1; i <= m; i++) b[i] = (a[i] - b[i] + 1);
}
struct pir {
  long long x, v;
} p[N];
bool cmp1(pir x, pir y) {
  if (!x.v) {
    if (!y.v) return x.x < y.x;
    return 1;
  }
  if (!y.v) return 0;
  if ((a[x.x] / x.v) != (a[y.x] / y.v)) return (a[x.x] / x.v) > (a[y.x] / y.v);
  if (a[x.x] % x.v == 0) {
    if (a[y.x] % y.v == 0) return x.x < y.x;
    return 0;
  }
  if (a[y.x] % y.v == 0) return 1;
  return x.x < y.x;
}
bool cmp2(pir x, pir y) {
  if ((a[x.x] / x.v) != (a[y.x] / y.v)) return (a[x.x] / x.v) < (a[y.x] / y.v);
  if (a[x.x] % x.v == 0) {
    if (a[y.x] % y.v == 0) return x.x > y.x;
    return 1;
  }
  if (a[y.x] % y.v == 0) return 0;
  return x.x > y.x;
}
int main() {
  scanf("%lld", &T);
  for (long long gg = 1; gg <= T; gg++) {
    scanf("%lld%lld%lld", &n, &m, &k);
    long long ls = 0, mx = 0, ans = INF;
    for (long long i = 1; i <= m; i++) {
      long long x;
      scanf("%lld", &x);
      a[i] = x - ls;
      ls = x;
    }
    a[++m] = n - ls;
    for (long long i = 1; i <= m; i++) mx = max(mx, a[i]);
    long long l = 2, r = mx * 2;
    while (l < r - 1) {
      long long mid = (l + r) >> 1;
      if (check1(mid, m + 1) < k + m)
        r = mid;
      else
        l = mid;
    }
    long long e = (check1(r, m + 1) >= k + m) ? r : l;
    cb1(e);
    for (long long i = 1; i <= m; i++) p[i] = (pir){i, b[i]};
    sort(p + 1, p + m + 1, cmp1);
    long long L = 1, R = m;
    while (L < R - 1) {
      long long Mid = (L + R) >> 1;
      if (!p[Mid].v)
        L = Mid;
      else {
        if (check1(a[p[Mid].x] / p[Mid].v * 2 + (a[p[Mid].x] % p[Mid].v != 0),
                   p[Mid].x) <= k)
          L = Mid;
        else
          R = Mid;
      }
    }
    if (p[R].v)
      if (check1(a[p[R].x] / p[R].v * 2 + (a[p[R].x] % p[R].v != 0), p[R].x) <=
          k)
        L = R;
    for (long long i = 1; i <= L; i++) {
      if (p[i].v)
        ans =
            min(ans, (a[p[i].x] / p[i].v + (a[p[i].x] % p[i].v != 0)) - e / 2);
      if (p[i].v == a[p[i].x]) break;
      e = min(e, ((a[p[i].x] % (p[i].v + 1) == 0)
                      ? (a[p[i].x] / (p[i].v + 1) * 2)
                      : ((a[p[i].x] / (p[i].v + 1)) * 2 + 1)));
    }
    l = 2, r = mx * 2;
    while (l < r - 1) {
      long long mid = (l + r) >> 1;
      if (check2(mid, 0) < n - k)
        l = mid;
      else
        r = mid;
    }
    e = (check2(l, 0) >= n - k) ? l : r;
    for (long long i = 1; i <= m; i++) b[i]++;
    for (long long i = 1; i <= m; i++) p[i] = (pir){i, b[i]};
    sort(p + 1, p + m + 1, cmp2);
    L = 1, R = m;
    while (L < R - 1) {
      long long Mid = (L + R) >> 1;
      if (p[Mid].v > a[p[Mid].x])
        L = Mid;
      else {
        if (check2(a[p[Mid].x] / p[Mid].v * 2 + (a[p[Mid].x] % p[Mid].v != 0),
                   p[Mid].x) <= n - k - m)
          L = Mid;
        else
          R = Mid;
      }
    }
    if (p[R].v <= a[p[R].x])
      if (check2(a[p[R].x] / p[R].v * 2 + (a[p[R].x] % p[R].v != 0), p[R].x) <=
          n - k - m)
        L = R;
    for (long long i = 1; i <= L; i++) {
      if (p[i].v <= a[p[i].x]) ans = min(ans, (e + 1) / 2 - a[p[i].x] / p[i].v);
      if (p[i].v == 1) break;
      e = max(e, ((a[p[i].x] % (p[i].v - 1) == 0)
                      ? (a[p[i].x] / (p[i].v - 1) * 2)
                      : ((a[p[i].x] / (p[i].v - 1)) * 2 + 1)));
    }
    printf("%lld\n", ans);
  }
  return 0;
}