码迷,mamicode.com
首页 > 其他好文 > 详细

4560 NOIP2015 D2T2 子串

时间:2016-07-09 19:19:18      阅读:185      评论:0      收藏:0      [点我收藏+]

标签:

4560 NOIP2015 D2T2 子串

 

 时间限制: 1 s
 空间限制: 128000 KB
 题目等级 : 黄金 Gold
 
 
题目描述 Description

有两个仅包含小写英文字母的字符串A和B。现在要从字符串A中取出k个互不重叠的非空子串,然后把这k个子串按照其在字符串A中出现的顺序依次连接起来得到一个新的字符串,请问有多少种方案可以使得这个新串与字符串B相等?注意:子串取出的位置不同也认为是不同的方案。

输入描述 Input Description

第一行是三个正整数n,m,k,分别表示字符串A的长度,字符串B的长度,以及问题描述中所提到的k,每两个整数之间用一个空格隔开。 

第二行包含一个长度为n的字符串,表示字符串A。 第三行包含一个长度为m的字符串,表示字符串B。

 

输出描述 Output Description

输出共一行,包含一个整数,表示所求方案数。由于答案可能很大,所以这里要求输出答案对1,000,000,007取模的结果。

样例输入 Sample Input

【Input1】

6 3 1 

aabaab 

aab

【Input2】

6 3 2 

aabaab 

aab

【Input3】

6 3 3 

aabaab 

aab

 

样例输出 Sample Output

【Output1】

2

【Output2】

7

【Output3】

7

数据范围及提示 Data Size & Hint

对于第1组数据:1≤n≤500,1≤m≤50,k=1; 

对于第2组至第3组数据:1≤n≤500,1≤m≤50,k=2; 

对于第4组至第5组数据:1≤n≤500,1≤m≤50,k=m; 

对于第1组至第7组数据:1≤n≤500,1≤m≤50,1≤k≤m; 

对于第1组至第9组数据:1≤n≤1000,1≤m≤100,1≤k≤m; 

对于所有10组数据:1≤n≤1000,1≤m≤200,1≤k≤m。

 

样例解释:

所有合法方案如下:(加下划线的部分表示取出的子串)

样例一:aab aab / aab aab

样例二:a ab aab / a aba ab / a a ba ab / aab a ab / aa b aab / aa baa b / aab aa b

样例三:a a b aab / a a baa b / a ab a a b / a aba a b / a a b a a b / a a ba a b / aab a a b

 

题解:(dp) 
f(k,i,j)表示分了k段,用了第一个串中的前i个数字,已经构成了第二个串的前j个的方案数 

f(k,i,j)=∑f(k1,l,j1)                     当s1[i]==s2[j]s1[i1]!=s2[j1]其0<l<i

           ∑f(k1,l,j1)+f(k,i1,j1)  当s1[i]==s2[j]s1[i1]==s2[j1]0<l<i

第一个转移是独立开辟出一个部分的可能数,第二个是算上之前就已经分了k个部分后把第k个部分扩大的方案数。还是比较好理解的,当然f数组要滚动第一维,不然会爆内存。 
我觉得这一道题的难度,都和当年乌龟棋在当时的难度差不多了。但这道题的70分算法O(nm2k)比较好写,所以拉不卡差距。 
相信大家在看了我的状态转移之后能想出O(nmk)的方法。其实我们需要维护的就只是那一个用一个tmp数组保存一下就好了(详见代码) 

AC代码:

#include<cstdio>
#include<iostream>
using namespace std;
#define N 1005
#define M 205
#define mod 1000000007
char a[N],b[M];
int f[2][N][M],n,m,K;
int main(){
    scanf("%d%d%d",&n,&m,&K);
    scanf("%s%s",a+1,b+1);
    for(int i=0;i<=n;i++) f[0][i][0]=1;
    int k=0;
    for(int h=1;h<=K;h++){
        k^=1;
        for(int i=0;i<=n;i++) f[k][i][h-1]=0;
        for(int j=h;j<=m;j++){
            for(int i=j;i<=n;i++){
                if(a[i]==b[j]){
                    f[k][i][j]=(f[k][i-1][j-1]+f[k][i-1][j])%mod;
                    f[k][i][j]=(f[k][i][j]+f[k^1][i-1][j-1])%mod;
                    if(i>=2) f[k][i][j]=(f[k][i][j]-f[k][i-2][j-1]+mod)%mod;
                }
                else f[k][i][j]=f[k][i-1][j];
            }
        }
    }
    printf("%d\n",f[k][n][m]);
    return 0;
}

 

4560 NOIP2015 D2T2 子串

标签:

原文地址:http://www.cnblogs.com/shenben/p/5656283.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!