CF427D

CF427D

SAの賢いが、役に立たない、実際にはボード。

質問の意味:

与えられた2つの文字列が、最短では国民は、連続文字列に1回だけ発生会います

解像度:

公共と最短のサブストリングは、この質問は何の関係もありません尋ねるかのように通常の状況下では、SAは、最長の共通のプレフィックスを見つけるために使用されています。
しかし、我々はまだ類推のアイデアで描くことができます:

あなたが最大の要素ZZを見つけたい理由を考えてみて?
最大未満であればあるため、最大値はこのシーケンスに含まれます。
だから、答えは一つの要素は、どの要素のz値がこれよりも大きくないということです、もちろん、それは、zの最大値を選択することで
、上記の考え方からすると小さい取得する方法への回答を見つける
+自然の二番目に大きいz値が十分ではありませんが、しかし、二番目に大きい値を見つけました1の条件を満足
一方は、それが最大値よりも小さい場合、最大値が一意に含むされ、一方、第二の最大値よりも大きい場合、最大値のみが含まれ
、+1だけ二番目に大きい値を許容することが可能ですA。
したがって、上記の方法、およびすべての$ S_1 $サフィックスのための二番目に大きい求めて最大値、二番目に大きい値に応じて、その後の答えを更新し
た後、我々はこの問題を解決するSAを楽しんでいることができます。

コード:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>

using namespace std;

#define LL long long
#define N 10010

string s1,s2,str;
int SA[N],rk[N],tp[N],cnt[N];
int len,tot,m,a[N],height[N];

void qsort() {
    for(int i = 1 ; i <= m ; i++) cnt[i] = 0;
    for(int i = 1 ; i <= tot; i++) cnt[rk[i]]++;
    for(int i = 1 ; i <= m ; i++) cnt[i] += cnt[i - 1];
    for(int i = tot ; i >= 1 ; i--) SA[cnt[rk[tp[i]]]--] = tp[i]; 
}
inline bool cmp(int *f,int x,int y,int w) {
    return f[x] == f[y] && f[x + w] == f[y + w];
}
void build_SA() {
    m = 127;
    for(int i = 1 ; i <= tot ; i++) {
        rk[i] = a[i];
        tp[i] = i;
    }
    qsort();
    for(int w = 1 , p = 0 ; p < tot ; w += w,m = p) {
        p = 0;
        for(int i = tot - w + 1 ; i <= tot ; i++) tp[++p] = i;
        for(int i = 1 ; i <= tot ; i++) {
            if(SA[i] > w) tp[++p] = SA[i] - w;
        } 
        qsort();
        swap(rk,tp);
        rk[SA[1]] = p = 1;
        for(int i = 2 ; i <= tot ; i++) 
            rk[SA[i]] = cmp(tp,SA[i],SA[i - 1],w) ? p : ++p;
    }
    int j = 0, k = 0; 
    for(int i = 1 ; i <= tot ; height[rk[i++]] = k) {
        for(k = k ? k - 1 : k, j = SA[rk[i] - 1] ; a[i + k] == a[j + k] ; k++); 
    }
}
inline bool check(int k,int div) {
    int cnt1 = 0,cnt2 = 0;
    for(int i = 1 ; i <= tot ; i++) {
        if(height[i] < k) {
            if(cnt1 == 1 && cnt2 == 1) return true;
            cnt1 = cnt2 = 0;
            if(SA[i] <= div) cnt1++;
            else if(SA[i] >= div) cnt2++;
            continue;
        }
        if(SA[i] <= div) cnt1++;
        else if(SA[i] >= div) cnt2++;
    }
    return cnt1 == 1 && cnt2 == 1;
}

int main() {
    cin>>s1>>s2;
    len = s1.length();
    str = s1 + '#' + s2;//加入'#'表示两个字符串的分界点。
    tot = len + s2.length() + 1;
    for(int i = 1 ; i <= tot ; i++) a[i] = str[i - 1];
    build_SA();
    int ans = -1;
    for(int i = 1 ; i <= len ; i++) {
        if(check(i,len)) {
            ans = i;
            break;
        }
    }
    printf("%d \n",ans);
    //system("pause");
    return 0;
}

おすすめ

転載: www.cnblogs.com/Repulser/p/11373550.html