Problem:
 Score : 1600 points    Problem Statement Joisino has a bar of length N , which has M marks on it.
The distance from the left end of the bar to the i -th mark is X_i .  She will place several squares on this bar.
Here, the following conditions must be met:   Only squares with integral length sides can be placed.  Each square must be placed so that its bottom side touches the bar.  The bar must be completely covered by squares.
That is, no square may stick out of the bar, and no part of the bar may be left uncovered.  The boundary line of two squares may not be directly above a mark.     Examples of arrangements that satisfy/violate the conditions   The beauty of an arrangement of squares is defined as the product of the areas of all the squares placed.
Joisino is interested in the sum of the beauty over all possible arrangements that satisfy the conditions.
Write a program to find it.
Since it can be extremely large, print the sum modulo 10^9+7 .      Constraints  All input values are integers.  1 \leq N \leq 10^9  0 \leq M \leq 10^5  1 \leq X_1 < X_2 < ... < X_{M-1} < X_M \leq N-1         Input Input is given from Standard Input in the following format:  N  M  X_1  X_2  ...  X_{M-1}  X_M       Output Print the sum of the beauty over all possible arrangements that satisfy the conditions, modulo 10^9+7 .        Sample Input 1 3 1
2      Sample Output 1 13  There are two possible arrangements:   Place a square of side length 1 to the left, and place another square of side length 2 to the right  Place a square of side length 3   The sum of the beauty of these arrangements is (1 \times 1 \times 2 \times 2) + (3 \times 3) = 13 .       Sample Input 2 5 2
2 3      Sample Output 2 66       Sample Input 3 10 9
1 2 3 4 5 6 7 8 9      Sample Output 3 100       Sample Input 4 1000000000 0      Sample Output 4 693316425

Solution (Python):
class Matrix():
    mod=10**9+7

    def set_mod(m):
        Matrix.mod=m

    def __init__(self,L):
        self.row=len(L)
        self.column=len(L[0])
        self._matrix=L
        for i in range(self.row):
            for j in range(self.column):
                self._matrix[i][j]%=Matrix.mod

    def __getitem__(self,item):
        if type(item)==int:
            raise IndexError("you must specific row and column")
        elif len(item)!=2:
            raise IndexError("you must specific row and column")

        i,j=item
        return self._matrix[i][j]

    def __setitem__(self,item,val):
        if type(item)==int:
            raise IndexError("you must specific row and column")
        elif len(item)!=2:
            raise IndexError("you must specific row and column")

        i,j=item
        self._matrix[i][j]=val

    def __add__(self,other):
        if (self.row,self.column)!=(other.row,other.column):
            raise SizeError("sizes of matrixes are different")

        res=[[0 for j in range(self.column)] for i in range(self.row)]
        for i in range(self.row):
            for j in range(self.column):
                res[i][j]=self._matrix[i][j]+other._matrix[i][j]
                res[i][j]%=Matrix.mod
        return Matrix(res)

    def __mul__(self,other):
        if type(other)!=int:
            if self.column!=other.row:
                raise SizeError("sizes of matrixes are different")

            res=[[0 for j in range(other.column)] for i in range(self.row)]
            for i in range(self.row):
                for j in range(other.column):
                    temp=0
                    for k in range(self.column):
                        temp+=self._matrix[i][k]*other._matrix[k][j]
                    res[i][j]=temp%Matrix.mod
            return Matrix(res)
        else:
            n=other
            res=[[(n*self._matrix[i][j])%Matrix.mod for j in range(self.column)] for i in range(self.row)]
            return Matrix(res)

    def __pow__(self,m):
        if self.column!=self.row:
            raise MatrixPowError("the size of row must be the same as that of column")

        n=self.row
        res=Matrix([[int(i==j) for i in range(n)] for j in range(n)])
        while m:
            if m%2==1:
                res=res*self
            self=self*self
            m//=2
        return res

    def __str__(self):
        res=[]
        for i in range(self.row):
            for j in range(self.column):
                res.append(str(self._matrix[i][j]))
                res.append(" ")
            res.append("\n")
        res=res[:len(res)-1]
        return "".join(res)

N,M=map(int,input().split())
X=list(map(int,input().split()))+[N]
dp=[0]*(M+2)
dp[0]=1
A=Matrix([[4,-2,1],[1,0,0],[0,1,0]])
DP=Matrix([[0],[0],[0]])
first=0
base=Matrix([[5],[1],[0]])

for i in range(1,M+2):
    DP=(A**(X[i-1]-first))*(DP+base*(-dp[i-1]))
    dp[i]=DP[2,0]
    first=X[i-1]

print((-dp[-1])%(10**9+7))

Solution (Python):
import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines

N,M,*X = map(int,read().split())

MOD = 10 ** 9 + 7

def mult(a,b,c,d,e,f):
    
    a,b,c,d,e = a*d,a*e+b*d,a*f+b*e+c*d,b*f+c*e,c*f
    b += e; c -= 4*e; d += 2*e; e = 0
    a += d; b -= 4*d; c += 2*d; d = 0
    a %= MOD; b %= MOD; c %= MOD
    return a,b,c


M = 10 ** 5
A1 = [0] * (M+1)
a,b,c = 1,0,0
for i in range(M+1):
    A1[i] = (a,b,c)
    a,b,c = b+4*a,c-2*a,a
    a %= MOD; b %= MOD; c %= MOD


A2 = [0] * (M+1)
a,b,c = 1,0,0
d,e,f = A1[M]
for i in range(M+1):
    A2[i] = (a,b,c)
    a,b,c = mult(a,b,c,d,e,f)

def power(n):
    
    q,r = divmod(n,M)
    a,b,c = A1[r]
    d,e,f = A2[q]
    return mult(a,b,c,d,e,f)

X.append(N)

a,b,c = 0,1,1
prev_x = 0
for x in X:
    a,b,c = mult(a,b,c,*power(x - prev_x))
    b -= a; c -= a
    prev_x = x

answer = a
print(answer)