Out of curiosity, I was trying to generate an identity involving $\zeta(s)$ and $\zeta(s+1)$, for integers $s>1$, and after a lot of scribbling ended up with the following:
$(\zeta(s)-1)*(\zeta(s+1)-1)=\sum_{c=4}^\infty((\sum \frac 1f)*\frac 1{c^s})$, where the outer summation runs for all composites $c\geq4$ and the inner summation is for all the factors $f$ of $c$ (excluding $1$ and $c$ themselves).
Simply put, the entire RHS would look like: $(\frac 12)*\frac 1{4^s}+(\frac 12+\frac 13)*\frac 1{6^s}+(\frac 12+\frac 14)*\frac 1{8^s}+(\frac 13)*\frac 1{9^s}+(\frac 12 + \frac 15)*\frac 1{10^s}+(\frac 12 + \frac 13 + \frac 14 + \frac 16)*\frac 1{12^s}+...$
Now, before I ask my question, I must tell you that I don't have formal training or an advanced degree in mathematics. My question:
Can the above identity be derived starting only with either LHS or RHS? I ask this because my derivation involves construction of a variation of $\zeta$, doing it for all the numbers upto a number $N$, adding those constructions, and then taking $\lim_{N\to\infty}$ (I can share the entire derivation if asked). I also want to know if there is something obviously intuitive about this identity that I am missing (that led me to deriving it the long way).
Also if somebody could point me to some literature that deals with similar identities, I would be highly grateful. Thanks!
PS: I tested the identity with $s=2$ (therefore using $\zeta(2)=\frac {\pi^2}6$ and the value of $\zeta(3)$ from https://oeis.org/A002117 upto around a hundred decimal places), and with the composites $c$ upto 10,000 for the outer summation. The LHS and RHS values do match to 4 decimal places. I am sure adding more composites to the summation will give more accurate results.