Problem: A championship is held in Berland, in which n players participate. The player with the number i has a_i (a_i ≥ 1) tokens.

The championship consists of n-1 games, which are played according to the following rules:

  * in each game, two random players with non-zero tokens are selected; 
  * the player with more tokens is considered the winner of the game (in case of a tie, the winner is chosen randomly); 
  * the winning player takes all of the loser's tokens; 



The last player with non-zero tokens is the winner of the championship.

All random decisions that are made during the championship are made equally probable and independently.

For example, if n=4, a = [1, 2, 4, 3], then one of the options for the game (there could be other options) is: 

  * during the first game, the first and fourth players were selected. The fourth player has more tokens, so he takes the first player's tokens. Now a = [0, 2, 4, 4]; 
  * during the second game, the fourth and third players were selected. They have the same number of tokens, but in a random way, the third player is the winner. Now a = [0, 2, 8, 0]; 
  * during the third game, the second and third players were selected. The third player has more tokens, so he takes the second player's tokens. Now a = [0, 0, 10, 0]; 
  * the third player is declared the winner of the championship. 



Championship winners will receive personalized prizes. Therefore, the judges want to know in advance which players have a chance of winning, i.e have a non-zero probability of winning the championship. You have been asked to find all such players. 

Input

The first line contains one integer t (1 ≤ t ≤ 10^4) — the number of test cases. Then t test cases follow.

The first line of each test case consists of one positive integer n (1 ≤ n ≤ 2 ⋅ 10^5) — the number of players in the championship.

The second line of each test case contains n positive integers a_1, a_2, …, a_n (1 ≤ a_i ≤ 10^9) — the number of tokens the players have.

It is guaranteed that the sum of n over all test cases does not exceed 2 ⋅ 10^5. 

Output

For each test case, print the number of players who have a nonzero probability of winning the championship. On the next line print the numbers of these players in increasing order. Players are numbered starting from one in the order in which they appear in the input. 

Example

Input


2
4
1 2 4 3
5
1 1 1 1 1


Output


3
2 3 4 
5
1 2 3 4 5 

Solution (in CPP):

Solution:
//by aman kumar
#include <cassert>
#include <cctype>
#include <cerrno>
#include <cfloat>
#include <ciso646>
#include <climits>
#include <clocale>
#include <cmath>
#include <csetjmp>
#include <csignal>
#include <cstdarg>
#include <cstddef>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <ccomplex>
#include <cfenv>
#include <cinttypes>
#include <cstdbool>
#include <cstdint>
#include <ctgmath>
#include <cwchar>
#include <cwctype>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <iterator>
#include <limits>
#include <list>
#include <locale>
#include <map>
#include <memory>
#include <new>
#include <numeric>
#include <ostream>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <typeinfo>
#include <utility>
#include <valarray>
#include <vector>
#include <array>
#include <atomic>
#include <chrono>
#include <condition_variable>
#include <forward_list>
#include <future>
#include <initializer_list>
#include <mutex>
#include <random>
#include <ratio>
#include <regex>
#include <scoped_allocator>
#include <system_error>
#include <thread>
#include <tuple>
#include <typeindex>
#include <type_traits>
#include <unordered_map>
#include <unordered_set>
using namespace std;
#define int           long long int
#define pb           push_back
#define pob          pop_back
#define pr           pair<int,int>
#define all(a)       a.begin(),a.end()
#define allr(a)      a.rbegin(),a.rend()
#define F            first
#define S            second
#define mst          multiset<int,int>
#define rep(i,a,b)   for(int i=a;i<b;i++)
#define pii          pair<int,int>
#define endl         "\n"
#define vi          vector<int>
#define elasped_time 1.0 * clock() / CLOCKS_PER_SEC
void optimise() {ios_base::sync_with_stdio(false); cin.tie(NULL);}
int mod = (int)1e9 + 7;
void normalize(int &a) {a = (a + mod) % mod;}
int modmul(int a, int b) {a = a % mod; b = b % mod; normalize(a); normalize(b); return (a * b) % mod;}
int modadd(int a, int b) {a = a % mod; b = b % mod; normalize(a); normalize(b); return (a + b) % mod;}
int modsub(int a, int b) {a = a % mod; b = b % mod; normalize(a); normalize(b); return (a - b + mod) % mod;}
int me(int x, int n) {x %= mod; if (n == 0)return 1; int u = me(x, n / 2) % mod; u = (u * u) % mod; if (n % 2)u = (u * x) % mod; return u;}
int me1(int x, int n ) {if (n == 0)return 1; int u = me1(x, n / 2); u = u * u; if (n % 2)u = u * x; return u;}
inline  int modInv(int a) { return me(a, mod - 2); }
inline  int modDiv(int a, int b) { return modmul(a, modInv(b)); }
float power(float x, int y)
{
  float temp;
  if (y == 0)
    return 1;
  temp = power(x, y / 2);
  if (y % 2 == 0)
    return temp * temp;
  else
  {
    if (y > 0)
      return x * temp * temp;
    else
      return (temp * temp) / x;
  }
}
int __gcd(int a, int b) {if (b == 0)return a; else return __gcd(b, a % b);}
pii operator+(pii a, pii b) {return {a.F + b.F, a.S + b.S};}
pii operator-(pii a, pii b) {return {a.F - b.F, a.S - b.S};}
pii operator*(pii a, pii b) {return {a.F * b.F, a.S * b.S};}
template <typename T>ostream &operator<<(ostream &os, const set<T>& v) {for (auto it : v) os << it << " "; return os;}
template <typename T>ostream &operator<<(ostream &os, const vector<T>& v) {for (int i = 0; i < v.size(); ++i)os << v[i] << " "; return os;}




void solve()
{

  int n;
  cin >> n;
  vector<pii>p(n);
  set<int>s;
  for (int i = 0; i < n; i++)
  {
    cin >> p[i].F;
    p[i].S = i;
    s.insert(i);
  }

  sort(p.begin(), p.end());
  int sum = 0;
  int j = 0;
  for (int i = 0; i < n; i++)
  {
    if (sum < p[i].F)
    {
      for (; j < i; j++)
      {
        s.erase(p[j].S);
      }
    }
    sum += p[i].F;

  }
  cout << s.size() << endl;
  for (auto y : s)
    cout << y + 1 << " ";





}

signed main()
{
  optimise();
#ifndef ONLINE_JUDGE
  freopen("input.txt",  "r",   stdin);
  freopen("output.txt", "w",   stdout);
  freopen("error.txt",  "w",   stderr);
#endif
  int t ;
  t = 1;
  cin >> t;

  cout << fixed << setprecision(3);

  while (t--)

  {
    solve();
    cout << endl;


  }
  cerr << endl << "Time Elasped : " << elasped_time << endl;
}