#include<bits/stdc++.h>
using namespace std;

const int n=4;
int A[4][4];
bool vis[4][4];
const int dy[]={0, 1, 0,-1};
const int dx[]={1, 0,-1, 0};


bool dfs(int i, int j){
	if(vis[i][j]) return false;
	vis[i][j]=true;
	
	bool ok=true;
	for(int k=0;k<n;k++) for(int l=0;l<n;l++){
		ok &= A[k][l] == (k*4 + l+1) %16;
	}
	if(ok) return true;
	
	for(int d=0;d<4;d++){
		int yy= i + dy[d], xx= j + dx[d];
		if(yy<0 || yy>=n || xx<0 || xx>=n) continue;
		swap(A[i][j], A[yy][xx]);
		if(dfs(yy, xx)) return true;
		swap(A[i][j], A[yy][xx]);
	}
	
	vis[i][j]=false;
	return false;
}

int main(){
	while(1){
		for(int i=0;i<n;i++) for(int j=0;j<n;j++){
			if(!(cin >> A[i][j])) return 0;
		}
		memset(vis,0,sizeof(vis));
		for(int i=0;i<n;i++) for(int j=0;j<n;j++){
			if(A[i][j]==0){
				bool ans=dfs(i,j);
				cout << (ans? "Yes":"No") << endl;
				goto aa;
			}
		}
		aa:;
	}
	return 0;
}