CF 314B: Sereja and Periods

最新推荐文章于 2020-04-25 07:59:31 发布

frog1902

最新推荐文章于 2020-04-25 07:59:31 发布

阅读量1.7k

点赞数

分类专栏： ACM_动态规划

本文链接：https://blog.csdn.net/liguan1/article/details/9060445

版权

ACM_动态规划专栏收录该内容

25 篇文章 0 订阅

订阅专栏

博客介绍了Codeforces竞赛中的一道题目314B - Sereja and Periods，讨论如何找到最大的整数p，使得[q, p]是字符串w的子序列。文章提到了两种解法，一种是预处理匹配策略，另一种利用了倍增思想，并提供了相关参考代码链接。" 132704311,19673616,使用OpenCV在图像上添加边框的Python教程,"['计算机视觉', '图像处理', 'OpenCV', 'Python']

摘要由CSDN通过智能技术生成

题目链接：http://codeforces.com/contest/314/problem/B

题意：

对于字符串s和正整数n，规定[s,n]=s+s+……+s;

即 ["abc", 2] ="abcabc"

现给出w = [a, b] , q = [c, d]

求最大的p，使得 [q, p]是w的子序列

如无合法p，则输出0

由于CF可以看代码的缘故，一般不写CF题解。

但是觉得这道题确实比较巧妙，于是记录一下。

做这道题时参考了两位神犇的代码：

http://codeforces.com/contest/314/submission/3832474

http://codeforces.com/contest/314/submission/3845338

谢谢！

算法：

解法一：

预处理若从c串的第i位开始匹配，

那么在a串中最多可以循环匹配c串的多少位。

解法二：

借用倍增的思想。

go[dep][get_id(i,j)]表示从字符串a的第i位，字符串c的第j位开始，匹配字符串c的(2^dep)位后

.first表示两个字符串匹配后到达的位置

.second表示“消耗”了多少个a串

注意，很显然在存在合法解的情况下

消耗a串的数量和c串的数量都不会超过int，

但是在不存在合法解的情况下，

.second初始化为INF，经过多次倍增可能会超过int

代码如下：

解法一：

#include<cstdio>
#include<iostream>
#include<sstream>
#include<cstdlib>
#include<cstring>
#include<string>
#include<climits>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
#include<stack>
#include<set>
#include<map>
#define INF 0x3f3f3f3f
#define eps 1e-8
using namespace std;

const int MAXN=110;
char a[MAXN],c[MAXN];
int nxt[2][MAXN];

int main() {
    int b,d;
    scanf("%d%d",&b,&d);
    scanf("%s%s",a,c);
    int len1=strlen(a);
    int len2=strlen(c);
    for(int j=0; j<len2; j++) {
        int cur=j;
        nxt[0][j]=0;
        for(int i=0; i<len1; i++) {
            if(c[cur]==a[i]) {
                cur++;
                if(cur==len2) {
                    nxt[0][j]++;
                    cur=0;
                }
            }
        }
        nxt[1][j]=cur;
    }
    long long ans=0LL;
    int cur=0;
    for(int i=0; i<b; i++) {
        ans+=nxt[0][cur];
        cur=nxt[1][cur];
    }
    printf("%I64d\n",ans/d);
    return 0;
}

解法二：

#include<cstdio>
#include<iostream>
#include<sstream>
#include<cstdlib>
#include<cstring>
#include<string>
#include<climits>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
#include<stack>
#include<set>
#include<map>
#define INF 0x3f3f3f3f
#define eps 1e-8
using namespace std;

const int MAXN=110;
const int MAXM=30;
char a[MAXN],c[MAXN];
int len1,len2;
pair<int,long long> go[MAXM][MAXN*MAXN];

int get_id(int i, int j) {
    return i*len2+j;
}

long long solve(int dis) {
    int cur=get_id(0,0);
    long long cot=0;
    for (int i=0; i<MAXM; i++) {
        if(dis>>i&1) {
            cot+=go[i][cur].second;
            cur=go[i][cur].first;
        }
    }
    return cot;
}

int main() {
    int b,d;
    scanf("%d%d",&b,&d);
    scanf("%s%s",a,c);
    len1=strlen(a);
    len2=strlen(c);
    for(int i=0; i<len1; i++) {
        for(int j=0; j<len2; j++) {
            int id=get_id(i,j);
            go[0][id]=make_pair(id,INF);
            for(int k=0; k<len1; k++) {
                if(a[(i+k)%len1]==c[j]) {
                    go[0][id]=make_pair((get_id((i+k+1)%len1,(j+1)%len2)),(i==0)||((i+k)/len1));
                    break;
                }
            }
        }
    }
    for(int i=1; i<MAXM; i++) {
        for(int j=0; j<len1*len2; j++) {
            int nxt=go[i-1][j].first;
            go[i][j]=make_pair(go[i-1][nxt].first, go[i-1][j].second+go[i-1][nxt].second);
        }
    }
    int l=0;
    int r=INF;
    while(l<r) {
        int mid=(l+r+1)>>1;
        if(solve(mid)<=b) {
            l=mid;
        } else {
            r=mid-1;
        }
    }
    printf("%d\n",l/(len2*d));
    return 0;
}