program main implicit none integer(16) :: i, j, m, k, mat0(0:49, 0:49), mat(0:49, 0:49) read *, m, k do i = 0, m - 1 mat0(i, i) = 1 do j = 0, m - 1 mat(mod(i+j, m), i) = mat(mod(i+j, m), i) + 1 mat(mod(i*j, m), i) = mat(mod(i*j, m), i) + 1 end do end do do while (k > 0) if (mod(k, 2) == 1) then mat0 = mod(matmul(mat0, mat), 1000000007) end if mat = mod(matmul(mat, mat), 1000000007) k = k / 2 end do print *, mat0(0, 0) end program main