real_t median(const std::initializer_list<real_t> vars) {
    real_t tmp[15];
    const unsigned x = vars.size() / 2;
    if (x & 1) {
      std::partial_sort_copy(vars.begin(),vars.end(),&tmp[0],&tmp[x]);
      return tmp[x];
    }
    const unsigned y = x + 1;
    std::partial_sort_copy(vars.begin(),&tmp[y]);
    return (tmp[x] + tmp[y]) / 2;
  }

我正在使用部分排序来降低复杂性,因为我只需要排序列表的一半.

此外,我假设std :: partial_sort_copy比std :: partial_sort或std :: nth_element更快,因为排序算法(It1！= It2)中没有必要的混洗.我的假设是否正确？

注意：假设real_t可能是双倍的,所以请不要批评使用除法.

NBB：我使用-pedantic和vars已知不超过15个元素.