//Let's join Kaede Takagaki Fan Club !! #pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define int long long //#define L __int128 typedef long long ll; typedef pair P; typedef pair P1; typedef pair P2; #define pu push #define pb push_back #define eb emplace_back #define mp make_pair #define eps 1e-7 #define INF 1000000000 #define a first #define b second #define fi first #define sc second //#define rng(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define rep(i,x) for(int i=0;i bool chmax(t&a,u b){if(a bool chmin(t&a,u b){if(b using vc=vector; template ostream& operator<<(ostream& os,const pair& p){ return os<<"{"< ostream& operator<<(ostream& os,const vc& v){ os<<"{"; for(auto e:v)os<> 30)) * 0xbf58476d1ce4e5b9; x = (x ^ (x >> 27)) * 0x94d049bb133111eb; return x ^ (x >> 31); } size_t operator()(uint64_t x) const { static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count(); return splitmix64(x + FIXED_RANDOM); } //don't make x negative! size_t operator()(pair x)const{ return operator()(uint64_t(x.first)<<32|x.second); } }; //unordered_set -> dtype, null_type //unordered_map -> dtype(key), dtype(value) using namespace __gnu_pbds; template using hash_table=gp_hash_table; template void g(T &a){ cin >> a; } template void o(const T &a,bool space=false){ cout << a << (space?' ':'\n'); } //ios::sync_with_stdio(false); const ll mod = 1000000007; //const ll mod = 1000000007; mt19937_64 mt(chrono::steady_clock::now().time_since_epoch().count()); template void add(T&a,T b){ a+=b; if(a >= mod) a-=mod; } ll modpow(ll x,ll n){ ll res=1; while(n>0){ if(n&1) res=res*x%mod; x=x*x%mod; n>>=1; } return res; } #define _sz 1 ll F[_sz],R[_sz]; void make(){ F[0] = 1; for(int i=1;i<_sz;i++) F[i] = F[i-1]*i%mod; R[_sz-1] = modpow(F[_sz-1], mod-2); for(int i=_sz-2;i>=0;i--) R[i] = R[i+1] * (i+1) % mod; } ll C(int a,int b){ if(b < 0 || a < b) return 0; return F[a]*R[b]%mod*R[a-b]%mod; } template using mat = vc>; template mat make_e(int n, int m, T v){ matret(n, vc(m)); rep(i, min(n, m)) ret[i][i] = v; return ret; } template mat operator+(mat a, mat b){ assert(a.size() == b.size()); if(a.size()) assert(a[0].size() == b[0].size()); rep(i, a.size()) rep(j, a[0].size()){ a[i][j] += b[i][j]; if(a[i][j] >= mod) a[i][j] -= mod; } return a; } template mat operator-(mat a, mat b){ assert(a.size() == b.size()); if(a.size()) assert(a[0].size() == b[0].size()); rep(i, a.size()) rep(j, a[0].size()){ a[i][j] += mod - b[i][j]; if(a[i][j] >= mod) a[i][j] -= mod; } return a; } template mat operator*(mat a, mat b){ if(a.empty() || b.empty()) return mat(); int n = a.size(); int m = b[0].size(); int len = a[0].size(); assert(len == b.size()); matret(n); rep(i, n) ret[i] = vc(m); rep(i, n) rep(j, len) rep(k, m){ ll add = (ll)(a[i][j]) * b[j][k]; add %= mod; ret[i][k] += add; if(ret[i][k] >= mod) ret[i][k] -= mod; } return ret; } template ll det(matm){ int n = m.size(); assert(n > 0 && m[0].size() == n); ll ret = 1; for(int j=0;j>bs; template vcAxB(matA, vcb){ int r = A.size(), c = A[0].size(); assert(b.size() == r); int nxt = 0, free = 0; vcret(c, 0); rep(j, c){ int pivot = -1; for(int i=nxt;ifix(c, 0); vctop(nxt); rep(i, nxt){ rep(j, c){ if(A[i][j] != 0){ top[i] = j; fix[j] = 1; ret[j] = b[i]; break; } } } //error for(int i=nxt;i(); //get basis /*vc>basis(free); int id = 0; rep(j, c){ if(fix[j] == 1) continue; vctmp(c, 0); tmp[j] = 1; for(int i=0;i matinv(matA){ int r = A.size(), c = A[0].size(); if(c == r){ rep(i, r){ A[i].resize(r+r, 0); rep(j, r) A[i][r+j] = (i == j); } c += r; } assert(c == r+r); rep(j, r){ int pivot = -1; for(int i=j;i(); } if(j != pivot){ rep(k, c) swap(A[j][k],A[pivot][k]); } ll rev = modpow(A[j][j],mod-2); for(int g=j;gret = make_e(r, r, 0LL); rep(i, r) rep(j, r){ assert(A[i][j] == (i==j)); ret[i][j] = A[i][r+j]; } return ret; } int a, b, n, t[1005]; void solve(){ cin>>a>>b>>n;repn(i,n)cin>>t[i]; matp = make_e(3, 3, 1LL); auto q = make_e(3, 3, 0LL); p[2][0] = a; p[2][1] = b; q[0][2] = 1; q[1][0] = 1; auto pq = q * p; matbeg = make_e(3, 1, 0LL); repn(i,n){ beg[0][0] = 1; beg[1][0] = 1; beg[2][0] = 0; matbin[65]; bin[0] = pq; rep(i, 63) bin[i+1] = bin[i] * bin[i]; matzan = make_e(3, 3, 1LL); bool f = t[i]%2; t[i] /= 2; rep(ii, 63){ if(((t[i]>>ii)&1LL)){ zan = zan * bin[ii]; } } if(f){ zan = p * zan; } beg = zan * beg; int ans = 0; rep(i, 3) add(ans, beg[i][0]); o(ans); } } signed main(){ cin.tie(0); ios::sync_with_stdio(0); cout<> t; while(t--) solve(); }