pext.cpp

template <typename T>
constexpr T pext(T x, T mask) {
  T ret = 0;
  for (T bb = 1; mask != 0; bb += bb) {
    if (x & mask & -mask) {
      ret |= bb;
    }
    mask &= (mask - 1);
  }
  return ret;
}