import sys,random,bisect from collections import deque,defaultdict from heapq import heapify,heappop,heappush from itertools import permutations from math import log,gcd input = lambda :sys.stdin.readline().rstrip() mi = lambda :map(int,input().split()) li = lambda :list(mi()) mod = 10**9 + 7 def mat_mul(X,Y): n,m = len(X),len(Y[0]) res = [[0 for j in range(m)] for i in range(n)] for i in range(n): for j in range(m): for k in range(len(Y)): res[i][j] += X[i][k] * Y[k][j] res[i][j] %= mod return res L,N,M = mi() K = li() cond = [] cond_to_i = {} t = 0 for i in range(1,N+1): if i in K: for j in range(1,i+2): cond.append((i,j)) cond_to_i[i,j] = t t += 1 else: cond.append((i,1)) cond_to_i[i,1] = t t += 1 cond.append((-1,-1)) cond_to_i[-1,-1] = t NN = len(cond) A = [[0 for j in range(NN)] for i in range(NN)] for i,j in cond: for n in range(1,N+1): if (i,j)==(-1,-1): ni,nj = (-1,-1) else: if i in K: if j==i+1: if n==i: ni,nj = i,j else: ni,nj = n,1 elif j==i: if n==i: ni,nj = i,j+1 else: ni,nj = -1,-1 else: if n==i: ni,nj = i,j+1 else: ni,nj = n,1 else: if n==i: ni,nj = i,j else: ni,nj = n,1 x,y = cond_to_i[i,j],cond_to_i[ni,nj] A[y][x] += 1 E = [[0]] * NN for i in range(1,N+1): t = cond_to_i[i,1] E[t] = [1] L -= 1 while L: if L & 1: E = mat_mul(A,E) A = mat_mul(A,A) L >>= 1 res = 0 for i,j in cond: t = cond_to_i[i,j] if (i,j)==(-1,-1): res += E[t][0] elif i in K and j==i: res += E[t][0] res %= mod #print(i,j,E[t][0]) print(res)