合并回文子串

n = int(input())
for _ in range(n):
    a, b = input(), input()
    s, t = len(a), len(b)
    dp = [
        [[[0 for l in range(t + 2)] for k in range(t + 2)] for j in range(s + 2)]
        for i in range(s + 2)
    ]
    ans = 0

    for x in range(s + 1):  # a串中的子串长度
        for y in range(t + 1):  # b串中的子串长度
            #             空串：i=j，x=j-i=0,子串为[i:j] 范围：0<=i=j-x<=s-x
            for i in range(s - x + 1):  # a串右端点
                j = i + x  # a串左端点
                for k in range(t - y + 1):  # b串右端点
                    l = k + y  # b串左端点
                    if x == 0 and y == 0:
                        dp[i][j][k][l] = 1
                    elif x == 0:
                        if y == 1:
                            dp[i][j][k][l] = 1
                        elif b[k] == b[l - 1]:
                            dp[i][j][k][l] = dp[i][j][k + 1][l - 1]
                    elif y == 0:
                        if x == 1:
                            dp[i][j][k][l] = 1
                        elif a[i] == a[j - 1]:
                            dp[i][j][k][l] = dp[i + 1][j - 1][k][l]

                    elif (i < s and k < t) and (
                        (a[i] == a[j - 1] and x > 1 and dp[i + 1][j - 1][k][l])
                       &nbs***bsp;(b[k] == b[l - 1] and y > 1 and dp[i][j][k + 1][l - 1])
                       &nbs***bsp;(l > 0 and a[i] == b[l - 1] and dp[i + 1][j][k][l - 1])
                       &nbs***bsp;(j > 0 and a[j - 1] == b[k] and dp[i][j - 1][k + 1][l])
                    ):
                        dp[i][j][k][l] = 1
                        ans = max(ans, x + y)
    print(ans)

#include<iostream>
#include<string>
#include<cstring>
using namespace std;

const int INF=0x3f3f3f3f;

string a,b;
int cache[51][51][51][51];

int dp(int al,int ar,int bl,int br){
  if(al==ar&&bl==br) return 0;
  if((al+1==ar && bl==br)||(al==ar && bl+1==br)) return 1; 

  if(cache[al][ar][bl][br]) return cache[al][ar][bl][br];
  int maxl=-INF;
  if(ar-al>=2 && a[al] == a[ar-1]) maxl= max(maxl,dp(al+1,ar-1,bl,br)+2);
  if(br-bl>=2 && b[bl] == b[br-1]) maxl= max(maxl,dp(al,ar,bl+1,br-1)+2);
  if(br-bl>=1 && ar-al>=1){
    if(a[al] == b[br-1]) maxl= max(maxl,dp(al+1,ar,bl,br-1)+2);
    if(a[ar-1] == b[bl]) maxl= max(maxl,dp(al,ar-1,bl+1,br)+2);
  }
  return cache[al][ar][bl][br] =maxl;
}

int main(){
  int n;
  cin>>n;
  for(int i=0;i<n;i++){
    cin>>a>>b;
    
    memset(cache,0,sizeof(cache));
    int res=-INF;
    
    for(int i=0;i<a.size();i++)
      for(int j=i+1;j<=a.size();j++)
        for(int k=0;k<b.size();k++)
          for(int l=k+1;l<=b.size();l++)
            res=max(res,dp(i,j,k,l));

    cout<<res<<endl;
  }
}