#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using namespace __gnu_pbds; typedef tree, rb_tree_tag, tree_order_statistics_node_update> ordered_set; typedef long long LL; typedef double D; #define all(v) (v).begin(), (v).end() mt19937 gene(chrono::system_clock::now().time_since_epoch().count()); typedef complex Complex; #define fi first #define se second #define ins insert #define pb push_back inline char GET_CHAR(){ const int maxn = 131072; static char buf[maxn],*p1=buf,*p2=buf; return p1==p2&&(p2=(p1=buf)+fread(buf,1,maxn,stdin),p1==p2)?EOF:*p1++; } inline int getInt() { int res(0); char c = getchar(); while(c < '0') c = getchar(); while(c >= '0') { res = res * 10 + (c - '0'); c = getchar(); } return res; } inline LL fastpo(LL x, LL n, LL mod) { LL res(1); while(n) { if(n & 1) { res = res * (LL)x % mod; } x = x * (LL) x % mod; n /= 2; } return res; } template struct Num { LL a; Num operator + (const Num & b) { return Num{(a + b.a) % mod}; } Num operator - (const Num & b) { return Num{(a - b.a + mod) % mod}; } Num operator * (const Num & b) { return Num{a * b.a % mod}; } Num operator / (const Num & b) { return Num{a * fastpo(b.a, mod - 2, mod) % mod}; } void operator += (const Num & b) {if((a += b.a) >= mod) a -= mod;} void operator -= (const Num & b) {if((a -= b.a) < 0) a += mod;} void operator *= (const Num & b) { a = a * b.a % mod; } void operator /= (const Num & b) { a = a * fastpo(b.a, mod - 2, mod) % mod; } void operator = (const Num & b) { a = b.a; } void operator = (const LL & b) { a = b; } }; template ostream & operator << (ostream & os, const Num & a) { os << a.a; return os; } LL gcd(LL a, LL b) { return b ? gcd(b, a % b) : a; } inline string itoa(LL x, int width = 0) { string res; if(x == 0) res.push_back('0'); while(x) { res.push_back('0' + x % 10); x /= 10; } while((int)res.size() < width) res.push_back('0'); reverse(res.begin(), res.end()); return res; } const int _B = 131072; char buf[_B]; int _bl = 0; inline void flush() { fwrite(buf, 1, _bl, stdout); _bl = 0; } __inline void _putchar(char c) { if(_bl == _B) flush(); buf[_bl++] = c; } inline void print(LL x, char c) { static char tmp[20]; int l = 0; if(!x) tmp[l++] = '0'; else { while(x) { tmp[l++] = x % 10 + '0'; x /= 10; } } for(int i = l - 1; i >= 0; i--) _putchar(tmp[i]); _putchar(c); } typedef double C; struct P { C x, y; void scan() { double _x, _y; scanf("%lf%lf", &_x, &_y); x = _x; y = _y; } void print() { cout << '(' << x << ' ' << y << ')' << endl; } P operator + (const P & b) const { return P{x + b.x, y + b.y}; } P operator - (const P & b) const { return P{x - b.x, y - b.y}; } C operator * (const P & b) const { return x * b.y - y * b.x; } C operator % (const P & b) const { return x * b.x + y * b.y; } }; P operator * (const C & x, const P & b) { return P{x * b.x, x * b.y}; } const int N = 1555; const int LOG = 20; const int mod = 1e9 + 7; const int inf = 1e9 + 7; int n, m; int dx[4] = {1, 0, -1, 0}; int dy[4] = {0, 1, 0, -1}; int rela[N]; int getr(int x) { int p = x; while(rela[p] != p) p = rela[p]; int p1 = p; p = x; while(rela[p] != p) { int p2 = rela[p]; rela[p] = p1; p = p2; } return p1; } char s[N]; int a[N]; int dp[N][N]; int sgt[N * 4]; void build(int k, int le, int ri) { memset(sgt, 0, sizeof(sgt)); } int getsum(int k, int le, int ri, int ql, int qr) { if(ql > qr) return 0; if(ql <= le && ri <= qr) { return sgt[k]; }else { int res = 0; int mid = (le + ri) / 2; if(ql <= mid) res = (res + getsum(k * 2, le, mid, ql, qr)) % mod; if(qr >= mid + 1) res = (res + getsum(k * 2 + 1, mid + 1, ri, ql, qr)) % mod; return res; } } void inc(int k, int le, int ri, int pos, int val) { if(le == ri) { sgt[k] = (sgt[k] + val) % mod; }else { int mid = (le + ri) / 2; if(pos <= mid) inc(k * 2, le, mid, pos, val); else inc(k * 2 + 1, mid + 1, ri, pos, val); sgt[k] = (sgt[k * 2] + sgt[k * 2 + 1]) % mod; } } int main() { int n, k; scanf("%d%d", &n, &k); vector vec; for(int i = 1; i <= n; i++) { scanf("%d", &a[i]); vec.pb(a[i]); } sort(all(vec)); for(int i = 1; i <= n; i++) { a[i] = lower_bound(all(vec), a[i]) - vec.begin() + 1; } scanf("%s", s + 1); for(int j = 1; j <= n; j++) { dp[0][j] = 1; } for(int i = 1; i <= k; i++) { build(1, 1, n); for(int j = 1; j <= n; j++) { if(s[i] == '>') { dp[i][j] = getsum(1, 1, n, a[j] + 1, n); }else { dp[i][j] = getsum(1, 1, n, 1, a[j] - 1); } inc(1, 1, n, a[j], dp[i - 1][j]); } } int ans = 0; for(int i = 1; i <= n; i++) ans = (ans + dp[k][i]) % mod; cout << ans << endl; }