Automation of next-to-leading order computations in QCD: the FKS subtraction