「BZOJ3864」Hero meet devil 题解

简要题意

给你一个只由 A G C T AGCT AGCT 组成的字符串 S S S ,对于每个 0 ≤ i ≤ ∣ S ∣ 0 \leq i \leq |S| 0iS,问有多少个只由 A G C T AGCT AGCT 组成的长度为 m m m 的字符串 T T T,使得 L C S ( S , T ) = i LCS(S,T)=i LCS(S,T)=i

SOLUTION

刚开始想的时候,怎么定义状态都会造成重复等各种问题,于是搜题解(emmm)。。。

考虑求 L C S LCS LCS时的 d p dp dp:定义 l c s [ i ] [ j ] lcs[i][j] lcs[i][j] 表示 A A A串的前 i i i位与B串的前 j j j位的 L C S LCS LCS长度

l c s [ i ] [ j ] = m a x ( l c s [ i − 1 ] [ j − 1 ] + ( A [ i ] = = B [ i ] ) , l c s [ i ] [ j − 1 ] , l c s [ i − 1 ] [ j ] ) lcs[i][j]=max({lcs[i−1][j−1]+ (A[i]==B[i]),lcs[i][j−1],lcs[i−1][j]}) lcs[i][j]=max(lcs[i1][j1]+(A[i]==B[i])lcs[i][j1]lcs[i1][j])

可以发现,当 i i i 相同时, l c s [ i ] [ j ] lcs[i][j] lcs[i][j] l c s [ i ] [ j + 1 ] lcs[i][j+1] lcs[i][j+1] 最多相差1

把每一位之间的差状压,表示为一个状态进行计数,就可以了

i i i位是 1 1 1意味着第 i i i上的字符对 l c s lcs lcs有贡献,即是1的位们就是公共子串部分

定义 d p [ i ] [ s ] dp[i][s] dp[i][s] 表示:长度为 i i i 的字符串,与 S S S l c s lcs lcs 状态为 s s s 的方案数

t r a n s ( s , k ) trans(s,k) trans(s,k) 表示在 l c s lcs lcs 状态为 s s s 的情况下,在后面添加字符 k k k 转移到的状态

在不同串但是 l c s lcs lcs同种状态的时候,在后面加上同一个字符,转移到的状态是一样的

转移方程: d p [ i ] [ t r a n s ( s , k ) ] + = d p [ i − 1 ] [ s ] dp[i][trans(s,k)]+=dp[i−1][s] dp[i][trans(s,k)]+=dp[i1][s]

注意:先想清楚求 l c s lcs lcs的过程在来看这个trans的方程转移

CODE

#include<bits/stdc++.h> 
using namespace std;
const int M=1e3+2,L=20,N=(1<<15)+2,P=1e9+7;
char ch[]="ACGT";
int T,n,m,dp[2][N],tr[N][4],ans[M],f[L],g[L];
char S[L];
void transwork(){
	for(int s=0;s<(1<<n);s++){
		memset(f,0,sizeof f);
		memset(g,0,sizeof g);
		for(int j=1;j<=n;j++){
			f[j]=f[j-1]+((s>>(j-1))&1);
		}//f[j]=lcs[n][j]
		for(int k=0;k<4;k++){//要加的字符 
			for(int j=1;j<=n;j++){
				g[j]=max(g[j-1],f[j]);
				if(ch[k]==S[j]){//可以有贡献 
					g[j]=max(g[j],f[j-1]+1);
				}//lcs方程过程 
			}
			tr[s][k]=0;//trans清0 
			for(int j=1;j<=n;j++){
				if(g[j]-g[j-1]){//有贡献 
					tr[s][k]|=1<<(j-1);
				}
			}		
		}
	}
}
int main(){
	scanf("%d",&T);
	while(T--){
		scanf("%s%d",S+1,&m);
		n=strlen(S+1);
		transwork(); 
		memset(dp,0,sizeof dp);
		dp[0][0]=1;
		for(int i=1;i<=m;i++){
			memset(dp[i&1],0,sizeof dp[i&1]);
			for(int s=0;s<(1<<n);s++){
				for(int k=0;k<4;k++){
					(dp[i&1][tr[s][k]]+=dp[(i&1)^1][s])%=P;
				}
			}
		}
		memset(ans,0,sizeof ans);
		for(int s=0;s<(1<<n);s++){
			(ans[__builtin_popcount(s)]+=dp[m&1][s])%=P;
			//有多少位为1
		}
		for(int i=0;i<=n;i++){
			printf("%d\n",ans[i]);
		}
	}
	return 0;
}

完结撒花❀

★,°:.☆( ̄▽ ̄)/$:.°★

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值