#include <bits/stdc++.h>
using namespace std;
const int MAXN = 5e3 + 10;
template <class T>
inline void _read(T &x)
{
	x = 0;
	char t = getchar();
	while (!isdigit(t)) t = getchar();
	while (isdigit(t))
	{
		x = x * 10 + t - '0';
		t = getchar();
	}
}
int n, m;
char ch[MAXN];
int sum[MAXN], f[MAXN][MAXN], las[MAXN], tot;
const int mo = 1e9 + 7;
int add(int x, int y)
{
	return x + y >= mo ? x + y - mo : x + y;
}
int main()
{
	_read(n), _read(m);
	scanf("%s", ch + 1);
	for (int i = 1; i <= n + 1; ++i)
	{
		sum[i] = sum[i - 1], las[i] = i;
		if (ch[i] == '1') ++tot, ++sum[i];
	}
	for (int i = 1, x, y; i <= m; ++i) _read(x), _read(y), las[x] = max(las[x], y);
	for(int i = 1; i <= n + 1; ++i)las[i] = max(las[i], las[i - 1]);
	f[1][sum[las[1]]] = 1;
	for(int i = 1; i <= n; ++i)
		for(int j = 0; j <= tot; ++j)
			if(f[i][j])
			{
				int a = las[i], b = las[i + 1];
				if(j)f[i + 1][j + sum[b] - sum[a] - 1] = add(f[i + 1][j + sum[b] - sum[a] - 1], f[i][j]);
				if(las[i] + 1 - i - j) f[i + 1][j + sum[b] - sum[a]] = (f[i + 1][j + sum[b] - sum[a]], f[i][j]);
			}
	printf("%d\n", f[n + 1][0]);
}