[BZOJ 3145][Feyat cup 1.5]Str 解题报告

最新推荐文章于 2020-07-05 12:08:12 发布
最新推荐文章于 2020-07-05 12:08:12 发布
阅读量162 收藏
点赞数
原文链接：http://www.cnblogs.com/Candyouth/p/5596918.html
版权

[Feyat cup 1.5]Str

Description
Arcueid，白姬，真祖的公主。在和推倒贵看电影时突然对一个问题产生了兴趣：
我们都知道真祖和死徒是有类似的地方。那么从现代科学的角度如何解释
呢？自然就得研究遗传密码了。Arcueid得知了两者的DNA片段，想寻求一个
DNA片段，使得其在两者的DNA中都出现过。
我们知道公主的脑袋有点不太灵活，如果两个DNA片段只有一个位置不
同，她也会将其认为是相同的。所以请您找出这样的最长的DNA片段吧。
Input
两行，每行一个字符串。
Output
一个整数，表示最长的 DNA 片段的长度。
Sample Input
aabbe
acbbc
Sample Output
4
HINT
100% 的数据 n<=10^5;m<=10^5 。

 

 

 

Sol：

2015年集训队论文里有讲到这个题目

Sam +　Sa

将a，b拼接起来，中间加一个特殊字符

考虑后缀自动机上的一个节点，对应A，B串的Right集合，已知这些点的最长公共后缀为这个节点的len，暴力的做法可以枚举两个集合的点对(a,b)，对于a+2,b+2求lcp更新

但是复杂度很高，考虑将A，B集合中(a+2)(b+2)这些点按照后缀排序，要求的只有相邻的所属的集合不同的两个后缀的lcp来更新答案

由于空间的限制，要在parent树自底向上维护set，启发式合并维护信息

要注意的地方是跨越两个串的地方不能直接丢进set里，会影响rank的比较以及ans的更新

所以记录一些奇奇怪怪的东西。。这个节点在A，B串后缀出现的位置的最大，最小值(语死早说不清楚)

 

代码异常丑陋，还是不要看的好。。

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <map>
#include <set>
  
#define maxn 500010
  
using namespace std;
  
int mark[maxn];
struct Node{int len, link; map<int, int> nxt;}st[maxn];
int root, size, last;
void init(){
    root = size = last = 0;
    st[root].len = 0;
    st[root].link = -1;
}
void Extend(int c){
    int p = last, cur = ++ size;
    st[cur].len = st[p].len + 1;
    for(; ~p && st[p].nxt[c] == 0; p = st[p].link)
        st[p].nxt[c] = cur;
    if(p == -1)
        st[cur].link = root;
    else{
        int q = st[p].nxt[c];
        if(st[q].len == st[p].len + 1)
            st[cur].link = q;
        else{
            int clone = ++ size;
            st[clone] = st[q];
            st[clone].len = st[p].len + 1;
            for(; ~p && st[p].nxt[c] == q; p = st[p].link)
                st[p].nxt[c] = clone;
            st[q].link = st[cur].link = clone;
        }
    }
    last = cur;
}
  
int n, m, N, checker;
char a[maxn], b[maxn];
int str[maxn];
  
int sa[maxn], t1[maxn], t2[maxn], c[maxn], rk[maxn], ht[maxn], ST[maxn][20], lg[maxn], stsize;
int lcp(int x, int y){
    if(x > N || y > N)return 0;
    x = rk[x], y = rk[y];
    if(x > y)swap(x, y);
    y --; int k = lg[y-x+1];
    return min(ST[x][k], ST[y-(1<<k)+1][k]);
}
  
/*
void Get_sa(int n, int m) {
    int *x = wa, *y = wb, *t, i, j, p;
    for(i = 0 ; i < m ; ++ i) ws[i] = 0;
    for(i = 0 ; i < n ; ++ i) ++ ws[x[i] = w[i]];
    for(i = 1 ; i < m ; ++ i) ws[i] += ws[i-1];
    for(i = n-1 ; i >= 0 ; -- i) sa[-- ws[x[i]]] = i;
    for(j = 1, p = 1 ; p < n ; j <<= 1, m = p) {
    //  printf("j = %d, p = %d \n", j, p);
        for(i = n-j, p = -1 ; i < n ; ++ i) y[++ p] = i;//puts("NO");
        for(i = 0 ; i < n ; ++ i) if(sa[i]>=j) y[++ p] = sa[i]-j;
        for(i = 0 ; i < n ; ++ i) wv[i] = x[y[i]];
        for(i = 0 ; i < m ; ++ i) ws[i] = 0;
        for(i = 0 ; i < n ; ++ i) ++ ws[wv[i]];
        for(i = 1 ; i < m ; ++ i) ws[i] += ws[i-1];
        for(i = n-1 ; i >= 0 ; -- i) sa[-- ws[wv[i]]] = y[i];
    //  puts("Oh NO");
        for(t = x, x = y, y = t, i = 1, x[sa[0]] = 0, p = 1 ; i < n ; ++ i) {
            x[sa[i]] = cmp(t, sa[i-1], sa[i], j)?p-1:p ++;
        }
    }
    return;
}
*/
  
void getsa(int n, int m){
    int *x = t1, *y = t2;
    for(int i = 0; i < m; i ++)c[i] = 0;
    for(int i = 0; i < n; i ++)c[x[i] = str[i]] ++;
    for(int i = 1; i < m; i ++)c[i] += c[i-1];
    for(int i = n-1; ~ i; i --)sa[-- c[x[i]]] = i;
  
    for(int k = 1; k <= n; k <<= 1){
        int p = 0;
        for(int i = n-k; i < n; i ++)y[p ++] = i;
        for(int i = 0; i < n; i ++)if(sa[i] >= k)y[p ++] = sa[i] - k;
          
        for(int i = 0; i < m; i ++)c[i] = 0;
        for(int i = 0; i < n; i ++)c[x[y[i]]] ++;
        for(int i = 1; i < m; i ++)c[i] += c[i-1];
        for(int i = n-1; ~ i; i --)sa[-- c[x[y[i]]]] = y[i];
          
        swap(x, y); x[sa[0]] = 0, p = 1;
        for(int i = 1; i < n; i ++)
            x[sa[i]] = y[sa[i]] == y[sa[i-1]] && y[sa[i]+k] == y[sa[i-1]+k] ? p-1 : p++;
        if(p >= n)break;
        m = p; 
    }
      
    int k = 0;
    for(int i = 0; i < n; i ++)rk[sa[i]] = i;
  
    for(int i = 0; i < n; i ++){
        if(rk[i] == 0){ht[0] = 0; continue;}
        if(k) k --; int j = sa[rk[i]-1];
        while(str[i+k] == str[j+k]) k ++;
        ht[rk[i]] = k;
    }
    for(int i = n; i; i --)sa[i] = sa[i-1] + 1;
    for(int i = 1; i <= n; i ++)rk[sa[i]] = i;
    for(int i = 1; i < n; i ++)ST[i][0] = ht[i];
    stsize = n - 1, lg[0] = -1;
    for(int i = 1; i < n; i ++)lg[i] = lg[i>>1] + 1;
    for(int j = 1; 1<<j <= stsize; j ++)
        for(int i = 1; i+(1<<j)-1 <= stsize; i ++)
            ST[i][j] = min(ST[i][j-1], ST[i+(1<<j-1)][j-1]);
}
  
int h[maxn], cnt, ans, Nw, fg;
struct Edge{int to, nxt;} edge[maxn];
void addedge(int u, int v){
    edge[++ cnt] = (Edge){v, h[u]}; h[u] = cnt;
}
  
struct cmp{
    bool operator ()(const int& x, const int& y){
        return rk[x] < rk[y];
    }
};
  
typedef set<int, cmp> se;
typedef set<int, cmp>::iterator iter;
vector<int> V[maxn];
se s[maxn];
  
inline void cmax(int x, int y){
    if(x > y)swap(x, y);
    if(y <= n || x > n)return;
    if(x <= n ^ y <= n)
        ans = max(ans, Nw + 1 + lcp(x, y));
}
  
bool ha[maxn], hb[maxn], fg1, FG1, fg2, FG2, FG3, FG4;
  
int mx1[maxn], mx2[maxn];
  
void merge(se& x, se& y, int u, int v){
    if(x.size() < y.size())swap(x, y), swap(V[u], V[v]);
    for(iter it = y.begin(); it != y.end(); it ++){
        iter p = x.upper_bound(*it);
        if(p != x.end())cmax(*p, *it);
        if(p != x.begin())p --, cmax(*p, *it);
    }
      
    FG1 = fg1 = ha[v], FG2 = fg2 = hb[v], FG3 = mx1[v] - Nw >= 1, FG4 = mx2[v] - Nw >= n + 2;
    for(int i = 0; i < V[v].size(); i ++){
        int nw = V[v][i];
        fg1 |= nw <= n, FG1 |= nw < n;
        fg2 |= nw > n && nw <= checker, FG2 |= nw > n && nw < checker;
        if(nw <= n)FG3 |= nw - Nw >= 1; else FG4 |= nw - Nw >= n + 2;
    }
      
    for(int i = 0; i < V[u].size(); i ++){
        int w = V[u][i];
        if(w <= n && fg2){
            ans = max(ans, Nw + (FG2 && ((w == n - 1 || w == n + m))));
            ans = max(ans, Nw + (FG4 && w - Nw != 0 && w - Nw != n + 1));
        }
        if(w > n && fg1){
            ans = max(ans, Nw + (FG1 && (w == n - 1 || w == n + m)));
            ans = max(ans, Nw + (FG3 && w - Nw != 0 && w - Nw != n + 1));
        }
    }
      
    FG1 = fg1 = ha[u], FG2 = fg2 = hb[u], FG3 = mx1[u] - Nw >= 1, FG4 = mx2[u] - Nw >= n + 2;
    for(int i = 0; i < V[u].size(); i ++){
        int nw = V[u][i];
        fg1 |= nw <= n, FG1 |= nw < n;
        fg2 |= nw > n && nw <= checker, FG2 |= nw > n && nw < checker;
        if(nw <= n)FG3 |= nw - Nw >= 1; else FG4 |= nw - Nw >= n + 2;
    }
  
    for(int i = 0; i < V[v].size(); i ++){
        int w = V[v][i];
        if(w <= n && fg2){
            ans = max(ans, Nw + (FG2 && ((w == n - 1 || w == n + m))));
            ans = max(ans, Nw + (FG4 && w - Nw != 0 && w - Nw != n + 1));
        }
        if(w > n && fg1){
            ans = max(ans, Nw + (FG1 && (w == n - 1 || w == n + m)));
            ans = max(ans, Nw + (FG3 && w - Nw != 0 && w - Nw != n + 1));
        }
    }
      
    for(int i = 0; i < V[v].size(); i ++)
        V[u].push_back(V[v][i]);
      
    for(iter it = y.begin(); it != y.end(); it ++)
        x.insert(*it);
    V[v].clear(), y.clear();
}
  
void dfs(int u){
    if(~mark[u]){
        if((mark[u] > n+1 && mark[u]+2 <= checker) || mark[u]+2 <= n){
            s[u].insert(mark[u]+2);
            ha[u] |= mark[u] <= n;
            hb[u] |= mark[u] > n;
        }
        else V[u].push_back(mark[u]);
        if(mark[u] <= n)mx1[u] = max(mx1[u], mark[u]);
        else mx2[u] = max(mx2[u], mark[u]);
    }
    for(int i = h[u]; i; i = edge[i].nxt){
        int v = edge[i].to;
        dfs(v);
        Nw = st[u].len;
        merge(s[u], s[v], u, v);
        ha[u] |= ha[v];
        hb[u] |= hb[v];
        mx1[u] = max(mx1[u], mx1[v]);
        mx2[u] = max(mx2[u], mx2[v]);
    }
}
  
int nd[10];
  
int main(){
    init();ans = 1;
    scanf("%s%s", a + 1, b + 1);
    n = strlen(a + 1), m = strlen(b + 1);
      
    for(int i = 1; i <= n; i ++)Extend(a[i] - 'a'); Extend(27);
    for(int i = 1; i <= m; i ++)Extend(b[i] - 'a');
      
    for(int i = 1; i <= size; i ++)addedge(st[i].link, i);
    memset(mark, -1, sizeof mark); int cur = root; mark[root] = 0;
    for(int i = 1; i <= n; i ++)cur = st[cur].nxt[a[i] - 'a'], mark[cur] = i;
    cur = st[cur].nxt[27];
    for(int i = 1; i <= m; i ++)cur = st[cur].nxt[b[i] - 'a'], mark[cur] = i+n+1;
      
    N = 0, checker = n + m + 1;
    for(int i = 1; i <= n; i ++)str[N ++] = a[i] - 'a' + 1; str[N ++] = 27;
    for(int i = 1; i <= m; i ++)str[N ++] = b[i] - 'a' + 1;
    getsa(n+m+2, 30);
      
    dfs(0);
    printf("%d\n", ans);
    return 0;
}

　　

转载于:https://www.cnblogs.com/Candyouth/p/5596918.html

dcf56843
关注
【动态规划】B君的多边形
qq_42397914的博客
10-29 252
描述 有一个正n多边形，我们要连接一些对角线，把这个多边形分成若干 个区域。要求连接的对角线不能相交，每个点可以连出也可以不连出对角 线。（即最终不要求所有区域均为三角形） 问总方案数mod (10^9+7)的结果。 输入 一行一个整数n 输出 一行一个整数表示答案。 样例输入 5 样例输出 11 提示 对于所有的的数据，满足n≤10^6 代码 这是一道找规律的题，没研究明白，先记着吧~ #include<iostream> using namespace std; long long inv
3145: [Feyat cup 1.5]Str
CRZbulabula的博客
01-29 703
3145: [Feyat cup 1.5]Str Time Limit: 40 Sec  Memory Limit: 512 MB Submit: 142  Solved: 25 [Submit][Status][Discuss] Description Arcueid，白姬，真祖的公主。在和推倒贵看电影时突然对一个问题产生了兴趣： 我们都知道真祖和死徒是有类似的地方。那么从现代科
OpenJudge超详细题解，动画图文题解
热门推荐
07-05 3万+
内容来源于微信公众号：大神编程。已经过原文作者授权。 更新时间：2020-7-5 因为OpenJudge题库和一本通题库很多题都是一样，所以很多链接都是直接使用一本通的题解，完全不影响。 为什么要转载： 动画、图文结合。 我看了部分文章，一些我觉得很简单的题，但从中学到了以前没有学过的知识点和细节，重点是从中学习别人思维方式。 让更多的同学能学到更多的知识。 官方的QQ群：893157498 这是我建的QQ群：795233394 欢迎各位志同道合的同学们^_^ 已更新的题目颜色是蓝色，点击相
BZOJ-4706 B君的多边形 OEIS
weixin_34362875的博客
09-20 137
题面 题意:有一个正n多边形,我们要连接一些对角线,把这个多边形分成若干个区域,要求连接的对角线不能相交,每个点可以连出也可以不连出对角线,即最终不要求所有区域均为三角形,问总方案数mod (10^9+7)的结果。n≤10^6 题解:手推1,1,3,11,45,利用OEIS知道这是 超级卡特兰数（又称大施罗德数） (n+1) * f(n) = (6*n-3) * f(n-1) - (n-2)...
超级卡特兰数（bzoj 4706: B君的多边形）
加载中...
07-16 1292
4706: B君的多边形 Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 159  Solved: 92 [Submit][Status][Discuss] Description 有一个正n多边形，我们要连接一些对角线，把这个多边形分成若干个区域。要求连接的对角线不能相交，每个点可以连出也可以不连出对角线。（即最终不要求所有区域均为
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
BZOJ题目镜像
09-14
【BZOJ题目镜像】是为了解决BZOJ（Blue Zoo Online Judge）系统因服务器问题时常出现访问不稳定或无法访问的情况而设立的一个备份资源。BZOJ是中国早期的在线编程竞赛平台，提供了大量的算法题目供编程爱好者进行...
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
bzoj problem
08-21
【标题】"bzoj problem" 指的是“八中OJ（Online ...总之，"bzoj problem" 提供的资源对于想要提升编程技能和算法理解的用户来说是宝贵的资料，通过解题实践，不仅能巩固理论知识，还能提高实际编程解决问题的能力。
BZOJ3230相似子串测试数据
08-19
标题"BZOJ3230相似子串测试数据"指的是一个与在线判题平台BZOJ（BestCoder ZOJ，中国大学算法竞赛在线评测系统）相关的编程问题，具体问题编号为3230，涉及的主要知识点是"相似子串"的搜索和处理。在编程竞赛中，...
bzoj4706: B君的多边形
mengbi_er的博客
07-24 724
递推式： f[n]=(f[n−1]∗(6n−3)−f[n−2]∗(n−2))/(n+1) f[n]=(f[n-1]*(6n-3)-f[n-2]*(n-2))/(n+1) 预处理逆元即可。#include<iostream> using namespace std; int mod=1000000007; long long f[1000005],inv[1000005]; int main()
NOI题库 B君的多边形
水墨苍穹的博客
10-24 1362
9282:B君的多边形 总时间限制: 20000ms 单个测试点时间限制: 1000ms 内存限制: 131072kB 描述 有一个正n多边形，我们要连接一些对角线，把这个多边形分成若干  个区域。要求连接的对角线不能相交，每个点可以连出也可以不连出对角  线。（即最终不要求所有区域均为三角形）  问总方案数mod (10^9+7)的结果。 
BZOJ 4706: B君的多边形 找规律
weixin_33889665的博客
10-03 196
4706: B君的多边形 题目连接： http://www.lydsy.com/JudgeOnline/problem.php?id=4706 Description 有一个正n多边形，我们要连接一些对角线，把这个多边形分成若干个区域。要求连接的对角线不能相交，每个点 可以连出也可以不连出对角线。（即最终不要求所有区域均为三角形）问总方案数mod (10^9+7)的结果。 Input 一行一个整数...
[2018.07.14 T1] B君的第四题
ShadyPi的IT笔记
07-14 431
暂无链接 B君的第四题 【问题描述】 已有的事后必再有，已行的事后必再行。 大家还记得Luogo P3959 [NOIP2017]宝藏 吗？ B君看到很多人用很多弱智的方法水过这个题，感到非常气愤。 这个题目就是把原题改成了一个计数题。 考虑到有些同学没有看过这个题目，我们简述这个题目如下。 输入一个图，我们定义一个有根生成树的权值是 ∑e={x,y}we×m...
BZOJ4706 B君的多边形 (超级卡特兰数/施罗德数)
Ark
10-30 298
题目 权限题QAQ 题解 超级卡特兰数/施罗德数 CODE #include <bits/stdc++.h> using namespace std; const int MAXN = 1000005; const int mod = 1e9 + 7; int inv[MAXN], f[MAXN], n; int main () { scanf("%d", &n); ...
B君的教育【省选模拟赛】
Welcome to yjjr's blog!
03-12 268
标签：数学 题目 分析 类似二进制拆分，短除法 直接看代码就好 复杂度O(log n) code #include&lt;iostream&gt; #include&lt;cstdio&gt; #include&lt;cstdlib&gt; #include&lt;cmath&gt; #include&lt;cstring&gt; #include&lt;algori...
日常训练 20170602 B君的教育
MintGreen
06-02 534
题意简述：有一个性质是设 p=−1+ip=−1+i ，对于每个高斯整数 x+yix+yi，我们都可以找到一个非负整数集合 SS 满足 ∑i∈Spi=n\sum_{i\in S} p^i = n，在输入 xx 和 yy 的情况下，输出 SS 。 对于一个复数（Complex Number），如果它的实部和虚部都是整数，我们称之为高斯整数。 很多大佬们都想到了可以做类似二进制拆分，就用了短除法，复杂
有序分数求和
sdzh2012的专栏
11-03 815
/* * Copyright (c) 2012, 烟台大学计算机学院 * All rights reserved. * 作 者：张浩 * 完成日期：2012 年11月3 日 * 版 本 号：v1.0 * 输入描述： * 问题描述：求1-1/3+3/5-5/7...19/21的和 * 程序输出:要用到double型 * 问题分析：找出每一项的规律来求和 * 算法
BZOJ4707 : B君的技巧
weixin_34112181的博客
09-19 101
建立线段树，设$f[x][l][r]$表示当前考虑$x$点，最左端是$l$，最右端是$r$的最少代价。 如果$a$在$x&lt;&lt;1$，$d$在$x&lt;&lt;1|1$， 设$g[a][c]=\min(f[x&lt;&lt;1][a][b]+w[b][c])$， 则$f[x][a][d]=\min(g[a][c]+f[x&lt;&lt;1|1][c][d])$。 对于$a$在$x...
bzoj做题计划 excel
最新发布
07-16
bzoj网站上有很多高质量的题解和优秀的解题代码，可以参考和学习。还有相关的讨论区，可以与其他程序员交流和讨论，共同进步。 最后，要坚持并保持思考。做题不是单纯为了刷数量，更重要的是学会思考和总结。遇到...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值