The system is not time invariant. This is plainly seen from the coefficient of the $ \delta $ functional in the $ Y_k[n] $ column. As the input changes, the output changes as the square of the subscript in the input sequence.
