for(b = e, c = f; d = b % c; b = c, c = d);

ef的最大公约数  c即答案