#include<iostream>
#include<vector>
#include<cmath>
#include<iomanip>
 
 
std::vector<std::vector<int>> mat_voigt_ind(const int N)
{
    /** this function will create pointer vector of matrix mat
     * in voigt notation
     */
 
    /** for every element A(i,j) of given
     * it will stores it index in the voigt vetor
     * given that matrix is square matrix of size N x N*/
    std::vector<std::vector<int>> mat = std::vector<std::vector<int>>(N, std::vector<int>(N, 0)) ;
 
 
    // starting and end index of the column to be filled
    std::size_t row_s = 0 ;
    std::size_t row_e = N - 1 ;
 
    std::size_t col_s = 0 ;
    std::size_t col_e = N - 1 ;
 
    /**we will fill diagonal first
     * then the last column
     * then first row
     * 
     * now next iteration
     */
 
    // starting index for matrix
    int i = row_s ;
    int j = col_s ;
 
    // k upto (N*N + N)/2 means we will run loop only for upper triangulat matrix
    for(std::size_t k = 0; k < (N*N + N)/2; k++)
    {
        mat[i][j] = k ;
 
        /** if non diagonal el, then add the el of lower triangle also */
        if(i != j)
            mat[j][i] = k + (N*N - N)/2 ;
 
        /** now update the index for matrix */
 
        // last column
        if(j == col_e and i > row_s)
            i-- ;
 
        // first row
        else if(i == row_s and j > col_s + 1)
            j-- ;
 
        // condition for next cycle
        else if(i == row_s and j == col_s + 1)
        {
            row_s += 1 ;
            row_e -= 2 ;
 
            col_s += 2 ;
            col_e -= 1 ;
 
            i = row_s ;
            j = col_s ;
        }
 
        // along the diagonal
        else
        {
            i++ ;
            j++ ;
        }
    }
 
    return mat ;
}
 
 
template <class T>
std::vector<std::vector<T>> dyadic(const std::vector<std::vector<T>>& A, const std::vector<std::vector<T>>& B)
{
    /** it will retun C as A dyadic B
     * given dimensions of A and B are equal
     */
 
    int N = A.size() ;
 
    std::vector<std::vector<T>> C(N*N, std::vector<float>(N*N,0)) ;  // N is the size of matrix A, i.e. which is calling the function
 
    /** if size of A and B are not equal */
    if(A.size() != B.size())
        return C ;
 
    /** generate a matrix of same size as A and B
     * which will store the index of voigt's matrix for that index of matrix
     */
    std::vector<std::vector<int>> voigt_mat = mat_voigt_ind(N) ;
 
    /** Cijkl = Aij * Bkl */
    for(int i= 0; i<N; i++)
        for(int j= 0; j<N; j++)
            for(int k= 0; k<N; k++)
                for(int l= 0; l<N; l++)
                    C[voigt_mat[i][j]][voigt_mat[k][l]] =  A[i][j] * B[k][l] ;
                    // C ij kl = A ij * B kl
 
    return C ;
}
 
template <class T>
std::vector<std::vector<T>> dyadic_up(const std::vector<std::vector<T>>& A, const std::vector<std::vector<T>>& B)
{
    /** it will retun C as A dyadic up B
     * given dimensions of A and B are equal
     */
 
    int N = A.size() ;
 
    std::vector<std::vector<T>> C(N*N, std::vector<float>(N*N,0)) ;  // N is the size of matrix A, i.e. which is calling the function
 
    /** if size of A and B are not equal */
    if(A.size() != B.size())
        return C ;
 
    /** generate a matrix of same size as A and B
     * which will store the index of voigt's matrix for that index of matrix
     */
    std::vector<std::vector<int>> voigt_mat = mat_voigt_ind(N) ;
 
    /** Cijkl = Aik * Bjl */
    for(int i= 0; i<N; i++)
        for(int j= 0; j<N; j++)
            for(int k= 0; k<N; k++)
                for(int l= 0; l<N; l++)
                    C[voigt_mat[i][j]][voigt_mat[k][l]] =  A[i][k] * B[j][l] ;
                    // C ij kl = A ik * B jl
 
    return C ;
}
 
template <class T>
std::vector<std::vector<T>> dyadic_down(const std::vector<std::vector<T>>& A, const std::vector<std::vector<T>>& B)
{
    /** it will retun C as A dyadic down B
     * given dimensions of A and B are equal
     */
 
    int N = A.size() ;
 
    std::vector<std::vector<T>> C(N*N, std::vector<float>(N*N,0)) ;  // N is the size of matrix A, i.e. which is calling the function
 
    /** if size of A and B are not equal */
    if(A.size() != B.size())
        return C ;
 
    /** generate a matrix of same size as A and B
     * which will store the index of voigt's matrix for that index of matrix
     */
    std::vector<std::vector<int>> voigt_mat = mat_voigt_ind(N) ;
 
    /** Cijkl = Aik * Bjl */
    for(int i= 0; i<N; i++)
        for(int j= 0; j<N; j++)
            for(int k= 0; k<N; k++)
                for(int l= 0; l<N; l++)
                    C[voigt_mat[i][j]][voigt_mat[k][l]] =  A[i][l] * B[j][k] ;
                    // C ij kl = A il * B jk
 
    return C ;
}
 
template <class T>
void printMatrix(const std::vector<std::vector<T>>& C, int N = 3)
{
    for(int i=0; i<C.size(); i++)
    {
        for(int j= 0; j<C[i].size(); j++)
        {
            std::cout<<std::setw(10) << std::fixed << std::setprecision(4)<<C[i][j] ;
 
            if((j+1)%N == 0)
                std::cout<<"\t" ;
        }
 
 
        if((i+1)%N == 0)
            std::cout<<std::endl ; 
 
        std::cout<<std::endl ; 
    }
}
 
 
int main()
{
    constexpr int N = 3 ;
 
    std::vector<std::vector<float>> A(N, std::vector<float>(N,0)) ;
    std::vector<std::vector<float>> B(N, std::vector<float>(N,0)) ;
 
 
    /** Problem no 5.3
     * considering alpha as 0
    */
    float alp = M_PI/2 ;
 
    // Rotation Matrix about V1
    A = { {1 , 0 , 0}, {0 , cos(alp), sin(alp)}, {0, -sin(alp), cos(alp)}} ;
 
    std::cout<<"\nTransformation Matrix for second order tensor for rotation about e1:"<<std::endl ;
    printMatrix(dyadic_up(A, A)) ;
 
    // Rotation Matrix about V2
    B = { {cos(alp), 0, -sin(alp)}, {0 , 1, 0}, {sin(alp), 0, cos(alp)}} ;
 
    std::cout<<"\nTransformation Matrix for second order tensor for rotation about e2:"<<std::endl ;
    printMatrix(dyadic_up(B, B)) ;
 
 
    /** for problem 5.4 */
    A = { {5, 3, 0}, {2, 4, 6}, {0, 6, 2}} ;
    B = { {4, 5, 2}, {0, 3, 4}, {0, 2, 5}} ;
 
    std::cout<<"\nA dyadic B:"<<std::endl ;
    printMatrix(dyadic(A, B)) ;
 
    std::cout<<"\nA dyadic up B:"<<std::endl ;
    printMatrix(dyadic_up(A, B)) ;
 
    return 0 ;
}
 
 

I2luY2x1ZGU8aW9zdHJlYW0+CiNpbmNsdWRlPHZlY3Rvcj4KI2luY2x1ZGU8Y21hdGg+CiNpbmNsdWRlPGlvbWFuaXA+CgoKc3RkOjp2ZWN0b3I8c3RkOjp2ZWN0b3I8aW50Pj4gbWF0X3ZvaWd0X2luZChjb25zdCBpbnQgTikKewogICAgLyoqIHRoaXMgZnVuY3Rpb24gd2lsbCBjcmVhdGUgcG9pbnRlciB2ZWN0b3Igb2YgbWF0cml4IG1hdAogICAgICogaW4gdm9pZ3Qgbm90YXRpb24KICAgICAqLwoKICAgIC8qKiBmb3IgZXZlcnkgZWxlbWVudCBBKGksaikgb2YgZ2l2ZW4KICAgICAqIGl0IHdpbGwgc3RvcmVzIGl0IGluZGV4IGluIHRoZSB2b2lndCB2ZXRvcgogICAgICogZ2l2ZW4gdGhhdCBtYXRyaXggaXMgc3F1YXJlIG1hdHJpeCBvZiBzaXplIE4geCBOKi8KICAgIHN0ZDo6dmVjdG9yPHN0ZDo6dmVjdG9yPGludD4+IG1hdCA9IHN0ZDo6dmVjdG9yPHN0ZDo6dmVjdG9yPGludD4+KE4sIHN0ZDo6dmVjdG9yPGludD4oTiwgMCkpIDsKCgogICAgLy8gc3RhcnRpbmcgYW5kIGVuZCBpbmRleCBvZiB0aGUgY29sdW1uIHRvIGJlIGZpbGxlZAogICAgc3RkOjpzaXplX3Qgcm93X3MgPSAwIDsKICAgIHN0ZDo6c2l6ZV90IHJvd19lID0gTiAtIDEgOwoKICAgIHN0ZDo6c2l6ZV90IGNvbF9zID0gMCA7CiAgICBzdGQ6OnNpemVfdCBjb2xfZSA9IE4gLSAxIDsKCiAgICAvKip3ZSB3aWxsIGZpbGwgZGlhZ29uYWwgZmlyc3QKICAgICAqIHRoZW4gdGhlIGxhc3QgY29sdW1uCiAgICAgKiB0aGVuIGZpcnN0IHJvdwogICAgICogCiAgICAgKiBub3cgbmV4dCBpdGVyYXRpb24KICAgICAqLwoKICAgIC8vIHN0YXJ0aW5nIGluZGV4IGZvciBtYXRyaXgKICAgIGludCBpID0gcm93X3MgOwogICAgaW50IGogPSBjb2xfcyA7CgogICAgLy8gayB1cHRvIChOKk4gKyBOKS8yIG1lYW5zIHdlIHdpbGwgcnVuIGxvb3Agb25seSBmb3IgdXBwZXIgdHJpYW5ndWxhdCBtYXRyaXgKICAgIGZvcihzdGQ6OnNpemVfdCBrID0gMDsgayA8IChOKk4gKyBOKS8yOyBrKyspCiAgICB7CiAgICAgICAgbWF0W2ldW2pdID0gayA7CgogICAgICAgIC8qKiBpZiBub24gZGlhZ29uYWwgZWwsIHRoZW4gYWRkIHRoZSBlbCBvZiBsb3dlciB0cmlhbmdsZSBhbHNvICovCiAgICAgICAgaWYoaSAhPSBqKQogICAgICAgICAgICBtYXRbal1baV0gPSBrICsgKE4qTiAtIE4pLzIgOwoKICAgICAgICAvKiogbm93IHVwZGF0ZSB0aGUgaW5kZXggZm9yIG1hdHJpeCAqLwoKICAgICAgICAvLyBsYXN0IGNvbHVtbgogICAgICAgIGlmKGogPT0gY29sX2UgYW5kIGkgPiByb3dfcykKICAgICAgICAgICAgaS0tIDsKICAgICAgICAKICAgICAgICAvLyBmaXJzdCByb3cKICAgICAgICBlbHNlIGlmKGkgPT0gcm93X3MgYW5kIGogPiBjb2xfcyArIDEpCiAgICAgICAgICAgIGotLSA7CgogICAgICAgIC8vIGNvbmRpdGlvbiBmb3IgbmV4dCBjeWNsZQogICAgICAgIGVsc2UgaWYoaSA9PSByb3dfcyBhbmQgaiA9PSBjb2xfcyArIDEpCiAgICAgICAgewogICAgICAgICAgICByb3dfcyArPSAxIDsKICAgICAgICAgICAgcm93X2UgLT0gMiA7CgogICAgICAgICAgICBjb2xfcyArPSAyIDsKICAgICAgICAgICAgY29sX2UgLT0gMSA7CgogICAgICAgICAgICBpID0gcm93X3MgOwogICAgICAgICAgICBqID0gY29sX3MgOwogICAgICAgIH0KCiAgICAgICAgLy8gYWxvbmcgdGhlIGRpYWdvbmFsCiAgICAgICAgZWxzZQogICAgICAgIHsKICAgICAgICAgICAgaSsrIDsKICAgICAgICAgICAgaisrIDsKICAgICAgICB9CiAgICB9CgogICAgcmV0dXJuIG1hdCA7Cn0KCgp0ZW1wbGF0ZSA8Y2xhc3MgVD4Kc3RkOjp2ZWN0b3I8c3RkOjp2ZWN0b3I8VD4+IGR5YWRpYyhjb25zdCBzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxUPj4mIEEsIGNvbnN0IHN0ZDo6dmVjdG9yPHN0ZDo6dmVjdG9yPFQ+PiYgQikKewogICAgLyoqIGl0IHdpbGwgcmV0dW4gQyBhcyBBIGR5YWRpYyBCCiAgICAgKiBnaXZlbiBkaW1lbnNpb25zIG9mIEEgYW5kIEIgYXJlIGVxdWFsCiAgICAgKi8KCiAgICBpbnQgTiA9IEEuc2l6ZSgpIDsKCiAgICBzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxUPj4gQyhOKk4sIHN0ZDo6dmVjdG9yPGZsb2F0PihOKk4sMCkpIDsgIC8vIE4gaXMgdGhlIHNpemUgb2YgbWF0cml4IEEsIGkuZS4gd2hpY2ggaXMgY2FsbGluZyB0aGUgZnVuY3Rpb24KCiAgICAvKiogaWYgc2l6ZSBvZiBBIGFuZCBCIGFyZSBub3QgZXF1YWwgKi8KICAgIGlmKEEuc2l6ZSgpICE9IEIuc2l6ZSgpKQogICAgICAgIHJldHVybiBDIDsKCiAgICAvKiogZ2VuZXJhdGUgYSBtYXRyaXggb2Ygc2FtZSBzaXplIGFzIEEgYW5kIEIKICAgICAqIHdoaWNoIHdpbGwgc3RvcmUgdGhlIGluZGV4IG9mIHZvaWd0J3MgbWF0cml4IGZvciB0aGF0IGluZGV4IG9mIG1hdHJpeAogICAgICovCiAgICBzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxpbnQ+PiB2b2lndF9tYXQgPSBtYXRfdm9pZ3RfaW5kKE4pIDsKCiAgICAvKiogQ2lqa2wgPSBBaWogKiBCa2wgKi8KICAgIGZvcihpbnQgaT0gMDsgaTxOOyBpKyspCiAgICAgICAgZm9yKGludCBqPSAwOyBqPE47IGorKykKICAgICAgICAgICAgZm9yKGludCBrPSAwOyBrPE47IGsrKykKICAgICAgICAgICAgICAgIGZvcihpbnQgbD0gMDsgbDxOOyBsKyspCiAgICAgICAgICAgICAgICAgICAgQ1t2b2lndF9tYXRbaV1bal1dW3ZvaWd0X21hdFtrXVtsXV0gPSAgQVtpXVtqXSAqIEJba11bbF0gOwogICAgICAgICAgICAgICAgICAgIC8vIEMgaWoga2wgPSBBIGlqICogQiBrbAoKICAgIHJldHVybiBDIDsKfQoKdGVtcGxhdGUgPGNsYXNzIFQ+CnN0ZDo6dmVjdG9yPHN0ZDo6dmVjdG9yPFQ+PiBkeWFkaWNfdXAoY29uc3Qgc3RkOjp2ZWN0b3I8c3RkOjp2ZWN0b3I8VD4+JiBBLCBjb25zdCBzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxUPj4mIEIpCnsKICAgIC8qKiBpdCB3aWxsIHJldHVuIEMgYXMgQSBkeWFkaWMgdXAgQgogICAgICogZ2l2ZW4gZGltZW5zaW9ucyBvZiBBIGFuZCBCIGFyZSBlcXVhbAogICAgICovCgogICAgaW50IE4gPSBBLnNpemUoKSA7CgogICAgc3RkOjp2ZWN0b3I8c3RkOjp2ZWN0b3I8VD4+IEMoTipOLCBzdGQ6OnZlY3RvcjxmbG9hdD4oTipOLDApKSA7ICAvLyBOIGlzIHRoZSBzaXplIG9mIG1hdHJpeCBBLCBpLmUuIHdoaWNoIGlzIGNhbGxpbmcgdGhlIGZ1bmN0aW9uCgogICAgLyoqIGlmIHNpemUgb2YgQSBhbmQgQiBhcmUgbm90IGVxdWFsICovCiAgICBpZihBLnNpemUoKSAhPSBCLnNpemUoKSkKICAgICAgICByZXR1cm4gQyA7CgogICAgLyoqIGdlbmVyYXRlIGEgbWF0cml4IG9mIHNhbWUgc2l6ZSBhcyBBIGFuZCBCCiAgICAgKiB3aGljaCB3aWxsIHN0b3JlIHRoZSBpbmRleCBvZiB2b2lndCdzIG1hdHJpeCBmb3IgdGhhdCBpbmRleCBvZiBtYXRyaXgKICAgICAqLwogICAgc3RkOjp2ZWN0b3I8c3RkOjp2ZWN0b3I8aW50Pj4gdm9pZ3RfbWF0ID0gbWF0X3ZvaWd0X2luZChOKSA7CgogICAgLyoqIENpamtsID0gQWlrICogQmpsICovCiAgICBmb3IoaW50IGk9IDA7IGk8TjsgaSsrKQogICAgICAgIGZvcihpbnQgaj0gMDsgajxOOyBqKyspCiAgICAgICAgICAgIGZvcihpbnQgaz0gMDsgazxOOyBrKyspCiAgICAgICAgICAgICAgICBmb3IoaW50IGw9IDA7IGw8TjsgbCsrKQogICAgICAgICAgICAgICAgICAgIENbdm9pZ3RfbWF0W2ldW2pdXVt2b2lndF9tYXRba11bbF1dID0gIEFbaV1ba10gKiBCW2pdW2xdIDsKICAgICAgICAgICAgICAgICAgICAvLyBDIGlqIGtsID0gQSBpayAqIEIgamwKCiAgICByZXR1cm4gQyA7Cn0KCnRlbXBsYXRlIDxjbGFzcyBUPgpzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxUPj4gZHlhZGljX2Rvd24oY29uc3Qgc3RkOjp2ZWN0b3I8c3RkOjp2ZWN0b3I8VD4+JiBBLCBjb25zdCBzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxUPj4mIEIpCnsKICAgIC8qKiBpdCB3aWxsIHJldHVuIEMgYXMgQSBkeWFkaWMgZG93biBCCiAgICAgKiBnaXZlbiBkaW1lbnNpb25zIG9mIEEgYW5kIEIgYXJlIGVxdWFsCiAgICAgKi8KCiAgICBpbnQgTiA9IEEuc2l6ZSgpIDsKCiAgICBzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxUPj4gQyhOKk4sIHN0ZDo6dmVjdG9yPGZsb2F0PihOKk4sMCkpIDsgIC8vIE4gaXMgdGhlIHNpemUgb2YgbWF0cml4IEEsIGkuZS4gd2hpY2ggaXMgY2FsbGluZyB0aGUgZnVuY3Rpb24KCiAgICAvKiogaWYgc2l6ZSBvZiBBIGFuZCBCIGFyZSBub3QgZXF1YWwgKi8KICAgIGlmKEEuc2l6ZSgpICE9IEIuc2l6ZSgpKQogICAgICAgIHJldHVybiBDIDsKCiAgICAvKiogZ2VuZXJhdGUgYSBtYXRyaXggb2Ygc2FtZSBzaXplIGFzIEEgYW5kIEIKICAgICAqIHdoaWNoIHdpbGwgc3RvcmUgdGhlIGluZGV4IG9mIHZvaWd0J3MgbWF0cml4IGZvciB0aGF0IGluZGV4IG9mIG1hdHJpeAogICAgICovCiAgICBzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxpbnQ+PiB2b2lndF9tYXQgPSBtYXRfdm9pZ3RfaW5kKE4pIDsKCiAgICAvKiogQ2lqa2wgPSBBaWsgKiBCamwgKi8KICAgIGZvcihpbnQgaT0gMDsgaTxOOyBpKyspCiAgICAgICAgZm9yKGludCBqPSAwOyBqPE47IGorKykKICAgICAgICAgICAgZm9yKGludCBrPSAwOyBrPE47IGsrKykKICAgICAgICAgICAgICAgIGZvcihpbnQgbD0gMDsgbDxOOyBsKyspCiAgICAgICAgICAgICAgICAgICAgQ1t2b2lndF9tYXRbaV1bal1dW3ZvaWd0X21hdFtrXVtsXV0gPSAgQVtpXVtsXSAqIEJbal1ba10gOwogICAgICAgICAgICAgICAgICAgIC8vIEMgaWoga2wgPSBBIGlsICogQiBqawoKICAgIHJldHVybiBDIDsKfQoKdGVtcGxhdGUgPGNsYXNzIFQ+CnZvaWQgcHJpbnRNYXRyaXgoY29uc3Qgc3RkOjp2ZWN0b3I8c3RkOjp2ZWN0b3I8VD4+JiBDLCBpbnQgTiA9IDMpCnsKICAgIGZvcihpbnQgaT0wOyBpPEMuc2l6ZSgpOyBpKyspCiAgICB7CiAgICAgICAgZm9yKGludCBqPSAwOyBqPENbaV0uc2l6ZSgpOyBqKyspCiAgICAgICAgewogICAgICAgICAgICBzdGQ6OmNvdXQ8PHN0ZDo6c2V0dygxMCkgPDwgc3RkOjpmaXhlZCA8PCBzdGQ6OnNldHByZWNpc2lvbig0KTw8Q1tpXVtqXSA7CgogICAgICAgICAgICBpZigoaisxKSVOID09IDApCiAgICAgICAgICAgICAgICBzdGQ6OmNvdXQ8PCJcdCIgOwogICAgICAgIH0KCiAgICAgICAgCiAgICAgICAgaWYoKGkrMSklTiA9PSAwKQogICAgICAgICAgICBzdGQ6OmNvdXQ8PHN0ZDo6ZW5kbCA7IAoKICAgICAgICBzdGQ6OmNvdXQ8PHN0ZDo6ZW5kbCA7IAogICAgfQp9CgoKaW50IG1haW4oKQp7CiAgICBjb25zdGV4cHIgaW50IE4gPSAzIDsKCiAgICBzdGQ6OnZlY3RvcjxzdGQ6OnZlY3RvcjxmbG9hdD4+IEEoTiwgc3RkOjp2ZWN0b3I8ZmxvYXQ+KE4sMCkpIDsKICAgIHN0ZDo6dmVjdG9yPHN0ZDo6dmVjdG9yPGZsb2F0Pj4gQihOLCBzdGQ6OnZlY3RvcjxmbG9hdD4oTiwwKSkgOwoKCiAgICAvKiogUHJvYmxlbSBubyA1LjMKICAgICAqIGNvbnNpZGVyaW5nIGFscGhhIGFzIDAKICAgICovCiAgICBmbG9hdCBhbHAgPSBNX1BJLzIgOwoKICAgIC8vIFJvdGF0aW9uIE1hdHJpeCBhYm91dCBWMQogICAgQSA9IHsgezEgLCAwICwgMH0sIHswICwgY29zKGFscCksIHNpbihhbHApfSwgezAsIC1zaW4oYWxwKSwgY29zKGFscCl9fSA7CgogICAgc3RkOjpjb3V0PDwiXG5UcmFuc2Zvcm1hdGlvbiBNYXRyaXggZm9yIHNlY29uZCBvcmRlciB0ZW5zb3IgZm9yIHJvdGF0aW9uIGFib3V0IGUxOiI8PHN0ZDo6ZW5kbCA7CiAgICBwcmludE1hdHJpeChkeWFkaWNfdXAoQSwgQSkpIDsKCiAgICAvLyBSb3RhdGlvbiBNYXRyaXggYWJvdXQgVjIKICAgIEIgPSB7IHtjb3MoYWxwKSwgMCwgLXNpbihhbHApfSwgezAgLCAxLCAwfSwge3NpbihhbHApLCAwLCBjb3MoYWxwKX19IDsKCiAgICBzdGQ6OmNvdXQ8PCJcblRyYW5zZm9ybWF0aW9uIE1hdHJpeCBmb3Igc2Vjb25kIG9yZGVyIHRlbnNvciBmb3Igcm90YXRpb24gYWJvdXQgZTI6Ijw8c3RkOjplbmRsIDsKICAgIHByaW50TWF0cml4KGR5YWRpY191cChCLCBCKSkgOwoKCiAgICAvKiogZm9yIHByb2JsZW0gNS40ICovCiAgICBBID0geyB7NSwgMywgMH0sIHsyLCA0LCA2fSwgezAsIDYsIDJ9fSA7CiAgICBCID0geyB7NCwgNSwgMn0sIHswLCAzLCA0fSwgezAsIDIsIDV9fSA7CgogICAgc3RkOjpjb3V0PDwiXG5BIGR5YWRpYyBCOiI8PHN0ZDo6ZW5kbCA7CiAgICBwcmludE1hdHJpeChkeWFkaWMoQSwgQikpIDsKCiAgICBzdGQ6OmNvdXQ8PCJcbkEgZHlhZGljIHVwIEI6Ijw8c3RkOjplbmRsIDsKICAgIHByaW50TWF0cml4KGR5YWRpY191cChBLCBCKSkgOwogICAgCiAgICByZXR1cm4gMCA7Cn0KCg==