A general matrix formula is proposed for the weight‐average molecular weights of the polymer systems formed through simultaneous scission, branching and crosslinking of N types of chains, assuming the chain connection statistics are Markovian. For the polymerization systems in which chains are generated consecutively, such as for free‐radical polymerization, the present theory can be applied by increasing the number of chain types N to infinity, by considering the chains formed at different times as different types of chains. The gel point determination reduces to the eigenvalue problem and the present theory extends the classical gelation theory to non‐random, history‐dependent reaction systems. From the mathematical point of view, this theory is capable of describing complex molecular build‐up processes through end‐linking, T‐ and H‐shaped chain connections, irrespective of reaction/reactor types used.