#include <iostream>
#include <vector>
#include <numeric>
#include <cmath>

using VecD = std::vector<double>;

double getMean(const VecD& data) {
return std::accumulate(data.cbegin(), data.cend(), 0.0) / data.size();
}

double getVariance(const VecD& data, double meanVal) {
double sum = 0.0;
for (double x : data) {
sum += std::pow(meanVal - x, 2);
}
return sum / data.size();
}

double getStdDev(const VecD& data, double meanVal) {
return sqrt(getVariance(data, meanVal));
}

VecD zScoreNormalize(const VecD& data) {
double meanVal = getMean(data);
double stdDev = getStdDev(data, meanVal);
VecD normalizedData(data.size());
for (size_t i = 0; i < data.size(); i++) {
normalizedData[i] = (data[i] - meanVal) / stdDev;
}
return normalizedData;
}

int main() {
VecD v = {50, 60, 70, 80, 90, 600};
auto ans = zScoreNormalize(v);
for (double x : ans) {
std::cout << x << ' ';
}
std::cout << std::endl;
return 0;
}