#include "bits/stdc++.h"
#define rep(i, a, n) for (auto i = a; i <= (n); ++i)
#define revrep(i, a, n) for (auto i = n; i >= (a); --i)
#define all(a) a.begin(), a.end()
#define sz(a) (int)(a).size()
template<class T> inline bool chmax(T &a, T b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> inline bool chmin(T &a, T b) { if (b < a) { a = b; return 1; } return 0; }
using namespace std;

template<class A, class B> string to_string(const pair<A, B> &p);
template<class A, class B, class C> string to_string(const tuple<A, B, C> &p);
template<class A, class B, class C, class D> string to_string(const tuple<A, B, C, D> &p);
string to_string(const string &s) { return '"' + s + '"'; }
string to_string(const char *s) { return to_string((string) s); }
string to_string(char c) { return "'" + string(1, c) + "'"; }
string to_string(bool x) { return x ? "true" : "false"; }
template<class A, class T = typename A::value_type> string to_string(const A &v) {
	bool first = 1;
	string res = "{";
	for (const auto &x: v) {
		if (!first) res += ", ";
		first = 0;
		res += to_string(x);
	}
	res += "}";
	return res;
}
template<class A, class B> string to_string(const pair<A, B> &p) {
	return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";
}
template<class A, class B, class C> string to_string(const tuple<A, B, C> &p) {
	return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " + to_string(get<2>(p)) + ")";
}
template<class A, class B, class C, class D> string to_string(const tuple<A, B, C, D> &p) {
	return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " + to_string(get<2>(p)) + ", " + to_string(get<3>(p)) + ")";
}

void debug_out() { cerr << endl; }
template<class H, class... T> void debug_out(const H& h, const T&... t) {
	cerr << " " << to_string(h);
	debug_out(t...);
}
#ifndef ONLINE_JUDGE
	#define debug(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)
#else
	#define debug(...) if (0) puts("No effect.")
#endif

using ll = long long;
using pii = pair<int, int>;
using vi = vector<int>;
using vvi = vector<vi>;

/**
 * Author: Yuhao Yao
 * Date: 22-12-03
 * Description: Kuhn Matching algorithm for \textbf{bipartite} graph $G = (L \cup R, E)$. Edges $E$ should be described as pairs such that pair $(x, y)$ means that there is an edge between the $x$-th vertex in $L$ and the $y$-th vertex in $R$. Returns a vector $lm$, where $lm[i]$ denotes the vertex in $R$ matched to the $i$-th vertex in $R$.
 * Time: O((|L| + |R|) |E|).
 * Status: tested on https://www.luogu.com.cn/problem/P2764, https://www.luogu.com.cn/problem/P3386.
 */
vi Kuhn(int n, int m, const vector<pii> &es) {
	vector<vi> g(n);
	for (auto [x, y]: es) g[x].push_back(y);
	vi rm(m, -1);
	rep(i, 0, n - 1) {
		vi vis(m);
		auto dfs = [&](auto &dfs, int x) -> int {
			for (auto y: g[x]) if (vis[y] == 0) {
				vis[y] = 1;
				if (rm[y] == -1 || dfs(dfs, rm[y])) {
					rm[y] = x;
					return 1;
				}
			}
			return 0;
		};
		dfs(dfs, i);
	}
	vi lm(n, -1);
	rep(i, 0, m - 1) if (rm[i] != -1) lm[rm[i]] = i;
	return lm;
}


int main() {
	ios::sync_with_stdio(0); cin.tie(0);
	
	int n; cin >> n;
	vector as(n - 1, vi(n));
	for (auto &vec : as) for (auto &x : vec) {
		cin >> x;
		x--;
	}
	auto Solve = [&](int s) -> void {
		vector<pii> es;
		vector<vector<int>> lsts(n - 1);
		rep(i, 0, n - 2) {
			auto &vec = as[i];
			int ind = find(all(vec), s) - vec.begin();
			rep(j, ind + 1, n - 1) {
				int v = vec[j];
				es.emplace_back(i, v);
			}
			lsts[i] = vector<int>(vec.begin() + ind + 1, vec.end());
		}
		auto lm = Kuhn(n - 1, n, es);
		if (count(all(lm), -1) > 0) {
			puts("No");
			return;
		} else {
			puts("Yes");
			vector<int> rm(n, -1);
			rep(i, 0, n - 2) rm[lm[i]] = i;
			vector<int> ords;
			vector<int> vis(n - 1);
			vector<int> used(n);
			vi sta;

			rep(i, 0, n - 2) if (vis[i] == 0) {
				auto dfs = [&](auto &dfs, int now) {
					vis[now] = 2;
					sta.push_back(now);
					while (1) {
						assert(sz(lsts[now]) > 0);
						int v = lsts[now].back();
						lsts[now].pop_back();
						if (used[v]) {
							continue;
						}
						int want = lm[now];
						if (want == v) {
							vis[now] = 1;
							used[v] = 1;
							sta.pop_back();
							ords.push_back(now);
							return;
						} else if (vis[rm[v]] == 2) {
							int u = rm[v];
							while (1) {
								int cur = sta.back();
								sta.pop_back();
								vis[cur] = 1;
								ords.push_back(cur);
								if (cur == u) break;
							}
							used[v] = 1;
							return;
						} else {
							assert(vis[rm[v]] == 0);
							dfs(dfs, rm[v]);
							if (vis[now] == 1) {
								used[v] = 1;
								return;
							}
						}
					}
				};
				dfs(dfs, i);
			}
			rep(i, 0, sz(ords) - 1) printf("%d%c", ords[i] + 1, " \n"[i == sz(ords) - 1]);
		}
	};
	rep(i, 0, n - 1) Solve(i);
	return 0;
}