Solving gapped Hamiltonians locally

We show that any short-range Hamiltonian with a gap between the ground and excited states can be written as a sum of local operators, such that the ground state is an approximate eigenvector of each operator separately. We then show that the ground state of any such Hamiltonian is close to a generalized matrix product state. The range of the given operators needed to obtain a good approximation to the ground state is proportional to the square of the logarithm of the system size times a characteristic ``factorization length.'' Applications to many-body quantum simulation are discussed. We also consider density matrices of systems at non zero temperature.