  We introduce the notion of "generalized bosons" whose exchange statistics
resemble those of bosons, but the local bosonic commutator
$[a_i,a_i^\dagger]=1$ is replaced by an arbitrary single-mode operator that is
diagonal in the generalized Fock basis. Examples of generalized bosons include
boson pairs and spins. We consider the analogue of the boson sampling task for
these particles and observe that its output probabilities are still given by
permanents, so that the results regarding hardness of sampling directly carry
over. Finally, we propose implementations of generalized boson sampling in
circuit-QED and ion-trap platforms.
