HDU 5769 접미사 배열 + 2 점

제목:
제목 링크:http://acm.hdu.edu.cn/showproblem.php?pid=5769 문자열 에 문자 x 를 포함 하 는 다른 하위 문자열 이 몇 개 있 는 지 찾 습 니 다.
생각:
height 배열 을 처리 합 니 다. 모든 접미사 i 에 대해 2 분 동안 i 이후 의 첫 번 째 문자 x 의 위 치 를 얻 을 수 있 습 니 다. max (x, sa [i] + height [i]) 이후 의 위치 끝 에 얻 은 문자열 은 모두 유일 합 니 다. 만 남 코드 입 니 다.
코드:

#include 
using namespace std;
typedef long long LL;
const int MAXN = 1e5 + 10;
const int INF = 0x3f3f3f3f;

int t1[MAXN], t2[MAXN], c[MAXN];

bool cmp(int *r, int a, int b, int l) {
    return r[a] == r[b] && r[a + l] == r[b + l];
}

void build(int a[],int sa[],int rk[],int height[],int n,int m) {
    n++;
    int i, j, p, *x = t1, *y = t2;
    //       ，  s      ，       
    for(i = 0; i < m; i++) c[i] = 0;
    for(i = 0; i < n; i++) c[x[i] = a[i]]++;
    for(i = 1; i < m; i++) c[i] += c[i-1];
    for(i = n-1; i >= 0; i--) sa[--c[x[i]]] = i;
    for(j = 1; j <= n; j <<= 1) {
        p = 0;
        //    sa         
        for(i = n-j; i < n; i++)y[p++] = i;//   j            
        for(i = 0; i < n; i++)if(sa[i] >= j)y[p++] = sa[i] - j;
        //    y                 
        //         
        for(i = 0; i < m; i++) c[i] = 0;
        for(i = 0; i < n; i++) c[x[y[i]]]++;
        for(i = 1; i < m; i++) c[i] += c[i-1];
        for(i = n-1; i >= 0; i--) sa[--c[x[y[i]]]] = y[i];
        //  sa x      x  
        swap(x,y);
        p = 1;
        x[sa[0]] = 0;
        for(i = 1; i < n; i++)
            x[sa[i]] = cmp(y,sa[i-1],sa[i],j)?p-1:p++;
        if(p >= n)break;
        m = p;//          
    }
    int k = 0;
    n--;
    for(i = 0; i <= n; i++)rk[sa[i]] = i;
    for(i = 0; i < n; i++)
    {
        if(k)   k--;
        j = sa[rk[i]-1];
        while(a[i+k] == a[j+k])
            k++;
        height[rk[i]] = k;
    }
}

int sa[MAXN], height[MAXN], rk[MAXN], a[MAXN];
char s[MAXN], t[2];
vector <int> vec;
vector <int> :: iterator it;

int main(){
    //freopen("in.txt", "r", stdin);
    int T, cs = 0;
    scanf("%d", &T);
    while (T--) {
        scanf("%s%s", t, s);
        int n = strlen(s);
        vec.clear();
        for (int i = 0; i < n; i++) {
            a[i] = s[i];
            if (s[i] == t[0]) vec.push_back(i);
        }
        a[n] = 0;
        build(a, sa, rk, height, n, 128);
        LL ans = 0;
        for (int i = 1; i <= n; i++) {
            it = lower_bound(vec.begin(), vec.end(), sa[i]);
            if (it == vec.end()) continue;
            int tmp = n - sa[i] - max(height[i], *it - sa[i]);
            if (tmp > 0) ans += tmp;
        }
        printf("Case #%d: %I64d
", ++cs, ans);
    }
    return 0;
}