// 行列の累乗をすればいい
// DPというより普通の再帰
// 引数: exponent l, O(log K)個しかみない
// 結果: G^l mod p
// 遷移: O(N^3)

#pragma region Header
#ifdef DEBUG
#include "debug.hpp"
#else
#define debug(...)
#define IFD(...)
#ifndef NDEBUG
#define NDEBUG
#endif
#endif
#include <utility>
#include <cstdint>
#include <cmath>
#include <cassert>
#include <vector>
#include <array>
#include <deque>
#include <queue>
#include <stack>
#include <tuple>
#include <set>
#include <map>
#include <unordered_map>
#include <bitset>
#include <string>
#include <functional>
#include <algorithm>
#include <numeric>
#include <stdexcept>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <boost/hana/functional/fix.hpp>
using boost::hana::fix;
using namespace std;
using i32 = int32_t;
using i64 = int64_t;
using u32 = uint32_t;
using u64 = uint64_t;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define repd(i, n) for (int i = (n)-1; i >= 0; --i)
#define FOR(i, a, b) for (int i = (a); i <= (b); ++i)
#define FORD(i, b, a) for (int i = (b); i >= (a); --i)
template <typename T>
int SZ(T const &x)
{
    return static_cast<int>(x.size());
}
template <class T>
bool chmax(T &a, T const &b)
{
    if (a < b)
    {
        a = b;
        return true;
    }
    return false;
}
template <class T>
bool chmin(T &a, T const &b)
{
    if (b < a)
    {
        a = b;
        return true;
    }
    return false;
}
#pragma endregion

template <u32 p>
struct modint
{
private:
    u32 data; // in [0,p)

    constexpr static inline i64 p_i64 = static_cast<i64>(p);

public:
    constexpr modint() noexcept
        : data(0) {}

    constexpr modint(i64 const x) noexcept
        : data(static_cast<u32>((x % p_i64 >= 0) ? x % p_i64 : x % p_i64 + p_i64)) {}

    constexpr modint(u32 const x, bool const) noexcept
        : data(x) {}

    // convert to u32
    constexpr u32 const &val() const noexcept
    {
        return data;
    }

    constexpr modint &operator+=(modint const x) noexcept
    {
        data += x.data;
        if (data >= p)
        {
            data -= p;
        }
        return *this;
    }

    constexpr modint &operator-=(modint const x) noexcept
    {
        if (data < x.data)
        {
            data += p;
        }
        data -= x.data;
        return *this;
    }

    constexpr modint &operator*=(modint const x) noexcept
    {
        data = static_cast<u32>(static_cast<u64>(data) * x.data % p);
        return *this;
    }

    constexpr modint &operator/=(modint x)
    {
        if (x.data == 0)
        {
            throw domain_error("Division by modint zero");
        }
        u32 exp = p - 2;
        // (*this) * (x ** exp) : invariant
        while (exp > 0)
        {
            if (exp & 1)
            {
                *this *= x;
            }
            x *= x;
            exp /= 2;
        }
        return *this;
    }

    constexpr modint &operator++() noexcept
    {
        return *this += modint(1, true);
    }

    constexpr modint operator++(int) noexcept
    {
        modint const ret(*this);
        ++(*this);
        return ret;
    }

    constexpr modint &operator--() noexcept
    {
        return *this -= modint(1, true);
    }

    constexpr modint operator--(int) noexcept
    {
        modint const ret(*this);
        --(*this);
        return ret;
    }

    constexpr modint operator-() const noexcept
    {
        return modint(0, true) -= *this;
    }

    constexpr modint operator+(modint const x) const noexcept
    {
        return modint(*this) += x;
    }

    constexpr modint operator-(modint const x) const noexcept
    {
        return modint(*this) -= x;
    }

    constexpr modint operator*(modint const x) const noexcept
    {
        return modint(*this) *= x;
    }

    constexpr modint operator/(modint const x) const
    {
        return modint(*this) /= x;
    }

    // x >= 0 のときべき乗を計算
    constexpr modint pow(u64 exp) const noexcept
    {
        modint ret(1, true), base(*this);
        // ret * (base ** exp) : invariant, == (*this) ** x
        while (exp > 0)
        {
            if (exp & 1)
            {
                ret *= base;
            }
            base *= base;
            exp /= 2;
        }
        return ret;
    }

    constexpr modint inv() const noexcept
    {
        return modint(1, true) /= *this;
    }

    constexpr bool operator==(modint const x) const noexcept
    {
        return data == x.data;
    }
};

constexpr u32 mod = 1000'000'007;
using mint = modint<mod>;

int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);
    cout << fixed << setprecision(15);

    int n;
    cin >> n;
    i64 k;
    cin >> k;
    vector a_s(n, vector<mint>(n));
    rep(y, n)
    {
        rep(x, n)
        {
            int a; cin>> a;
            a_s[y][x] = a;
        }
    }

    auto mul = [n](vector<vector<mint>> const &m1, vector<vector<mint>> const &m2) {
        vector m3(n, vector<mint>(n, 0));
        rep(y, n)
        {
            rep(x, n)
            {
                rep(k, n)
                {
                    m3[y][x] += m1[y][k] * m2[k][x];
                }
            }
        }
        return m3;
    };

    vector dp(n, vector<mint>(n, 0));
    rep(i, n)
    {
        dp[i][i] = 1;
    }
    // 計算結果をdpに入れる
    auto loop = fix([&](auto self, i64 e) {
        if (e == 0)
        {
            return;
        }
        if (e & 1)
        {
            self(e - 1);
            dp = mul(dp, a_s);
        }
        else
        {
            self(e / 2);
            dp = mul(dp, dp);
        }
    });
    loop(k);

    mint ans = accumulate(dp.begin(), dp.end(), mint{}, [](mint acc, auto v) {
        return acc + accumulate(v.begin(), v.end(), mint{});
    });
    cout << ans.val() << '\n';

    return 0;
}


Measured memory_bytes: 3716