#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, " " ) );
}