$$M=y(4lnx + 4lny +1)=y(4ln(xy)+1),N=x(lnx+lny+1)=x(ln(xy)+1)$$
$$M_y =4ln(xy) +5, N_x=ln(xy)+2$$
Since $M_y\ne N_x$, it is not exact.
$$R=\frac{M_y-N_x}{N}=\frac{3ln(xy)+3}{x(ln(xy)+1)} = \frac{3(ln(xy)+1)}{x(ln(xy)+1)}=\frac{3}{x}$$
R is a function of x only, so there is an integrating factor of the form
$$\mu=e^{\int R dx}=e^{\int \frac{3}{x} dx}=e^{3lnx }=x^3$$
multiply the differential equation by $\mu$
$$x^3y(4ln(xy) +1)+x^4(ln(xy)+1)y'=0$$
Since $M_y = N_x$, It is exact now. There exist a $\psi(x,y)$ such that $\psi_x=M, \psi_y=N$ .
$$\psi=\int M dx=\int x^3y(4ln(xy) +1)dx=x^4yln(xy)+h(y)$$
$$\psi_y=x^4ln(xy)+x^4+h'(y)=x^4(ln(xy)+1)=N$$
$$∴h'(y)=0, h(y)=C$$
$$Therefore, \psi(x,y)= x^4yln(xy)=C$$
as x=1, y=1, we get C=0
$$∴x^4yln(xy)=0$$