#include <bits/stdc++.h>

#define MOD 1000000007
#define INF (1<<30)
#define INFL (1<<62)
#define pe(str) return cout<<(str)<<endl,0
#define px(str) cout<<(str)<<endl,exit(0)
#define accu(ac,v) cout<<setprecision(ac)<<v<<endl
#define bpc(a) __builtin_popcount(a)
#define pr(str) cout<<(str)<<endl
#define prpr(str) cout<<str<<endl
#define db(str) cerr<<"debug:"<<(str)<<endl
#define dbdb(str) cerr<<"debugs:"<<str<<endl
#define re(i,n) for(int i=0;i<(n);i++)
#define rre(i,n) for(int i=(n);i>=0;i--)
#define rep(i,a,b) for(int i=(a);i<(b);i++)
#define rrep(i,a,b) for(int i=(b);i>=(a);i--)
#define bw(a,b,c) (((a)<=(b))&&((b)<=(c)))
#define hello cout<<"hello"<<endl
#define spre(a) setprecision(a)
#define all(a) (a).begin(),(a).end()
#define rall(a) (a).rbegin(),(a).rend()
#define IF(C,T,F) ((C)?(T):(F))
#define max3(a,b,c) max((a),max((b),(c)))
#define min3(a,b,c) min((a),min((b),(c)))
#define max4(a,b,c,d) max((a),max((b),max((c),(d))))
#define min4(a,b,c,d) min((a),min((b),min((c),(d))))
#define enl cout<<endl
#define declare(t,n) t n;cin>>n
#define EPS 1e-6
#define EPSIN(a,b) ((b)-EPS<=(a)&&(a)<=(b)+EPS)

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef pair<int,int> pii;
typedef complex<int> point;

template<class InputIterator>
void dump(InputIterator first,InputIterator last,char delim=' '){
	for(InputIterator it=first;it!=last;it++){
		if(it!=first)cout<<delim;
		cout<<*it;
	}
}

ull p(ull x,ull y){
	if(y==0) return 1;
	ull A=p(x,y/2)%MOD;
	return (y&1?((A*A)%MOD*x)%MOD: (A*A)%MOD);
	// return (y&1?((A*A)*x):(A*A));
}

template<class T>
void swp(T &a,T &b){
	T t=a; a=b; b=t;
}

int c[4][4];
int s=0;
int x,y;

void dfs(int i,int j){
	if(!(bw(0,i,4)&&bw(0,j,4)))return;
	if(c[i][j]!=1)return;
	c[i][j]=0;
	dfs(i+1,j);
	dfs(i-1,j);
	dfs(i,j-1);
	dfs(i,j+1);
}

int main(void){
	int a[4][4];
	int b[4][4]={ {1,2,3,4},{5,6,7,8},{9,10,11,12},{13,14,15,0} };

	re(i,4) re(j,4) cin>>a[i][j];
	re(i,4){
		re(j,4){
			re(k,4){
				re(l,4){
					if(a[i][j]==b[k][l]){
						c[i][j]=abs(i-k)+abs(j-l);
						if(a[i][j]!=0&&abs(i-k)+abs(j-l)>=2)pe("No");
						if(c[i][j]==1)s++;
						if(a[i][j]==0){x=i;y=j;c[i][j]=1;}
						goto B;
					}
				}
			}
B:;
		}
	}
	// re(i,4){
		// dump(c[i],c[i]+4);
		// cout<<endl;
	// }
	if(s==0)pe("Yes");
	if(s!=0&&c[x][y]==0)pe("No");
	dfs(x,y);
	re(i,4) re(j,4) if(c[i][j]==1)pe("No");
	pe("Yes");
	return 0;
}