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.
Comment: 24 pages, 2 figures