Программа сумм подмножеств

Я новичок в C ++ и пишу программу суммы подмножеств, которая принимает определенный пользователем набор чисел, первое число которого считается общим. Я пытался использовать DDD для отладки этой программы, однако продолжаю получать ошибку за пределами допустимого. Я не могу понять, почему это происходит. Есть какие-нибудь подсказки? Благодарю.
Вот ошибка:

terminate called after throwing an instance of 'std::out_of_range'
what(): vector::_M_range_check

Код:

#include <iostream>
#include <vector>
#include <cassert>
#include <iomanip>
#include <climits>
#include <math.h>
#include <algorithm>

typedef unsigned int uint;
using namespace std;

//////// Function declerations //////////
void sum_of_subsets( uint index,
uint weight,
uint total,
vector<bool> &include,
vector<uint> &w,
uint W );

bool promising ( uint index,
uint weight,
uint W,
vector<uint> w,
uint total );

.

/////////////// Main //////////////////
int main()
{
//string sortingCode = "-as";
vector<uint> w;             // vector of weights
vector<bool> include;
uint W;                     // the total
uint index = 0;
uint weight = 0;            // the current weight of subsets
uint total = 0;             // the superset total weight
while( ! cin.eof() )
{
uint value;
if( cin >> value && ! cin.eof() )
w.push_back( value );
}

W = w.front();
w.erase( w.begin() );
// instantiate the include vector to false
for( uint k = 0; k <= w.size(); k++ )
include.push_back(0);
// calculate the superset total
for( uint k = 0; k <= w.size()-1; k++ )
total += w.at(k);
// calculate the sum of subsets accordig to CL argument
sum_of_subsets( index, weight, total, include, w, W );

// report success
return 0;
}

.

////////// Function Bodies ///////////
void sum_of_subsets( uint index,
uint weight,
uint total,
vector<bool> &include,
vector<uint> &w,
uint W )
{
cout << "inside sumb_of_subsets" << endl;
if( promising(index, weight, W, w, total) )
{
cout << "promising is true, continue" << endl;
if( weight == W )
{
for( uint k = 0; k <= index; k++ )
{
if(include.at(k))
cout << w.at(k) << " ";;
}
cout << endl;
}
else
{
include.at(index + 1) = 1;
cout << "index1 = " << index << endl;
sum_of_subsets( index + 1,
weight + w.at(index + 1 ),
total - w.at(index + 1),
include, w, W ) ;
include.at(index + 1) = 0;
cout << "index2 = " << index << endl;
sum_of_subsets( index + 1,
weight,
total - w.at(index + 1),
include, w, W );
}
}
}

.

bool promising ( uint index,
uint weight,
uint W,
vector<uint> w,
uint total )
{
cout << "inside promising" << endl;
cout << "W = " << W << endl;
cout << "weight = " << weight << endl;
return ( weight + total >= W )
&& ( (weight == W) || (weight + w.at(index+1) <= W));
}