0
  typedef pair<double, double> dd; 

  const double epsilon = 1e-6;

  struct sort_by_polar_angle { 
    dd center; 
    // Constuctor of any type 
    // Just find and store the center 
    template<typename T> sort_by_polar_angle(T b, T e) { 
        int count = 0;
        center = dd(0,0); 
        while(b != e) { 
                    center.first += b->first;
                    center.second += b->second;
               b++; 
            count++;
        } 
               double k = count ? (1.0/count) : 0;
        center.first *= k;
               center.second *= k;
   } 
   // Compare two points, return true if the first one is earlier 
   // than the second one looking by polar angle 
   // Remember, that when writing comparator, you should 
   // override not ‘operator <’ but ‘operator ()’ 
   bool operator () (const dd& a, const dd& b) const { 
        double p1 = atan2(a.second-center.second, a.first-center.first); 
        double p2 = atan2(b.second-center.second, b.first-center.first); 
        return p1 + epsilon < p2; 
   } 
  }; 

// ... 

vector < dd >  points; 

sort(all(points), sort_by_polar_angle(all(points))); 

当 sort_by_polar_angle() 被调用时,它是否起到了构造函数的作用?重载的operator()如何正确使用?

4

3 回答 3

4

当你调用函数sort_by_polar_angle()sort(),你正在创建一个临时类型的对象sort_by_polar_angle(即调用它的构造函数)。在排序算法中,您传递的函子对象使用了类似的东西functor(),它将被翻译成functor.operator().

于 2009-04-23T07:07:59.930 回答
0

它不适用于从中心直线上的点。这些的“角度”是相等的,所以它们的顺序是不确定的。对这些点进行排序将返回未确定的结果。

这是因为任何“排序”操作都应该是“严格”的:如果 a < b,则 b > a。(实际上定义要复杂一些。)

dd a0 = { 0, 1 };
dd b0 = { 0, 2 };

assert( sort_by_polar_angle ()( a0, b0 ) && ! sort_by_polar_angle () ( b0, a0 ) );
于 2009-04-23T09:16:49.837 回答
0

尝试在这里了解函子是什么的基本介绍,或在谷歌上搜索 c++ 函子

于 2009-04-24T02:13:45.047 回答