#include <bits/stdc++.h>
#define all(vec) vec.begin(), vec.end()
using namespace std;
using ll = long long;
using P = pair<int, int>;
constexpr ll INF = (1LL << 30) - 1LL;
constexpr ll LINF = (1LL << 60) - 1LL;
constexpr ll MOD = 1e9 + 7;
template <typename T> void chmin(T &a, T b) { a = min(a, b); }
template <typename T> void chmax(T &a, T b) { a = max(a, b); }
template <std::uint_fast64_t Modulus> class modint {
    using u64 = std::uint_fast64_t;

  public:
    u64 a;
    constexpr modint(const u64 x = 0) noexcept : a(x % Modulus) {}
    constexpr u64 &value() noexcept { return a; }
    constexpr const u64 &value() const noexcept { return a; }
    constexpr modint operator+(const modint rhs) const noexcept {
        return modint(*this) += rhs;
    }
    constexpr modint operator-(const modint rhs) const noexcept {
        return modint(*this) -= rhs;
    }
    constexpr modint operator*(const modint rhs) const noexcept {
        return modint(*this) *= rhs;
    }
    constexpr modint operator/(const modint rhs) const noexcept {
        return modint(*this) /= rhs;
    }
    constexpr modint operator^(const u64 rhs) const noexcept {
        return modint(*this) ^= rhs;
    }
    constexpr modint &operator+=(const modint rhs) noexcept {
        a += rhs.a;
        if (a >= Modulus) {
            a -= Modulus;
        }
        return *this;
    }
    constexpr modint &operator-=(const modint rhs) noexcept {
        if (a < rhs.a) {
            a += Modulus;
        }
        a -= rhs.a;
        return *this;
    }
    constexpr modint &operator*=(const modint rhs) noexcept {
        a = a * rhs.a % Modulus;
        return *this;
    }
    constexpr modint &operator/=(modint rhs) noexcept {
        u64 exp = Modulus - 2;
        while (exp) {
            if (exp % 2) {
                *this *= rhs;
            }
            rhs *= rhs;
            exp /= 2;
        }
        return *this;
    }
    constexpr modint &operator^=(u64 exp) {
        modint rhs = modint(*this);
        a = 1;
        while (exp) {
            if (exp % 2) {
                *this *= rhs;
            }
            rhs *= rhs;
            exp /= 2;
        }
        return *this;
    }
};
using mint = modint<MOD>;
vector<mint> f, fi;
void comb(int n) {
    f.resize(n + 1);
    fi.resize(n + 1);
    f[0] = 1;
    for (ll i = 1; i <= n; i++) {
        f[i] = f[i - 1] * mint(i);
    }
    fi[n] = mint(1) / f[n];
    for (ll i = n - 1; i >= 0; i--) {
        fi[i] = fi[i + 1] * mint(i + 1LL);
    }
}
mint ncr(ll n, ll r) {
    if (n < r) {
        return mint(0);
    }
    if (r == 0) {
        return mint(1);
    }
    return f[n] * fi[n - r] * fi[r];
}
int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    int n, m;
    cin >> n >> m;
    string s;
    cin >> s;
    vector<int> co(n + 1);
    for (int i = 0; i < n; i++) {
        co[i + 1] = co[i] + (s[i] == '1');
    }
    int mr = 0;
    vector<P> v;
    vector<int> rs(n, -1);
    for (int i = 0; i < m; i++) {
        int l, r;
        cin >> l >> r;
        --l;
        --r;
        chmax(rs[l], r);
    }
    for (int i = 0; i < n; i++) {
        if (rs[i] != -1 && rs[i] > mr) {
            mr = rs[i];
            v.push_back(P(i, rs[i]));
        }
    }
    v.push_back(P(n, n));
    vector<vector<mint>> dp(m + 1, vector<mint>(n + 1));
    dp[0][co[v[0].first]] = 1;
    comb(n + 10);
    for (int i = 0; i + 1 < v.size(); i++) {
        int l = v[i].first, r = v[i].second;
        int a = r - l + 1, c = a;
        chmin(c, v[i + 1].first - l);
        for (int j = 0; j <= n; j++) {
            if (dp[i][j].a == 0) {
                continue;
            }
            int b = co[r + 1] - j;
            for (int d = max(0, b - a + c); d <= min(b, c); d++) {
                int to = j + d;
                if (c == a) {
                    to = co[v[i + 1].first];
                }
                dp[i + 1][to] += dp[i][j] * ncr(c, d);
            }
        }
    }
    mint res = 0;
    for (int i = 0; i <= n; i++) {
        res += dp[v.size() - 1][i];
    }
    cout << res.a << endl;
}

Measured memory_bytes: 70912