int binom_iter(int a, int b)
{
  int i;
  float e = 1.0f;
  for(i=b; i >= 0; i--)
      e *= (a-i)/(b-i);
  return (int) e;
}