#include <bits/stdc++.h>
#define syosu(x) fixed<<setprecision(x)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> P;
typedef pair<double,double> pdd;
typedef pair<ll,ll> pll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<double> vd;
typedef vector<vd> vvd;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<string> vs;
typedef vector<P> vp;
typedef vector<vp> vvp;
typedef vector<pll> vpll;
typedef pair<int,P> pip;
typedef vector<pip> vip;
const int inf=1<<30;
const ll INF=1ll<<60;
const double pi=acos(-1);
const double eps=1e-8;
const ull mod=1e9+7;
const int dx[4]={-1,0,1,0},dy[4]={0,-1,0,1};

const int M=1005;
ll n,m,dp[M][M];
vi a,b;
vip e;
vvi g;

void dfs(int v,int t){
	b[v]=t;
	for(auto u:g[v]) if(b[u]==inf) dfs(u,t);
}

int main(){
	cin>>n>>m;
	a=vi(n);
	e=vip(m);
	for(int i=0;i<n;i++) cin>>a[i];
	for(int i=0;i<m;i++){
		int u,v;
		cin>>u>>v;
		u--;v--;
		e[i]={a[v],{u,v}};
	}
	sort(e.rbegin(),e.rend());
	dp[1][1]++;
	for(int i=2;i<M;i++) for(int j=1;j<=i;j++) dp[i][j]=(dp[i-1][j]+(j?dp[i-1][j-1]:0))*j%mod;
	ll res=0;
	for(int i=0;i<n;i++){
		int A=a[i];
		b=vi(n,inf);
		b[i]=A;
		g=vvi(n);
		for(int j=0;j<m;j++){
			P p=e[j].second;
			int u=p.first,v=p.second;
			g[u].push_back(v);
			if(b[u]<inf&&b[v]==inf) dfs(v,min(A,a[v]));
		}
		for(int j=0;j<n;j++) if(a[j]<=b[j]&&b[j]<inf) (res+=dp[A][a[j]])%=mod;
	}
	cout<<res<<endl;
}