I was toying around with the Riemann Zeta function recently and noticed that I could get to a particular representation (valid for $Re(s)>1$) in a couple of different odd ways.
The first was by looking at perfect powers and non-perfect powers (by which I mean positive integers that are $>1$ and not equal to any perfect power). Here perfect powers are defined to be any integers of the form $a^b$ where $a,b\in \mathbb Z$ and $a,b>2$ (as a consequence, $1$ is neither perfect nor non-perfect).
It is clear that any positive integer $>1$ can be represented uniquely as some $a^k$ where $a$ is a non-perfect power and $k\in \mathbb Z^+$.
Using this fact, we can rewrite the traditional series representation of the Zeta function to get: $$\zeta(s)=\sum_{m=1}^{\infty}\frac{1}{m^s}=1+\sum_{r=1}^{\infty}\sum_{n=1}^{\infty}\frac{1}{a_r^{ns}}$$
where $a_r$ denotes the $r$th non-perfect number. One notes that the inner sum is a simple geometric sum- it may be evaluated to give $$\zeta(s)=1+\sum_{r=1}^{\infty}\bigg(\frac{1}{a_r^s-1}\bigg)$$
Now, we use the fact that each integer $>1$ is either a perfect power or a non-perfect power to give $$\zeta(s)=1+\sum_{m=2}^{\infty}\bigg(\frac{1}{m^s-1}\bigg)-\sum_{r=1}^{\infty}\bigg(\frac{1}{p_r^s-1}\bigg)$$
where $p_r$ denotes the $r$th perfect power. Finally, we note that every perfect power may be expressed uniquely as some $a^k$ where $a$ is a non-perfect power and $k$ is an integer $\geq 2$. Using this gives us $$\zeta(s)=1+\sum_{m=2}^{\infty}\bigg(\frac{1}{m^s-1}\bigg)-\sum_{n=2}^{\infty}\sum_{r=1}^{\infty}\bigg(\frac{1}{a_r^{ns}-1}\bigg)$$
But now we can apply the same tricks again on the inner sum, i.e. $$\zeta(s)=1+\sum_{m=2}^{\infty}\bigg(\frac{1}{m^s-1}\bigg)-\sum_{n=2}^{\infty}\Bigg(\sum_{m=2}^{\infty}\bigg(\frac{1}{m^{ns}-1}\bigg)-\sum_{r=1}^{\infty}\bigg(\frac{1}{p_r^{ns}-1}\bigg)\Bigg)\\=1+\sum_{m=2}^{\infty}\bigg(\frac{1}{m^s-1}\bigg)-\sum_{n=2}^{\infty}\sum_{m=2}^{\infty}\bigg(\frac{1}{m^{ns}-1}\bigg)+\sum_{n=2}^{\infty}\sum_{r=1}^{\infty}\bigg(\frac{1}{p_r^{ns}-1}\bigg)\\=1+\sum_{m=2}^{\infty}\bigg(\frac{1}{m^s-1}\bigg)-\sum_{n=2}^{\infty}\sum_{m=2}^{\infty}\bigg(\frac{1}{m^{ns}-1}\bigg)+\sum_{n_1=2}^{\infty}\sum_{n_2=2}^{\infty}\sum_{r=1}^{\infty}\bigg(\frac{1}{a_r^{n_1n_2s}-1}\bigg)$$
We can apply these couple of tricks repeatedly to get $$\zeta(s)=1+\sum_{m=2}^{\infty}\bigg(\frac{1}{m^s-1}\bigg)-\sum_{n_1=2}^{\infty}\sum_{m=2}^{\infty}\bigg(\frac{1}{m^{n_1s}-1}\bigg)+\sum_{n_1=2}^{\infty}\sum_{n_2=2}^{\infty}\sum_{m=2}^{\infty}\bigg(\frac{1}{m^{n_1n_2s}-1}\bigg)-...$$
But, there's another way we can get to this expansion- by using the Mobius function.
Since $\sum_{d:d|n}\mu(d)=1$ if $n=1$ and is $0$ otherwise, for any function $f$, we have $$f(s)=\sum_{i=1}^{\infty}\bigg(f(is)\sum_{d:d|i}\mu(d)\bigg)$$
If we rearrange the sum by collecting all the like $\mu(d)$ terms we get
$$f(s)=\sum_{d=1}^{\infty}\sum_{j=1}^{\infty}\bigg(\mu(d)f((dj)s)\bigg)$$
When we let $f(s)=\zeta(s)-1$, we get $$\zeta(s)-1=\sum_{d=1}^{\infty}\sum_{j=1}^{\infty}\bigg(\mu(d)(\zeta(djs)-1)\bigg)=\sum_{d=1}^{\infty}\mu(d)\sum_{j=1}^{\infty}\sum_{m=2}^{\infty}\frac{1}{m^{djs}}=\sum_{d=1}^{\infty}\mu(d)\sum_{m=2}^{\infty}\frac{1}{m^{ds}-1}$$
Then, using this expansion of the Mobius function as the sum of the number of integer points on an n-dimensional hyperboloid and adding $1$ to both sides of the equation the same series may be obtained: $$\zeta(s)=1+\sum_{m=2}^{\infty}\bigg(\frac{1}{m^s-1}\bigg)-\sum_{n_1=2}^{\infty}\sum_{m=2}^{\infty}\bigg(\frac{1}{m^{n_1s}-1}\bigg)+\sum_{n_1=2}^{\infty}\sum_{n_2=2}^{\infty}\sum_{m=2}^{\infty}\bigg(\frac{1}{m^{n_1n_2s}-1}\bigg)-...$$
My question is, why do these different approaches (one, involving the Mobius function and, the other, using the basic properties of perfect powers) converge to the same result? Is there any way to remove the $\zeta$ function from all of this and get a clear-cut relationship between the Mobius function and the perfect/non-perfect numbers (in particular, could examining the error term of the expansion yield any results) or is this equality of approaches fundamentally dependent on the $\zeta$ function's properties? Or perhaps both approaches distract from a simpler possible derivation?
(please comment or edit for any corrections)