Algorithm C++

#include 
#include 
#include 
#include 
#include 
#include 
using namespace std;
int main( )
{
   const int last_number = 50;
   vector v( last_number - 1, 1 );
   v[0] = 2;
   partial_sum( v.begin(), v.end(), v.begin() );
   vector::iterator stop = v.end();
   for( vector::iterator start = v.begin(); start != stop;++start )
   // remove all subsequent numbers that are not divisible by the current one
   stop = remove_if( start+1, stop,not1( bind2nd( modulus(), *start ) ) );
   cout << "Prime numbers: ";
   copy( v.begin(), stop, ostream_iterator( cout, " " ) );
}