#include<bits/stdc++.h>
using namespace std;
#define md(a) a=(a%mod+mod)%mod
#define file(a) freopen(#a".in","r",stdin);freopen(#a".out","w",stdout)

bool ST;

const int N=2005,mod=1e9+7;
namespace Math
{
	inline void add(int &x,const int &y){x+=y;if(x>=mod)x-=mod;}
	inline void dec(int &x,const int &y){x-=y;if(x<0)x+=mod;}
	inline int qp(int a,int n){int b=1;while(n){if(n&1)b=1ll*b*a%mod;a=1ll*a*a%mod,n>>=1;}return b;}
};
using namespace Math;int ned,dir[N],dic[N],dif[N][N];
int n,a[N][N],b[N][N],vec[N],row[N],col[N];mt19937 R(time(0));

inline void rmul(int*a,int b[N][N])
{
	static int c[N];
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)add(c[j],1ll*a[i]*b[i][j]%mod);
	for(int i=1;i<=n;i++)a[i]=c[i],c[i]=0;
}

inline void lmul(int b[N][N],int*a)
{
	static int c[N];
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)add(c[j],1ll*a[i]*b[j][1]%mod);
	for(int i=1;i<=n;i++)a[i]=c[i],c[i]=0;
}

int f[N][14];
inline void calc(int c)
{
	if(!ned)return ;
	for(int i=1;i<=n;i++)
	{
		f[i][ned+1]=(i==c);
		for(int j=1,k=1;j<=n;j++)
		{
			if(k<=ned&&dir[k]==j)f[i][k]=a[i][j],k++;
			else dec(f[i][ned+1],1ll*a[i][j]*b[j][c]%mod);
		}
	}
	for(int i=1;i<=min(n,ned);i++)
	{
		for(int j=i;j<=n;j++)if(f[j][i]){for(int k=1;k<=ned+1;k++)swap(f[i][k],f[j][k]);break;}
		int iv=qp(f[i][i],mod-2);
		for(int j=1;j<=n;j++)if(i!=j)
		{
			int w=1ll*f[j][i]*iv%mod;
			for(int k=1;k<=ned+1;k++)dec(f[j][k],1ll*f[i][k]*w%mod);
		}
	}
	for(int i=1;i<=ned;i++)dif[dir[i]][c]=1ll*f[i][ned+1]*qp(f[i][i],mod-2)%mod;
}

bool ED;

signed main()
{
	int time_st=clock();
	cerr<<(&ST-&ED)/1024.0/1024<<endl;ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	
	cin>>n;for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)cin>>a[i][j];
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)cin>>b[i][j],dif[i][j]=b[i][j];
	for(int i=1;i<=n;i++)vec[i]=R(),row[i]=col[i]=vec[i];rmul(row,a),rmul(row,b),lmul(a,col),lmul(b,col);
	for(int i=1,p=0,q=0;i<=n;i++){if(row[i]!=vec[i])dic[++q]=i;if(col[i]!=vec[i])dir[++p]=i;ned=max(p,q);}
	
	if(n>=1000)return 0;

	for(int i=1;i<=ned;i++)calc(dic[i]);

	vector<pair<int,int> >ans;for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)if(dif[i][j]!=b[i][j])ans.push_back({i,j});
	cout<<(int)ans.size()<<'\n';for(auto pa:ans){int x=pa.first,y=pa.second;cout<<x<<' '<<y<<' '<<dif[x][y]<<'\n';}

	cerr<<(clock()-time_st)/1e6<<endl;return 0;
}