// Instructions:
//    compile with g++ pretty2dick.cc -o pretty2dick
// Running:
//   pretty2dick inputfile1 [inputfile2 inputfile3 ... inputfileN] outputfilename 

#include <iostream>
#include <string>
#include <vector>
#include <cstdlib>
#include <fstream>
using namespace std;

void usage();

int main(int ac, char **av)
{
	if(ac < 3) {usage(); return 1;}
	
	int samsize, numpolys;
	vector<int> positions;
	vector<string> indivs;
	vector< vector<string> > polys;
	ifstream istr;

	for(int i=1; i<ac-1; i++)
	{
		istr.open(av[i], ios::in);
		string line;
		getline(istr, line);
		while(line.length()!=0)
		{
			if(positions.size() == 0 || atoi(line.c_str()) != positions[positions.size()-1])
			{
				positions.push_back(atoi(line.c_str()));
				if(positions.size() == 1)
				{
					line.erase(0, line.find("\t")+1);
					indivs.push_back(line);
					indivs[indivs.size()-1].erase(line.find("\t"), line.length());
					line.erase(0, line.find("\t")+1);
				}
				else
				{
					line.erase(0, line.find("\t")+1);
					line.erase(0, line.find("\t")+1);
				}

				vector<string> placeholder;
				placeholder.push_back(line);
				placeholder[0].erase(line.find("\t"), line.length());
				line.erase(0, line.find("\t")+1);
				placeholder.push_back(line);
				polys.push_back(placeholder);

/*				cout << positions[positions.size()-1] << "\t"; //debugging
				cout <<indivs.size()<< indivs[indivs.size()-1] << "\t" << polys[polys.size()-1][0]
					<< "\t" << polys[polys.size()-1][1] << endl; //debugging
*/			}
			else
			{
				if(positions.size() == 1)
				{
					line.erase(0, line.find("\t")+1);
					indivs.push_back(line);
					indivs[indivs.size()-1].erase(line.find("\t"), line.length());
					line.erase(0, line.find("\t")+1);
				}
				else
				{
					line.erase(0, line.find("\t")+1);
					line.erase(0, line.find("\t")+1);
				}
				
				polys[polys.size()-1].push_back(line);
				polys[polys.size()-1][polys[polys.size()-1].size()-1].erase(line.find("\t"), line.length());
				line.erase(0, line.find("\t")+1);
				polys[polys.size()-1].push_back(line);

/*				cout << positions[positions.size()-1] << "\t"; //debugging
				cout <<indivs.size()<< indivs[indivs.size()-1] << "\t" << polys[polys.size()-1][polys[polys.size()-1].size()-2]
					<< "\t" << polys[polys.size()-1][polys[polys.size()-1].size()-1] << endl; //debugging
*/			}

			getline(istr, line);
		}
		istr.close();
	}

	ofstream outfile(av[ac-1], ios::out | ios::trunc);
	outfile << "//"<<endl;
	outfile << "\t" << indivs.size() << "\t" << polys.size() << endl;
	outfile << "\t";
	for(int a=0; a < positions.size(); a++)
	{
		outfile << positions[a];
		if(a != positions.size() - 1) outfile << "\t";
		else outfile << endl;
	}

	for(int j=0; j<indivs.size(); j++)
	{
		outfile << indivs[j] << "B\t";
		for(int k=0; k < positions.size(); k++)
		{
			outfile << polys[k][2*j];
			if(k != positions.size() - 1) outfile << "\t";
			else outfile << endl;
		}
		outfile << indivs[j] << "D\t";
		for(int l=0; l < positions.size(); l++)
		{
			outfile << polys[l][2*j+1];
			if(l != positions.size() - 1) outfile << "\t";
			else outfile << endl;
		}
	}

	return 0;
}

void usage()
{
	cout << "pretty2dick infile1 [infile2] ... targetfile" << endl;
	cout << "An unlimited number of infiles may be passed" << endl;
	cout << "All will be concatinated into the target file" << endl;
}