以后遇见题在加上去。
一、字符串线性DP
1.Cheapest Palindrome(区间dp)
题目链接
d[i][j]
:从i
到j
为回文串花费的最小
状态转移方程:
if(s[i]==s[j])
d[i][j]=d[i+1][j-1]
else
d[i][j]=min(d[i][j-1]+add[s[j]],d[i+1][j]+add[s[i]],del[s[i]]+d[i+1][j],del[s[j]]+d[i][j-1])
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N = 2010;
int f[N][N];
int n,m;
int add[N],del[N];
int main()
{
cin>>n>>m;
string s;
cin>>s;s=" "+s;
char op[2];
int x,y;
for(int i=0;i<n;i++)
{
cin>>op;
cin>>x>>y;
add[op[0]-'a']=x;
del[op[0]-'a']=y;
}
memset(f,0x3f,sizeof f);
f[0][0]=0;
for(int len=1;len<=m;len++)
{
for(int i=1;i+len-1<=m;i++)
{
int j=i+len-1;
if(s[i]==s[j])
{
if(len<=2)f[i][j]=0;
else f[i][j]=f[i+1][j-1];
}
else
{
int x=min(f[i+1][j]+add[s[i]-'a'],f[i][j-1]+add[s[j]-'a']);
int y=min(f[i+1][j]+del[s[i]-'a'],f[i][j-1]+del[s[j]-'a']);
f[i][j]=min(x,y);
}
}
}
cout<<f[1][m]<<endl;
}
2.Palindrome
题目链接
MLE做法:(类似于1的区间dp)
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=5010;
int f[N][N];
int main()
{
int n;
string s;
cin>>n;
cin>>s;
s=" "+s;
for(int len=1;len<=n;len++)
{
for(int i=1;i+len-1<=n;i++)
{
int j=i+len-1;
if(len==1)f[i][j]=0;
else if(s[i]==s[j])
{
if(len==2)f[i][j]=0;
else f[i][j]=f[i+1][j-1];
}
else
{
f[i][j]=min(f[i+1][j],f[i][j-1])+1;
}
}
}
cout<<f[1][n]<<endl;
}
AC做法
res=n-LCS
因为字符串和它的反向字符串的最大公共子序列的长度LCS,就等同于在原串中有LCS个字符是对应相等的。答案就是那些对应不相等的在对应位置添加相同的字符。所以res=n-LCS;
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=5010;
int f[2][N]={0};
int main()
{
int n;
cin>>n;
string a,b;
cin>>b;a=b;
reverse(b.begin(),b.end());
a=" "+a;
b=" "+b;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
f[i&1][j]=max(f[i-1&1][j],f[i&1][j-1]);
if(a[i]==b[j])
{
f[i&1][j]=max(f[i&1][j],f[i-1&1][j-1]+1);
}
}
}
cout<<n-f[n&1][n]<<endl;
}
3.Nun Heh Heh Aaaaaaaaaaa----子序列个数问题
题目链接
动态规划+快速幂+后缀和
复杂的O(n)
详细内容请看代码注解
先不考虑a,先算芳香字符串非a个数,最后个数乘以2的a个数次方-1,相加取模。
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 1e5+5;
typedef long long LL;
int dp[10][N];
int a[N];//a的后缀数量
const int mod = 998244353;
//预处理后缀中a的个数
void init(string s)
{
memset(a,0,sizeof a);
for(int i=s.size();i>=1;i--)
{
a[i]=a[i+1]+(s[i]=='a');
}
}
// 这里别不小心用了龟速乘,我刚才不小心用了龟速乘一直TLE
//后面有几个a,每个a可以选也可以不选所以是2的a次方,但a个数不能为0所以-1
LL qmi(int b)
{
LL res=1%mod;
LL a=2;
while(b){
if(b&1)res=res*a%mod;
a=a*a%mod;
b>>=1;
}
return (res-1)%mod;
}
int main()
{
std::ios::sync_with_stdio(false);
std::cin.tie(nullptr);
int T;
cin>>T;
while(T--)
{
string s;
cin>>s;
//使下标从1开始
s=" "+s;
init(s);
//string ss=" nunhehheh";
//dp过程 nunhehheh
//dp[j][i]表示s前i个字符中ss[1~j]字符串出现次数
LL res=0;
for(int i=1;i<=s.size();i++)
{
dp[1][i]=(dp[1][i-1]+(s[i]=='n'))%mod;
dp[2][i]=(dp[2][i-1]+dp[1][i-1]*(s[i]=='u'))%mod;
dp[3][i]=(dp[3][i-1]+dp[2][i-1]*(s[i]=='n'))%mod;
dp[4][i]=(dp[4][i-1]+dp[3][i-1]*(s[i]=='h'))%mod;
dp[5][i]=(dp[5][i-1]+dp[4][i-1]*(s[i]=='e'))%mod;
dp[6][i]=(dp[6][i-1]+dp[5][i-1]*(s[i]=='h'))%mod;
dp[7][i]=(dp[7][i-1]+dp[6][i-1]*(s[i]=='h'))%mod;
dp[8][i]=(dp[8][i-1]+dp[7][i-1]*(s[i]=='e'))%mod;
dp[9][i]=(dp[8][i-1]*(s[i]=='h'))%mod;
if(dp[9][i])res=res%mod+(dp[9][i]*qmi(a[i+1]))%mod;
res%=mod;
}
cout<<res<<endl;
}
}
分享leetcode相似题目
4.不同的子序列
题目链接
(动态规划) O(nm)
可以换一种考虑问题的方式:用S
中的字符,按顺序匹配T
中的字符,问有多少种方式可以匹配完T
中的所有字符。
可以用动态规划来做:
f[i][j]
表示用 S
的前 i
个字符,能匹配完T
的前j
个字符的方案数。
初始化:因为 S
可以从任意一个字符开始匹配,所以 f[i][0]=1,∀i∈[0,len(S)]。
状态转移:
如果 S[i−1]≠T[j−1]
,则 S[i−1]
不能匹配 T[j−1]
,所以f[i][j]=f[i−1][j];
如果 S[i−1]=T[j−1]
,则 S[i−1]
既可以匹配 T[j−1]
,也可以不匹配 T[j−1]
,所以 f[i][j]=f[i−1][j]+f[i−1][j−1]
;
时间复杂度分析:假设 S
的长度是n
,T
的长度是 m
,则共有 nm
个状态,状态转移的复杂度是 O(1)
,所以总时间复杂度是 O(nm)
。
题解 摘自yxc
作者:yxc
来源:AcWing
class Solution {
public:
int numDistinct(string s, string t) {
int n=s.size(),m=t.size();
vector<vector<unsigned long long>>f(n+1,vector<unsigned long long>(m+1));
for(int i=0;i<=n;i++)f[i][0]=1;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
f[i][j]=f[i-1][j];
if(s[i-1]==t[j-1])
f[i][j]+=f[i-1][j-1];
}
}
return f[n][m];
}
};
二、字符串区间DP
1.字符串的最长回文子序列
题目链接
class Solution {
public:
int longestPalindromeSubseq(string s) {
int n=s.size();
s=" "+s;
vector<vector<int>>f(n+1,vector<int>(n+1));
for(int len=1;len<=n;len++)
{
for(int i=1;i+len-1<=n;i++)
{
int j=i+len-1;
if(i==j)f[i][i]=1;
else
{
f[i][j]=max(f[i+1][j],f[i][j-1]);
if(s[i]==s[j])f[i][j]=max(f[i+1][j-1]+2,f[i][j]);
}
}
}
return f[1][n];
}
};
2.字符串的最长回文子串
题目链接
最优解法为:O(n)的manacher算法
本题用区间dp
class Solution {
public:
int getLongestPalindrome(string s) {
int n=s.size();
s=" "+s;
int res=0;
vector<vector<int>>f(n+1,vector<int>(n+1));
for(int len=1;len<=n;len++)
{
for(int i=1;i+len-1<=n;i++)
{
int j=i+len-1;
if(s[i]==s[j])
{
if(j-i<3)f[i][j]=true;
else f[i][j]=f[i+1][j-1];
}
else f[i][j]=false;
if(f[i][j])res=max(res,j-i+1);
}
}
return res;
}
};
3.合并回文子串
题目链接
相当于在最长回文子串的基础上扩展二维
qcjj详细题解入口:传送门
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=110;
void solve()
{
bool f[N][N][N][N];
int res=0;
string a,b;
cin>>a>>b;
int n=a.size(),m=b.size();
a=" "+a,b=" "+b;
for(int len1=0;len1<=n;len1++)
{
for(int len2=0;len2<=m;len2++)
{
for(int i=1;i+len1-1<=n;i++)
{
for(int k=1;k+len2-1<=m;k++)
{
int j=i+len1-1,l=k+len2-1;
if(len1+len2<=1)f[i][j][k][l]=1;
else
{
bool &v=f[i][j][k][l];
v=0;
if(len1>1)v|=(f[i+1][j-1][k][l]&&(a[i]==a[j]));
if(len1&&len2)v|=(f[i+1][j][k][l-1]&&(a[i]==b[l]));
if(len1&&len2)v|=(f[i][j-1][k+1][l]&&(a[j]==b[k]));
if(len2>1)v|=(f[i][j][k+1][l-1]&&(b[k]==b[l]));
}
if(f[i][j][k][l])res=max(len1+len2,res);
}
}
}
}
cout<<res<<endl;
}
int main()
{
int T;
cin>>T;
while(T--)
{
solve();
}
}
4.小A的回文串
题目链接
这题很简单,区间dp常见的循环问题,扩展到2n就行了
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
const int N=5005;
bool f[2*N][2*N];
void solve()
{
string s;
cin>>s;
int n=s.size();
s+=s;
s=" "+s;
int res=0;
for(int len=1;len<=n;len++)
{
for(int i=1;i+len-1<=2*n;i++)
{
int j=i+len-1;
if(s[i]==s[j])
{
if(j-i<3)f[i][j]=true;
else f[i][j]=f[i+1][j-1];
}
else f[i][j]=false;
if(f[i][j])res=max(res,j-i+1);
}
}
cout<<res<<endl;
}
int main()
{
solve();
}