#include <bits/stdc++.h>


using namespace std;
typedef long long ll;

const int N = 5e3 + 2, MOD = (int)1e9 + 7;

int n, h[N], k;

void add(int &x, int y) {
    x += y;
    if(x >= MOD) x -= MOD;
}
int dp[N];
map<int, int> pd[N][N];
void test() {
    cin >> n;
    for(int i = 1; i <= n; i++) {
        cin >> h[i];
    }
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= n; j++) {
            pd[i][j].clear();
        }
    }
    dp[0] = 1;
    for(int i = 1; i <= n; i++) {
        if(i == 1) pd[i][1][h[i]] = dp[i - 1];
        else {
            for(int j = h[i - 1] - 1; j <= h[i - 1] + 1; j++) {
                if(j != h[i]) {
                    for(int f = 1; f <= i; f++) {
                        add(pd[i][1][h[i]], pd[i - 1][f][j]);
                    }
                }
            }
        }
        for(int j = 1; j <= i; j++) {
            if(pd[i - 1][j].count(h[i] - 1)) pd[i][j][h[i] - 1] = pd[i - 1][j][h[i] - 1];
            if(pd[i - 1][j].count(h[i] + 1)) pd[i][j][h[i] + 1] = pd[i - 1][j][h[i] + 1];
            if(j > 1) {
                if(pd[i - 1][j - 1].count(h[i])) pd[i][j][h[i]] = pd[i - 1][j - 1][h[i]];
            }
        }
        int j = i - 1;
        dp[i] = dp[i - 1];
        for(; j >= 1 && abs(h[j] - h[i]) == 1; j--) {
            add(dp[i], dp[j - 1]);
        }
        for(int f = 1; f <= i; f++) {
            if(pd[j][f].count(h[i] - 1)) add(dp[i], f * 1ll * pd[j][f][h[i] - 1] % MOD);
            if(pd[j][f].count(h[i] + 1)) add(dp[i], f * 1ll * pd[j][f][h[i] + 1] % MOD);
        }
    }

    cout << dp[n] << '\n';
}

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

    int t = 1;
    cin >> t;

    while(t--) {
        test();
    }
}