  We consider the link average of the half-BPS Wilson loop operators in N = 6
superconformal Chern-Simons-matter theory, which is called ABJM theory. We show
that this loop average is reduced to a (super)matrix integral by the
localization method, in a similar way to the bosonic U(N) Chern-Simons theory.
Using this matrix integral, we compute the two- and three-link averages with an
operator formalism inspired by a three-dimensional topological field theory. We
obtain a factorization of the link average, and the Verlinde formula in a
sector of supergroup representations. We also propose a refined version of ABJM
theory, and compute some refined link averages.
