#include <bits/stdc++.h>

using namespace std;

const int N = 1e5 + 5;

int n;
string S, T;
vector<int> res;
set<int> z[2][2];

void add(int x)
{
    assert(S[x] == '1');
    res.push_back(x);
    z[S[(x + n - 1) % n] - '0'][T[(x + n - 1) % n] - '0'].erase((x + n - 1) % n);
    z[S[x] - '0'][T[x] - '0'].erase(x);
    z[S[(x + 1) % n] - '0'][T[(x + 1) % n] - '0'].erase((x + 1) % n);
    S[(x + n - 1) % n] ^= 1;
    S[x] ^= 1;
    S[(x + 1) % n] ^= 1;
    T[x] ^= 1;
    z[S[(x + n - 1) % n] - '0'][T[(x + n - 1) % n] - '0'].insert((x + n - 1) % n);
    z[S[x] - '0'][T[x] - '0'].insert(x);
    z[S[(x + 1) % n] - '0'][T[(x + 1) % n] - '0'].insert((x + 1) % n);
    return;
}

signed main()
{
    int tests;
    scanf("%d", &tests);
    while (tests--) {
        scanf("%d", &n);
        T = "";
        cin >> S;
        for (int i = 0; i < n; ++i) {
            S[i] ^= 1;
        }
        for (int i = 0; i < 2; ++i) {
            for (int j = 0; j < 2; ++j) {
                T += (char)(i + '0');
                T += (char)(j + '0');
                for (int k = 1; k < n - 1; ++k) {
                    T += (char)(T[k - 1] ^ T[k] ^ S[k]);
                }
                if ((T[0] ^ T[1] ^ T[n - 1]) == S[0] && (T[n - 2] ^ T[n - 1] ^ T[0]) == S[n - 1]) {
                    break;
                }
                T = "";
            }
            if (!T.empty()) {
                break;
            }
        }
        if (T.empty()) {
            puts("0");
            continue;
        }
        z[0][0].clear(), z[0][1].clear();
        z[1][0].clear(), z[1][1].clear();
        for (int i = 0; i < n; ++i) {
            z[S[i] - '0'][T[i] - '0'].insert(i);
        }
        res.clear();
        while (z[1][0].size() || z[1][1].size()) {
            if (z[1][1].size()) {
                add(*z[1][1].begin());
            } else {
                int x = *z[1][0].begin(), y, z;
                if (T[(x + 1) % n] == '1') {
                    y = (x + 1) % n, z = (y + 1) % n;
                } else {
                    y = (x + n - 1) % n, z = (y + n - 1) % n;
                }
                add(x), add(y), add(x), add(z);
            }
        }
        printf("%d\n", (int)res.size());
        for (int i = 0; i < (int)res.size(); ++i) {
            printf("%d%c", res[i] + 1, " \n"[i + 1 == (int)res.size()]);
        }
    }
    return 0;
}