Noncommutative multi-indices are noncommutative monomials in a $\mathbb{N}$-indexed family of indeterminates. We define on them a $\mathbb{Z}$-graded operadic structure, with the help of a shifting derivation. Multi-indices of degree 0 are called populated: they form a suboperad, isomorphic to the operad of Novikov algebras. This operadic structure, and the relation between pre-Lie and Novikov algebras, induces two bialgebraic structure in cointeraction on commutative multi-indices. We show how to combinatorially embed this double bialgebra into the Connes-Kreimer Hopf algebra of rooted trees, with its two coproducts based, firstly on cuts, secondly, on contraction of edges, and how this embedding can be characterized by a Dyson-Schwinger equation. We also study the unique polynomial invariant compatible with the two bialgebraic structures on multi-indices and use to describe the antipode for the first coproduct.