LOJ 2172 「FJOI2016」所有公共子序列问题——序列自动机

题目：https://loj.ac/problem/2172

在两个序列自动机上同时走，这样暴搜。

先走字典序小的字符，一边搜一边输出，就是按字典序排序的。

方案数很多，需要高精度？空间很小，要压位。1e9的20位恰好够。

不开 n*n 的DP数组，给出现的状态分配一个位置，开 3e6 的DP数组，空间就能了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=,K=;
int n,m,c[N][K],d[N][K],lst[K];
char a[N],b[N];
int Id(char ch)
{
  if(ch<='Z')return ch-'A';
  return ch-'a'+;
}
char Id2(int k)
{
  if(k<)return k+'A';
  return k-+'a';
}
namespace S1{
  int ans,top; char s[N];
  void dfs(int p0,int p1)
  {
    printf("%s\n",s+); ans++;
    for(int i=;i<;i++)
      if(c[p0][i]<=n&&d[p1][i]<=m)
    {
      s[++top]=Id2(i);
      dfs(c[p0][i],d[p1][i]);
      s[top]=s[top+];top--;
    }
  }
  void solve()
  {dfs(,); printf("%d\n",ans);}
}
namespace S2{
  const int M=3e6,bs=1e9;
  int tot; int dy[N][N];
  struct Node{
    int a[];//
    void Inc(int k)
    {
      a[]+=k;
      for(int i=;i<=a[];i++)
    if(a[i]>=bs)a[i]-=bs,a[i+]++;
      while(a[a[]+])a[]++;
    }
    void Inc(Node k)
    {
      int lm=max(a[],k.a[]);
      for(int i=;i<=lm;i++)
    {
      a[i]+=k.a[i];
      if(a[i]>=bs)a[i+]+=a[i]/bs,a[i]%=bs;
    }
      while(a[a[]+])a[]++;
    }
    void print()
    {
      printf("%d",a[a[]]);
      for(int i=a[]-;i;i--)
    printf("%09d",a[i]);puts("");
    }
  }dp[M];
  int dfs(int p0,int p1)
  {
    if(dy[p0][p1])return dy[p0][p1];
    int cr=++tot; dy[p0][p1]=cr; dp[cr].Inc();
    for(int i=,x,y;i<;i++)
      if((x=c[p0][i])<=n&&(y=d[p1][i])<=m)
    {
      int v=dfs(x,y); dp[cr].Inc(dp[v]);
    }
    return cr;
  }
  void solve()
  { int v=dfs(,); dp[v].print();}
}
int main()
{
  scanf("%d%d",&n,&m);
  scanf("%s",a+); scanf("%s",b+);
  for(int i=;i<;i++)lst[i]=n+;
  for(int i=n;i>=;i--)
    {
      for(int j=;j<;j++)
    c[i][j]=lst[j];
      lst[Id(a[i])]=i;
    }
  for(int i=;i<;i++)lst[i]=m+;
  for(int i=m;i>=;i--)
    {
      for(int j=;j<;j++)
    d[i][j]=lst[j];
      lst[Id(b[i])]=i;
    }
  int op;scanf("%d",&op);
  if(op==)S1::solve();
  else S2::solve();
  return ;
}