#include #include using namespace std; const long long mod=1000000007; const vector>a{ {1,0,0,0,1,1}, {0,1,0,1,0,1}, {0,0,1,1,1,0}, {3,0,0,0,0,0}, {0,3,0,0,0,0}, {0,0,3,0,0,0} }; const vector> b={ {1},{0},{0},{3},{0},{0} }; vector> matrix_mul(const vector>& A,const vector>&B){ vector>C(A.size(),vector(B.at(0).size(),0)); for(int i=0;i> matrix_exp(vector> X,long long n){ vector>Y(X.size(),vector(X.size(),0)); for(int i=0;i0){ if(n&1){ Y=matrix_mul(Y,X); } X=matrix_mul(X,X); n>>=1; } return Y; } long long calculate_nth(long long n,vector> A){ A=matrix_exp(A,n); vector> res; res=matrix_mul(A,b); return res.at(3).at(0); } long long mod_pow(long long x,long long n,long long m){ long long res=1; while(n>0){ if(n&1){ res=res*x%m; } x=x*x%m; n>>=1; } return res; } void solve(){ long long N; cin>>N; long long ans,c,An; ans=calculate_nth(N,a); c=mod_pow(mod_pow(3,mod-2,mod),N+1,mod); An=ans*c%mod; cout<>T; for(int i=0;i