#include <bits/stdc++.h>
using namespace std;
const long long N = 1e5 + 10;
const long long M = 1e7 + 19;
long long Max(long long a, long long b) {
  if (a > b) return a;
  return b;
}
long long Min(long long a, long long b) {
  if (a > b) return b;
  return a;
}
long long big_mul(long long a, long long b, long long m) {
  long long x = 0, y = a % m;
  while (b > 0) {
    if (b & 1) x = (x + y) % m;
    y = (2 * y) % m;
    b >>= 1;
  }
  return x % m;
}
long long big_mod(long long a, long long b, long long m) {
  long long x = 1, y = a % m;
  while (b > 0) {
    if (b & 1) x = (x * y) % m;
    y = (y * y) % m;
    b >>= 1;
  }
  return x % m;
}
long long Big_mod(long long a, long long b, long long m) {
  long long x = 1, y = a % m;
  while (b > 0) {
    if (b & 1) x = big_mul(x, y, m);
    y = big_mul(y, y, m);
    b >>= 1;
  }
  return x % m;
}
long long mod_inverse(long long a, long long m) { return big_mod(a, m - 2, m); }
int ara[500][500];
vector<int> vec;
vector<int> parent(500);
int find_parent(int x) {
  if (parent[x] == x) return x;
  return parent[x] = find_parent(parent[x]);
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  cout.tie(0);
  ;
  int n;
  cin >> n;
  vec.push_back(0);
  for (int i = 0; i < n; i++) {
    int a;
    cin >> a;
    vec.push_back(a);
  }
  cin.ignore();
  for (int i = 1; i < n + 1; i++) {
    string str;
    cin >> str;
    cin.ignore();
    for (int j = 0; j < str.size(); j++) ara[i][j + 1] = str[j] - '0';
  }
  for (int i = 0; i < n + 1; i++) parent[i] = i;
  for (int i = 1; i < n + 1; i++) {
    for (int j = 1; j < n + 1; j++) {
      if (ara[i][j]) {
        int p = find_parent(i);
        int q = find_parent(j);
        if (p != q) parent[p] = parent[q];
      }
    }
  }
  vector<int> ans;
  for (int i = 1; i < n + 1; i++) {
    ans.clear();
    for (int j = 1; j < n + 1; j++) {
      if (find_parent(j) == i) ans.push_back(vec[j]);
    }
    sort(ans.begin(), ans.end());
    for (int j = n; j > 0; j--) {
      if (find_parent(j) == i) {
        vec[j] = ans[ans.size() - 1];
        ans.pop_back();
      }
    }
  }
  for (int i = 1; i < n + 1; i++) cout << vec[i] << " ";
  return 0;
}
