So, n is prime iff M_1(n)=n+1. That's much simpler than the first equation listed there!
Indeed, looking things up, it seems that in general the functions M_a can be written as a linear combination (note: with polynomial coefficients, not constant) of the sigma_k (sigma_k is the sum of the k'th power of the divisors). So this result becomes a lot less surprising once you know that...