#include<stdio.h>

#include<string.h>

#include<algorithm>

#include<map>

#include<math.h>

using namespace std;

#define mod 1000000007

long long int a[2];

long long int b[2][2];

void init()

{

    a[0]=1;

    a[1]=0;

    b[0][0]=3;

    b[0][1]=1;

    b[1][0]=1;

    b[1][1]=3;

}

void cheng()

{

    long long q=((a[0]*b[0][0])%mod+(a[1]*b[1][0])%mod)%mod;

    long long w=((a[0]*b[0][1])%mod+(a[1]*b[1][1])%mod)%mod;

    a[0]=q;

    a[1]=w;

    return ;

}

void ch()

{

    long long int z[2][2];

    for(int i=0;i<2;i++)

    {

        for(int k=0;k<2;k++)

        {

            z[i][k]=0;

            for(int j=0;j<2;j++)

            {

                z[i][k]+=(b[i][j]*b[j][k])%mod;

            }

            z[i][k]%=mod;

        }

    }

    for(int i=0;i<2;i++)

    for(int k=0;k<2;k++)

    {

        b[i][k]=z[i][k];

    }

    return ;

}

void cal(long long int n)

{

    if(n==0)

    return ;

    if(n%2==0)

    {

        n/=2;

        ch();

        cal(n);

    }

    else if(n%2!=0)

    {

        n--;

        cheng();

        cal(n);

    }

}

long long n;

int main(){

while(~scanf("%lld",&n))

{

    init();

    cal(n);

    printf("%lld\n",a[0]);

}

}