[SDOI2008]Sandy的卡片
题目描述
Sandy和Sue的热衷于收集干脆面中的卡片。
然而,Sue收集卡片是因为卡片上漂亮的人物形象,而Sandy则是为了积攒卡片兑换超炫的人物模型。
每一张卡片都由一些数字进行标记,第i张卡片的序列长度为Mi,要想兑换人物模型,首先必须要集够N张卡片,对于这N张卡片,如果他们都有一个相同的子串长度为k,则可以兑换一个等级为k的人物模型。相同的定义为:两个子串长度相同且一个串的全部元素加上一个数就会变成另一个串。
Sandy的卡片数远远小于要求的N,于是Sue决定在Sandy的生日将自己的卡片送给Sandy,在Sue的帮助下,Sandy终于集够了N张卡片,但是,Sandy并不清楚他可以兑换到哪个等级的人物模型,现在,请你帮助Sandy和Sue,看看他们最高能够得到哪个等级的人物模型。
输入格式:
第一行为一个数N,表示可以兑换人物模型最少需要的卡片数,即Sandy现在有的卡片数
第i+1行到第i+N行每行第一个数为第i张卡片序列的长度Mi,之后j+1到j+1+Mi个数,用空格分隔,分别表示序列中的第j个数
输出格式:
一个数k,表示可以获得的最高等级。
输入样例#1:
22 1 23 4 5 9
输出样例#1:
2
说明
数据范围:
30%的数据保证n<=50
100%的数据保证n<=1000
题解:
连数据范围都不给是要闹鬼咯,垃圾题目,吃枣药丸。
将字符串差分后,就是最长公共子串。
这样就可以把所有串连起来,求height数组,然后对答案进行二分,找到一段height值全部大于mid且包含所有串的一段区间,这个mid值就合法。
代码:
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int max_n = 300010;
struct point
{
int wh;
int height;
}p[max_n];
char s[max_n];
int t1[max_n],t2[max_n],sa[max_n],rank[max_n],cc[max_n],height[max_n];
LL val[max_n],maxn[max_n],minn[max_n],ans[max_n],num[max_n],size[max_n];
int fa[max_n];
int len;
LL inf;
inline bool cmp1(point a,point b)
{
return a.height>b.height;
}
inline bool cmp(int *y,int a,int b,int k)
{
int ranka1=y[a];
int rankb1=y[b];
int ranka2=a+k<len ? y[a+k] : -1;
int rankb2=b+k<len ? y[b+k] : -1;
return ranka1==rankb1 && ranka2==rankb2;
}
inline void make_sa()
{
int *x=t1,*y=t2,m=26;
for(int i=0; i<m; ++i) cc[i]=0;
for(int i=0; i<len; ++i) cc[x[i]=s[i]-'a']++;
for(int i=1; i<m; ++i) cc[i]+=cc[i-1];
for(int i=len-1; i>=0; --i) sa[--cc[x[i]]]=i;
for(int k=1; k<=len; k<<=1)
{
int p=0;
for(int i=len-k; i<len; ++i) y[p++]=i;
for(int i=0; i<len; ++i)
if(sa[i]>=k) y[p++]=sa[i]-k;
for(int i=0; i<m; ++i) cc[i]=0;
for(int i=0; i<len; ++i) cc[x[y[i]]]++;
for(int i=1; i<m; ++i) cc[i]+=cc[i-1];
for(int i=len-1; i>=0; --i) sa[--cc[x[y[i]]]]=y[i];
swap(x,y);
m=1; x[sa[0]]=0;
for(int i=1; i<len; ++i)
x[sa[i]]=cmp(y,sa[i],sa[i-1],k) ? m-1 : m++;
}
for(int i=0; i<len; ++i) rank[sa[i]]=i;
}
inline void make_height()
{
int k=0;
for(int i=0; i<len; ++i)
{
if(!rank[i]) continue;
int j=sa[rank[i]-1];
if(k) k--;
while(s[i+k]==s[j+k]) k++;
height[rank[i]]=k;
}
}
inline int find(int a)
{
if(fa[a]==a) return a;
else return fa[a]=find(fa[a]);
}
inline int unionn(int a,int b)
{
int a1=find(a),b1=find(b);
fa[b1]=a1;
size[a1]+=size[b1];
maxn[a1]=max(maxn[a1],maxn[b1]);
minn[a1]=min(minn[a1],minn[b1]);
return a1;
}
inline void get_ans()
{
for(int i=0; i<len; ++i)
p[i].wh=i,p[i].height=height[i];
for(int i=0; i<len; ++i) fa[i]=i,size[i]=1,maxn[i]=val[i],minn[i]=val[i];
sort(p,p+len,cmp1);
for(int i=0; i<len; ++i)
{
if(!p[i].wh) continue;
int a1=find(sa[p[i].wh]),b1=find(sa[p[i].wh-1]);
num[p[i].height]+=size[a1]*size[b1];
ans[p[i].height]=max(ans[p[i].height],max(maxn[a1]*maxn[b1],minn[a1]*minn[b1]));
unionn(a1,b1);
}
for(int i=len-2; i>=0; --i) num[i]+=num[i+1],ans[i]=max(ans[i],ans[i+1]);
for(int i=len-1; i>=0; --i)
if(ans[i]<-inf)
ans[i]=0;
for(int i=0; i<len; ++i)
printf("%lld %lld\n",num[i],ans[i]);
}
int main()
{
freopen("savour.in","r",stdin);
freopen("savour.out","w",stdout);
scanf("%d%s",&len,s);
memset(ans,-0x7f,sizeof(ans));
inf=-ans[0]-1;
for(int i=0; i<len; ++i)
scanf("%lld",&val[i]);
make_sa();
make_height();
get_ans();
}