Weights hoj 2635

#include <stdio.h>
#include <cstring>
#define mod 9999997
long long tmp[3][3],ans[3][3],a[3][3];
void mul(long long a[][3],long long b[][3])
{
    for(int i=1; i<=2; i++)
        for(int j=1; j<=2; j++)
        {
            tmp[i][j]=0;
            for(int k=1; k<=2; k++)
                tmp[i][j]+=a[i][k]*b[k][j];
            if(tmp[i][j]>=mod) tmp[i][j]%=mod;
        }
    memcpy(a,tmp,sizeof(tmp));
}
void pow(long long ans[][3],long long b[][3],int n)
{
    ans[1][1]=ans[2][2]=1;
    ans[1][2]=ans[2][1]=0;
    while(n)
    {
        if(n&1) mul(ans,b);
        mul(b,b);
        n>>=1;
    }
}
int main()
{
    int n;
    while(scanf("%d",&n)==1)
    {
        memset(ans,0,sizeof(ans));
        a[1][1]=3;
        a[1][2]=a[2][2]=1;
        a[2][1]=0;
        pow(ans,a,n);
        printf("%d\n",ans[1][2]%mod);
    }
    return 0;
}