Eu tenho tentado encontrar a interseção entre dois std :: set em C ++, mas continuo recebendo um erro.

Eu criei um pequeno teste de amostra para este

#include <iostream>
#include <vector>
#include <algorithm>
#include <set>
using namespace std;

int main() {
  set<int> s1;
  set<int> s2;

  s1.insert(1);
  s1.insert(2);
  s1.insert(3);
  s1.insert(4);

  s2.insert(1);
  s2.insert(6);
  s2.insert(3);
  s2.insert(0);

  set_intersection(s1.begin(),s1.end(),s2.begin(),s2.end());
  return 0;
}

O último programa não gera nenhuma saída, mas espero ter um novo conjunto (vamos chamá-los3) com os seguintes valores:

s3 = [ 1 , 3 ]

Em vez disso, estou recebendo o erro:

test.cpp: In function ‘int main()’:
test.cpp:19: error: no matching function for call to ‘set_intersection(std::_Rb_tree_const_iterator<int>, std::_Rb_tree_const_iterator<int>, std::_Rb_tree_const_iterator<int>, std::_Rb_tree_const_iterator<int>)’

O que eu entendo deste erro, é que não há definição emset_intersection que aceitaRb_tree_const_iterator<int> como parâmetro.

Além disso, suponho que ostd::set.begin() método retorna um objeto desse tipo,

existe uma maneira melhor de encontrar a interseção de doisstd::set em C ++? De preferência, uma função embutida?

Muito obrigado!