#include <cstdio>
#include <map>
#include <vector>
#include <algorithm>

using namespace std;

using ll = long long int;

int N;
const int MAXN=205;
int A[MAXN];
bool dp[MAXN][MAXN][MAXN];
int  ans[MAXN][MAXN][MAXN];
bool dp2[MAXN][MAXN];
int ans2[MAXN][MAXN];

int main(){

    scanf("%d",&N);
    for(int i=1;i<=N;i++)
        scanf("%d",A+i);

    dp2[0][0]=true;
    for(int i=0;i<N;i++){
        for(int j=0;j<200;j++){
            if (!dp2[i][j]) continue;
            int v=A[i+1]%200;
            dp2[i+1][j]=true;
            dp2[i+1][(j+v)%200]=true;
            ans2[i+1][(j+v)%200]=1;
        }
    }

    vector<int> x,y;
    int curv=0;
    for(int i=N;i>0;i--){
        
        if (ans2[i][curv]==1){
            x.push_back(i);
            curv=(curv+200-(A[i]%200))%200;
        }else x.push_back(i),y.push_back(i);
    }

    if (x.size()!=y.size()&&y.size()!=0&&x.size()!=0){
        puts("Yes");
        printf("%d ",(int)x.size());
        while(!x.empty()) printf("%d ",x.back()),x.pop_back();
        putchar('\n');
        printf("%d ",(int)y.size());
        while(!y.empty()) printf("%d ",y.back()),y.pop_back();
        putchar('\n');
        return 0;
    }

    dp[0][0][0]=true;
    for(int i=0;i<N;i++){
        for(int j=0;j<200;j++){
            for(int k=0;k<200;k++){
                const auto& cur=dp[i][j][k];
                if (!cur) continue;
                int v=A[i+1]%200;

                if (!dp[i+1][j][k]){
                    dp[i+1][j][k]=true;
                    ans[i+1][j][k]=0;
                    // printf("%d %d\n",j,k);
                }
                if (!dp[i+1][(j+v)%200][k]||ans[i][j][k]!=1||ans[i+1][(j+v)%200][k]==0){
                    dp[i+1][(j+v)%200][k]=true;
                    ans[i+1][(j+v)%200][k]=1;
                    // printf("%d %d\n",(j+v)%200,k);
                }
                if (!dp[i+1][j][(k+v)%200]||ans[i][j][k]!=2||ans[i+1][j][(k+v)%200]==0){
                    dp[i+1][j][(k+v)%200]=true;
                    ans[i+1][j][(k+v)%200]=2;
                    // printf("%d %d\n",j,(k+v)%200);
                }
            }
        }
    }

    vector<int> X,Y;
    int ANS=0;
    for(int i=199;i>=0;i--){
        if (dp[N][i][i]){
            // printf("i=%d\n",i);
            int a=i,b=i;
            ANS=i;

            for(int i=N;i>0;i--){
                // printf("a=%d b=%d\n",a,b);
                // printf("ans=%d\n",ans[i][a][b]);
                int v=A[i]%200;
                if (ans[i][a][b]==1){
                    X.push_back(i);
                    a=(a+200-v)%200;
                }else if (ans[i][a][b]==2){
                    Y.push_back(i);
                    b=(b+200-v)%200;
                }
            }

            break;
        }
    }

    // printf("%d %d\n",(int)X.size(),(int)Y.size());
    if (X.size()!=0&&Y.size()!=0){
        puts("Yes");
        printf("%d ",(int)X.size());
        while(!X.empty()) printf("%d ",X.back()),X.pop_back();
        putchar('\n');
        printf("%d ",(int)Y.size());
        while(!Y.empty()) printf("%d ",Y.back()),Y.pop_back();
        putchar('\n');
    }
    // }else if (ANS==0&&(X.size()+Y.size())!=N&&(X.size()+Y.size())!=0){
    //     puts("Yes");

    //     for(int i=1;i<=N;i++) 
    //         if (find(X.begin(),X.end(),i)==X.end()&&find(Y.begin(),Y.end(),i)==Y.end()){
    //             X.push_back(i),Y.push_back(i);
    //         }
    //     sort(X.begin(),X.end());
    //     sort(Y.begin(),Y.end());
    //     printf("%d ",(int)X.size());
    //     while(!X.empty()) printf("%d ",X.back()),X.pop_back();
    //     putchar('\n');
    //     printf("%d ",(int)Y.size());
    //     while(!Y.empty()) printf("%d ",Y.back()),Y.pop_back();
    //     putchar('\n');
    // }
    else puts("No");
    

    return 0;
}
#pragma region head
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using vi = vector<int>;
using vll = vector<ll>;
using pi = pair<int, int>;
using pll = pair<ll, ll>;
template <class T>
using vv = vector<vector<T>>;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define repi(i, a, b) for (int i = (int)(a); i < (int)(b); i++)
#define rrep(i, n) for (int i = (int)(n)-1; i >= 0; i--)
#define rrepi(i, a, b) for (int i = (int)(b)-1; i >= (int)(a); i--)
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
#define bit(n) (1LL << (n))
template <class T>
inline bool chmax(T &a, const T &b) {
    if (a < b) {
        a = b;
        return 1;
    }
    return 0;
}
template <class T>
inline bool chmin(T &a, const T &b) {
    if (a > b) {
        a = b;
        return 1;
    }
    return 0;
}
const int INF = 1002003004;
const ll LINF = 1002003004005006007ll;
struct preprocess {
    preprocess() {
        cin.tie(nullptr);
        ios::sync_with_stdio(false);
        cout << fixed << setprecision(20);
    }
} ____;
#pragma endregion head

#pragma region library
#pragma endregion library

int main() {
    int n;
    cin >> n;
    string s, t;
    cin >> s >> t;
    ll ans = 0;
    int er = -1;
    queue<int> wait;
    int cnt = 0;
    rep(i, n) {
        if (s[i] == '1') cnt++;
        if (t[i] == '1') cnt--;
    }
    if (cnt < 0 || cnt % 2) {
        cout << -1 << '\n';
        return 0;
    }
    rep(i, n) {
        if (er != -1) {
            if (t[i] == '1') {
                wait.push(i);
            }
            if (s[i] == '1') {
                ans += i - er;
                er = -1;
            }
        } else {
            if (t[i] == '1') {
                wait.push(i);
            }
            if (s[i] == '1') {
                if (wait.empty()) {
                    er = i;
                } else {
                    ans += i - wait.front();
                    wait.pop();
                }
            }
        }
    }
    if (er != -1 || !wait.empty()) {
        cout << -1 << '\n';

    } else {
        cout << ans << '\n';
    }
}
