#include<iostream>
#include<vector>

using std::cout;
using std::vector;

#define MATCH 1
#define MISMATCH -1
#define GAP -2

int max(int i, int j) {
  if (i>=j)
    return i;
  else
    return j;
}

int computeAlign(const char * x,
		 const char * y,
		 vector<vector<int> >& scoremat, 
		 vector<vector<int> >& pointermat) {
  int m=strlen(x)+1;
  int n=strlen(y)+1;
  
  //init scoremat
  for (int i=0; i<m; i++) 
    scoremat[i][0]=GAP*i;
  for (int j=0; j<n; j++)
    scoremat[0][j]=GAP*j;
  
  //init pointermat
  for(int i=1; i<m; i++) 
    pointermat[i][0]=1;
  for(int j=1; j<n; j++) 
    pointermat[0][j]=3;

  //dynamic programming
  for (int i=1; i<m; i++) 
    for (int j=1; j<n; j++) { 
      int a=scoremat[i-1][j]+GAP;
      int b=scoremat[i-1][j-1]+(x[i-1]==y[j-1]? MATCH:MISMATCH);
      int c=scoremat[i][j-1]+GAP;
      scoremat[i][j]=max(a,max(b,c)); //favor counterclockwise order
      if (scoremat[i][j]==a) 
	pointermat[i][j]=1;
      else if (scoremat[i][j]==b)
	pointermat[i][j]=2;
      else pointermat[i][j]=3;
    }
  return scoremat[m-1][n-1];
}

void trackAlign(const vector<vector<int> >& pointermat,
		int i, int j,
		vector<int>& v) {
  if (i!=0 || j!=0) {
    if (pointermat[i][j]==1) {
      trackAlign(pointermat, i-1, j, v);
      v.push_back(1);
    }
    else if (pointermat[i][j]==2) {
      trackAlign(pointermat, i-1, j-1, v);
      v.push_back(2);
    }
    else { 
      trackAlign(pointermat, i, j-1, v);
      v.push_back(3);
    }
  }
}

void printAlign(const char * x, 
		const char * y,
		const vector<int>& v) {
  for (int i=0; i<v.size(); i++) {
    if (v[i]==3)
      cout<<"-";
    else
      cout<<*x++;
  }
  cout<<"\n";
  for (int i=0; i<v.size(); i++) {
    if (v[i]==2)
      cout<<"|";
    else
      cout<<" ";
  }
  cout<<"\n";
  for (int i=0; i<v.size(); i++) {
    if (v[i]==1)
      cout<<"-";
    else
      cout<<*y++;
  }
  cout<<"\n\n";
}

  
int main(int ac, char ** av) {
  
  if (ac<3) {
    cout<<"expecting two strings\n";
    exit(0);
  }
  
  //the strings
  char * x=av[1];
  char * y=av[2];
  int m=strlen(x)+1;
  int n=strlen(y)+1;

  vector<vector<int> > scoremat(m, vector<int>(n, 0));
  vector<vector<int> > pointermat(m, vector<int>(n, 0));
  cout<<"\n\nthe score of the alignment is "<<computeAlign(x, y, scoremat, pointermat)<<"\n\n";
  vector<int> v;
  trackAlign(pointermat, m-1, n-1, v);
  printAlign(x, y, v);
}