1. 1
2. 2
3. 3
4. 4
5. 5
6. 6
7. 7
8. 8
9. 9
10. 10
11. 11
12. 12
13. 13
14. 14
15. 15
16. 16
17. 17
18. 18
19. 19
20. 20
21. 21
22. 22
23. 23
24. 24
25. 25
26. 26
27. 27
28. 28
29. 29
30. 30
31. 31
32. 32
33. 33
34. 34
35. 35
36. 36
37. 37
38. 38
39. 39
40. 40
41. 41
42. 42
43. 43
44. 44
45. 45
46. 46
47. 47
48. 48
49. 49
50. 50
51. 51
52. 52
53. 53
54. 54
55. 55
56. 56
57. 57
58. 58
59. 59
60. 60
61. 61
62. 62
63. 63
64. 64
65. 65
66. 66
67. 67
68. 68
69. 69
70. 70
71. 71
72. 72
73. 73
74. 74
75. 75
76. 76
77. 77
78. 78
79. 79
80. 80
81. 81
82. 82
83. 83
84. 84
85. 85
86. 86
87. 87
88. 88
89. 89
90. 90
91. 91
92. 92
93. 93
94. 94
95. 95
96. 96
97. 97
98. 98
99. 99
100. 100
101. 101
102. 102
103. 103
104. 104
105. 105
106. 106
107. 107
108. 108

#include <bits/stdc++.h>
#include <ext/pb_ds/tree_policy.hpp>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/rope>
#define FOR(i, a, b, in) for (int i = (a); i < (b); i += in)
#define REP(i, a, b) for (int i = (a); i < (b); i++)
#define RFOR(i, a, b, in) for (int i = (a) - 1; i >= (b); i -= in)
#define RREP(i, a, b) for (int i = (a) - 1; i >= (b); i--)
#define trav(a, x) for (auto& a : x)
#define mp make_pair
#define pb push_back
#define eb emplace_back
#define f first
#define s second
#define lb lower_bound
#define ub upper_bound
#define sz(x) (int)x.size()
#define all(x) begin(x), end(x)
#define rsz resize
#define mem(a, b) memset(a, (b), sizeof(a))
using namespace std;
using namespace __gnu_pbds;
using namespace __gnu_cxx;
typedef string str;
typedef long long ll;
typedef long double ld;
typedef complex<ld> cd;
typedef pair<int, int> ii; typedef pair<ll, ll> pl; typedef pair<ld, ld> pd;
typedef vector<int> vi; typedef vector<ll> vl; typedef vector<ld> vd;
typedef vector<ii> vii; typedef vector<pl> vpl; typedef vector<pd> vpd;
typedef tree<int, null_type, less<int>, rb_tree_tag, tree_order_statistics_node_update> ordered_set;
typedef tree<int, null_type, less_equal<int>, rb_tree_tag, tree_order_statistics_node_update> ordered_multiset; // WARNING: May be broken
constexpr int INF = 1e9;
constexpr ll LINF = 1e18;
const ll MOD = 1e9+7;
const ll INV = 5e8+4;
const int S = 200;
const ld PI = 4*atan((ld)1);
void io(str s) {
    if (fopen((s+".in").c_str(), "r")) {
        freopen((s+".in").c_str(), "r", stdin);
        freopen((s+".out").c_str(), "w", stdout);
    } 
    ios_base::sync_with_stdio(false);
    cin.tie(NULL), cout.tie(NULL);
}


int N, K, A[50001], cl[50001][21], cr[50001][21], ans[200001];

void solve(int l, int r, vector<pair<ii, int>>& q) {
    int m = (l + r) >> 1;
    for (int i = l; i <= r; i++) memset(cl[i], 0, sizeof(cl[i])), memset(cr[i], 0, sizeof(cr[i]));
    int c[21][21] = { 0 };
    for (int i = m; i >= l; i--) {
        if (i < m) memcpy(cl[i], cl[i + 1], sizeof(cl[i]));
        for (int j = A[i]; j <= K; j++) {
            for (int k = A[i]; k <= K; k++) {
                cl[i][j] = (c[k][j] + cl[i][j]) % MOD;
                c[A[i]][j] = (c[k][j] + c[A[i]][j]) % MOD;
            }
        }
        cl[i][A[i]]++, c[A[i]][A[i]]++;
    }
    memset(c, 0, sizeof(c));
    for (int i = m + 1; i <= r; i++) {
        if (i > m + 1) memcpy(cr[i], cr[i - 1], sizeof(cr[i]));
        for (int j = A[i]; j > 0; j--) {
            for (int k = A[i]; k > 0; k--) {
                cr[i][j] = (c[j][k] + cr[i][j]) % MOD;
                c[j][A[i]] = (c[j][k] + c[j][A[i]]) % MOD;
            }
        }
        cr[i][A[i]]++, c[A[i]][A[i]]++;
    }
    vector<pair<ii, int>> ql, qr;
    for (auto& x : q) {
        if (x.f.f <= m && x.f.s > m) {
            ll cnt = 0;
            for (int i = 1; i <= K; i++) {
                cnt = (cl[x.f.f][i] + cnt) % MOD;
                ans[x.s] = (cnt * cr[x.f.s][i] + cl[x.f.f][i] + cr[x.f.s][i] + ans[x.s]) % MOD;
            }
        }
        else if (x.f.s <= m) ql.pb(x);
        else qr.pb(x);
    }
    if (l == r) return;
    solve(l, m, ql), solve(m + 1, r, qr);
}

int main() {
    io("nondec");

    cin >> N >> K;
    for (int i = 0; i < N; i++) cin >> A[i];

    int Q;
    cin >> Q;
    vector<pair<ii, int>> q(Q);
    for (int i = 0; i < Q; i++) {
        cin >> q[i].f.f >> q[i].f.s;
        q[i].f.f--, q[i].f.s--;
        q[i].s = i;
    }
    solve(0, N - 1, q);
    for (int i = 0; i < Q; i++) cout << ans[i] + (q[i].f.f == q[i].f.s) + 1 << '\n';    
}