Taking advantage of the simplicity of the Hubbard or Pariser Parr Pople Hamiltonian, it is possible to define a multi-reference few-state coupled cluster expansion. The reference functions are all the most neutral valence-bond (VB) determinants and the targets are the lowest eigenstates. Test calculations show that the method gives reliable results even when the \t/U\ ratio is as large as 1, and provides excellent additional roots.