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
#include <bits/stdc++.h>
using namespace std;
const long long mod=1000000000+7;
long long QuickPow(long long a, long long b)
{
    long long wynik=1;
    while(b>0)
    {
        if(b&1){wynik=(wynik*a)%mod;}
        a=(a*a)%mod;
        b>>=1;
    }
    return wynik;
}
vector<int>D[200009];
long long silnia[400009],silnia2[400009];
bool tab[200009];
bool zaj[200009];
bool cyc[200009];
bool CzyDwu=1;
int ile1=0,ile2=0;
int jed1=0,jed2=0;
void DFS(int v)
{
    zaj[v]=1;
    if(cyc[v]==0)
    {
        ile1++;
        if(tab[v]){jed1++;}
    }
    else
    {
        ile2++;
        if(tab[v]){jed2++;}
    }
    for(int som : D[v])
    {
        if(zaj[som]==0)
        {
            cyc[som]=(cyc[v]^1);
            DFS(som);
        }
        else if(cyc[v]==cyc[som]){CzyDwu=0;}
    }
}
long long N(int a, int b)
{
    return (((silnia[a]*silnia2[b])%mod)*silnia2[a-b])%mod;
}
int main()
{ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    silnia[0]=1;
    for(int i=1;i<=400000;i++){silnia[i]=(silnia[i-1]*i)%mod;}
    silnia2[400000]=QuickPow(silnia[400000],mod-2);
    for(int i=400000;i>=1;i--){silnia2[i-1]=(silnia2[i]*i)%mod;}
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=n;i++){cin>>tab[i];}
    for(int i=1;i<=m;i++)
    {
        int a,b;
        cin>>a>>b;
        D[a].push_back(b);
        D[b].push_back(a);
    }
    long long WYNIK=1;
    for(int i=1;i<=n;i++)
    {
        ile1=0;ile2=0;
        jed1=0;jed2=0;
        CzyDwu=1;
        if(zaj[i]==0)
        {
            DFS(i);
            long long wynik=0;
            if(CzyDwu)
            {
                int r=jed2-jed1;
                for(int i=0;i<=ile1;i++)
                {
                    int j=i+r;
                    if(0<=j&&j<=ile2)
                    {
                        wynik=(wynik+N(ile1,i)*N(ile2,j))%mod;
                    }
                }
            }
            else
            {
                int ile=ile1+ile2;
                int bit=(jed1+jed2)%2;
                for(int i=bit;i<=ile;i+=2)
                {
                    wynik=(wynik+N(ile,i))%mod;
                }
            }
            WYNIK=(WYNIK*wynik)%mod;
        }
    }
    cout<<WYNIK<<'\n';
    return 0;
}