How confinement of quarks is implemented for multi-hadron systems in the quark cluster model is reviewed. In order to learn the nature of the confining interaction for fermions we first study 1+1 dimensional QED and QCD, in which the gauge field can be eliminated exactly and generates linear interaction of fermions. Then, we compare the two-body potential model, the flip-flop model and the Born-Oppenheimer approach in the strong coupling lattice QCD for the meson-meson system. Having shown how the long-range attraction between hadrons, van der Waals interaction, shows up in the two-body potential model, we discuss two distinct attempts beyond the two-body potential model: one is a many-body potential model, the flip-flop model, and the other is the Born-Oppenheimer approach in the strong coupling lattice QCD. We explain how the emergence of the long-range attraction is avoided in these attempts. Finally, we present the results of the application of the flip-flop model to the baryon-baryon scattering in the quark cluster model.