Codeforces645E Intellectual Inquiry(dp)

题目大意
给定一个长度为n的只包含前K个小写字母的字符串,现在要在它后面补上m个前K个小写字母,求能得到的最多的不同的子序列个数(空串也算),答案对109+7取模. 
n,m≤105,K≤26.
题解
我们先考虑怎么求一个给定字符串的不同子序列个数. 
关键是dp时怎么保证不重复计数? 
呃…说实话我只是自己瞎定义了一个状态,然后发现可以通过奇妙的转移避免重复计算,所以我也无法很好地解释思路究竟是怎样的(=.=||). 
设dp[i][j]表示从前i个字符中取出的结尾为字符j的方案数.可以写出转移方程: 
显然若str[i]≠j,那么dp[i][j]=dp[i−1][j]. 
而对str[i]=j的情况,有:
dp[i][j]=∑kdp[i−1][k]+1.


即在前面所有的方案后面都加上一个str[i],还算上一种只取str[i]的方案. 
这样转移就保证了每个dp值表示的所有方案中没有重复的,而在结尾不同的方案后面接上str[i]也不会有重复,从而避免了重复计数.加上空串后最终答案为∑dp[n][i]+1. 
于是维护一个dp值之和就可以O(n)求出给定串的不同子序列个数.当然dp的第一维可以省掉. 
然而接下来的m个字符需要我们自己构造,只要写出上面的dp方程就不难了. 
容易发现在前面的dp值确定的情况下不管怎么构造得到的dp值都是不变的,而改变的只是dp值之和.所以可以贪心地选dp值最小的进行构造. 
注意还有一点是这里的dp值是取模过的,所以不能直接比较大小,可以发现dp值显然是单调递增的,所以可以机智地用下标进行比较. 
复杂度是O(n+m⋅K).
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e6+5,MAX_K=26,P=1e9+7;
int last[MAX_K],dp[MAX_K];
char str[N];
inline void mod_add(int &a,int b){
    if((a+=b)>=P)a-=P;
}
int main(){
    int n,m,K,sum=0,ans=1;
    scanf("%d%d%s",&m,&K,str);
    n=strlen(str);
    for(int i=0;i<K;++i){
        dp[i]=0;
        last[i]=-1;
    }
    for(int i=0;i<n;++i){
        last[str[i]-='a']=i;
        int tmp=dp[str[i]];
        dp[str[i]]=(sum+1)%P;
        sum=(((sum<<1)+1-tmp)%P+P)%P;
    }
    for(int i=0;i<m;++i){
        int ptr=-1;
        for(int j=0;j<K;++j)
            if(ptr==-1||last[j]<last[ptr])ptr=j;
        int tmp=dp[ptr];
        dp[ptr]=(sum+1)%P;
        sum=(((sum<<1)+1-tmp)%P+P)%P;
        last[ptr]=n+i;
    }
    for(int i=0;i<K;++i)
        mod_add(ans,dp[i]);
    printf("%d\n",ans);
    return 0;
}
/*





*/



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
区间DP是一种动态规划的方法,用于解决区间范围内的问题。在Codeforces竞赛中,区间DP经常被用于解决一些复杂的字符串或序列相关的问题。 在区间DP中,dp[i][j]表示第一个序列前i个元素和第二个序列前j个元素的最优解。具体的转移方程会根据具体的问题而变化,但是通常会涉及到比较两个序列的元素是否相等,然后根据不同的情况进行状态转移。 对于区间长度为1的情况,可以先进行初始化,然后再通过枚举区间长度和区间左端点,计算出dp[i][j]的值。 以下是一个示例代码,展示了如何使用区间DP来解决一个字符串匹配的问题: #include <cstdio> #include <cstring> #include <string> #include <iostream> #include <algorithm> using namespace std; const int maxn=510; const int inf=0x3f3f3f3f; int n,dp[maxn][maxn]; char s[maxn]; int main() { scanf("%d", &n); scanf("%s", s + 1); for(int i = 1; i <= n; i++) dp[i][i] = 1; for(int i = 1; i <= n; i++) { if(s[i] == s[i - 1]) dp[i][i - 1] = 1; else dp[i][i - 1] = 2; } for(int len = 3; len <= n; len++) { int r; for(int l = 1; l + len - 1 <= n; l++) { r = l + len - 1; dp[l][r] = inf; if(s[l] == s[r]) dp[l][r] = min(dp[l + 1][r], dp[l][r - 1]); else { for(int k = l; k <= r; k++) { dp[l][r] = min(dp[l][r], dp[l][k] + dp[k + 1][r]); } } } } printf("%d\n", dp[n]); return 0; } 希望这个例子能帮助你理解区间DP的基本思想和应用方法。如果你还有其他问题,请随时提问。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值