Let $A, B$ be given $n\times n$ matrices with positive entries. Does the matrix equation $XAY=B$ always has a solution $X, Y$ with nonnegative entries?
I tried to use Kronecker product to transform the equation as $(Y^\top\otimes X)\mathrm{vec}\,A=\mathrm{vec}\,B$... how to proceed then? Thanks.