  We examine the renormalization of the first order formulation of Yang-Mills
theory, by using the BRST idenities. These preserve the gauge invariance of the
theory and enable a recursive proof of renormalizability to higher orders in
perturbation theory. The renormalization involves non-linear mixings as well as
re-scalings of the fields and sources, which lead to a renormalized action at
all orders.
