  We prove an explicit integral formula for computing the product of two
shifted Riemann zeta functions everywhere in the complex plane. We show that
this formula implies the existence of infinite families of exact exponential
sum identities involving the divisor functions, and we provide examples of
these identities. We conjecturally propose a method to compute divisor
functions by matrix inversion, without employing arithmetic techniques.
