P2463 [SDOI2008]Sandy的卡片
题目链接
题意分析
首先对于一个可以匹配的字符串
我们发现 差分之后出除了最后一位外是相等的
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;
}
