#include <bits/stdc++.h>
#define sort_(a) sort(a.begin(), a.end())
#define rsort(a) sort(a.rbegin(), a.rend())
#define rep(i, n) for (int i = 0; i < n; i++)
#define rrep(i, n) for (int i = n - 1; 0 <= i; i--)
#define for_(i, start, end) for (int i = start; i < end; i++)
#define rfor(i, start, end) for (int i = start - 1; 0 <= i; i--)
#define all(a) a.begin(), a.end()
using ll = long long;
using ld = long double;
using namespace std;
constexpr int INF32       = 1'050'000'000;
constexpr long long INF64 = 4'000'000'000'000'000'000;
constexpr int MOD7        = 1'000'000'007;
constexpr int MOD53       = 998'244'353;
constexpr double PI       = 3.14159265358979323846;
template <class T> inline bool chmin(T &a, T b) {
    if (a > b) {
        a = b;
        return true;
    }
    return false;
}
template <class T> inline bool chmax(T &a, T b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}
void print() { std::cout << '\n'; }
template <class H, class... T> void print(H &&head, T &&... args) {
    std::cout << head;
    sizeof...(args) == 0 ? std::cout << "" : std::cout << ' ';
    print(std::forward<T>(args)...);
}
template <class T> void print(std::vector<T> &v) {
    for (int i = 0; i < v.size(); i++) {
        std::cout << v[i];
        i == v.size() - 1 ? std::cout << '\n' : std::cout << ' ';
    }
}
template <class T> void print(std::vector<std::vector<T>> &v) {
    for (int i = 0; i < v.size(); i++) {
        for (int j = 0; j < v[i].size(); j++) {
            std::cout << v[i][j];
            j == v[i].size() - 1 ? std::cout << '\n' : std::cout << ' ';
        }
    }
}
void scan() {}
template <class H, class... T> void scan(H &&head, T &&... args) {
    std::cin >> head;
    scan(std::forward<T>(args)...);
}
template <class T> void scan(std::vector<T> &v) {
    for (auto &&i : v) {
        scan(i);
    }
}

template <int mod> class modint {
  private:
    long long val;

  public:
    constexpr modint(long long x = 0) noexcept : val(x % mod) {
        if (val < 0)
            val += mod;
    };
    constexpr modint operator-() const noexcept { return val ? mod - val : 0; }
    constexpr modint operator+(const modint &x) const noexcept {
        return modint(*this) += x;
    }
    constexpr modint &operator++() noexcept {
        val += 1;
        if (mod <= val)
            val -= mod;
        return *this;
    }
    constexpr modint operator++(int) noexcept {
        modint<mod> ret = *this;
        val += 1;
        if (mod <= val)
            val -= mod;
        return ret;
    }
    constexpr modint operator-(const modint &x) const noexcept {
        return modint(*this) -= x;
    }
    constexpr modint &operator--() noexcept {
        val -= 1;
        if (val < 0)
            val += mod;
        return *this;
    }
    constexpr modint operator--(int) noexcept {
        modint<mod> ret = *this;
        val -= 1;
        if (val < 0)
            val += mod;
        return ret;
    }
    constexpr modint operator*(const modint &x) const noexcept {
        return modint(*this) *= x;
    }
    constexpr modint operator/(const modint &x) const noexcept {
        return modint(*this) /= x;
    }
    constexpr modint &operator+=(const modint &x) noexcept {
        val += x.val;
        if (mod <= val)
            val -= mod;
        return *this;
    }
    constexpr modint &operator-=(const modint &x) noexcept {
        val -= x.val;
        if (val < 0)
            val += mod;
        return *this;
    }
    constexpr modint &operator*=(const modint &x) noexcept {
        val = val * x.val % mod;
        return *this;
    }
    constexpr modint &operator/=(const modint &x) noexcept {
        val *= x.inv();
        return *this;
    }
    constexpr modint inv() const noexcept {
        long long a = val, b = mod, u = 1, v = 0, t;
        while (0 < b) {
            t = a / b;
            a -= t * b;
            swap(a, b);
            u -= t * v;
            swap(u, v);
        }
        return modint(u);
    }
    constexpr bool operator==(const modint &x) const noexcept {
        return this->val == x.val;
    }
    constexpr bool operator!=(const modint &x) const noexcept {
        return this->val != x.val;
    }
    constexpr bool operator<(const modint &x) const noexcept {
        return val < x.val;
    }
    constexpr bool operator<=(const modint &x) const noexcept {
        return val <= x.val;
    }
    constexpr bool operator>(const modint &x) const noexcept {
        return x.val < val;
    }
    constexpr bool operator>=(const modint &x) const noexcept {
        return x.val <= val;
    }
    friend constexpr ostream &operator<<(ostream &os, modint<mod> &x) noexcept {
        return os << x.val;
    }
    friend constexpr istream &operator>>(istream &is, modint<mod> &x) noexcept {
        return is >> x.val;
    }
    friend constexpr modint<mod> power(const modint<mod> &radix,
                                       long long x) noexcept {
        if (x == 0)
            return 1;
        auto t = power(radix, x / 2);
        t      = t * t;
        if (x & 1)
            t = t * radix;
        return t;
    }
};

using mint = modint<MOD7>;

int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);

    ll N, K;
    scan(N, K);
    vector<int> a(N);
    scan(a);

    ll cnt  = 0;
    ll cnt2 = 0;
    rep(i, N) {
        rep(j, N) {
            if (a[j] < a[i])
                cnt2++;
            if (i < j && a[j] < a[i])
                cnt++;
        }
    }
    mint ans = cnt2;
    ans *= K * (K - 1) / 2;
    ans += cnt * K;
    print(ans);

    return 0;
}

Measured memory_bytes: 892