#include #include #include using namespace atcoder; using mint = modint1000000007; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf 1000000001 template struct matrix{ F0 func0; F1 func1; vector> value; int height,width; T init_value0,init_value1; matrix(vector> X,F0 f0,F1 f1,T iv0,T iv1):func0(f0),func1(f1){ height = X.size(); width = X[0].size(); value = X; init_value0 = iv0,init_value1 = iv1; } matrix(int h,int w,F0 f0,F1 f1,T iv0,T iv1):func0(f0),func1(f1){ vector> d(h,vector(w,iv0)); height = h,width = w; value = d; init_value0 = iv0,init_value1 = iv1; } void set_1(){ for(int i=0;i R(height,another.width,func0,func1,init_value0,init_value1); for(int i=0;i R(height,width,func0,func1,init_value0,init_value1); R.set_1(); auto temp = *this; while(cnt!=0LL){ if(cnt%2==1){ R *= temp; } temp *= temp; cnt/=2; } return R; } }; int main(){ int t; cin>>t; auto f0 = [](mint a,mint b){ return a+b; }; auto f1 = [](mint a,mint b){ return a*b; }; rep(_,t){ int N,M,X; cin>>N>>M>>X; vector> temp = {{1,M},{M,1}}; matrix A(temp,f0,f1,0,1); A = A.beki(N); temp = {{1,0}}; matrix B(temp,f0,f1,0,1); B *= A; cout<