CF49E Common ancestor（dp+dp+dp）

纪念卡常把自己卡死的一次自闭模拟赛

QWQ
一开始看这个题，以为是个图论，仔细一想，貌似可以直接dp啊。

首先，因为规则只有从两个变为1个，貌似可以用类似区间(dp)的方式来(check)一段区间能不能合成某一个字母！

那我们定义(f[i][j][k])表示第一个串，([l,r])区间，是否可以合成(k)这个字母

然后转移的时候，枚举区间，枚举规则，枚举断点，满足(f[l][k][p1]==1)且(f[k+1][r][p2]==1) 才能使当前状态合法。
其中(p1,p2)表示当前规则的两个字母

for (int i=1;i<=n;i++) f[i][i][cc(s[i])]=1; 
  for (register int i=2;i<=n;++i)
    for (register int l=1;l<=n-i+1;++l)
      {
      	 int r = l+i-1;
      	 for (register int j=1;j<=26;++j)
      	 {
      	 	for (register int p=1;p<=num[j];++p)
      	 	{
      	 	  for (register int k=l;k<=r;++k) 
      	 	  {
      	 	    f[l][r][j]=max(f[l][r][j],f[l][k][a[j][p].a]&f[k+1][r][a[j][p].b]);
      	 	    if (f[l][r][j]) break;
      	 	  }
      	 	  if (f[l][r][j]) break;
      	    }
		 }
	  }

同时定义(g[l][r][k])数组表示第二个串区间([l,r])能否合成k。处理和f类似。

统计答案的时候呢
还需要一个(dp[i][j])表示第一个串的前i个字符和第二个串的前j个字符的最短公共祖先
那么，考虑枚举两个断点，两个串的后面两段能合成同一个字母，那么就可以从那个断点之前的状态转移过来

QWQ
详细直接看代码吧

memset(dp,127/3,sizeof(dp));
  dp[0][0]=0;
  for (register int i=1;i<=nn;++i)
  {
  	 for (register int k=1;k<=n;++k)
  	 {
  	   for (register int j=1;j<=i;++j)
  	     for (register int p=1;p<=k;++p)
  	     {
  	     	if (dp[j-1][p-1]==dp[maxn-3][maxn-3]) continue;
  	 	    bool flag=false;
  	 	    for (register int o=1;o<=26;o++)
  	 	      if (g[j][i][o] && f[p][k][o]) flag=true;
  	        if (flag) dp[i][k]=min(dp[i][k],dp[j-1][p-1]+1);
  	     }
     }
  }

最后复杂度就是(O(n^4*26))

我也不知道为啥能跑过啊
qwqwqwq

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<set>
#define mk makr_pair
#define ll long long
#include<ctime>
using namespace std;
inline int read()
{
  int x=0,f=1;char ch=getchar();
  while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
  while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
  return x*f;
}
const int maxn = 110;
struct Node{
	int a,b;
};
Node a[maxn][maxn];
int num[maxn];
int f[maxn][maxn][maxn];
int g[maxn][maxn][maxn];
int n,m;
char s[maxn];
char ss[maxn];
int nn;
string ans;
int dp[maxn][maxn];
inline int cc(char c)
{
	return c-'a'+1;
}
int main()
{
  scanf("%s",s+1);
  scanf("%s",ss+1);
  s[0]=ss[0]='*';
  n=strlen(s+1);
  nn=strlen(ss+1);  
  m=read();
  for (register int i=1;i<=m;++i)
  {
  	 char ymh[10];
  	 scanf("%s",ymh+1);
  	 int now = ymh[1]-'a'+1;
  	 num[now]++;
  	 a[now][num[now]].a = ymh[4]-'a'+1;
	 a[now][num[now]].b = ymh[5]-'a'+1; 
  }
  for (int i=1;i<=n;i++) f[i][i][cc(s[i])]=1; 
  for (register int i=2;i<=n;++i)
    for (register int l=1;l<=n-i+1;++l)
      {
      	 int r = l+i-1;
      	 for (register int j=1;j<=26;++j)
      	 {
      	 	for (register int p=1;p<=num[j];++p)
      	 	{
      	 	  for (register int k=l;k<=r;++k) 
      	 	  {
      	 	    f[l][r][j]=max(f[l][r][j],f[l][k][a[j][p].a]&f[k+1][r][a[j][p].b]);
      	 	    if (f[l][r][j]) break;
      	 	  }
      	 	  if (f[l][r][j]) break;
      	    }
		 }
	  }
  for (int i=1;i<=nn;i++) g[i][i][cc(ss[i])]=1;
  for (register int i=2;i<=nn;i++)
    for (register int l=1;l<=nn-i+1;++l)
      {
      	 int r = l+i-1;
      	 for (register int j=1;j<=26;++j)
      	 {
      	 	for (register int p=1;p<=num[j];++p)
      	 	{
      	 	  for (register int k=l;k<=r;++k)
			  { 
      	 	    g[l][r][j]=max(g[l][r][j],g[l][k][a[j][p].a]&g[k+1][r][a[j][p].b]);
      	 	    if (g[l][r][j]) break;
      	 	  }
			  if (g[l][r][j]) break;
      	    }
		 }
	  }
  memset(dp,127/3,sizeof(dp));
  dp[0][0]=0;
  for (register int i=1;i<=nn;++i)
  {
  	 for (register int k=1;k<=n;++k)
  	 {
  	   for (register int j=1;j<=i;++j)
  	     for (register int p=1;p<=k;++p)
  	     {
  	     	if (dp[j-1][p-1]==dp[maxn-3][maxn-3]) continue;
  	 	    bool flag=false;
  	 	    for (register int o=1;o<=26;o++)
  	 	      if (g[j][i][o] && f[p][k][o]) flag=true;
  	        if (flag) dp[i][k]=min(dp[i][k],dp[j-1][p-1]+1);
  	     }
     }
  } 
  if(dp[nn][n]==dp[maxn-3][maxn-3]) dp[nn][n]=-1; 
  cout<<dp[nn][n]<<endl;
  return 0;
}