import itertools as iter import collections as coll import heapq as hq import bisect as bis from decimal import Decimal as dec from copy import deepcopy as dcopy import math import sys sys.setrecursionlimit(10**6) def input(): return sys.stdin.readline().rstrip() def getN(): return int(sys.stdin.readline()) def getNs(): return map(int,sys.stdin.readline().split()) def getList(): return list(map(int,sys.stdin.readline().split())) def strinps(n): return [sys.stdin.readline().rstrip() for _ in range(n)] pi = 3.141592653589793 mod = 10**9+7 MOD = 998244353 INF = math.inf dx = [1,0,-1,0]; dy = [0,1,0,-1] def prod_mat(A, B): res = [[0, 0], [0, 0]] for i in range(2): for j in range(2): for k in range(2): res[i][j] += A[i][k] * B[k][j] res[i][j] %= mod return res def exp_mat(A, n): res = [[1, 0], [0, 1]] while(n): if(n & 1): res = prod_mat(A, res) n >>= 1 A = prod_mat(A, A) return res """ Main Code """ a,b = getNs() n = getN() t = [getN() for _ in [0] * n] p = [[a, b], [1, 0]] for i in t: k = i >> 1 r = i & 1 q = exp_mat(p, k) v = [q[0][0] + q[0][1], q[1][0] + q[1][1]] if(r): print((v[0] + v[1] + v[0] * a % mod + v[1] * b % mod) % mod) else: print((v[0] + v[1]) % mod)