#include <bits/stdc++.h>
using namespace std;

const int mod = 1000000007;
const int maxn = 1000025;

int madd(int a, int b) {
    return a + b >= mod ? a + b - mod : a + b;
}

struct BI {
    vector<int> evt;
    int b[maxn];
    void revert() {
        for (int p: evt)
            for (++p; p < maxn; p += p & -p)
                b[p] = 0;
    }
    void add(int p, int v) {
        evt.emplace_back(p);
        for (++p; p < maxn; p += p & -p)
            b[p] = madd(b[p], v);
    }
    int query(int p) {
        int r = 0;
        for (++p; p > 0; p -= p & -p)
            r = madd(r, b[p]);
        return r;
    }
} BIT;

int main() {
    cin.tie(nullptr)->sync_with_stdio(false);

    int n, k;
    cin >> n >> k;

    vector<int> u;
    vector<vector<int>> v(n, vector<int>(k));
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < k; j++) {
            cin >> v[i][j];
            u.push_back(v[i][j]);
        }
    }

    sort(u.begin(), u.end());
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < k; j++) {
            v[i][j] = lower_bound(u.begin(), u.end(), v[i][j])
                - u.begin();
        }
    }
    BIT.add(0, 1);

    vector<int> dp(k);
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < k; j++) {
            dp[j] = BIT.query(v[i][j]);
        }
        BIT.revert();
        for (int j = 0; j < k; j++) {
            BIT.add(v[i][j], dp[j]);
        }
    }

    int z = 0;
    for (int j = 0; j < k; j++)
        z = madd(z, dp[j]);
    cout << z << '\n';
    return 0;
}