  We consider the problem of calculating the large number of Wick contractions
necessary to compute states with the quantum numbers of many baryons in lattice
QCD. We consider a constructive approach and a determinant-based approach and
show that these methods allow the required contractions to be performed in
computationally manageable amount of time for certain choices of interpolating
operators. Examples of correlation functions computed using these techniques
are shown for the quantum numbers of the light nuclei, He, Be, C, O and Si.
