洛谷 P1140 相似基因 题解

每日一题 day23 打卡

Analysis

dp[i][j]表示序列A中前i个与序列B中前j个匹配的相似度最大值

所以，dp方程很容易想到：

1.让a[i]与b[j]匹配

2.让a[i]与B序列中一个空位匹配

3.让b[j]与A序列中一个空位匹配

dp[i][j]=max(dp[i][j],dp[i-1][j-1]+form[a[i]][b[j]],dp[i-1][j]+form[a[i]][5],dp[i][j-1]+form[5][b[j]]);

对于初始化，就是分别让每个a[i]和每个b[j]与一个空串匹配

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define int long long
 #define maxn 100+10
 #define INF 2147483647
 using namespace std;
 inline int read()
 {
     int x=;
     bool f=;
     char c=getchar();
     for(; !isdigit(c); c=getchar()) if(c=='-') f=;
     for(; isdigit(c); c=getchar()) x=(x<<)+(x<<)+c-'';
     if(f) return x;
     return -x;
 }
 inline void write(int x)
 {
     if(x<){putchar('-');x=-x;}
     if(x>)write(x/);
     putchar(x%+'');
 }
 int n,m;
 char s1[maxn],s2[maxn];
 int str1[maxn],str2[maxn],dp[maxn][maxn];
 int form[maxn][maxn]=
 {
     {,,,,,,},
     {,,-,-,-,-},
     {,-,,-,-,-},
     {,-,-,,-,-},
     {,-,-,-,,-},
     {,-,-,-,-,}
 };
 inline int max_four(int x,int y,int z,int o)
 {
     return max(max(x,y),max(z,o));
 }
 signed main()
 {
     n=read();
     for(int i=;i<=n;i++) cin>>s1[i];
     m=read();
     for(int i=;i<=m;i++) cin>>s2[i];
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
             dp[i][j]=-INF;
     for(int i=;i<=n;i++)
     {
         if(s1[i]=='A') str1[i]=;
         else if(s1[i]=='C') str1[i]=;
         else if(s1[i]=='G') str1[i]=;
         else if(s1[i]=='T') str1[i]=;
     }
     for(int i=;i<=m;i++)
     {
         if(s2[i]=='A') str2[i]=;
         else if(s2[i]=='C') str2[i]=;
         else if(s2[i]=='G') str2[i]=;
         else if(s2[i]=='T') str2[i]=;
     }
     for(int i=;i<=n;i++) dp[i][]=dp[i-][]+form[str1[i]][];
     for(int i=;i<=m;i++) dp[][i]=dp[][i-]+form[][str2[i]];
     for(int i=;i<=n;i++)
         for(int j=;j<=m;j++)
         {
             dp[i][j]=max_four(dp[i][j],dp[i-][j-]+form[str1[i]][str2[j]],dp[i-][j]+form[str1[i]][],dp[i][j-]+form[][str2[j]]);
         }
     write(dp[n][m]);
     return ;
 }

请各位大佬斧正（反正我不认识斧正是什么意思）