import sequtils,strutils

var
    M, K : int
    mat : array[51,array[51, int64]]
    modnum : int64 = 1000_000_007
(M, K) = stdin.readline.split.map(parseInt)


for n in 0..<M:
    var t : int
    for i in 0..<M:
        t = (n + i) mod M
        mat[t][n] += 1
        t = (n * i) mod M
        mat[t][n] += 1
proc multi(A, B : array[0..50, array[0..50, int64]]): array[51, array[51,int64]]=
    var L = M
    for i in 0..<L:
        for j in 0..<L:
            for k in 0..<L:
                result[i][j] += A[i][k] * B[k][j]
                result[i][j] = result[i][j] mod modnum

proc powmat(A : array[0..50, array[0..50, int64]], m : int) : array[51, array[51,int64]]=
    if m == 1:
        return A

    if m == 0:
        for i in 0..M:
            result[i][i] = 1

    if m mod 2 == 0:
        return powmat(multi(A,A), m div 2)

    return multi(powmat(multi(A, A) , m div 2), A)

mat = powmat(mat, K)


echo mat[0][0]