#include<bits/stdc++.h>
using namespace std;
const int N=100;
int n;
char s[N+5][N+5];
char t[N+5][N+5];
struct oper{
    int t,x;
    oper(int t=0,int x=0):t(t),x(x){}
};
vector<oper> op;
vector<pair<int,int>> dif[2];
void work(int tp,int x,int c){
    for(int i=0;i<c;i++)
        op.push_back({tp,x});
    if(tp==1){
        static char tmp[N+5];
        for(int i=0;i<n;i++) tmp[i]=s[x][i];
        for(int i=0;i<n;i++) s[x][(i+c)%n]=tmp[i];
    }
    else{
        static char tmp[N+5];
        for(int i=0;i<n;i++) tmp[i]=s[i][x];
        for(int i=0;i<n;i++) s[(i+c)%n][x]=tmp[i];
    }
}
void rotate(int x1,int y1,int x2,int y2,int tp){
    if(tp==0){
        work(1,x1,(y2-y1+n)%n);
        work(2,y2,(x1-x2+n)%n);
        work(1,x1,(y1-y2+n)%n);
        work(2,y2,(x2-x1+n)%n);
    }
    else{
        work(2,y2,(x1-x2+n)%n);
        work(1,x1,(y2-y1+n)%n);
        work(2,y2,(x2-x1+n)%n);
        work(1,x1,(y1-y2+n)%n);
    }
}
int main(){
    scanf("%d",&n);
    for(int i=0;i<n;i++)
        scanf("%s",s[i]);
    for(int j=0;j<n;j++)
        scanf("%s",t[j]);
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            if(s[i][j]!=t[i][j]){
                dif[s[i][j]=='#'].push_back({i,j});
            }
    if(dif[0].size()!=dif[1].size()){
        puts("No");
        return 0;
    }
    for(int i=0;i<dif[0].size();i++){
        auto [x1,y1]=dif[0][i];
        auto [x2,y2]=dif[1][i];
        // printf("%d %d %d %d\n",x1,y1,x2,y2);
        if(x1!=x2&&y1!=y2){
            rotate(x1,y1,x2,y2,s[x1][y1]!=s[x1][y2]);
        }
        else if(x1!=x2){
            rotate(x1,(y1-1+n)%n,x2,y1,s[x1][(y1-1+n)%n]==s[x1][y1]);
        }
        else if(y1!=y2){
            rotate(x1,y1,(x1+1)%n,y2,s[(x2+1)%n][y1]!=s[x1][y1]);
        }
    }
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            assert(s[i][j]==t[i][j]);
    puts("Yes");
    printf("%d\n",op.size());
    for(auto [t,x]:op){
        printf("%d %d\n",t,x+1);
    }
    return 0;
}