#include <bits/stdc++.h>
#define ff first
#define ss second
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

const int MOD = 1e9 + 7;
struct mint {
	int x;
	mint() : x(0) {}
	mint(int x) : x(x) {}
	mint operator+(int ot) const { return x + ot >= MOD ? x + ot - MOD : x + ot; }
	mint operator-(int ot) const { return x - ot < 0 ? x - ot + MOD : x - ot; }
	mint operator*(int ot) const { return 1ll * x * ot % MOD; }
	mint operator+=(int ot) { return *this = *this + ot; }
	mint operator-=(int ot) { return *this = *this - ot; }
	mint operator*=(int ot) { return *this = *this * ot; }
	operator int() const { return x; }
};
mint mpow(mint a, ll x) {
	mint r = 1;
	while (x) {
		if (x & 1) r *= a;
		a *= a;
		x /= 2;
	}
	return r;
}

mt19937 rnd(1557);
int rng() {
	return uniform_int_distribution<int>(0, MOD - 1)(rnd);
}
const int MAXN = 2020;
int n; mint A[MAXN][MAXN], C[MAXN][MAXN];
mint x[MAXN], y[MAXN], z[MAXN];

mint M[MAXN][MAXN], V[MAXN];
void gaussian(int n, int k) {
	// (n, k) * (k, 1) = (n, 1)
	for (int i = 0; i < k; i++) { // there are at most k pivots
		int p = -1;
		for (int j = i; j < n; j++) {
			if (M[j][i]) p = j;
		}
		for (int j = 0; j <= k; j++) swap(M[i][j], M[p][j]);

		mint inv = mpow(M[i][i], MOD - 2);
		for (int j = 0; j <= k; j++) M[i][j] *= inv;
		for (int j = 0; j < n; j++) if (i != j) {
			if (M[j][i]) {
				for (int l = k; l >= i; l--) M[j][l] -= M[j][i] * M[i][l];
			}
		}
	}
}
int main() {
	cin.tie(0); ios_base::sync_with_stdio(0);
	cin >> n;
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			int x; cin >> x;
			A[i][j] = x;
		}
	}
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			int x; cin >> x;
			C[i][j] = x;
		}
	}

	// find 12 "bad" rows and columns
	vector<int> BR, BC;
	for (int i = 0; i < n; i++) x[i] = rng();
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			y[i] += A[j][i] * x[j];
		}
	}
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			z[i] += C[j][i] * y[j];
		}
	}
	for (int i = 0; i < n; i++) if (x[i] != z[i]) BC.push_back(i);
	for (int i = 0; i < n; i++) y[i] = z[i] = 0;

	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			y[i] += A[i][j] * x[j];
		}
	}
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			z[i] += C[i][j] * y[j];
		}
	}
	for (int i = 0; i < n; i++) if (x[i] != z[i]) BR.push_back(i);

	vector<tuple<int, int, mint>> R;
	for (int y : BC) {
		for (int i = 0; i < n; i++) {
			for (int j = 0; j <= BR.size(); j++) M[i][j] = 0;
		}
		M[y][BR.size()] = 1;

		for (int i = 0; i < BR.size(); i++) {
			int x = BR[i];
			for (int j = 0; j < n; j++) {
				M[j][i] = A[j][x];
			}
		}
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				M[i][BR.size()] -= A[i][j] * C[j][y];
			}
		}

		// cout << "before\n";
		// for (int i = 0; i < n; i++) {
		// 	for (int j = 0; j <= BR.size(); j++) {
		// 		cout << M[i][j] << ' ';
		// 	}
		// 	cout << '\n';
		// }

		gaussian(n, BR.size());

		// cout << "after\n";
		// for (int i = 0; i < n; i++) {
		// 	for (int j = 0; j <= BR.size(); j++) {
		// 		cout << M[i][j] << ' ';
		// 	}
		// 	cout << '\n';
		// }

		for (int i = 0; i < BR.size(); i++) {
			if (M[i][BR.size()] != 0) {
				R.push_back({BR[i] + 1, y + 1, C[BR[i]][y] + M[i][BR.size()]});
			}
		}
	}
	cout << R.size() << '\n';
	for (auto [x, y, z] : R) {
		cout << x << ' ' << y << ' ' << z << '\n';
	}
}