I=lambda:map(int,raw_input().split());R=xrange;p=10**9+7 
N,M=I();w=I();W=max(w)+1;E=[N*[0]for i in R(N)]
D=[W*[0]for i in R(W)];D[0][0]=1
for i in R(1,W):
	for j in R(i+1):D[i][j]=j*(D[i-1][j]+D[i-1][j-1])%p
for i in R(N):E[i][i]=w[i]
for t in R(M):i,j=I();E[i-1][j-1]=w[i-1]
for k in R(N):
	for i in R(N):
		for j in R(N):E[i][j]=max(E[i][j],min(E[i][k],E[k][j]))
A=0
for i in R(N):
	for j in R(N):
		if w[j]<=E[i][j]:A=(A+D[w[i]][w[j]])%p
print A