#include <bits/stdc++.h>
#define INIT std::ios::sync_with_stdio(false);std::cin.tie(0);
// VAR(int, x);
#define VAR(type, ...)type __VA_ARGS__;Scan(__VA_ARGS__);
template<typename T> void Scan(T& t){std::cin >> t;}
template<typename First,typename...Rest>void Scan(First& first,Rest&...rest){std::cin>>first;Scan(rest...);}
#define OUT(d) std::cout<<(d);
#define FOUT(n, d) std::cout<<std::fixed<<std::setprecision(n)<<(d);
#define SOUT(n, c, d) std::cout<<std::setw(n)<<std::setfill(c)<<(d);
#define SP std::cout<<" ";
#define TAB std::cout<<"\t";
#define BR std::cout<<"\n";
#define ENDL std::cout<<std::endl;
#define FLUSH std::cout<<std::flush;
#define VEC(type, c, n) std::vector<type> c(n);for(auto& i:c)std::cin>>i;
#define MAT(type, c, m, n) std::vector<std::vector<type>> c(m, std::vector<type>(n));for(auto& r:c)for(auto& i:r)std::cin>>i;
#define ALL(a) (a).begin(),(a).end()
#define FOR(i, a, b) for(int i=(a);i<(b);++i)
#define RFOR(i, a, b) for(int i=(b)-1;i>=(a);--i)
#define REP(i, n) for(int i=0;i<int(n);++i)
#define RREP(i, n) for(int i=(n)-1;i>=0;--i)
#define PAIR std::pair<int, int>
#define IN(a0, y, a1, b0, x, b1) (a0<=y && y<a1 && b0<=x && x<b1)
#define SHOW(d) {std::cout << #d << "\t:" << d << "\n";}
#define SHOWVECTOR(v) {std::cout << #v << "\t:";for(const auto& i : v){std::cout << i << " ";}std::cout << "\n";}
#define SHOWVECTOR2(v) {std::cout << #v << "\t:\n";for(const auto& i : v){for(const auto& j : i){std::cout << j << " ";}std::cout << "\n";}}
#define SHOWPAIRVECTOR(v) {for(const auto&& i:v){std::cout<<'('<<i.first<<", "<<i.second<<") ";}std::cout<<"\n";}
#define CHECKTIME(state) {auto start=std::chrono::system_clock::now();state;auto end=std::chrono::system_clock::now();auto res=std::chrono::duration_cast<std::chrono::nanoseconds>(end-start).count();std::cerr<<"[Time:"<<res<<"ns  ("<<res/(1.0e9)<<"s)]\n";}
 
//#define int ll
using ll = long long;
constexpr ll INF = 1LL<<60;
constexpr double EPS = 0.0000001;

signed main(){
    INIT;
    CHECKTIME(([&]
    {
        MAT(int, a, 4, 4);
        std::vector<std::vector<int>> base(4, std::vector<int>(4));
        REP(i, 4)REP(j, 4){
            base[i][j] = i*4+j+1;
        }
        base[3][3] = 0;
        std::vector<std::vector<bool>> t(4, std::vector<bool>(4, false));
        int x = 3, y = 3;
        bool conti;
        while(true){
            conti = false;
            if(a[y][x] == 0) break;
            if(y!=0 && base[y-1][x] == a[y][x]){
                if(t[y-1][x]) break;
                std::swap(base[y-1][x], base[y][x]);
                t[y][x] = true;
                --y;
                conti = true;
            }
            if(y!=3 && base[y+1][x] == a[y][x]){
                if(t[y+1][x]) break;
                std::swap(base[y+1][x], base[y][x]);
                t[y][x] = true;
                ++y;
                conti = true;
            }
            if(x!=0 && base[y][x-1] == a[y][x]){
                if(t[y][x-1]) break;
                std::swap(base[y][x-1], base[y][x]);
                t[y][x] = true;
                --x;
                conti = true;
            }
            if(x!=3 && base[y][x+1] == a[y][x]){
                if(t[y][x+1]) break;
                std::swap(base[y][x+1], base[y][x]);
                t[y][x] = true;
                ++x;
                conti = true;
            }
            if(!conti) break;
        }
        OUT(a==base?"Yes":"No")BR;
    }
    ()));
    return 0;
}