1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
//Jakub Trela
#include <bits/stdc++.h>
using namespace std;

typedef long long LL;

vector<vector<int>> g;
vector<int> tab;
vector<bool> visited;
vector<int> kolor;
vector<LL> silnia;

LL MOD=1e9+7;

LL pot(LL a,LL b){
    if(b==0)return 1;
    if(b%2==0){
        LL w=pot(a,b/2);
        return (w*w)%MOD;
    }else{
        LL w=pot(a,b-1);
        return (w*a)%MOD;
    }
}

LL newton(LL n,LL k){
    return (silnia[n]*pot((silnia[k]*silnia[n-k])%MOD,MOD-2))%MOD;
}

int s,s0,s1,p0,p1;
bool czyDwu;

void dfs(int w,int k){
    visited[w]=1;
    kolor[w]=k;
    s++;
    if(k)s1++;
    else s0++;
    if(k)p1+=tab[w];
    else p0+=tab[w];
    for(int v:g[w]){
        if(visited[v] && kolor[v]==k)czyDwu=0;
        if(visited[v])continue;
        dfs(v,k^1);
    }
}

int main(){
    ios_base::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);

    int n,m;
    cin>>n>>m;

    silnia.resize(n+1);
    silnia[0]=1;
    for(int i=1;i<=n;i++){
        silnia[i]=(silnia[i-1]*(LL)i)%MOD;
    }


    tab.resize(n+1);
    g.resize(n+1);
    visited.resize(n+1);
    kolor.resize(n+1);

    for(int i=1;i<=n;i++)cin>>tab[i];

    for(int i=0;i<m;i++){
        int a,b;
        cin>>a>>b;

        g[a].push_back(b);
        g[b].push_back(a);
    }

    LL ans=1;

    for(int i=1;i<=n;i++){
        s=0;s0=0;s1=0;
        p0=0;p1=0;
        czyDwu=1;
        dfs(i,0);
        if(czyDwu){
            if(p0>p1){
                swap(s0,s1);
                swap(p0,p1);
            }
            LL res=0;
            for(int j=0;j<=s;j++){
                if(j<=s0 && j+p1-p0<=s1){
                    res+=(newton(s0,j)*newton(s1,j+p1-p0))%MOD;
                    res%=MOD;
                }
            }
            ans*=res;
            ans%=MOD;
        }else{
            ans*=pot(2,s-1);
            ans%=MOD;
        }
    }

    cout<<ans<<endl;

    return 0;
}