POJ 2250解题报告

求两个字符串的公共子序列。因为是最基本的算法，所以一次就过了。并且用的是C++的vector和string。。。table数组记录最长公共子序列的长度，如果当前比较的两个字符（本题中为字符串）相当，则table[i][j] = table[i - 1][j - 1] + 1，否则，table[i][j] = max(table[i - 1][j], table[i][j - 1])。table数组的空间肯定可以改进，但想来就挺麻烦的。。。path数组记录路径，0表示当前点包含在最长公共子序列中，1表示向左找，2表示向上找。通过从头往前查path数组就可以构造出最长公共子序列。

注：代码中注释的部分是把一个string按照空格分隔成单词，保存在vector<string>里面，相当于Python或Java里面的split，只是不能指定分隔符。

#include <iostream>
#include <string>
#include <vector>
#include <sstream>
#include <iterator>

using namespace std;

string longestcommonsubsequence(const vector<string> &stra, const vector<string> &strb)
{
int **table = new int*[stra.size()];
for(int i = 0; i < stra.size(); ++i)
{
table[i] = new int[strb.size()];
}

int **path = new int*[stra.size()];
for(int i = 0; i < stra.size(); ++i)
{
path[i] = new int[strb.size()];
}

for(int i = 0; i < stra.size(); ++i)
{
for(int j = 0; j < strb.size(); ++j)
{
if(stra[i] == strb[j])
{
if(i == 0 || j == 0)
{
table[i][j] = 1;

}
else
{
table[i][j] = table[i - 1][j - 1] + 1;
}
path[i][j] = 0;
}
else
{
if(i == 0 && j == 0)
{
table[i][j] = 0;
path[i][j] = 1;
}
else if(i != 0 && j == 0)
{
table[i][j] = table[i - 1][j];
path[i][j] = 1;
}
else if(i == 0 && j != 0)
{
table[i][j] = table[i][j - 1];
path[i][j] = 2;
}
else
{
int len1 = table[i - 1][j];
int len2 = table[i][j - 1];
if(len1 > len2)
{
table[i][j] = len1;
path[i][j] = 1;
}
else
{
table[i][j] = len2;
path[i][j] = 2;
}
}
}
}
}
string lcs;
bool first = true;
for(int i = stra.size() - 1, j = strb.size() - 1; i >= 0 && j >= 0;)
{
if(path[i][j] == 0)
{
if(first)
{
lcs.insert(0, stra[i]);
first = false;
}
else
lcs.insert(0, stra[i] + " ");
i--;
j--;
}
else if(path[i][j] == 1)
{
i--;
}
else
{
j--;
}
}
return lcs;
}

int main()
{
/*
string sentence = "die einkommen der landwirte sind fuer die abgeordneten ein buch mit sieben siegeln um dem abzuhelfen muessen dringend alle subventionsgesetze verbessert werden";
istringstream iss(sentence);
vector<string> stra;
copy(istream_iterator<string>(iss),
istream_iterator<string>(),
back_inserter<vector<string> >(stra));
sentence = "die steuern auf vermoegen und einkommen sollten nach meinung der abgeordneten nachdruecklich erhoben werden dazu muessen die kontrollbefugnisse der finanzbehoerden dringend verbessert werden";
//iss.str(sentence);
istringstream iss2(sentence);
vector<string> strb;
copy(istream_iterator<string>(iss2),
istream_iterator<string>(),
back_inserter<vector<string> >(strb));
*/
vector<string> stra, strb;
string str;
while(true)
{
stra.clear();
strb.clear();
while(true)
{
cin>>str;
if(cin.eof())
return 0;
if(str == "#")
break;
stra.push_back(str);
}
while(true)
{
cin>>str;
if(str == "#")
break;
strb.push_back(str);
}
cout<<longestcommonsubsequence(stra, strb)<<endl;
}
return 0;
}