#include <iostream>
using namespace std;

bool inside(int a, int b, int alim, int blim){
	return 0 <= a && a < alim && 0 <= b && b < blim;
}

int main(){
	int a[4][4], dd[] = {0, 1, 0, -1, 0}, d;
	for(int i = 0; i < 4; i++) for(int j = 0; j < 4; j++){
		cin >> a[i][j];
		bool f = a[i][j] == (4 * i + j + 1) % 16;
		for(int k = 0; k < 4; k++){
			int ni = i + dd[k], nj = j + dd[k + 1];
			if(inside(ni, nj, 4, 4)) f |= a[i][j] == (4 * ni + nj + 1) % 16;
		}
		if(!f){
			cerr << i << ' ' << j << endl;
			cout << "No" << endl;
			return 0;
		}
		if(a[i][j] == 0) d = 3 - i + 3 - j;
	}
	int count = 0;
	for(int i = 0; i < 4; i++) for(int j = 0; j < 4; j++){
		int p = (4 * i + j + 1) % 16;
		if(a[i][j] == p) continue;
		for(int k = 0; k < 4; k++){
			int ni = i + dd[k], nj = j + dd[k + 1];
			if(inside(ni, nj, 4, 4) && a[ni][nj] == p){
				swap(a[i][j], a[ni][nj]);
				count++;
				break;
			}
		}
	}
	cout << (d % 2 == count % 2 ? "Yes" : "No") << endl;
}