I am always struggling when it comes to dealing with maximum functions. I am trying to find the solution to this integral
$$\int_{\mathbb{Z}_p^3}||xy,xz,yz||_p^sd\mu(x,y,z),$$
where $||xy,xz,yz||_p^s=\sup\{|xy|_p^s,|xz|_p^s,|yz|_p^s\}$ and $d\mu$ is the normalized Haar measure.
Note: this is the same $p$-adic integration in Igusa local zeta function.
Ex: $\int_{\mathbb{z}_p}|x|_p^sd\mu(x)=\frac{1-p^{-1}}{1-p^{-s-1}}$
I'll solve this by a similar method I use in your other question: problem in $p$-adic integration
$$\int \max(|xy|, |xz|, |yz|)^s d\mu$$ Expand over each first digit,
$$\frac{1}{p^3} \sum_{a=0}^{p-1}\sum_{b=0}^{p-1}\sum_{c=0}^{p-1}\int \max(|(a+px)(b+py)|, |(a+px)(c+pz)|, |(b+py)(c+pz)|)^s d\mu$$
Now we split the sums between $0$ and $1$ to $p-1$, since $0$ as a first digit causes the p-adic absolute value to be smaller, the other digits cause it to be 1.
When there are no 0s or only one 0 as a digit, then the max is simply 1.
$$\frac{1}{p^3}\left((p-1)^3 + 3*(p-1)^2 + 3*\sum_{a=1}^{p-1}\int \max(|(a+px)py|, |(a+px)pz|, |p^2yz|)^s d\mu + \int\max(|p^2xy|, |p^2xz|, |p^2yz|)^s d\mu\right)$$
If we focus on just on this first integral in the sum, it simplifies substantially because $|a+px|=1$ $$\int \max(|py|, |pz|, |p^2yz|)^s d\mu$$ The 3rd term in the max is always smaller than the first two, and we can factor out a $p^{-s}$
$$p^{-s}\int \max(|y|, |z|)^s d\mu$$
This integral is exactly the one solved in the question linked at the start, so I'll put it here,
$$p^{-s}\frac{1-p^{-2}}{1-p^{-s-2}}$$
Now plugging that back in we have,
$$\frac{1}{p^3}\left((p-1)^3 + 3(p-1)^2 + 3p^{-s}\frac{1-p^{-2}}{1-p^{-s-2}} + p^{-2s}\int\max(|xy|, |xz|, |yz|)^s d\mu\right)$$
At this point we've got the original integral we started with! So all that's left is using the "snake eating it's own tail" trick. Call the integral J to make it cleaner to solve,
$$J = \frac{ (p-1)^3 + 3(p-1)^2}{p^3} + 3p^{-s-3}\frac{1-p^{-2}}{1-p^{-s-2}} + p^{-2s-3}J $$
$$J(1-p^{-2s-3}) = \frac{ (p-1)^3 + 3(p-1)^2}{p^3} + 3p^{-s-3}\frac{1-p^{-2}}{1-p^{-s-2}} $$
$$J= \frac{ (p-1)^3 + 3(p-1)^2}{p^3(1-p^{-2s-3}) } + 3p^{-s-3}\frac{1-p^{-2}}{(1-p^{-s-2})(1-p^{-2s-3}) } $$
It can be simplified more but I'll leave it here.