Go: package main

import (
	"bufio"
	"os"
	"strconv"
	"strings"
	"fmt"
)

var scanner = bufio.NewScanner(os.Stdin)

type Point struct {
	x, y int
}

func main() {
	dx := []int{1, 0, -1, 0}
	dy := []int{0, -1, 0, 1}
	maze := make([][]rune, 30)
	for i := range maze {
		maze[i] = make([]rune, 30)
	}
	next := make([][]rune, 30)
	for i := range next {
		next[i] = make([]rune, 30)
	}
	check := make([][]int, 30)
	for i := range check {
		check[i] = make([]int, 30)
	}
	for scanner.Scan() {
		buf := strings.Split(scanner.Text(), " ")
		W, _ := strconv.Atoi(buf[0])
		H, _ := strconv.Atoi(buf[1])
		if W == 0 && H == 0 {
			break
		}
		for i := 0; i < H; i++ {
			scanner.Scan()
			buf_ := scanner.Text()
			for j, v := range buf_ {
				maze[i][j] = v
			}
		}
		man := 0
		for i := 0; i < H; i++ {
			for j := 0; j < H; j++ {
				if maze[i][j] != '.' && maze[i][j] != 'X' && maze[i][j] != '#' {
					switch maze[i][j] {
					case 'E':
						maze[i][j] = 0
						man++
					case 'N':
						maze[i][j] = 1
						man++
					case 'W':
						maze[i][j] = 2
						man++
					case 'S':
						maze[i][j] = 3
						man++
					}
				}
			}
		}
		ans := -1
		for t := 0; t <= 180; t++ {
			if man == 0 {
				ans = t
				break
			}
			for i := 0; i < H; i++ {
				for j := 0; j < W; j++ {
					check[i][j] = 0x7FFFFFFF
				}
			}
			for i := 0; i < H; i++ {
				for j := 0; j < W; j++ {
					if maze[i][j] > 3 {
						continue
					}
					l := (maze[i][j] + 3) % 4
					for k := 0; k < 4; k++ {
						y := i + dy[l]
						x := j + dx[l]
						if 0 <= y && y < H && 0 <= x && x < W && (maze[y][x] == '.' || maze[y][x] == 'X') {
							maze[i][j] = l
							check[y][x] = min(check[y][x], int((l+2)%4))
							break
						}
						l = (l + 1) % 4
					}
				}
			}
			for i := 0; i < H; i++ {
				for j := 0; j < W; j++ {
					next[i][j] = maze[i][j]
				}
			}
			for i := 0; i < H; i++ {
				for j := 0; j < W; j++ {
					if maze[i][j] > 3 {
						continue
					}
					y := i + dy[maze[i][j]]
					x := j + dx[maze[i][j]]
					if 0 <= y && y < H && 0 <= x && x < W && check[y][x] == int((maze[i][j]+2)%4) {
						if maze[y][x] == '.' {
							next[y][x] = maze[i][j]
						} else if maze[y][x] == 'X' {
							man--
						}
						next[i][j] = '.'
					}
				}
			}
			for i := 0; i < H; i++ {
				for j := 0; j < W; j++ {
					maze[i][j]=next[i][j]
				}
			}
		}
		if ^ans != 0 {
			fmt.Println(ans)
		} else {
			fmt.Println("NA")
		}
	}
}

func min(a, b int) int {
	if a < b {
		return a
	}
	return b
}

C++: #include <iostream>
#include <complex>
#include <sstream>
#include <string>
#include <algorithm>
#include <deque>
#include <list>
#include <map>
#include <numeric>
#include <queue>
#include <vector>
#include <set>
#include <limits>
#include <cstdio>
#include <cctype>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <ctime>
using namespace std;
#define REP(i, j) for(int i = 0; i < (int)(j); ++i)
#define FOR(i, j, k) for(int i = (int)(j); i < (int)(k); ++i)
#define SORT(v) sort((v).begin(), (v).end())
#define REVERSE(v) reverse((v).begin(), (v).end())
typedef pair<int, int> P;
const int L = 4;

int W, H;
int my[] = {0, -1, 0, 1};
int mx[] = {1, 0, -1, 0};
map<char, int> ctoi;

class C{
  public:
    int y, x, d, f, g;
    C(){ f = true; }
    C(int yy, int xx, int dd) { y = yy; x = xx; d = dd; f = true; g = false; }
    P front() { return P(y + my[d], x + mx[d]); }
    void move() { y = y + my[d]; x = x + mx[d]; }
};

bool isPerson(int y, int x, vector<string> &v) { return (v[y][x] != '#' && v[y][x] != '.' && v[y][x] != 'X'); }

bool canMoveBeforeTurn(int tar, vector<string> &v, vector<C> &ps){
  if(!ps[tar].f) return false;
  P next = ps[tar].front();
  if(v[next.first][next.second] != '.' && v[next.first][next.second] != 'X') return false;
  REP(i, ps.size()) if(i != tar && ps[i].f && next == P(ps[i].y, ps[i].x)) return false;
  return true;
}

bool check(int tar, int n, vector<C> &ps){
  vector<int> v(4);
  v[0] = 2; v[1] = 3; v[2] = 0; v[3] = 1;
  return (v[tar] < v[n]);
}

bool canMoveAfterTurn(int tar, vector<C> &ps){
  REP(i, ps.size()){
    if(!ps[i].f || tar == i) continue;
    if(ps[i].front() == ps[tar].front() && !check(ps[tar].d, ps[i].d, ps)) return false;
  }
  return true;
}

bool isEnd(vector<C> &ps){
  REP(i, ps.size()) if(ps[i].f) return false;
  return true;
}

void disp(vector<string> &v, vector<C> &ps){
  vector<string> vv = v;
  REP(i, ps.size()){
    if(!ps[i].f) continue;
    if(ps[i].d == 0) vv[ps[i].y][ps[i].x] = 'E';
    if(ps[i].d == 1) vv[ps[i].y][ps[i].x] = 'N';
    if(ps[i].d == 2) vv[ps[i].y][ps[i].x] = 'W';
    if(ps[i].d == 3) vv[ps[i].y][ps[i].x] = 'S';
  }
  cout <<endl;
  REP(i, v.size()) cout <<vv[i] <<endl;
}

int main() {
  ctoi['E'] = 0; ctoi['N'] = 1; ctoi['W'] = 2; ctoi['S'] = 3;
  while(cin >>W >>H && W && H){
    vector<string> v(H);
    REP(i, H) cin >>v[i];
    vector<C> ps;
    REP(i, H) REP(j, W) if(isPerson(i, j, v)) { ps.push_back(C(i, j, ctoi[v[i][j]])); v[i][j] = '.'; }

    int res = 0;
    while(!isEnd(ps) && res++ < 180){
      //disp(v, ps);
      REP(i, ps.size()){
        if(!ps[i].f) continue;
        int nds[] = {(ps[i].d + 3) % L, ps[i].d, (ps[i].d + 1) % L, (ps[i].d + 2) % L};
        for(int nd : nds){
          ps[i].d = nd;
          if(canMoveBeforeTurn(i, v, ps)) break;
        }
      }
      REP(i, ps.size()) if(ps[i].f && canMoveAfterTurn(i, ps)) ps[i].g = 1;
      REP(i, ps.size()){
        if(!ps[i].f || ps[i].g == 0) continue;
        ps[i].g = 0;
        ps[i].move();
        if(v[ps[i].y][ps[i].x] == 'X') ps[i].f = false;
      }
    }
    if(res > 180) cout <<"NA" <<endl;
    else cout <<res <<endl;
  }
  return 0;
}

C++: #include <iostream>
#include <vector>
#include <algorithm>
#include <complex>

using namespace std;

typedef complex<int> P;

struct human {
	P p;
	int dir;

	human(P p, int dir) : p(p), dir(dir) { }
	human() { this -> p = P(0, 0); this -> dir = 0; }
};

bool cp(human a, human b) {
	return a.dir > b.dir;
}

int abs(int n) {
	return n < 0 ? -n : n;
}

P dir(int n) {
	n %= 4;
	return n == 0 ? P(0, -1) : n == 1 ? P(1, 0) : n == 2 ? P(0, 1) : P(-1, 0);
}

int W, H;
int m[30][30];
int state(P p) {
//	if (p.real() < 0 || p.real() >= W || p.imag() < 0 || p.imag() >= H) cout << "ERROR" << endl;
//	cout << p.real() << " " << p.imag() << endl;
	return m[p.imag()][p.real()];
}

int main() {
	while (cin >> W >> H, W || H) {
		vector<human> v;

		for (int y = 0; y < H; ++y) {
			string str; cin >> str;
			for (int x = 0; x < W; ++x)
				m[y][x] =
					  str[x] == '#' ? 1
					: str[x] == 'X' ? -1
					: str[x] == 'N' ? (v.push_back( human( P(x, y), 0) ), 2)
					: str[x] == 'E' ? (v.push_back( human( P(x, y), 1) ), 2)
					: str[x] == 'S' ? (v.push_back( human( P(x, y), 2) ), 2)
					: str[x] == 'W' ? (v.push_back( human( P(x, y), 3) ), 2)
					: 0
					;
		}

		int t = 0;
		while (v.size() && ++t <= 180) {
			// 向き変更
			for (int i = 0; i < v.size(); ++i)
				v[i].dir =
					  state(v[i].p + dir(v[i].dir + 1)) <= 0 ? (v[i].dir + 1) % 4
					: state(v[i].p + dir(v[i].dir))     <= 0 ? v[i].dir
					: state(v[i].p + dir(v[i].dir + 3)) <= 0 ? (v[i].dir + 3) % 4
					: state(v[i].p + dir(v[i].dir + 2)) <= 0 ? (v[i].dir + 2) % 4
					: v[i].dir
					;

			// 向きでソート
			sort( v.begin(), v.end(), cp );

			// 移動
			for (int i = 0; i < v.size(); ++i) {
				int s = state(v[i].p + dir(v[i].dir));
				if (s > 0) continue;

				m[v[i].p.imag()][v[i].p.real()] = 0;
				v[i].p += dir(v[i].dir);

				if (s == 0) m[v[i].p.imag()][v[i].p.real()] = 2;
				else v.erase( v.begin() + i-- );
			}
		}

		if (t > 180) cout << "NA" << endl;
		else cout << t << endl;
	}
}