牛客2021暑期训练5-D-Double Strings

牛客2021暑期训练5-D-Double Strings

题目链接

题意

对所给的两个字符串 A、B 求好的子序列个数,好的子序列为:一段相同的前缀(可以无)+ 一个不同字符(a比b小)+ 长度相同的任意后缀(可以无)

题解

f [ i ] [ j ] f[i][j] f[i][j] A [ 1... i ] A[1...i] A[1...i] B [ 1... j ] B[1...j] B[1...j] 可用前缀个数,可得递推式
f [ i ] [ j ] = f [ i − 1 ] [ j ] + f [ i ] [ j − 1 ] − f [ i − 1 ] [ j − 1 ] f[i][j]=f[i−1][j]+f[i][j−1]−f[i−1][j−1] f[i][j]=f[i1][j]+f[i][j1]f[i1][j1]
当 A[i]==B[i] 时,额外有 f [ i ] [ j ] + = f [ i − 1 ] [ j − 1 ] + 1 f[i][j]+=f[i−1][j−1]+1 f[i][j]+=f[i1][j1]+1

g [ i ] [ j ] g[i][j] g[i][j] 为从 {1, 2, …, i} 和 {1, 2, …, j} 中取出长度相同的序列的个数,等价于所有字符相等的 f f f,即
g [ i ] [ j ] = g [ i − 1 ] [ j ] + g [ i ] [ j − 1 ] + 1 g[i][j]=g[i−1][j]+g[i][j−1]+1 g[i][j]=g[i1][j]+g[i][j1]+1

因此可以先处理 f f f g g g 数组,再枚举 i , j i,j ij,当 A [ i ] < B [ j ] A[i]<B[j] A[i]<B[j] 时, a n s + = ( f [ i − 1 ] [ j − 1 ] + 1 ) × ( g [ l a − i ] [ l b − j ] + 1 ) ans+=(f[i−1][j−1]+1)×(g[la−i][lb−j]+1) ans+=(f[i1][j1]+1)×(g[lai][lbj]+1) ,加1是“可以无”

代码
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int N=5e3+9,mod=1e9+7;
int la,lb;
LL ans;
int f[N][N],g[N][N];
char a[N],b[N];
int main()
{
	cin>>a+1>>b+1;
	la=strlen(a+1);
	lb=strlen(b+1);
	for(int i=1;i<=la;i++)
		for(int j=1;j<=lb;j++)
		{
			f[i][j]=(f[i-1][j]+f[i][j-1]-f[i-1][j-1])%mod;
			if(a[i]==b[j]) f[i][j]=(f[i][j]+f[i-1][j-1]+1)%mod;
			g[i][j]=(g[i-1][j]+g[i][j-1]+1)%mod;
		}
	for(int i=1;i<=la;i++)
		for(int j=1;j<=lb;j++)
			if(a[i]<b[j])
				ans=(ans+1LL*(f[i-1][j-1]+1)*(g[la-i][lb-j]+1)%mod)%mod;
	printf("%lld",(ans+mod)%mod);
	return 0;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值