#![allow(non_snake_case,unused_imports,dead_code)]
use std::cmp::max;
use std::cmp::min;
use std::str::FromStr;
use std::collections::{HashMap,HashSet,VecDeque,BinaryHeap};
use std::cmp::Ordering;

use proconio::input;
use proconio::marker::{Chars,Usize1,Isize1};

fn comp_f64(a: &f64, b: &f64) -> Ordering {
    if a < b {
        return Ordering::Less;
    } else if a > b {
        return Ordering::Greater;
    }
    Ordering::Equal
}

fn main() {
    input!{
        N:usize,
        S1:Chars,
        S2:Chars
    }

    let MOD = 1_000_000_007;

    let mut ans = ModInt::new(3,MOD);

    let mut p = S1[0];
    let mut pn;
    if S1[0]==S2[0] {
        pn = 1;
    } else {
        pn = 2;
        ans *= 2;
    };
    for i in 1..N{
        if S1[i] == p{continue;}
        let n;
        if S1[i] == S2[i]{
            n = 1;
        } else {
            n = 2;
        }
        if n==1 && pn==2{
            ans *= 1;
        } else if n==1 && pn==1{
            ans *= 2;
        } else if n==2 && pn==1{
            ans *= 2;
        } else {
            ans *= 3;
        }
        p = S1[i];
        pn = n;
    }
    println!("{}",ans);
}

use std::ops::{Add, AddAssign, Sub, SubAssign, Mul, MulAssign, Div, DivAssign};
use std::fmt;
use std::fmt::Display;

#[derive(Copy, Clone, Debug)]
struct ModInt{
    n: isize,
    MOD: isize
}
impl ModInt {
    fn new(n:isize,MOD:isize) -> ModInt{
        ModInt {n:(n%MOD+MOD)%MOD, MOD:MOD}
    }
    fn pow_mod(a:isize,b:isize,MOD:isize) -> isize{
        if b==0{ return 1;}
        if b%2==1{ return a*ModInt::pow_mod(a,b-1,MOD)%MOD}
        ModInt::pow_mod(a*a%MOD,b/2,MOD)
    }
    fn inv(self) -> ModInt{
        let n = ModInt::pow_mod(self.n,self.MOD-2,self.MOD);
        ModInt::new(n,self.MOD)
    }
}
impl Display for ModInt{
    fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result{
        write!(f,"{}",self.n)
    }
}

impl Add for ModInt{
    type Output = ModInt;
    fn add(self, other:ModInt) -> ModInt
    {
        ModInt::new(self.n+ other.n,self.MOD)
    }
}
impl Add<isize> for ModInt{
    type Output = ModInt;
    fn add(self, other:isize) -> ModInt
    {
        ModInt::new(self.n+ other,self.MOD)
    }
}

impl AddAssign for ModInt{
    fn add_assign(&mut self, other:ModInt)
    {
        self.n = (self.n + other.n)%self.MOD;
    }
}
impl AddAssign<isize> for ModInt{
    fn add_assign(&mut self, other:isize)
    {
        self.n = (self.n + other)%self.MOD;
    }
}

impl Sub for ModInt{
    type Output = ModInt;
    fn sub(self, other:ModInt) -> ModInt
    {
        ModInt::new(self.n-other.n,self.MOD)
    }
}
impl Sub<isize> for ModInt{
    type Output = ModInt;
    fn sub(self, other:isize) -> ModInt
    {
        ModInt::new(self.n-other,self.MOD)
    }
}
impl SubAssign for ModInt{
    fn sub_assign(&mut self, other:ModInt)
    {
        self.n = (self.n + self.MOD - other.n%self.MOD)%self.MOD;
    }
}
impl SubAssign<isize> for ModInt{
    fn sub_assign(&mut self, other:isize)
    {
        self.n = (self.n + self.MOD - other%self.MOD)%self.MOD;
    }
}

impl Mul for ModInt{
    type Output = ModInt;
    fn mul(self, other:ModInt) -> ModInt
    {
        ModInt::new(self.n*other.n,self.MOD)
    }
}
impl Mul<isize> for ModInt{
    type Output = ModInt;
    fn mul(self, other:isize) -> ModInt
    {
        ModInt::new(self.n*other,self.MOD)
    }
}
impl MulAssign for ModInt{
    fn mul_assign(&mut self, other:ModInt)
    {
        self.n = self.n*(other.n%self.MOD)%self.MOD;
    }
}
impl MulAssign<isize> for ModInt{
    fn mul_assign(&mut self, other:isize)
    {
        let o = (other%self.MOD+self.MOD)%self.MOD;
        self.n = self.n*o%self.MOD;
    }
}

impl Div for ModInt{
    type Output = ModInt;
    fn div(self, other:ModInt) -> ModInt
    {
        other.inv() * self
    }
}
impl Div<isize> for ModInt{
    type Output = ModInt;
    fn div(self, other:isize) -> ModInt
    {
        ModInt::new(other,self.MOD).inv() * self
    }
}
impl DivAssign for ModInt{
    fn div_assign(&mut self, other:ModInt)
    {
        *self *= other.inv()
    }
}
impl DivAssign<isize> for ModInt{
    fn div_assign(&mut self, other:isize)
    {
        *self *= ModInt::pow_mod(other, self.MOD-2, self.MOD);
    }
}

Measured memory_bytes: 2136