using namespace std;
#include <iostream>
#include <iomanip>
#include <fstream>
#include <sstream>
#include <math.h>
#include <string>


int main(int argc, char **argv) {

	double avg, dev;
	double tmp;
	int i, j, col, zaehler;
	int ecut;
	int flag;
	string zeile;
	stringstream line;
	
	if (argc < 4) {
		cout << "Usage: " << argv[0] << " <Spalte> <X> <File1> <File2> ...\n";
		return(1);
	} 
  line.str(argv[1]);
  line >> col;
  line.clear();
  line.str(argv[2]);
  line >> ecut;
  cout << "Using " << col << " with Ecut of " << ecut << "." << endl;

	// get average
	avg = 0.;
	zaehler=0;
	for (i=3;i<argc;i++) {
		//cout << "Reading column " << col << " from file " << argv[i] << "\n";
		ifstream test(argv[i]);
		if (test == NULL) { cout << "Can't open File " << argv[i] << "\n"; return(255); }
		flag=1;
		//cout << "Looking for " << ecut << " \n";
		while (getline(test, zeile, '\n')) {
			//cout << zeile;
			istringstream input(zeile);
			input >> ws >> tmp;			
			if (tmp==(double)ecut) {	// found correct line!
				for (j=col;j>0;j--)
					if (!input.eof())
						input >> ws >> tmp;
					else
						flag=0;
				//cout << tmp << "\n";
				avg += tmp;
			} 
			//else cout << " skipping\n";
		}
		zaehler+=flag;
		test.close();
	}
	if (zaehler != 0) avg /= zaehler;	

	// get deviation
	dev = 0.;
	zaehler=0;
	for (i=3;i<argc;i++) {
		//cout << "Reading column " << col << " from file " << argv[i] << "\n";
		ifstream test(argv[i]);
		if (test == NULL) { cout << "Can't open File " << argv[i] << "\n"; return(255); }
		flag=1;
		while (getline(test, zeile, '\n')) {
			istringstream input(zeile);
			input >> ws >> tmp;			
			if (tmp==(double)ecut) {	// found correct line!
				for (j=col;j>0;j--)
					if (!input.eof())
						input >> ws >> tmp;
					else
						flag=0;
				//cout << tmp << "\n";
				dev += pow(tmp - avg, 2);
			}
		}
		zaehler+=flag;
		test.close();
	}
	//dev = 1/(zaehler)*dev;
	dev /= zaehler;
			
	cout << setprecision(8) << avg << "\t" << setprecision(8) << sqrt(dev) << "\n";
	return(0);
}