#include <bits/stdc++.h>
using namespace std;
const int N = 3010;
const int MOD = 1e9 + 7;
int n, m;
int li[N], ri[N];
int C[N][N], f[N][N];
string st;
void ADD(int &t, int d) {t += d; if (t >= MOD) t -= MOD;}
int main()
{
    cin >> n >> m;
    cin >> st;
    for (int i = 0; i <= n; ++ i)
    {
        C[i][0] = C[i][i] = 1;
        for (int j = 1; j < i; ++ j)
            C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % MOD;
    }
    for (int i = 1; i <= m; ++ i)
    {
        cin >> li[i] >> ri[i];
        if (ri[i] < ri[i - 1]) i --, m --;
    }
    li[m + 1] = ri[m + 1] = n + 1;
    f[0][0] = 1;
    for (int i = 1; i <= m; ++ i)
    {
        int c = 0;
        for (int j = max(ri[i - 1] + 1, li[i]); j <= ri[i]; ++ j)
            if (st[j - 1] == '1') c ++;
        for (int j = 0; j <= max(ri[i - 1] - li[i] + 1, 0); ++ j)
            for (int k = 0; k <= min(j + c, min(li[i + 1] - 1, ri[i]) - li[i] + 1); ++ k)
                if (j + c - k <= max(ri[i] - li[i + 1] + 1, 0))
                    ADD(f[i][j + c - k], 1ll * f[i - 1][j] * C[min(li[i + 1] - 1, ri[i]) - li[i] + 1][k] % MOD);
    }
    cout << f[m][0] << endl;
}
