cin>>n>>m;
 
        vector<vll> rev(n,vll(n,0));
        rep(i,m)
        {
            cin>>l>>r;
            l--;r--;
            rev[r][l]++;
        }
 
        vector<vll> dp(n,vll(1<<n,0));
        dp[0][1]=1;
 
        rep(state,(1<<n))
        {
            rep(curr,n)
            {
                if(state & (1<<curr))
                {
                    rep(ch,n)
                    {
                        if(state & (1<<curr))
                        {
                        dp[curr][state]+=(dp[ch][state^(1<<curr)]*rev[curr][ch]);
                        dp[curr][state]%=MOD;
                        }
                        
                    }
                }
            }
                
        }
 
        out(dp[n-1][(1<<n)-1]);