虽然不能用__gcd()

但是可以去stl_algo头文件里找这个函数怎么写的

这个头文件在<algorithm>里

template<typename _EuclideanRingElement>
    _EuclideanRingElement
    __gcd(_EuclideanRingElement __m, _EuclideanRingElement __n)
    {
      while (__n != 0)
	{
	  _EuclideanRingElement __t = __m % __n;
	  __m = __n;
	  __n = __t;
	}
      return __m;
    }

虽然看起来很恐怖而且不是递归的写法 但肯定是对的

把里面的东西抄抄改改就好了