【机试练习】【C++】【codeup 2432求最长公共子串（串）】

/**
 *2432: 求最长公共子串（串）
 *题目描述
 *求采用顺序结构存储的串s和串t的一个最长公共子串，若没有则输出false，若最长的有多个则输出最先出现的那一串。

 *输入
 *输入两个字符串

 *输出
 *输出公共子串

 *样例输入
 *abcdef
 *adbcef
 *样例输出
 *bc
 */
#include
#include
#include
#include
#include
using namespace std;
// 定义取模的数值
const int _P = 10000019;
const int _MOD = 1000000007;

const int MAXN = 10100;
string a, b;

long long Ha[MAXN], Hb[MAXN];

long long _POW[MAXN];
void _init(int len){
	_POW[0] = 1;
	for(int i = 1; i <= len; i++){ // 这里使用len，节省时间 
		_POW[i] = (_POW[i - 1] * _P) % _MOD;
	}
}


// 首先实现求字符串各个位置二十六进制取模哈希的函数
void calc_hash(long long H[], const string& a){
	H[0] = a[0];
	for(int i = 1; i < a.length(); i++){
		H[i] = (H[i - 1] * _P + a[i]) % _MOD;
	}
}

// 实现求子串[i ... j]哈希的函数
long long calc_sub_hash(long long H[], int i, int j){
	if(i == 0){
		return H[j];
	} 
	return ((H[j] - H[i - 1] * _POW[j - i + 1]) % _MOD + _MOD) % _MOD; // 到这里才知道要先实现POW 
}

struct Store{
	int i;
	int j;
	long long hashvalue;
};

// 调用calc_sub_hash求子串哈希值，子串长度的集合
void calc_all_sub_hash(long long H[], int len, vector<Store> &pr){ // len是字符串总长度 
	long long hashvalue = -1;
	for(int i = 0; i < len; i++){
		for(int j = i; j < len; j++){ // 不可能是反向的 
			hashvalue = calc_sub_hash(H, i, j);
			Store s;
			s.i = i;
			s.j = j;
			s.hashvalue = hashvalue;	
			pr.push_back(s);
		}
	}	
} 

string find_max_substr(vector<Store> pr_a, vector<Store> pr_b, const string &str_a){
	int maxi = -1, maxj = -1;
	int maxlen = -1;
	for(int i = 0; i < pr_a.size(); i++){
		for(int j = 0; j < pr_b.size(); j++){
			if(pr_a[i].hashvalue == pr_b[j].hashvalue){ // 这儿注意！ 不能都是i 否则答案错误 
				// 看长度一样不一样
				int lena = pr_a[i].j - pr_a[i].i + 1;
				int lenb = pr_b[j].j - pr_b[j].i + 1;
				if(lena == lenb){
					if(maxlen < lena){
						maxlen = lena;
						maxi = pr_a[i].i;
						maxj = pr_a[i].j;
					} else if(maxlen == lena){
						maxi = min(maxi, pr_a[i].i);
						maxj = min(maxj, pr_a[i].j);
					}
				}
			}
		}
	}
	if(maxi < 0) return "false";
	return str_a.substr(maxi, maxj);	
}




int main(){
	string str_a, str_b;
	getline(cin ,str_a);
	getline(cin ,str_b);
	_init(max(str_a.length(), str_b.length()));
	calc_hash(Ha, str_a);
	calc_hash(Hb, str_b);
	vector<Store> pr_a, pr_b;
	calc_all_sub_hash(Ha, str_a.length(),pr_a);
	calc_all_sub_hash(Hb, str_b.length(),pr_b);
//	for(int i = 0; i < pr_a.size(); i++){
//		cout << pr_a[i].hashvalue << " " << pr_a[i].i << " " << pr_a[i].j << endl;
//	}
//	cout << endl;
//	for(int i = 0; i < pr_b.size(); i++){
//		cout << pr_b[i].hashvalue << " " << pr_b[i].i << " " << pr_b[i].j << endl;
//	}
	cout<< find_max_substr(pr_a, pr_b, str_a);	
	return 0;
}