【LOJ】#2172. 「FJOI2016」所有公共子序列问题

题解

听说是什么序列自动机？
我们考虑对于每个位置的串，下面拼接相同的字符时，拼接最近的一个，这样可以保证不重不漏
为了实现这个我们需要什么呢，我们需要一个链表，记录一下每个位置的下一个字符会转移到哪里
例如
ABAB
ch[1]['A'] = 3 ch[1]['B'] = 2
听起来挺好建的，具体可以看代码

然后我们记录dp[a][b]是第一个字符串a位置的字符为开头，第二个字符串b位置的字符为开头，公共子序列的个数
可以记忆化搜索，每一次枚举下一个位置拼接上哪个字符
输出方案的话也是一样，枚举下一个位置拼上哪个字符

后五个点需要高精度啊= =

#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <map>
//#define ivorysi
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define mo 974711
#define MAXN 3055
#define RG register
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
    res = 0;char c = getchar();T f = 1;
        while(c < '0' || c > '9') {
            if(c == '-') f = -1;
            c = getchar();
        }
        while(c >= '0' && c <= '9') {
        res = res * 10 + c - '0';
        c = getchar();
    }
    res *= f;
}
template<class T>
void out(T x) {
    if(x < 0) {putchar('-');x = -x;}
    if(x >= 10) {
        out(x / 10);
    }
    putchar('0' + x % 10);
}
int N,M,MK,ci[205];
char x[MAXN],y[MAXN],ic[205],s[MAXN];
struct segAM {
    int head[62],next[MAXN],ch[MAXN][62],tot;
    segAM() {
        tot = 1;next[1] = 0;
        for(int i = 0 ; i <= 60; ++i) head[i] = 1;
    }
    void Insert(char c) {
        ++tot;
        next[tot] = head[ci[c]];
        for(int i = 0 ; i <= 51 ; ++i) {
            for(int j = head[i] ; j && !ch[j][ci[c]] ; j = next[j]) {
                ch[j][ci[c]] = tot;
            }
        }
        head[ci[c]] = tot;
    }
}S1,S2;
const int BASE = 100000000;
const int LEN = 8;
struct Bignum {
    vector<int> v;
    Bignum operator = (int64 x) {
        v.clear();
        do {
            v.pb(x % BASE);
            x /= BASE;
        }while(x);
        return *this;
    }
    friend Bignum operator + (const Bignum &a,const Bignum &b) {
        Bignum c;c.v.clear();
        int p = 0,g = 0;
        while(1) {
            int x = g;
            if(p < a.v.size()) x += a.v[p];
            if(p < b.v.size()) x += b.v[p];
            ++p;
            if(!x) break;
            c.v.pb(x % BASE);
            g = x / BASE;
        }
        return c;
    }
    void print() {
        int s = v.size() - 1;
        printf("%d",v[s]);
        --s;
        for(int i = s ; i >= 0 ; --i) {
            printf("%08d",v[i]);
        }
    }
}dp[MAXN][MAXN];
bool vis[MAXN][MAXN];
void get(int a,int b) {
    if(a == 0 || b == 0) return;
    if(vis[a][b]) return;
    dp[a][b] = 1;
    vis[a][b] = 1;
    for(int i = 0 ; i <= 51 ; ++i) {
        int ta = S1.ch[a][i],tb = S2.ch[b][i];
        get(ta,tb);
        dp[a][b] = dp[a][b] + dp[ta][tb];
    }
}
void get2(int a,int b,int dep) {
    if(a == 0 || b == 0) return;
    s[dep] = 0;if(dep != 0) printf("%s\n",s);
    for(int i = 0 ; i <= 51 ; ++i) {
        s[dep] = ic[i];
        get2(S1.ch[a][i],S2.ch[b][i],dep + 1);
    }
}
void Solve() {
    read(N);read(M);
    scanf("%s%s",x + 1,y + 1);
    read(MK);
    for(char i = 'A' ; i <= 'Z' ; ++i) 
        ci[i] = i - 'A',ic[i - 'A'] = i;
    for(char i = 'a' ; i <= 'z' ; ++i) 
        ci[i] = i - 'a' + 26,ic[i - 'a' + 26] = i;
    for(int i = 1 ; i <= N ; ++i) S1.Insert(x[i]);
    for(int i = 1 ; i <= M ; ++i) S2.Insert(y[i]); 
    if(MK == 1) get2(1,1,0);
    get(1,1);
    dp[1][1].print();enter;
}
int main() {
#ifdef ivorysi
    freopen("f1.in","r",stdin);
#endif
    Solve();
    return 0;
}

今天写的题解全是FJOI啊（其实还写了2道JOI的水题，太水了不写题解了）
zxg说让我们看看数学相关的集训队论文，然后做做FJ和SC的省选= =
FJ猫锟 SC YJQ
我们都表示毫不理解该省省选题和该省出题人有什么关系，真的不懂中老年人思维啊
但是反正都是要无脑提高刷题量嘛（大雾）
量变能带来质变
……吗？

39天NOI2018打卡

我想猫锟了QAQ