Dev - C++ - Calcular la función Z de Riemann en C++

#include <iostream>

#include <iomanip>

#include <vector>

#include <cmath>

using namespace std;

void numerosPrimos( vector<int> &v, int sz );

double terminoReimann( int primo, int s );

int main()

{

    vector<int> v;

    int terminos, s;

    double total = 0, temp;

    cout << "\nNumero de terminos: "; cin >> terminos;

    cout << "\nValor de S: "; cin >> s;

    numerosPrimos( v, 5 );

    cout << "\nTerminos:\n";

    cout << fixed;

    for( int i = 0; i < terminos; i++ ) {

        temp = terminoReimann( v[i], s );

        total += temp;

        cout << setprecision(6) << temp;

        if( i < terminos - 1 ) cout << "\n+ ";

        else cout << "\n= ";

    }

    cout << setprecision(6) << total << endl;

    return 0;

}

void numerosPrimos( vector<int> &v, int sz )

{

    int n = 0;

    int i = 1;

    while( n < sz ) {

        int k = 1;

        int a = 0;

        while( k <= i ) {

            if( i % k == 0 && i > 1 ) a++;

            k++;

        }

        if( a == 2 ) {

            v.push_back( i );

            n++;

        }

        i++;

    }

}

double terminoReimann( int primo, int s )

{

    return 1 / ( 1 - pow( primo, -s) );

}