2016 ChinaFinal F. Mr. Panda and Fantastic Beasts SAM

简略题意：问第一个串的子串中，不在其他n-1个串出现的长度最小且字典序最小的子串是什么。

将n-1个子串用不会出现的字符拼接起来，那么问题转化成
问A串的子串中，不在B串出现的长度最小且字典序最小的子串是什么。

将B串建立后缀自动机，用A串在自动机上跑。和湖南大学校赛的那题一样，动态维护A的子串匹配上B的子串的长度，当失配的时候更新答案即可。
需要注意的是，从一个节点走向下一个节点，x++的答案虽然是匹配上B的，但是对于这题，x需要比B中匹配的上的节点right集合中最大的那个串还要大才行..所以不能是简单的++。

#define others
#ifdef poj
#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <algorithm>
#include <vector>
#include <string>
#include <map>
#include <set>
#endif // poj
#ifdef others
#include <bits/stdc++.h>
#endif // others
//#define file
#define all(x) x.begin(), x.end()
using namespace std;
#define eps 1e-8
const double pi = acos(-1.0);

typedef long long LL;
typedef unsigned long long ULL;
void umax(int &a, int b) {
    a = max(a, b);
}
void umin(int &a, int b) {
    a = min(a, b);
}
int dcmp(double x) {
    return fabs(x) <= eps?0:(x > 0?1:-1);
}
void file() {
    freopen("data_in.txt", "r", stdin);
    freopen("data_out.txt", "w", stdout);
}

namespace solver {
    int res = 1e9, pos = 0;
    string s, t;
    namespace SAM{
        static const int maxnode = 2e6+10;//至少开两倍
        static const int maxn = 27;
        int p, q, np, nq;
        int cnt, ST_T, len;
        //int deg[maxnode], right[maxnode];
        vector<int>G;
        int trans[maxnode][maxn], l[maxnode], fa[maxnode];
        int newnode() {
            int x = ++cnt;
            memset(trans[x], 0, sizeof trans[x]);
            fa[x] = 0;
            l[x] = 0;
            return cnt;
        }
        void init() {
            cnt = 0;
            ST_T = newnode();
            /*
            memset(right, 0, sizeof right);
            memset(deg, 0, sizeof deg);
            G.clear();
            */
        }
        void add(int c) {
            //令当前串为T，新加的字符为x。
            p = ST_T; np = ST_T = newnode(); l[np] = l[p] + 1;//令p = ST(T)，新建np = ST(Tx)
            while(!trans[p][c]&&p) trans[p][c] = np, p = fa[p];//对于p的所有没有标号c的边的祖先v，trans[v][c] = np。
            if(!p) fa[np] = 1; //找到p的第一个祖先vp,他有标号c的边，如果没有这样的vp，那么fa[p]=root，结束该阶段。
            else {
                q = trans[p][c];//令q=trans[vp][c]
                if(l[p] + 1 == l[q]) fa[np] = q;//若l[q] = l[vp] + 1,令fa[np] = q,结束该阶段。
                else {
                    nq = newnode(); l[nq] = l[p] + 1;//否则建立新节点nq
                    memcpy(trans[nq],trans[q],sizeof(trans[q]));//trans(nq, *) = trans(q, *)
                    fa[nq] = fa[q];
                    fa[np] = fa[q] = nq;
                    while(trans[p][c] == q) trans[p][c] = nq, p = fa[p];//对于所有的trans(v, c) == q的p的祖先v, trans(v, c)改为nq。
                }
            }
        }

        void build(string &str) {
            len = str.size();
            for(int i = 0; i < len; i++)
                add(str[i]-'a');

            /*
            for(int i = 1; i <= cnt; i++)
                for(int j = 0; j < 26; j++)
                    if(trans[i][j]) deg[trans[i][j]]++;
            queue<int> q;
            q.push(1);

            while(!q.empty()) {
                int u = q.front(); q.pop();
                G.push_back(u);
                for(int i = 0; i < 26; i++)
                    if(trans[u][i]) {
                        --deg[trans[u][i]];
                        if(deg[trans[u][i]] == 0)
                            q.push(trans[u][i]);
                    }
            }

            int p = 1;
            for(int i = 0; i < len; i++) {
                 p = trans[p][str[i]-'a'];
                 right[p]++;
            }
            for(int i = G.size()-1; i >= 0; i--) {
                int val = G[i];
                right[fa[val]] += right[val];
            }
            求解right集合用
            */
        }
        void solve(string &str) {
            int x = 0, p = 1;
            for(int i = 0; i < str.size(); i++) {
                if(trans[p][str[i]-'a']) {
                    p = trans[p][str[i]-'a'];
                    x = l[fa[p]] + 1;
                } else {
                    while(p && !trans[p][str[i]-'a']) {
                        if(res > x + 1) {
                            res = x + 1;
                            pos = i;
                        } else if(res == x + 1) {
                            if(s.substr(i - res + 1, res) < s.substr(pos - res + 1, res)) {
                                pos = i;
                            }
                        }
                        p = fa[p];
                        x = l[fa[p]] + 1;
                    }
                    if(p == 0) {
                        x = 0, p = 1;
                    } else {
                        p = trans[p][str[i]-'a'];
                        x = l[fa[p]]+1;
                    }
                }
            }
        }
    };
    int tt;
    void solve() {
//        cin.tie(0);
//        ios::sync_with_stdio(0);
        cin >> tt;
        int f = 0;
        while(tt--) {
            res = 1e9;
            int n;
            cin >> n;
            cin >> s;
            cin >> t;
            for(int i = 2; i < n; i++) {
                string tmp;
                t += 'z'+1;
                cin >> tmp;
                t += tmp;
            }

            SAM::init();
            SAM::build(t);
            SAM::solve(s);
            cout<<"Case #"<<++f<<": ";
//            cout<<res<<" "<<pos<<endl;
            if(res == 1e9) {
                cout<<"Impossible"<<'\n';
            } else {
                cout<<s.substr(pos-res+1, res)<<'\n';
            }
        }
    }
}
/*
1. 字典序
1
2
abc
dddd
*/
int main() {
//    file();
    solver::solve();
    return 0;
}