标签：style blog http color 使用 os io for
一道字符串的题。
求两个字符串的长度大于K的子串中，相同的有多少对。
知道了要用到后缀数组，可是想了好久也没想好怎么使用。最后看了一下论文，提到了两个字符串可以合并，中间加上一个特殊字符就可以了。|
虽然这样还是写了好长时间啊。。。
首先是合并，然后建立后缀数组，求出sa[]，height[]，然后在进行加和，用到了一个题解称作单调栈的东西。不知是什么，反正就是合并嘛。
分别对于两个串进行加和。
#include <iostream> #include <cstdio> #include <cstring> using namespace std; #define maxn 200010 #define LL long long int rank[maxn], sa[maxn], t1[maxn], t2[maxn], sum[maxn], h[maxn]; char c[maxn]; void solve() { int *x = t1; int *y = t2; int m = 128; int l = strlen(c); for(int i = 0; i < m; i++) sum[i] = 0; for(int i = 0; i <= l; i++) sum[x[i]=c[i]]++; for(int i = 1; i < m; i++) sum[i] += sum[i-1]; for(int i = l; i >= 0; i--) sa[--sum[x[i]]] = i; for(int j = 1; j <= l; j <<= 1) { int num = 0; for(int i = l-j+1; i<=l; i++) y[num++] = i; for(int i = 0; i <= l; i++) if(sa[i] >= j) y[num++] = sa[i] - j; for(int i = 0; i < m; i++) sum[i] = 0; for(int i = 0; i <= l; i++) sum[x[i]]++; for(int i = 1; i < m; i++) sum[i] += sum[i-1]; for(int i = l; i >= 0; i--) sa[--sum[x[y[i]]]] = y[i]; swap(x,y); m = 1; x[sa[0]] = 0; for(int i = 1; i <= l; i++) x[sa[i]] = (y[sa[i]] == y[sa[i-1]] && y[sa[i]+j] == y[sa[i-1]+j]) ? m-1:m++; if(m > l) break; } m = 0; for(int i = 0; i <= l; i++) rank[sa[i]] = i; for(int i = 0; i < l; i++) { if(m) m--; int j = sa[rank[i]-1]; while(c[i+m] == c[j+m]) m++; h[rank[i]] = m; } } int main() { //freopen("in.txt", "r", stdin); //freopen("out(2).txt", "w", stdout); int n; while(scanf("%d", &n) && n) { scanf("%s", c); int l1 = strlen(c); c[l1] = 1; scanf("%s", c+l1+1); int l2 = strlen(c); c[l2] = 0; solve(); LL ans = 0; LL sum2 = 0; int l = 0; for(int i = 1; i <= l2; i++) { if(h[i] < n) { l = 0; sum2 = 0; } else { int tmp = 0; while(l && t1[l-1] >= h[i]) { l--; tmp += t2[l]; sum2 -= t2[l] * (t1[l] - h[i]); } t1[l] = h[i]; t2[l++] = tmp; if(sa[i-1] > l1) { t2[l-1]++; sum2 += h[i]- n + 1; } if(sa[i] < l1) ans += sum2; } } l = 0; sum2 = 0; for(int i = 1; i <= l2; i++) { if(h[i] < n) { l = 0; sum2 = 0; } else { int tmp = 0; while(l && t1[l-1] >= h[i]) { l--; tmp += t2[l]; sum2 -= t2[l] * (t1[l] - h[i]); } t1[l] = h[i]; t2[l++] = tmp; if(sa[i-1] < l1) { t2[l-1]++; sum2 += h[i] - n + 1; } if(sa[i] > l1) ans += sum2; } } cout << ans << endl; } }
这样的话对于求最长回文字串我也会了。。。
不过是将第二个串换成是它的逆序，然后找到最长的
标签：style blog http color 使用 os io for
原文地址：http://www.cnblogs.com/ye8370/p/3910851.html
