We study the reaction kinetics of end-functionalized polymer chains dispersed in an unreactive polymer melt. Starting from an infinite hierarchy of coupled equations for many-chain correlation functions, a closed equation is derived for the 2nd order rate constant k after postulating simple physical bounds. Our results generalize previous 2-chain treatments (valid in dilute reactants limit) by Doi [1], de Gennes [2], and Friedman and O'Shaughnessy [3], to arbitrary initial reactive group density n 0 and local chemical reactivity Q. Simple mean field (MF) kinetics apply at short times, k ∼ Q. For high Q, a transition occurs to diffusioncontrolled (DC) kinetics with k ≈ x 3 t /t (where x t is rms monomer displacement in time t) leading to a density decay n t ≈ n 0 − n 2 0 x 3 t . If n 0 exceeds the chain overlap threshold, this behavior is followed by a regime where n t ≈ 1/x 3 t during which k has the same power law dependence in time, k ≈ x 3 t /t, but possibly different numerical coefficient. For unentangled melts this gives n t ∼ t −3/4 while for entangled cases one or more of the successive regimes n t ∼ t −3/4 , t −3/8 and t −3/4 may be realized depending on the magnitudes of Q and n 0 . Kinetics at times longer than the longest polymer relaxation time τ are always MF. If a DC regime has developed before τ then the long time rate constant is k ≈ R 3 /τ where R is the coil radius. We propose measuring the above kinetics in a model experiment where radical end groups are generated by photolysis.