编辑距离问题

#include < iostream >
#include <fstream>
#include <string>



using namespace std;
//readfile

int readfile(string filename, string & str1,string &str2)
{
   ifstream ifile(filename.c_str(),ios_base::in);  // ios_base::binary使用该方式读取时会将\r\n分别读取，导致错误
   if (!ifile)
   {
       cerr << "cannot open file!\n";
       exit(1);
   }
   getline(ifile,str1);
   getline(ifile, str2);

   return 0;
}


//将数据写入output文件
template <typename T>
inline int writefile(string filename, T n)
{
   ofstream ofile(filename.c_str(), ios_base::binary | ios_base::out);
   if (!ofile)
   {
       cerr << "cannot open file!\n";
       exit(1);
   }
   ofile << n;
   cout << "write over." << endl;
   return 0;
}

inline int get_min(int x, int y,int z)
{
   int temp=(x > y ? y : x);
   return temp > z ? z : temp;
}

int show_2Dmatrix(int **ma, int row, int col)
{
   for (int i = 0; i < row; ++i){
       for (int j = 0; j < col; ++j){
           std::cout << "\t" << ma[i][j];
       }
       std::cout << std::endl;
   }
   return 0;
}

int editdist(string str1, string str2)
{
   int Len1 =str1.size();
   int Len2 = str2.size();

   //申请栈
   int ** dist = new int*[Len1+1];  
   int ** oper = new int*[Len1 + 1];
   for (int i = 0; i <= Len1; ++i){
       dist[i] = new int[Len2 + 1];
       oper[i] = new int[Len2 + 1];
   }

   //初始化
   for (int col = 0; col <= Len2; ++col){
       dist[0][col] = col;//只需增加
       oper[0][col] = 2;
   }
   for (int row = 0; row <= Len1; ++row){
       dist[row][0] = row;//只需删除
       oper[row][0] = 3;
   }
   oper[0][0] = 0;

   //递归
   for (int i = 1; i <= Len1; ++i){
       for (int j = 1; j <= Len2; ++j){
           
           int cost = str1[i - 1] == str2[j - 1] ? 0 : 1;

           int deletion = dist[i - 1][j] + 1;
           int insertion = dist[i][j-1] + 1;
           int substitution = dist[i - 1][j - 1] + cost;
           //dist[i][j] = get_min(deletion, insertion, substitution);
           if (deletion > insertion){
               if (substitution > insertion){
                   dist[i][j] = insertion;
                   oper[i][j] = 2;
               }
               else{
                   dist[i][j] = substitution;
                   oper[i][j] = cost;
               }
           }
           else{
               if (substitution > deletion){
                   dist[i][j] = deletion;
                   oper[i][j] = 3;
               }
               else{
                   dist[i][j] = substitution;
                   oper[i][j] = cost;
               }
           }
       }
   }
   
   show_2Dmatrix(dist,Len1+1,Len2+1);
   cout << endl;
   show_2Dmatrix(oper, Len1 + 1, Len2 + 1);
   
   int ans = dist[Len1][Len2];
   
   //释放栈
   for (int i = 0; i <= Len1; ++i){
       delete[] dist[i];
   }
   delete[] dist;

   return ans;
}

int main()
{
   string input, output;
   //cin >> input >> output;
   input = "input.txt"; output = "output.txt";
   string str1, str2;
   readfile(input, str1, str2);

   int ans=editdist(str1, str2);

   writefile(output,ans);

   system("pause");
   return 0;
}