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
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
#include <iostream>
#include <vector>
#include <set>
#include <algorithm>
#include <map>
int n,m,a,b;
std::vector<std::vector<int>>sa;
bool st[200000];
bool par[200000];
bool odw[200000];
int parz;
int nprz;
int p1;
int n1;
bool odd_cy;
void dfs(int w,bool p){
    odw[w]=true;
    par[w]=p;
    if(p){
        parz++;
        if(st[w])p1++;
    }else{
        nprz++;
        if(st[w])n1++;
    }
    for(auto i:sa[w]){
        if(odw[i]){
            if(par[i]==p)
                odd_cy=true;
        }else
            dfs(i,!p);
    }
}

long long x;
long long y;

void euk(long long a, long long b)
{
    if(b!=0)
    {
        euk(b, a%b);
        long long pom = y;
        y = x  - a/b*y;
        x = pom;
    }
}
long long euklides(long long a, long long b){
    x=1;
    y=0;
    euk(a,b);
    if(y<0)
        y+=a;
    return y;
}


long long res=1;
long long mod=1000000007LL;
void mn(long long m){
    res=(res*m)%mod;
    //std::cerr<<"mn"<<m<<std::endl;
}
long long mn(long long s1,long long s2){
    return (s1*s2)%mod;
}

long long su(long long s1,long long s2){
    return (s1+s2)%mod;
}
long long dz(long long ac,long long dz){
    //std::cerr<<ac<<" "<<dz<<std::endl;
    return (ac*euklides(mod,dz)%mod);
}
int main() {
    std::ios_base::sync_with_stdio(0);
    std::cin.tie(NULL);
    //std::cerr<<euklides(mod,2)<<std::endl;

    std::cin>>n>>m;
    for(int i=0;i<n;i++)
        std::cin>>st[i];
    sa.resize(n);
    for(int i=0;i<m;i++){
        std::cin>>a>>b;
        a--;b--;
        sa[a].emplace_back(b);
        sa[b].emplace_back(a);
    }
    for(int i=0;i<n;i++){
        //std::cerr<<i<<std::endl;
        parz=0;
        nprz=0;
        p1=0;
        n1=0;
        odd_cy=false;
        if(!odw[i]) {
            dfs(i,true);
            //std::cerr<<p1<<" "<<n1<<" "<<parz<<" "<<nprz<<std::endl;
            if(odd_cy)
                for(int i=1;i<parz+nprz;i++)
                    mn(2);
            else {
                if ((p1 == 0 && n1 == nprz) || (n1 == 0 && p1 == parz))
                    continue;
                a=std::min(p1,n1);
                p1-=a;
                n1-=a;
                long long acp=1;
                long long acn=1;
                long long acc=0;
                for(int i=0;i<p1;i++){
                    acp=mn(acp,parz-i);
                    acp=dz(acp,i+1);
                }
                for(int i=0;i<n1;i++){
                    acn=mn(acn,nprz-i);
                    acn=dz(acn,i+1);
                }
                acc+=mn(acp,acn);
                //std::cerr<<"acc"<<acp<<" "<<acn<<std::endl;
                while(p1<parz&&n1<nprz){
                    acp=mn(acp,parz-p1);
                    acp=dz(acp,p1+1);
                    acn=mn(acn,nprz-n1);
                    acn=dz(acn,n1+1);
                    p1++;
                    n1++;
                    acc+=mn(acp,acn);
                    //std::cerr<<"acc"<<acp<<" "<< acn<<std::endl;
                }
                mn(acc%mod);
            }


        }


    }
    std::cout<<res<<std::endl;
}