spoj LCS

最新推荐文章于 2021-01-28 05:07:33 发布
最新推荐文章于 2021-01-28 05:07:33 发布
阅读量105 收藏
点赞数
原文链接：http://www.cnblogs.com/gaudar/p/10805234.html
版权

解法1：后缀自动机

解法1.1用串1建后缀自动机，用串2匹配。并记录能匹配的长度。匹配成功+1，否则变为maxlen[最后fail到的位置]+1。

解法1.2串一先加，然后加入'#‘。再加串2.加串2的时候，如果新得到的节点的slink在一串，则说明出现了2次。

解法1.3.建立广义后缀自动机，用数组记录各节点出现在哪个串，在拓扑排序。然后出现2次的maxlen最大值就是。

解法2:2分+HASH。会超时。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#include <iomanip>
#include <cstring>
#include <map>
#include <queue>
#include <set>
#include <cassert>
#include <stack>
#include <bitset>
//#include <unordered_set>
#define mkp make_pair
#define err cout<<"here"<<endl
using namespace std;
const double EPS=1e-12;
typedef long long lon;
typedef unsigned long long ull;
typedef map<ull,int>::iterator IT;
const lon SZ=500010,SSZ=120,APB=26,mod=97,one=1;
const lon INF=0x7FFFFFFF;
int n,cnt,nex[SZ][APB],fail[SZ];
int maxlen[SZ],minlen[SZ],trans[SZ][APB],slink[SZ];
struct nd{
    int to,wt;
    nd(int a=0,int b=0):to(a),wt(b){}
};
char ch1[SZ],ch2[SZ];

int add(int pre,int cur)
{
    int z=++cnt;
    int u=pre;
    for(;u!=-1&&!trans[u][cur];u=slink[u])
    {
        trans[u][cur]=z;
    }
    if(u==-1)
    {
        maxlen[z]=maxlen[pre]+1;
        minlen[z]=1;
        slink[z]=0;
    }
    else
    {
        int x=trans[u][cur];
        if(maxlen[x]==maxlen[u]+1)
        {
            maxlen[z]=maxlen[pre]+1;
            minlen[z]=maxlen[x]+1;
            slink[z]=x;
        }
        else
        {
            int y=++cnt;
            memcpy(trans[y],trans[x],sizeof(trans[x]));
            maxlen[y]=maxlen[u]+1;
            minlen[y]=maxlen[slink[x]]+1;
            slink[y]=slink[x];
            maxlen[x]=maxlen[x];
            minlen[x]=maxlen[y]+1;
            slink[x]=y;
            maxlen[z]=maxlen[pre]+1;
            minlen[z]=maxlen[y]+1;
            slink[z]=y;
            for(;u!=-1&&trans[u][cur]==x;u=slink[u])
            {
                trans[u][cur]=y;
            }
        }
    }
    return z;
}

void init()
{
    cin>>ch1+1>>ch2+1;
    int pre=0;
    slink[0]=-1;
    for(int i=1;ch1[i];++i)
    {
        pre=add(pre,ch1[i]-'a');
    }
}

int get_nex(int sta,int c,int &len)
{
    int u=sta;
    if(trans[u][c])
    {++len;
        return trans[u][c];
    }
    for(;u!=-1&&!trans[u][c];)
    {
        u=slink[u];
        //if(u!=-1)len=min(len,maxlen[u]);
    }
    if(u==-1)u=0,len=0;
    else len=maxlen[u]+1,u=trans[u][c];
    return u;
}

void work()
{
    int cur=0,res=0,len=0;
    for(int i=1;ch2[i];++i)
    {
        int c=ch2[i]-'a';
        cur=get_nex(cur,c,len);
        res=max(res,len);
        //cout<<len<<endl;
    }
    cout<<res<<endl;
}

void release()
{
    
}

int main()
{
    std::ios::sync_with_stdio(0);
    //freopen("d:\\1.txt","r",stdin);
    int casenum;
    //cin>>casenum;
    //cout<<casenum<<endl;
    //for(int time=1;time<=casenum;++time)
    //for(int time=1;cin>>ch+1,ch[1]!='e';++time)
    {
        init();
        work();
        release();
    }
    return 0;
}
#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#include <iomanip>
#include <cstring>
#include <map>
#include <queue>
#include <set>
#include <cassert>
#include <stack>
#include <bitset>
//#include <unordered_set>
#define mkp make_pair
#define err cout<<"here"<<endl
using namespace std;
const double EPS=1e-12;
typedef long long lon;
typedef unsigned long long ull;
typedef map<ull,int>::iterator IT;
const lon SZ=1000010,SSZ=120,APB=27,mod=97,one=1;
const lon INF=0x7FFFFFFF;
int n,cnt,num,nex[SZ][APB],fail[SZ],ans;
int maxlen[SZ],minlen[SZ],trans[SZ][APB],slink[SZ];
struct nd{
    int to,wt;
    nd(int a=0,int b=0):to(a),wt(b){}
};
char ch1[SZ],ch2[SZ];
bool ok;

int add(int pre,int cur)
{
    int z=++cnt;
    int u=pre;
    for(;u!=-1&&!trans[u][cur];u=slink[u])
    {
        trans[u][cur]=z;
    }
    if(u==-1)
    {
        maxlen[z]=maxlen[pre]+1;
        minlen[z]=1;
        slink[z]=0;
    }
    else
    {
        int x=trans[u][cur];
        if(maxlen[x]==maxlen[u]+1)
        {
            maxlen[z]=maxlen[pre]+1;
            minlen[z]=maxlen[x]+1;
            slink[z]=x;
            if(ok&&x<=num)ans=max(ans,maxlen[x]);
        }
        else
        {
            int y=++cnt;
            memcpy(trans[y],trans[x],sizeof(trans[x]));
            maxlen[y]=maxlen[u]+1;
            minlen[y]=maxlen[slink[x]]+1;
            slink[y]=slink[x];
            maxlen[x]=maxlen[x];
            minlen[x]=maxlen[y]+1;
            slink[x]=y;
            maxlen[z]=maxlen[pre]+1;
            minlen[z]=maxlen[y]+1;
            slink[z]=y;
            for(;u!=-1&&trans[u][cur]==x;u=slink[u])
            {
                trans[u][cur]=y;
            }
            if(ok&&x<=num)ans=max(ans,maxlen[y]);
        }
    }
    return z;
}

void init()
{
    cin>>ch1+1>>ch2+1;
    int pre=0;
    slink[0]=-1;
    for(int i=1;ch1[i];++i)
    {
        pre=add(pre,ch1[i]-'a');
    }
    pre=add(pre,26);
    ok=1;
    num=cnt;
    for(int i=1;ch2[i];++i)
    {
        pre=add(pre,ch2[i]-'a');
    }
    cout<<ans<<endl;
}

void work()
{
    
}

void release()
{
    
}

int main()
{
    std::ios::sync_with_stdio(0);
    //freopen("d:\\1.txt","r",stdin);
    int casenum;
    //cin>>casenum;
    //cout<<casenum<<endl;
    //for(int time=1;time<=casenum;++time)
    //for(int time=1;cin>>ch+1,ch[1]!='e';++time)
    {
        init();
        work();
        release();
    }
    return 0;
}
#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#include <iomanip>
#include <cstring>
#include <map>
#include <queue>
#include <set>
#include <cassert>
#include <stack>
#include <bitset>
//#include <unordered_set>
#define mkp make_pair
#define err cout<<"here"<<endl
using namespace std;
const double EPS=1e-12;
typedef long long lon;
typedef unsigned long long ull;
typedef map<ull,int>::iterator IT;
const lon SZ=1000010,SSZ=120,APB=27,mod=97,one=1;
const lon INF=0x7FFFFFFF;
int n,cnt,nex[SZ][APB],fail[SZ];
int maxlen[SZ],minlen[SZ],trans[SZ][APB],slink[SZ];
struct nd{
    int to,wt;
    nd(int a=0,int b=0):to(a),wt(b){}
};
char ch1[SZ],ch2[SZ];
int sta[SZ],in[SZ];

int add(int pre,int cur,int num)
{
    if(trans[pre][cur])
    {
        int u=trans[pre][cur];
        if(maxlen[u]==maxlen[pre]+1)
        {
            sta[u]|=1<<num;
            return u;
        }
        else
        {
            int z=++cnt;
            memcpy(trans[z],trans[u],sizeof(trans[u]));
            maxlen[z]=maxlen[pre]+1;
            slink[z]=slink[u];
            slink[u]=z;
            minlen[z]=maxlen[slink[z]]+1;
            minlen[u]=maxlen[z]+1;
            sta[z]|=1<<num;
            return z;
        }
    }
    int z=++cnt;
    int u=pre;
    for(;u!=-1&&!trans[u][cur];u=slink[u])
    {
        trans[u][cur]=z;
    }
    if(u==-1)
    {
        maxlen[z]=maxlen[pre]+1;
        minlen[z]=1;
        slink[z]=0;
    }
    else
    {
        int x=trans[u][cur];
        if(maxlen[x]==maxlen[u]+1)
        {
            maxlen[z]=maxlen[pre]+1;
            minlen[z]=maxlen[x]+1;
            slink[z]=x;
        }
        else
        {
            int y=++cnt;
            memcpy(trans[y],trans[x],sizeof(trans[x]));
            sta[y]=sta[x];
            maxlen[y]=maxlen[u]+1;
            minlen[y]=maxlen[slink[x]]+1;
            slink[y]=slink[x];
            maxlen[x]=maxlen[x];
            minlen[x]=maxlen[y]+1;
            slink[x]=y;
            maxlen[z]=maxlen[pre]+1;
            minlen[z]=maxlen[y]+1;
            slink[z]=y;
            for(;u!=-1&&trans[u][cur]==x;u=slink[u])
            {
                trans[u][cur]=y;
            }
        }
    }
    sta[z]|=1<<num;
    return z;
}

void topo()
{
    for(int i=1;i<=cnt;++i)
    {
        ++in[slink[i]];
    }
    queue<int> q;
    for(int i=1;i<=cnt;++i)
    {
        if(in[i]==0)
        {
            q.push(i);
        }
    }
    for(;q.size();)
    {
        int fr=q.front();
        q.pop();
        int x=slink[fr];
        --in[x];
        sta[x]|=sta[fr];
        if(in[x]==0&&x)
        {
            q.push(x);
        }
    }
}

void init()
{
    cin>>ch1+1>>ch2+1;
    int pre=0;
    slink[0]=-1;
    for(int i=1;ch1[i];++i)
    {
        pre=add(pre,ch1[i]-'a',0);
    }
    pre=0;
    for(int i=1;ch2[i];++i)
    {
        pre=add(pre,ch2[i]-'a',1);
    }
    topo();
}

void work()
{
    int res=0;
    for(int i=1;i<=cnt;++i)
    {
        if(sta[i]==3)res=max(res,maxlen[i]);
    }
    cout<<res<<endl;
}

void release()
{
    
}

int main()
{
    std::ios::sync_with_stdio(0);
    //freopen("d:\\1.txt","r",stdin);
    int casenum;
    //cin>>casenum;
    //cout<<casenum<<endl;
    //for(int time=1;time<=casenum;++time)
    //for(int time=1;cin>>ch+1,ch[1]!='e';++time)
    {
        init();
        work();
        release();
    }
    return 0;
}
#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#include <iomanip>
#include <cstring>
#include <map>
#include <queue>
#include <set>
#include <cassert>
#include <stack>
#include <bitset>
//#include <unordered_set>
#define mkp make_pair
#define err cout<<"here"<<endl
using namespace std;
const double EPS=1e-12;
typedef long long lon;
typedef unsigned long long ull;
typedef map<ull,int>::iterator IT;
const lon SZ=250010,SSZ=120,APB=26,mod=97,one=1;
const lon INF=0x7FFFFFFF;
int n;
struct nd{
    int to,wt;
    nd(int a=0,int b=0):to(a),wt(b){}
};
char ch1[SZ],ch2[SZ];
ull arr[SZ];

bool chk(int x)
{
    vector<ull> st;
    ull cur=0;
    for(int i=1;ch1[i];++i)
    {
        if(i<=x)
        {
            cur=cur*mod+ch1[i];
        }
        else
        {
            cur=cur-arr[x-1]*ch1[i-x];
            cur=cur*mod+ch1[i];
        }
        if(i>=x)st.push_back(cur);
    }
    cur=0;
    sort(st.begin(),st.end());
    for(int i=1;ch2[i];++i)
    {
        if(i<=x)
        {
            cur=cur*mod+ch2[i];
        }
        else
        {
            cur=cur-arr[x-1]*ch2[i-x];
            cur=cur*mod+ch2[i];
        }
        if(i>=x&&*lower_bound(st.begin(),st.end(),cur)==cur)return 1;
    }
    return 0;
}

void init()
{
    //cin>>ch1+1>>ch2+1;
    scanf(" %s %s",ch1+1,ch2+1);
    int len=strlen(ch1+1);
    arr[1]=mod;
    for(int i=2;i<len+5;++i)
    {
        if(i&1)arr[i]=arr[i/2]*arr[i/2]*mod;
        else arr[i]=arr[i/2]*arr[i/2];
    }
    int lo=1,hi=strlen(ch1+1)+1;
    for(;lo<hi;)
    {
        int mid=(lo+hi)/2;
        bool ok=chk(mid);
        if(ok)lo=mid+1;
        else hi=mid;
    }
    cout<<lo-1<<endl;
}

void work()
{
    
}

void release()
{
    
}

int main()
{
    //std::ios::sync_with_stdio(0);
    //freopen("d:\\1.txt","r",stdin);
    int casenum;
    
    //cin>>casenum;
    //cout<<casenum<<endl;
    //for(int time=1;time<=casenum;++time)
    //for(int time=1;cin>>ch+1,ch[1]!='e';++time)
    {
        init();
        work();
        release();
    }
    return 0;
}

 

转载于:https://www.cnblogs.com/gaudar/p/10805234.html

weixin_30505485
关注
spoj 1811 LCS (后缀自动机 SAM)
ddyyxx的程序员之路
03-25 7824
http://www.spoj.com/problems/LCS/ 题意：求两个字符串A,B的最长公共子串。字符串长度不超过250000。   思路：这应该算是后缀自动机的经典应用了吧，我们先构造A的SAM，然后用A的SAM一次读入B的每一个字符，初始时状态在root处，此时最大匹配数为tmp=0，（这里的最大匹配数是指以当前读入的字符结尾，往前能匹配的最大长度），设当前到达的状态为p，最大
SPOJ 1812 LCS2 - Longest Common Substring II
Elijahqi
04-01 287
http://www.elijahqi.win/archives/2857 题意翻译 题面描述 给定一些字符串，求出它们的最长公共子串 输入格式 输入至多 10 10 行，每行包含不超过 100000 100000 个的小写字母，表示一个字符串 输出格式 一个数，最长公共子串的长度 若不存在最长公共子串，请输出 0 0 。 题目描述 A string is finite seq...
SPOJ LCS
lambda QAQ
09-02 361
求两个串的最长公共字串后缀自动机模板题#include<bits/stdc++.h> using namespace std;const int maxn = 251234,mlen = 26; int toid(char c){ return c - 'a'; }struct Node{ int len,fa,nex[mlen]; Node(int _len = 0){
SPOJ LCS （SAM）
დ♂Hany01's Blog Space♂ღ
04-01 217
Description 求两个串的lcs Solution 建出一个串的SAM，然后另一个串在上面匹配即可。 Source /************************************************ * Au: Hany01 * Date: Apr 1st, 2018 * Prob: SPOJ LCS * Email: hany01@foxmail...
SPOJ LCS (后缀自动机)
aekijekzn484879763的博客
07-21 110
SPOJ LCS Problem : 给两个串S、T，询问两个串的最长公共子串。 Solution :复(yu)习(xi)了一下后缀自动机。后缀自动机一定要结合后缀树来理解呀 ！！！后缀自动机的fail指针实现了原串的逆序的后缀树，后缀自动机的nt转移数组相当于后缀树上的后缀链。 #include <iostream> #include <string> usi...
SPOJ LCS（后缀自动机）
r_clover的博客
03-02 466
(http://acm.hust.edu.cn/vjudge/contest/view.action?cid=106071#problem/A) 题意：问串A，B的最长相同子串 解法：直接构造A的后缀自动机，搞一波lcp就好了。#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #includ
[SPOJ1812]LCS2 - Longest Common Substring II
srand(time(0))
04-04 1392
后缀自动机
字符串进阶——谈谈后缀全家桶(1) 后缀树
Myriad Dreamin 愿逐将属于我的奇迹
07-21 610
后缀是什么 然而如果单从字符串构造来讲任何S[i…len[S]]⊂S(1⩽i⩽len[S])S[i…len[S]]⊂S(1⩽i⩽len[S])S[i\dots len[S]]\subset S(1\leqslant i\leqslant len[S])均为字符串SSS的后缀。 后缀树 §§\S1 从Trie到后缀树的变化 在WikiEn上后缀树的定义如下Suffix Tree. I...
ACM中常用算法----字符串
码代码的猿猿
07-26 8058
ACM中常用算法—-字符串ACM中常用的字符串算法不多,主要有以下几种: Hash 字典树 KMP AC自动机 manacher 后缀数组 EX_KMP SAM(后缀自动机) 回文串自动机 下面来分别介绍一下:0. Hash字符串的hash是最简单也最常用的算法,通过某种hash函数将不同的字符串分别对应到不同的数字.进而配合其他数据结构或STL可以做到判重,统计,查询等操作. #### 字符串的h
mysql导入sql文件报错1067_导入数据库备份报错1067 – Invalid default value for ‘create_time’...
weixin_35873946的博客
01-28 1646
通过navicat工具导入psc数据库备份文件，报错如下，mysql版本5.7执行如下语句不通过DROP TABLE IF EXISTS `guard_user`;CREATE TABLE `guard_user` (`id` bigint(20) unsigned NOT NULL AUTO_INCREMENT COMMENT '操作员编号',`login_account` varchar(16...
spoj LCS 【后缀自动机】
beginendzrq
07-13 990
琦不会后缀自动机…… 是以前太浪了…… 所以所有东西都留到了noi前来学…… 马上狗牌退役了TAT（心塞qwq题目大意：给出两个串A,B，求A、B的最长公共子串对A建后缀自动机，然后用B去匹配，若能匹配上就转移到儿子，否则沿着parent树向上跳#include <iostream> #include <algorithm> #include <cstring> #include <cstdi
SPOJ - LCS Longest Common Substring （后缀自动机）
LP_Cong
10-17 242
求两个串的最长公共子串 保存我的模板！   #include&lt;bits/stdc++.h&gt; using namespace std; const int MAXN=500005; typedef long long ll; int step[MAXN]; int SZ=1; int ch[MAXN][30]; int fa[MAXN]; int last=1; void e...
spoj 1811 LCS - Longest Common Substring (后缀自动机)
whai的专栏
08-14 1719
spoj 1811 LCS - Longest Common Substring 题意： 给出两个串S, T, 求最长公共子串。 限制： |S|, |T| 思路： dp O(n^2) 铁定超时 后缀数组 O(nlog(n)) 在spoj上没试过，感觉也会被卡掉 后缀自动机 O(n) 我们考虑用SAM读入字符串B; 令当前状态为s，同时最大匹配长度为len;
SPOJ 1812 LCS2(SAM)
窝不是爱酱,喵~~~~
11-26 5547
转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents    by---cxlove   题目：求多个串的LCS http://www.spoj.pl/problems/LCS2/  先将一个串建SAM，然后用后面的串去匹配，对于每一个串，保存最大值，对于不同的串，更新最小值。 SAM结点多两个值，ml表示多个串的最小
SAM练习 spoj1811 LCS；spoj1812 LCS2；spoj8222
yxr0105
06-07 801
SAM的若干性质： 1、每个节点所代表的后缀一定是连续长度的。 2、x节点代表的[l1,r1]，pre[x]代表[l2,r2]，其中l1=r2+1。 3、pre指针的包含关系。 4、l排序满足拓扑排序。 spoj1811 LCS n 做的时候是在好久以前刚刚学完SAM的时候。。 当时naive啊。。好多SAM的性质都没理解上来照猫画虎硬搞。。 之后还给学弟口胡了几句。。。 按照
SPOJ LCS 最长公共子串 后缀自动机&后缀树(Ukkonen)
Butt Loem
09-17 2388
终于搞清楚了这两个恶心的数据结构。其实后缀树也不难写嘛。题目　　给定两个字符串a和b，求在a和b中都有出现的连续子串的最长长度。样例输入alsdfkjfjkdsal fdjskalajfkdsla 样例输出3 做法1　　使用后缀自动机。clj的课件讲得很详细了，这里不细说。主要说几件事： 后缀自动机的状态的本质是right集合（见课件），parent意味着right集合的最小扩充。时刻记着这一点可以
ssm智能家政保洁预约系统最新完整版程序+论文.zip
10-18
基于SSM框架的智能家政保洁预约系统，是一个旨在提高家政保洁服务预约效率和管理水平的平台。该系统通过集成现代信息技术，为家政公司、家政服务人员和消费者提供了一个便捷的在线预约和管理系统。 系统的主要功能包括： 1. **用户管理**：允许消费者注册、登录，并管理他们的个人资料和预约历史。 2. **家政人员管理**：家政服务人员可以注册并更新自己的个人信息、服务类别和服务时间。 3. **服务预约**：消费者可以浏览不同的家政服务选项，选择合适的服务人员，并在线预约服务。 4. **订单管理**：系统支持订单的创建、跟踪和管理，包括订单的确认、完成和评价。 5. **评价系统**：消费者可以在家政服务完成后对服务进行评价，帮助提高服务质量和透明度。 6. **后台管理**：管理员可以管理用户、家政人员信息、服务类别、预约订单以及处理用户反馈。 系统采用Java语言开发，使用MySQL数据库进行数据存储，通过B/S架构实现用户与服务的在线交互。系统设计考虑了不同用户角色的需求，包括管理员、家政服务人员和普通用户，每个角色都有相应的权限和功能。此外，系统还采用了软件组件化、精化体系结构、分离逻辑和数据等方法，以便于未来的系统升级和维护。 智能家政保洁预约系统通过提供一个集中的平台，不仅方便了消费者的预约和管理，也为家政服务人员提供了一个展示和推广自己服务的机会。同时，系统的后台管理功能为家政公司提供了强大的数据支持和决策辅助，有助于提高服务质量和管理效率。该系统的设计与实现，标志着家政保洁服务向现代化和网络化的转型，为管理决策和控制提供保障，是行业发展中的重要里程碑。
东方财富网抓取上市公司爬虫-EMSpider.zip
最新发布
10-18
东方财富网抓取上市公司爬虫-EMSpider
卫生健康系统-java-基于springBoot智能推荐的卫生健康系统（毕业论文+开题）
10-18
本基于智能推荐的卫生健康系统有管理员和用户两个角色。用户功能有个人中心，健康视频管理，视频类型管理，科室类型管理，医生信息管理，在线预约管理，我的收藏管理等。管理员功能有个人中心，用户管理，健康视频管理，视频类型管理，科室类型管理，医生信息管理，在线预约管理，健康论坛管理，我的收藏管理，留言板管理，系统管理等。因而具有一定的实用性。 本站是一个B/S模式系统，采用Spring Boot框架，MYSQL数据库设计开发，充分保证系统的稳定性。系统具有界面清晰、操作简单，功能齐全的特点，使得基于智能推荐的卫生健康系统管理工作系统化、规范化。
动态规划题解：45道经典问题解析
这些问题来源于不同的在线判题平台，如POJ、SPOJ、AOA和UVa。" 在动态规划（Dynamic Programming，简称DP）这个领域，这些题目提供了丰富的实践机会，帮助学习者深入理解如何应用动态规划方法解决实际问题。动态...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值