// Standard I/O
#include <iostream>
#include <sstream>
#include <cstdio>
// Standard Library
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <cmath>
// Template Class
#include <complex>
#include <string>
#include <vector>
#include <list>
#include <set>
#include <map>
#include <queue>
#include <stack>
// Container Control
#include <algorithm>

using namespace std;

#define rep( i, n ) for( int i = 0; i < n; ++i )
#define irep( i, n ) for( int i = n-1; i >= 0; --i )
#define reep( i, s, n ) for ( int i = s; i < n; ++i )
#define ireep( i, n, s ) for ( int i = n-1; i >= s; --i )
#define foreach(itr, x) for( typeof(x.begin()) itr = x.begin(); itr != x.end(); ++itr)

#define mp( a, b ) make_pair( a, b )
#define pb( a ) push_back( a )
#define all( v ) v.begin(), v.end()
#define fs first
#define sc second
#define vc vector

// for visualizer.html
double SCALE = 1.0;
double OFFSET_X = 0.0;
double OFFSET_Y = 0.0;
#define LINE(x,y,a,b) cerr << "line(" << SCALE*(x) + OFFSET_X << ","	\
	<< SCALE*(y) + OFFSET_Y << ","										\
	<< SCALE*(a) + OFFSET_X << ","										\
	<< SCALE*(b) + OFFSET_Y << ")" << endl;
#define CIRCLE(x,y,r) cerr << "circle(" << SCALE*(x) + OFFSET_X << ","	\
	<< SCALE*(y) + OFFSET_Y << ","										\
	<< SCALE*(r) << ")" << endl;

typedef long long ll;
typedef complex<double> Point;

typedef pair<int, int> pii;
typedef pair<int, pii> ipii;
typedef vector<int> vi;
typedef vector<double> vd;
typedef vector< vector<int> > vii;
typedef vector< vector<double> > vdd;

typedef vector<int>::iterator vi_itr;

const int IINF = 1 << 28;
const double INF = 1e30;
const double EPS = 1e-10;
const double PI = acos(-1.0);

// Direction : L U R D
const int dx[] = { -1, 0, 1, 0};
const int dy[] = { 0, -1, 0, 1 };

int a[4][4], field[4][4];
bool used[16];

bool check ()
{
	bool ret = true;
	rep(i, 4) rep(j, 4)
		if( a[i][j] != field[i][j] ){
			ret = false;
			break;
		}
	return ret;
}

bool solve ( int x, int y )
{
	if( check() ) return true;

	bool ret = false;
	rep(i, 4){
		int nx = x + dx[i];
		int ny = y + dy[i];
		if( nx < 0 || nx >= 4 || ny < 0 || ny >= 4 ) continue;
		if( field[ny][nx] == 0 || used[field[ny][nx]] ) continue;
		
		used[field[ny][nx]] = true;
		swap(field[y][x], field[ny][nx]);
		ret |= solve(nx, ny);
		swap(field[y][x], field[ny][nx]);
		used[field[ny][nx]] = false;

		if( ret ) break;
	}
	return ret;
}

int main()
{
	pii p;
	rep(i, 4) rep(j, 4){
		field[i][j] = (i*4+j+1)%16;
		cin >> a[i][j];
	}

	cout << (solve(3, 3) ? "Yes" : "No") << endl;
}