#include<bits/stdc++.h>
using namespace std;
int n,p[1010000],a[1010000];
void sol(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%d",&p[i]),a[i]=0;
    int s=0;
    vector<int>v1;
    for(int i=1;i<=n;i++)if(!a[i]){
        s++;
        int u=i;
        vector<int>vp;
        while(!a[u])vp.push_back(u),a[u]=1,u=p[u];
        if(s>1){
            printf("%d\n",vp.size());
            sort(vp.begin(),vp.end());
            for(int x:vp)printf("%d ",x);
            puts("");
            return;
        }
        v1=vp;
    }
    int e=-1,mi=n+1;
    for(int i=0;i<n;i++){
        int z=abs(v1[i]-v1[(i+1)%n]);
        if(mi>z)mi=z,e=i;
    }
    printf("%d\n",n-mi-1);
    puts("1"),printf("%d\n",v1[e]);
}
int main(){
    int T;scanf("%d",&T);while(T--)sol();
    return 0;
}