Method 1:
For non-negative $u,v$ let $f(u)=u^{p-1}$ and $g(v)=v^{q-1}.$ Since $(p-1)(q-1)=1,$ the graphs of $(u,f(u))$ and of $(g(v),v)$ are the same graph, which I will call $G.$
Now $x^p/p=\int_0^xf(u)du$ is the area of the figure bounded by $G$ and $\{x\}\times [0,f(x)]$ and the $u$-axis in the $(u,v)$-plane. And $y^q/q=\int_0^yg(y)dy$ is the area of the figure bounded by $G$ and $[0, g(y)]\times \{y\}$ and the $v$-axis.
If $x^p=y^q,$ the union of these two figures is the rectangle $R(x,y)$ with corners $(0,0),(0,y),(x,y),(x,0),$ which has area $xy,$ so $x^p/p+y^q/q=\int_0^xf(u)du+\int_0^yg(v)dv=xy.$
If $x>y^{q/p}$ then the union of these two figures is $R(x,y)$ plus a region bounded by $G$, the $u$-axis, and the segments $\{y^{q/p}\}\times [0,y]$ and $\{x\}\times [0,f(x)],$ so $x^p/p+y^q/q=xy+\int_{y^{q/p}}^x f(u)du>xy$.
Similarly if $x<y^{p/q}$ then $y>x^{q/p}$ and $x^p/p+y^q/q=xy+\int_{x^{q/p}}^y g(v)dv>xy.$
This is easy to see on a diagram.
Generally, if $f:[0,\infty)\to [0,\infty)$ is continuous and strictly monotonic, and $f(0)=0,$ then $xy\leq \int_0^xf(u)du+\int_0^yf^{-1}(v)dv.$
Method 2:
For non-negative $x,y$ let $H(x,y)=x^p/p+y^q/q-xy.$ Then $\frac {\partial H(x,y)}{\partial x}=\frac1px^{p-1}-y,$ which is negative if $0\leq x^{p-1}<y$ and positive if $x^{p-1}>y.$ So for a fixed $y,$ the function $H(x,y)$ has its minimum when $x^{p-1}=y.$ And we have $H(x, x^{p-1})=0.$