The multipartite quantum systems are of particular interest for the study of such phenomena as entanglement and non-local correlations. The symmetry group of the whole multipartite system is the wreath product of the group acting in the “local” Hilbert space and the group of permutations of the constituents. The dimension of the Hilbert space of a multipartite system depends exponentially on the number of constituents, which leads to computational difficulties. We describe an algorithm for decomposing representations of wreath products into irreducible components. The C implementation of the algorithm copes with representations of dimensions in quadrillions. The program, in particular, builds irreducible invariant projectors in the Hilbert space of a multipartite system. The expressions for these projectors are tensor product polynomials. This structure is convenient for efficient computation of quantum correlations in multipartite systems with a large number of constituents.