MOD = 10**9 + 7

def multiply_matrix(A, B):
    res = [[0] * 3 for _ in range(3)]
    for i in range(3):
        for j in range(3):
            for k in range(3):
                res[i][j] += A[i][k] * B[k][j]
                res[i][j] %= MOD
    return res

def matrix_power(mat, power):
    result = [[1 if i == j else 0 for j in range(3)] for i in range(3)]
    while power > 0:
        if power % 2 == 1:
            result = multiply_matrix(result, mat)
        mat = multiply_matrix(mat, mat)
        power = power // 2
    return result

def multiply_vector(mat, vec):
    res = [0, 0, 0]
    for i in range(3):
        res[i] = (mat[i][0] * vec[0] + mat[i][1] * vec[1] + mat[i][2] * vec[2]) % MOD
    return res

def main():
    import sys
    input = sys.stdin.read().split()
    N = int(input[0])
    a = int(input[1]) % MOD
    b = int(input[2]) % MOD
    c = int(input[3]) % MOD
    vec = [a, b, c]
    
    if N == 1:
        print(a, b, c)
        return
    
    M = [
        [1, (MOD - 1) % MOD, 0],
        [0, 1, (MOD - 1) % MOD],
        [(MOD - 1) % MOD, 0, 1]
    ]
    
    power = N - 1
    M_exp = matrix_power(M, power)
    final = multiply_vector(M_exp, vec)
    print(final[0], final[1], final[2])

if __name__ == '__main__':
    main()