整型求法
int gcd(int x,int y) { return b?gcd(b,a%b):a; }
浮点型求法
1 double gcd(double x,double y) { 2 while(fabs(x) > esp && fabs(y) > esp) { 3 if(x > y) 4 x -= floor(x / y) * y; 5 else 6 y -= floor(y / x) * x; 7 } 8 return x+y; 9 }
原文：http://www.cnblogs.com/zmin/p/7622116.html