#include<stdio.h>
#include<algorithm>
using namespace std;
long long mod=1000000007;
long long rj[1100][1100];
long long C[1100][1100];
int c[110000];
int g[210][210];
int main(){
	for(int i=0;i<1100;i++){
		rj[i][0]=1;
		for(int j=1;j<1100;j++)rj[i][j]=rj[i][j-1]*i%mod;
	}
	C[0][0]=1;
	for(int i=0;i<1050;i++){
		for(int j=0;j<=i;j++){
			C[i+1][j]=(C[i+1][j]+C[i][j])%mod;
			C[i+1][j+1]=(C[i+1][j+1]+C[i][j])%mod;
		}
	}
	int a,b;scanf("%d%d",&a,&b);
	for(int i=0;i<a;i++)scanf("%d",c+i);
	for(int i=0;i<a;i++)for(int j=0;j<a;j++)g[i][j]=-999999999;
	for(int i=0;i<a;i++)g[i][i]=999999999;
	
	for(int i=0;i<b;i++){
		int p,q;scanf("%d%d",&p,&q);p--;q--;
		g[p][q]=min(c[q],c[p]);
	}
	for(int k=0;k<a;k++)for(int i=0;i<a;i++)for(int j=0;j<a;j++)g[i][j]=max(g[i][j],min(g[i][k],g[k][j]));
	long long ret=0;
	for(int i=0;i<a;i++)for(int j=0;j<a;j++){
		if(g[i][j]>=c[j]){
			for(int k=0;k<c[j];k++){
				long long tmp=rj[c[j]-k][c[i]]*C[c[j]][k]%mod;
				if(k%2)ret=(ret+mod-tmp)%mod;
				else ret=(ret+tmp)%mod;
			}
		//	printf("%d %d: %lld\n",i,j,ret);
		}
	}
	printf("%lld\n",ret);
}