Problem: There are n distinct points on a coordinate line, the coordinate of i-th point equals to x_i. Choose a subset of the given set of points such that the distance between each pair of points in a subset is an integral power of two. It is necessary to consider each pair of points, not only adjacent. Note that any subset containing one element satisfies the condition above. Among all these subsets, choose a subset with maximum possible size.

In other words, you have to choose the maximum possible number of points x_{i_1}, x_{i_2}, ..., x_{i_m} such that for each pair x_{i_j}, x_{i_k} it is true that |x_{i_j} - x_{i_k}| = 2^d where d is some non-negative integer number (not necessarily the same for each pair of points).

Input

The first line contains one integer n (1 ≤ n ≤ 2 ⋅ 10^5) — the number of points.

The second line contains n pairwise distinct integers x_1, x_2, ..., x_n (-10^9 ≤ x_i ≤ 10^9) — the coordinates of points.

Output

In the first line print m — the maximum possible number of points in a subset that satisfies the conditions described above.

In the second line print m integers — the coordinates of points in the subset you have chosen.

If there are multiple answers, print any of them.

Examples

Input

6
3 5 4 7 10 12


Output

3
7 3 5

Input

5
-1 2 5 8 11


Output

1
8

Note

In the first example the answer is [7, 3, 5]. Note, that |7-3|=4=2^2, |7-5|=2=2^1 and |3-5|=2=2^1. You can't find a subset having more points satisfying the required property.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int N = 2e5 + 7;
int n;
long long a[N], fn[33];
int main() {
  fn[0] = 1;
  for (int i = 1; i <= 32; i++) {
    fn[i] = fn[i - 1] * 2;
  }
  scanf("%d", &n);
  for (int i = 1; i <= n; i++) {
    scanf("%lld", &a[i]);
  }
  int f = 0;
  sort(a + 1, a + n + 1);
  a[n + 1] = LLONG_MAX;
  long long x, y, z;
  for (int i = 1; i <= n; i++) {
    for (int j = 0; j <= 32; j++) {
      long long cnt = a[i] + fn[j];
      long long cnt1 = a[i] + fn[j] + fn[j];
      int p = lower_bound(a + 1, a + n + 1, cnt) - a;
      int p1 = lower_bound(a + 1, a + n + 1, cnt1) - a;
      if (a[p] == cnt && a[p1] == cnt1) {
        f = 3;
        x = a[i];
        y = cnt;
        z = cnt1;
        break;
      } else if (a[p] == cnt) {
        f = 2;
        x = a[i];
        y = cnt;
      } else if (a[p1] == cnt1) {
        f = 2;
        x = a[i];
        y = cnt1;
      }
    }
    if (f == 3) break;
  }
  if (f == 3) {
    cout << f << endl;
    cout << x << " " << y << " " << z << endl;
  } else if (f == 2) {
    cout << f << endl;
    cout << x << " " << y << endl;
  } else {
    cout << 1 << endl;
    cout << a[1] << endl;
  }
}