Problem: Professor GukiZ has two arrays of integers, a and b. Professor wants to make the sum of the elements in the array a sa as close as possible to the sum of the elements in the array b sb. So he wants to minimize the value v = |sa - sb|.

In one operation professor can swap some element from the array a and some element from the array b. For example if the array a is [5, 1, 3, 2, 4] and the array b is [3, 3, 2] professor can swap the element 5 from the array a and the element 2 from the array b and get the new array a [2, 1, 3, 2, 4] and the new array b [3, 3, 5].

Professor doesn't want to make more than two swaps. Find the minimal value v and some sequence of no more than two swaps that will lead to the such value v. Professor makes swaps one by one, each new swap he makes with the new arrays a and b.

Input

The first line contains integer n (1 ≤ n ≤ 2000) — the number of elements in the array a.

The second line contains n integers ai ( - 109 ≤ ai ≤ 109) — the elements of the array a.

The third line contains integer m (1 ≤ m ≤ 2000) — the number of elements in the array b.

The fourth line contains m integers bj ( - 109 ≤ bj ≤ 109) — the elements of the array b.

Output

In the first line print the minimal value v = |sa - sb| that can be got with no more than two swaps.

The second line should contain the number of swaps k (0 ≤ k ≤ 2).

Each of the next k lines should contain two integers xp, yp (1 ≤ xp ≤ n, 1 ≤ yp ≤ m) — the index of the element in the array a and the index of the element in the array b in the p-th swap.

If there are several optimal solutions print any of them. Print the swaps in order the professor did them.

Examples

Input

5
5 4 3 2 1
4
1 1 1 1


Output

1
2
1 1
4 2


Input

5
1 2 3 4 5
1
15


Output

0
0


Input

5
1 2 3 4 5
4
1 2 3 4


Output

1
1
3 1

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const long double pi = 3.1415926535897932384626433832795l;
template <typename T>
inline auto sqr(T x) -> decltype(x * x) {
  return x * x;
}
template <typename T1, typename T2>
inline bool umx(T1& a, T2 b) {
  if (a < b) {
    a = b;
    return 1;
  }
  return 0;
}
template <typename T1, typename T2>
inline bool umn(T1& a, T2 b) {
  if (b < a) {
    a = b;
    return 1;
  }
  return 0;
}
const int N = 2010;
struct Input {
  int n, m;
  int a[N];
  int b[N];
  bool read() {
    if (scanf("%d", &n) != 1) return 0;
    for (int i = int(0); i < int(n); ++i) scanf("%d", &a[i]);
    if (scanf("%d", &m) != 1) return 0;
    for (int i = int(0); i < int(m); ++i) scanf("%d", &b[i]);
    return 1;
  }
  void init(const Input& input) { *this = input; }
};
struct Data : Input {
  long long res;
  vector<pair<int, int> > ans;
  void write() {
    cout << res << endl;
    cout << ((int)(ans).size()) << endl;
    for (int i = int(0); i < int(((int)(ans).size())); ++i) {
      cout << ans[i].first + 1 << ' ' << ans[i].second + 1 << endl;
    }
  }
  virtual void solve() {}
  virtual void clear() { *this = Data(); }
};
const int inf = 1e9;
struct Solution : Data {
  set<pair<long long, pair<int, int> > > s;
  long long suma, sumb;
  void upd_ans(set<pair<long long, pair<int, int> > >::iterator it, int i,
               int j) {
    if (it == s.end()) return;
    long long tek = a[i] + a[j];
    long long tekres = abs(suma - sumb - 2 * tek + 2 * it->first);
    if (tekres < res) {
      res = tekres;
      ans.clear();
      ans.push_back(make_pair(i, it->second.first));
      ans.push_back(make_pair(j, it->second.second));
    }
  }
  void solve() {
    suma = 0, sumb = 0;
    for (int i = int(0); i < int(n); ++i) suma += a[i];
    for (int i = int(0); i < int(m); ++i) sumb += b[i];
    res = abs(suma - sumb);
    ans.clear();
    for (int i = int(0); i < int(m); ++i) {
      for (int j = int(0); j < int(i); ++j) {
        s.insert(make_pair(b[j] + b[i], make_pair(j, i)));
      }
    }
    for (int i = int(0); i < int(n); ++i) {
      for (int j = int(0); j < int(i); ++j) {
        long long tek = a[i] + a[j];
        long long x = (sumb - suma) / 2 + tek;
        auto tekit = s.lower_bound(make_pair(x, make_pair(-1, -1)));
        auto nextit = s.lower_bound(make_pair(x, make_pair(inf, inf)));
        upd_ans(tekit, i, j);
        upd_ans(nextit, i, j);
        if (tekit != s.begin()) {
          tekit--;
          upd_ans(tekit, i, j);
        }
      }
    }
    for (int i = int(0); i < int(n); ++i)
      for (int j = int(0); j < int(m); ++j) {
        if (abs(suma - sumb - 2 * a[i] + 2 * b[j]) < res) {
          res = abs(suma - sumb - 2 * a[i] + 2 * b[j]);
          ans.clear();
          ans.push_back(make_pair(i, j));
        }
      }
  }
  void clear() { *this = Solution(); }
};
Solution sol;
int main() {
  cout.setf(ios::showpoint | ios::fixed);
  cout.precision(20);
  sol.read();
  sol.solve();
  sol.write();
  return 0;
}