2018-09-12

原题传送门(洛谷)https://www.luogu.org/problemnew/show/P2679


模拟考试的时候完全没有想到 正确的DP方程呢

本来写了一个大致是对的转移方程 结果算了一下空间 大概不够 就放弃了(第一维可以滚动 掉的啊喂 傻孩子啊) 愣是写了50分暴力

下面是五十分的原代码

 #include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring> #define For(i,a,b) for(register int i=a;i<=b;++i)
#define Dwn(i,a,b) for(register int i=a;i>=b;--i)
#define Re register using namespace std;
const int md=1e9+; int na,nb,kx;
string As,Bs;
int sum=;
char c[];
int cf[]; int main(){
// freopen("ex2.in","r",stdin);
freopen("substring.in","r",stdin);
freopen("substring.out","w",stdout);
cin>>na>>nb>>kx;
cin>>As; cin>>Bs;
if(kx==){
int px=;
while(){
int p=As.find(Bs,px);
if(p==-)break;
sum++; px=p+;
if(sum>=md)sum%=md;
}
cout<<sum%md<<endl; return ;
}
if(kx==){
string s1,s2;
For(i,,nb-){
s1=""; s2="";
For(j,,i)s1=s1+Bs[j];
For(j,i+,nb-)s2=s2+Bs[j];
// cout<<s1<<" "<<s2<<" "<<sum<<endl; int px1,px2;
int sum2=;
px1=px2=; while(){
int p1=As.find(s1,px1);
if(p1==-)break;
sum2=;
px2=p1+s1.size();
if(p1>na-)break;
while(){
int p2=As.find(s2,px2);
if(p2==-)break;
sum2++; px2=p2+;
if(sum2>=md)sum2%=md;
}
sum+=sum2; if(sum>=md)sum%=md;
px1=p1+;
}
}
cout<<sum%md<<endl;
return ;
}
if(kx==nb){
For(i,,nb-)c[i+]=Bs[i],cf[i+]=;
For(i,,na-){
char ac=As[i];
Dwn(j,nb,){
if(ac==c[j]){
if(j==)cf[j]++;
else cf[j]+=cf[j-];
cf[j]%=md;
}
}
}
sum=cf[nb]%md;
cout<<sum<<endl;
return ;
}
}

当k=1,(即只能分一个字串) 直接s.find()去搜就好了

当k=2,(即分两个字串)还是暴力用find()去搜

然后 当k=m (即与B串等长)想了一个挺巧妙的方法。。 用 cf[i] 数组表示此时的B串前 i 位的方案数  按顺序 考虑 A串的一个字符

从后往前枚举 i  当 A[j]== B[i]  cf[i]+=cf[i-1]  (算了算了都是暴力 不看也罢


正确 方法

多维动归 !

用我们 奥赛教练的话来说就是

每次只要一考到 多维动归 我就替学生捏一把汗 这道题是能否拿到绝对高分的关键!!

然鹅 模拟测试时我没做出来呢

好了 下面正经写题解

DP数组  f[i][j][k][0/1] 表示A的前i位 B的前j位 用了k个字串 A的第i个字串是否使用的方案数

那么我们很容易就可以得到转移方程

A[i]==B[j]

1. 不取(0)f[i][j][k][0] <-- f[i-1][j][k][0] + f[i-1][j][k][1]

2. 取  (1)f[i][j][k][1] <-- f[i-1][j-1][k][1] (与前方字符连成一个字串)+ f[i-1][j-1][k-1][1](自己单独另起一串)+ f[i-1][j-1][k-1][0]

A[i]!=B[j]

1. 不取(0)f[i][j][k][0] <-- f[i-1][j][k][0] + f[i-1][j][k][1]

2. 取 显然只能等于 0

初始边界 所有的 f[i][0][0][0]=1

然后愉快地把第一维 i 用滚动数组滚掉即可

复杂度 O(nmk)

附上代码

 // luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm> #define For(i,a,b) for(register int i=a;i<=b;++i)
#define Re register
using namespace std;
const int Na=1e3+,Nb=;
const int md=1e9+;
long long f[][Nb][Nb][];
char A[Na],B[Nb];
int na,nb,kx;
int main(){
// freopen("substring9.in","r",stdin);
// freopen("substring.out","w",stdout);
scanf("%d%d%d",&na,&nb,&kx);
scanf("%s",A+); scanf("%s",B+);
f[][][][]=f[][][][]=;
For(i,,na){
int now=i%;
For(j,,nb) For(k,,kx){
if(A[i]==B[j]){
f[now][j][k][]=(f[now^][j-][k][]+f[now^][j-][k-][]+f[now^][j-][k-][])%md;
f[now][j][k][]=(f[now^][j][k][]+f[now^][j][k][])%md;
}else{
f[now][j][k][]=;
f[now][j][k][]=(f[now^][j][k][]+f[now^][j][k][])%md;
}
}
}
long long fn=f[na%][nb][kx][]+f[na%][nb][kx][];
cout<<fn%md<<endl;
fclose(stdin); fclose(stdout);
return ; }

【NOIP2015】 Day2 T2 字串 (多维动归)的更多相关文章

  1. 【NOIP2015】字串

    [NOIP2015]字串 标签: DP NOIP Description 有两个仅包含小写英文字母的字符串 A 和 B.现在要从字符串 A 中取出 k 个互不重叠的非空子串,然后把这 k 个子串按照其 ...

  2. [NOIP2002]字串变换 T2 双向BFS

    题目描述 已知有两个字串  A,B  及一组字串变换的规则(至多6个规则): A1−>B1 A2−>B2 规则的含义为:在  A$中的子串  A1可以变换为可以变换为B1.A2可以变换为可 ...

  3. 2019年华南理工大学程序设计竞赛(春季赛) K Parco_Love_String(后缀自动机)找两个串的相同字串有多少

    https://ac.nowcoder.com/acm/contest/625/K 题意: 给出Q 个询问 i , 求 s[0..i-1] 与 s[i...len-1] 有多少相同的字串 分析: 给出 ...

  4. 最长公共子序列(LCS)问题 Longest Common Subsequence 与最长公告字串 longest common substr

    问题描述:字符序列的子序列是指从给定字符序列中随意地(不一定连续)去掉若干个字符(可能一个也不去掉)后所形成的字符序列.令给定的字符序列X=“x0,x1,…,xm-1”,序列Y=“y0,y1,…,yk ...

  5. 04-String——课后作业1:字串加密

    题目:请编写一个程序,加密或解密用户输入的英文字串要求设计思想.程序流程图.源代码.结果截图. 程序设计思想:首先由用户选择是加密还是解密,利用String类中的charAt函数依次取出字串中的字符, ...

  6. 最大公共字串LCS问题(阿里巴巴)

    给定两个串,均由最小字母组成.求这两个串的最大公共字串LCS(Longest Common Substring). 使用动态规划解决. #include <iostream> #inclu ...

  7. 编程:使用递归方式判断某个字串是否回文(Palindrome)

    Answer: import java.util.Scanner; public class Palindrome { private static int len;//全局变量整型数据 privat ...

  8. NOIP2002字串变换[BFS]

    题目描述 已知有两个字串 A$, B$ 及一组字串变换的规则(至多6个规则): A1$ -> B1$ A2$ -> B2$ 规则的含义为:在 A$中的子串 A1$ 可以变换为 B1$.A2 ...

  9. 【BZOJ 4517】【SDOI 2016 Round1 Day2 T2】排列计数

    本蒟蒻第一次没看题解A的题竟然是省选$Round1$ $Day2$ $T2$ 这道组合数学题. 考试时一开始以为是莫队,后来想到自己不会组合数的一些公式,便弃疗了去做第三题,,, 做完第三题后再回来看 ...

随机推荐

  1. Avro Parquet

    行   支持数据追加 列  频繁进行小部分列查询

  2. 剑指Offer:树的子结构【26】

    剑指Offer:树的子结构[26] 题目描述 输入两棵二叉树A,B,判断B是不是A的子结构.(ps:我们约定空树不是任意一个树的子结构) 解题思路 分为两步: 第一步:在树A中找到和树B的根节点的值一 ...

  3. maven资料

    1.配置settings.xml:http://www.cnblogs.com/jingmoxukong/p/6050172.html?utm_source=gold_browser_extensio ...

  4. C++的逐过程和逐语句的区别

    1.逐语句是指在遇到函数调用语句的时候进入到函数内部执行. 2.逐过程是指在遇到函数调用语句时把函数当作一条语句执行.

  5. H3C-路由器密码恢复

    路由器密码恢复: 1.先关闭电源,重新启动路由器,注意终端上显示 press CTRL+B to enter extended boot menu 的时候,我们迅速按下ctrl+B,这样将进入扩展启动 ...

  6. 并不对劲的loj3048:p5283:[十二省联考]异或粽子

    题目大意 有\(n\)(\(n\leq5\times10^5\))个数\(a_1,a_2,...a_n\)(\(a_i\leq 2^{32}-1\)) 求区间异或和前\(k(k\leq2\times1 ...

  7. BZOJ_1004_[HNOI2008]Cards_burnside+DP

    BZOJ_1004_[HNOI2008]Cards_burnside+DP Description 小春现在很清闲,面对书桌上的N张牌,他决定给每张染色,目前小春只有3种颜色:红色,蓝色,绿色.他询问 ...

  8. gsoap开发webservice

    gSOAP编译工具提供了一个SOAP/XML 关于C/C++ 语言的实现,从而让C/C++语言开发web服务或客户端程序的工作变得轻松了很多.绝大多数的C++web服务工具包提供一组API函数类库来处 ...

  9. springboot开启远程调试

    远程调试maven设置 The run goal forks a process for the boot application. It is possible to specify jvm arg ...

  10. kafka之三:kafka java 生产消费程序demo示例

    kafka是吞吐量巨大的一个消息系统,它是用scala写的,和普通的消息的生产消费还有所不同,写了个demo程序供大家参考.kafka的安装请参考官方文档. 首先我们需要新建一个maven项目,然后在 ...