PS

BOJ 9249 : 최장 공통 부분 문자열

lickelon 2024. 11. 1. 23:28

코드

#include <bits/stdc++.h>

#define all(x) (x).begin(), (x).end()

#define INF 0x7FFFFFFF

using namespace std;

using ll = long long;
using ld = long double;
using pii = pair<int,int>;
using pll = pair<ll, ll>;

vector<int> getsa(string &s) {
    int n = s.size();
    vector<int> sa(n), g(n+1), ng(n+1);
    for(int i = 0; i < n; i++) {
        sa[i] = i;
        g[i] = s[i];
    }

    g[n] = -1;
    for(int t = 1; t < n; t <<= 1) {
        auto cmp = [&](int i, int j) {
            if(g[i] == g[j]) return g[i+t] < g[j+t];
            else return g[i] < g[j];
        };
        sort(all(sa), cmp);
        ng[sa[0]] = 0;
        ng[n] = -1;
        for(int i = 1; i < n; i++) {
            if(cmp(sa[i-1], sa[i])) ng[sa[i]] = ng[sa[i-1]] + 1;
            else ng[sa[i]] = ng[sa[i-1]];
        }
        g = ng;
    }
    return sa;
}

vector<int> getlcp(vector<int> &sa, string &s){
    int n = s.size();
    vector<int> lcp(n), isa(n);
    for(int i=0;i<n;++i) isa[sa[i]] = i;
    for(int k=0, i=0;i<n;++i) {
        if(isa[i]){
            for(int j=sa[isa[i]-1]; s[i+k]==s[j+k]; ++k);
            lcp[isa[i]] = (k ? k-- : 0);
        }
    }
    return lcp;
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0); cout.tie(0);

    string a, b;
    cin >> a >> b;
    string s = a+" "+b;
    auto sa = getsa(s);
    auto lcp = getlcp(sa, s);

    int l = 0, idx = -1;
    for(int i = 1; i < sa.size(); i++) {
        ll ta = sa[i-1] - a.length();
        ll tb = sa[i] - a.length();
        if(ta*tb < 0) {
            if(lcp[i] > l) {
                l = lcp[i];
                idx = sa[i];
            }
        }
    }

    cout << l << "\n" << s.substr(idx, l);

    return 0;
}

풀이

두 문자열을 구분자와 함께 연결해 한 문자열로 만든다.

접미사 배열을 구하고 접미사 배열에서 인접한 두 위치가 구분자를 기준으로 서로 다른 위치인지 확인한다.

다른 위치라면 두 위치의 lcp의 최댓값을 저장한다.

728x90