const DIVISOR: u64 = 1_000_000_007;

fn calc(a: Vec<Vec<u64>>, b: Vec<Vec<u64>>) -> Vec<Vec<u64>> {
    let length = 4usize;
    let mut result: Vec<Vec<u64>> = vec![vec![0; length]; length];
    for i in 0..length {
        for j in 0..length {
            for k in 0..length {
                result[i][j] += a[i][k] * b[k][j];
                result[i][j] %= DIVISOR;
            }
        }
    }
    result
}

fn power(val: Vec<Vec<u64>>, p: u64) -> Vec<Vec<u64>> {
    let I: Vec<Vec<u64>> = vec![vec![1, 0, 0, 0], vec![0, 1, 0, 0], vec![0, 0, 0, 1], vec![0, 0, 0, 1]];
    let F: Vec<Vec<u64>> = vec![vec![0, 1, 1, 1], vec![1, 0, 1, 1], vec![1, 1, 1, 0], vec![1, 1, 1, 0]];
    if p == 0 { return I; }
    if p == 1 { return F; }
    let temp = power(val, p/2);
    let temp2 = temp.clone();
    let mut ret = calc(temp, temp2);
    if p % 2 == 1 {
        ret = calc(ret, F);
    }
    ret
}

fn main() {
    let mut n = String::new();
    std::io::stdin().read_line(&mut n).ok();
    let n: u64 = n.trim().parse().unwrap();

    let F: Vec<Vec<u64>> = vec![vec![0, 1, 1, 1], vec![1, 0, 1, 1], vec![1, 1, 1, 0], vec![1, 1, 1, 0]];
    let op = power(F, n);
    println!("{}", op[0][0]);
}