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

P2463 [SDOI2008]Sandy的卡片

时间:2019-04-04 09:17:01      阅读:138      评论:0      收藏:0      [点我收藏+]

标签:temp   最长匹配   max   题意   and   math   lse   break   字符串   

题目链接

题意分析

首先对于一个可以匹配的字符串

我们发现 差分之后出除了最后一位外是相等的

A 1 4 6 7

B 3 6 8 2

差分之后

A 3 2 1 -3

B 3 2 -6 2

所以我们要求的就是拆分之后最长匹配长度+1

首先 我们将差分之后的拼成一个长串

然后建出\(SA\)

发现答案具有可二分性

然后我们再使用\(height\)数组 将\(lcp≥now\)的后缀分成一组

如果存在共计\(n\)个串后缀分成一组

说明答案存在

CODE:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<string>
#include<queue>
#include<map>
#include<stack>
#include<list>
#include<set>
#include<deque>
#include<vector>
#include<ctime>
#define ll long long
#define inf 0x7fffffff
#define N 100008
#define IL inline
#define M 2008611
#define D double
#define maxn 30000
#define R register
using namespace std;
template<typename T>IL void read(T &_)
{
    T __=0,___=1;char ____=getchar();
    while(!isdigit(____)) {if(____=='-') ___=0;____=getchar();}
    while(isdigit(____)) {__=(__<<1)+(__<<3)+____-'0';____=getchar();}
    _=___ ? __:-__;
}
/*-------------OI使我快乐-------------*/
int n,m,key,maxx,tot;
int s[M],tmp[M],vis[M],rnk[M],SA[M],cdy[M],wzy[M],hei[M];
int bel[M],sta[M];
IL void Rsort()
{
    for(R int i=0;i<=key;++i) vis[i]=0;
    for(R int i=1;i<=n;++i) vis[cdy[i]]++;
    for(R int i=1;i<=key;++i) vis[i]+=vis[i-1];
    for(R int i=n;i;--i) SA[vis[cdy[wzy[i]]]--]=wzy[i];
}
IL void get_SA()
{
    for(R int i=1;i<=n;++i) cdy[i]=s[i],wzy[i]=i;Rsort();
    for(R int x=1;x<=n;x<<=1)
    {
        int cnt=0;
        for(R int i=n-x+1;i<=n;++i) wzy[++cnt]=i;
        for(R int i=1;i<=n;++i) if(SA[i]>x) wzy[++cnt]=SA[i]-x;
        Rsort();swap(cdy,wzy);cdy[SA[cnt=1]]=1;
        for(R int i=2;i<=n;++i)
        cdy[SA[i]]=(wzy[SA[i]]==wzy[SA[i-1]]&&wzy[SA[i]+x]==wzy[SA[i-1]+x] ? cnt:++cnt);
        if(cnt==n) break;
        else key=cnt;
    }
}
IL void get_hei()
{
    for(R int i=1;i<=n;++i) rnk[SA[i]]=i;
    int lat,k=0;
    for(R int i=1;i<=n;++i)
    {
        if(k) --k;
        lat=SA[rnk[i]-1];
        while(s[lat+k]==s[i+k]) ++k;
        hei[rnk[i]]=k;
    }
}
IL bool check(int now)
{
    ++tot;int cnt=0;
    for(R int i=1;i<=n;++i)
    {
        if(hei[i]<now) ++tot,cnt=0;
        if(vis[bel[SA[i]]]!=tot) vis[bel[SA[i]]]=tot,++cnt;
        if(cnt==m) return 1;
    }
    return 0;
}
int main()
{
//  freopen("222.in","r",stdin);
//  freopen("222.out","w",stdout);
    read(m);key=m+N;
    for(R int x=1;x<=m;++x)
    {
        int len;read(len);maxx=max(maxx,len);
        for(R int i=1;i<=len+1;++i) tmp[i]=0;
        for(R int i=1;i<=len;++i) read(tmp[i]);
        for(R int i=1;i<=len;++i) s[++n]=tmp[i]-tmp[i+1]+maxn,bel[n]=x;
        s[++n]=x+N;
    }
//  for(R int i=1;i<=n;++i) printf("%d%c",s[i],(i==n ? '\n':' '));
//  for(R int i=1;i<=n;++i) printf("%d%c",bel[i],(i==n ? '\n':' '));
    get_SA();get_hei();
    int le=0,ri=n,ans=0;
    memset(vis,0,sizeof vis);
    while(le<=ri)
    {
        int mid=(le+ri)>>1;
        if(check(mid)) ans=mid,le=mid+1;
        else ri=mid-1;
    }
    printf("%d\n",ans+1);
//  for(R int i=1;i<=n;++i) printf("%d%c",SA[i],(i==n ? '\n':' '));
//  fclose(stdin);
//  fclose(stdout);
    return 0;
}

HEOI 2019 RP++

P2463 [SDOI2008]Sandy的卡片

标签:temp   最长匹配   max   题意   and   math   lse   break   字符串   

原文地址:https://www.cnblogs.com/LovToLZX/p/10652769.html

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