#include <bits/stdc++.h>
using namespace std;

#define li long long int
#define rep(i, to) for (li i = 0; i < ((li)(to)); i++)
#define repp(i, start, to) for (li i = (li)(start); i < ((li)(to)); i++)
#define pb push_back
#define sz(v) ((li)(v).size())
#define bgn(v) ((v).begin())
#define eend(v) ((v).end())
#define allof(v) (v).begin(), (v).end()
#define dodp(v, n) memset(v, (li)n, sizeof(v))
#define bit(n) (1ll << (li)(n))
#define mp(a, b) make_pair(a, b)
#define rin rep(i, n)
#define EPS 1e-12
#define ETOL 1e-8
#define MOD 1000000007
typedef pair<li, li> PI;

#define INF bit(60)

#define DBGP 1

#define idp if (DBGP)
#define F first
#define S second
#define p2(a, b) idp cout << a << "\t" << b << endl
#define p3(a, b, c) idp cout << a << "\t" << b << "\t" << c << endl
#define p4(a, b, c, d) idp cout << a << "\t" << b << "\t" << c << "\t" << d << endl
#define p5(a, b, c, d, e) idp cout << a << "\t" << b << "\t" << c << "\t" << d << "\t" << e << endl
#define p6(a, b, c, d, e, f) idp cout << a << "\t" << b << "\t" << c << "\t" << d << "\t" << e << "\t" << f << endl
#define p7(a, b, c, d, e, f, g) idp cout << a << "\t" << b << "\t" << c << "\t" << d << "\t" << e << "\t" << f << "\t" << g << endl
#define p8(a, b, c, d, e, f, g, h) idp cout << a << "\t" << b << "\t" << c << "\t" << d << "\t" << e << "\t" << f << "\t" << g << "\t" << h << endl
#define p9(a, b, c, d, e, f, g, h, i) idp cout << a << "\t" << b << "\t" << c << "\t" << d << "\t" << e << "\t" << f << "\t" << g << "\t" << h << "\t" << i << endl
#define p10(a, b, c, d, e, f, g, h, i, j) idp cout << a << "\t" << b << "\t" << c << "\t" << d << "\t" << e << "\t" << f << "\t" << g << "\t" << h << "\t" << i << "\t" << j << endl
#define foreach(it, v) for (__typeof((v).begin()) it = (v).begin(); it != (v).end(); ++it)
#define p2p(x) idp p2((x).F, (x).S)
#define dump(x, n)                         \
    idp                                    \
    {                                      \
        rep(i, n) { cout << x[i] << " "; } \
        puts("");                          \
    }
#define dump2(x, n)                                                     \
    idp                                                                 \
    {                                                                   \
        rep(i, n) { cout << "[" << x[i].F << " , " << x[i].S << "] "; } \
        puts("");                                                       \
    }
#define dumpi(x)                  \
    idp                           \
    {                             \
        foreach (it, x)           \
        {                         \
            cout << (*it) << " "; \
        }                         \
        puts("");                 \
    }
#define dumpi2(x)                                               \
    idp                                                         \
    {                                                           \
        foreach (it, x)                                         \
        {                                                       \
            cout << "[" << (it)->F << " , " << (it)->S << "] "; \
        }                                                       \
        puts("");                                               \
    }

#define read2d(a, w, h) rep(i, h) rep(j, w) cin >> a[i][j]
#define dump2d(a, w, h)                   \
    rep(i, h)                             \
    {                                     \
        rep(j, w) cout << a[i][j] << " "; \
        puts("");                         \
    }

typedef pair<li, li> PI;

li dp[111][3][3];
li n;
string s[2];

int main(void)
{
    cin >> n;
    cin >> s[0];
    cin >> s[1];
    if (s[0][0] == s[1][0]) {
        rep(i, 3) {
            dp[0][i][i] = 1;
        }
    } else {
        rep(i, 3) {
            rep(j, 3) {
                if (i == j)continue;
                dp[0][i][j] = 1;
            }
        }
    }

    repp(i, 1, n) {
        if (s[0][i] == s[1][i]) {
            rep(j, 3) {
                rep(pre1, 3) {
                    rep(pre2, 3) {
                        if (j == pre1 || j == pre2) {
                            continue;
                        }
                        dp[i][j][j] = (dp[i][j][j] + dp[i - 1][pre1][pre2]) % MOD;
                    }
                }
            }
        } else {
            rep(now1, 3) {
                rep(now2, 3) {
                    if (now1 == now2) {
                        continue;
                    }
                    rep(pre1, 3) {
                        if ((now1 == pre1) ^ (s[0][i] == s[0][i - 1])) {
                            continue;
                        }
                        rep(pre2, 3) {
                            if ((now2 == pre2) ^ (s[1][i] == s[1][i - 1])) {
                                continue;
                            }
                            dp[i][now1][now2] = (dp[i][now1][now2] + dp[i - 1][pre1][pre2]) % MOD;
                        }
                    }
                }
            }
        }
    }

    li res = 0;
    rep(i, 3) {
        rep(j, 3) {
            res = (res + dp[n - 1][i][j]) % MOD;
        }
    }
    cout << res << endl;
    /*
        rep(i, n) {
            rep(i1, 3) {
                rep(i2, 3) {
                    cout << dp[i][i1][i2] << " ";
                } cout << endl;
            } cout << endl;
        }
    */

    return 0;
}

Measured memory_bytes: 256