首先剔除$1$號心情不能到達的無用心情，然後采用分割法進行DFA的最小化。

每次遍歷所有集合，將集合中和集合中第一個心情行為或者轉移所在集合不同的心情放入新集合中。

最後按字典序依次給每個集合編號即可。

因為最多$O(N)$次分離，每次遍歷時間復雜度為$O(NP)$。

故總時間復雜度為$O(N^2P)$。

#include<cstdio>
const int N=1010;
int p,m,n,i,j,k,x,a[N][30],cnt,len[N],g[N][N],pos[N],id[N],q[N],all;bool vis[N],flag;
inline bool equal(int x,int y){
  if(a[x][0]!=a[y][0])return 0;
  for(int i=1;i<=p;i++)if(pos[a[x][i]]!=pos[a[y][i]])return 0;
  return 1;
}
void dfs(int x){
  if(vis[x])return;
  vis[x]=1;
  g[1][++len[1]]=x;
  for(int i=1;i<=p;i++)dfs(a[x][i]);
}
int main(){
  scanf("%d%d%d",&p,&m,&n);
  for(i=1;i<=n;i++)for(j=1;j<=p;j++)scanf("%d",&a[i][j]);
  for(i=1;i<=n;i++)scanf("%d",&a[i][0]);
  dfs(cnt=1);
  while(1){
    flag=0;
    for(i=1;i<=cnt;i++)for(j=1;j<=len[i];j++)pos[g[i][j]]=i;
    for(i=1;i<=cnt;i++)if(len[i]>1){
      for(j=2;j<=len[i];j++)if(!equal(g[i][1],g[i][j]))break;
      if(j>len[i])continue;
      cnt++,k=1;
      for(j=2;j<=len[i];j++)if(!equal(g[i][1],g[i][j]))g[cnt][++len[cnt]]=g[i][j];
      else g[i][++k]=g[i][j];
      len[i]=k;
      flag=1;
      break;
    }
    if(!flag)break;
  }
  printf("%d\n",cnt);
  id[1]=q[1]=all=1;
  for(i=1;i<=cnt;i++){
    x=g[q[i]][1];
    for(j=1;j<=p;j++){
      k=pos[a[x][j]];
      if(!id[k])q[id[k]=++all]=k;
      printf("%d%c",id[k],j<p?‘ ‘:‘\n‘);
    }
  }
  for(i=1;i<=cnt;i++)printf("%d%c",a[g[q[i]][1]][0],i<cnt?‘ ‘:‘\n‘);
  return 0;
}

BZOJ4267 : 小強的顏色