#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); }