#include<bits/stdc++.h>
using namespace std;

vector<double> readMFCC(const string &filename) {
// This function reads MFCC feature values from file
vector<double> features;
ifstream infile("mfcc_files/" + filename);
if (!infile) {
cerr << "Error opening file: " << filename << endl;
return features;
}
double val;
while (infile >> val) {
features.push_back(val);
}
infile.close();
return features;//
}

double euclideanDistance(const vector<double>& a, const vector<double>& b) {
double sum = 0;
for (size_t i = 0; i < a.size(); i++) {
sum += pow(a[i] - b[i], 2);
}
return sqrt(sum);
}

int main() {
vector<string> classes = {"AKAL", "AMAR", "ANEKO", "EDURO", "KANOK"};
vector<vector<string>> class_files = {
{"AKAL00.mfc", "AKAL01.mfc", "AKAL02.mfc", "AKAL03.mfc", "AKAL04.mfc"},
{"AMAR00.mfc", "AMAR01.mfc", "AMAR02.mfc", "AMAR03.mfc", "AMAR04.mfc"},
{"ANEK00.mfc", "ANEK01.mfc", "ANEK02.mfc", "ANEK03.mfc", "ANEK04.mfc"},
{"EDUR00.mfc", "EDUR01.mfc", "EDUR02.mfc", "EDUR03.mfc", "EDUR04.mfc"},
{"KANOK00.mfc", "KANOK01.mfc", "KANOK02.mfc", "KANOK03.mfc", "KANOK04.mfc"}
};

int numClasses = classes.size();
vector<vector<double>> distances(numClasses);
vector<double> avgError(numClasses, 0.0);

for (int c = 0; c < numClasses; c++) {
// Read seed (first file)
vector<double> seed = readMFCC(class_files[c][0]);
double sum = 0;

for (size_t f = 0; f < class_files[c].size(); f++) {
vector<double> sample = readMFCC(class_files[c][f]);
double dist = euclideanDistance(seed, sample);
distances[c].push_back(dist);
sum += dist;
}

avgError[c] = sum / class_files[c].size();
}

// Print Output
cout << left << setw(10) << "Classes"
<< setw(45) << "Class Members (Distances)"
<< "Average Error" << endl;
cout << string(80, '-') << endl;

for (int c = 0; c < numClasses; c++) {
cout << left << setw(10) << classes[c];
for (size_t d = 0; d < distances[c].size(); d++) {
cout << fixed << setprecision(2) << distances[c][d];
if (d != distances[c].size() - 1) cout << ", ";
}
cout << setw(10) << " " << avgError[c] << endl;
}

return 0;
}