POJ1743(二分+后缀数组)

题目

开始后缀数组啦,开心,什么都看不懂,但只要理解好几个数组的意思,抄起模板来随便过,呜呜呜,真的理解不了国家集训队的论文了,神仙高中生啊,照着看了几遍了还是只会抄代码
楼教主男人八题之一(我不配做男人啦
找序列的最长不重叠子串

题解

太经典了,论文中写的比我的思考清晰多了(什么你不知道是那篇论文?你还学什么后缀数组)

代码

主函数的格式有点小问题,Devc太垃圾了,没有一键缩进,一键校正代码的功能(但是用起来是真)

#include<iostream>
#include<algorithm>
#include<cstdio>
#define N 40005
using namespace std;
const int maxn=40005;
#define INF (1<<30)
int t1[N],t2[N],sa[N],h[N],rk[N],c[10*N],a[N];
char s[N];
int m,n;
void calcsa(int n,int m){
    int *x=t1,*y=t2,p=0,f=0;
    for(int i=1;i<=m;i++)c[i]=0;
    for(int i=1;i<=n;i++)c[x[i]=a[i]]++;
    for(int i=1;i<=m;i++)c[i]+=c[i-1];
    for(int i=n;i;i--)sa[c[x[i]]--]=i;
    for(int i=1;i<=n&&p<=n;i<<=1){
        p=0;
        for(int j=n-i+1;j<=n;j++)y[++p]=j;
        for(int j=1;j<=n;j++)if(sa[j]>i)y[++p]=sa[j]-i;
        for(int j=1;j<=m;j++)c[j]=0;
        for(int j=1;j<=n;j++)c[x[y[j]]]++;
        for(int i=1;i<=m;i++)c[i]+=c[i-1];
        for(int j=n;j;j--)sa[c[x[y[j]]]--]=y[j];
        swap(x,y);x[sa[1]]=1;p=2;
        for(int j=2;j<=n;j++)
        x[sa[j]]=y[sa[j]]==y[sa[j-1]]&&y[sa[j]+i]==y[sa[j-1]+i]?p-1:p++;
        m=p;
    }
    for(int i=1;i<=n;i++)rk[sa[i]]=i;
    for(int i=1;i<=n;i++){
        int j=sa[rk[i]-1];
        if(f)f--;while(a[i+f]==a[j+f])f++;
        h[rk[i]]=f;
    }
}
int jiancha(int k)
{
	int minm=INF,maxm=0;
	for(int i=2;i<=n;i++)
	{
		//cout<<h[i]<<endl;
		if(h[i]<k)
		{
			minm=INF;
			maxm=0;
		}
		else
		{
			minm=min(minm,min(sa[i],sa[i-1]));
			maxm=max(maxm,max(sa[i],sa[i-1]));
			if(maxm-minm>=k)
			{
				return 1;
			}
		}
	}
	return 0;
}
int b[maxn];
int main()
{
	int x1;
while(~scanf("%d",&n)&&n)
{
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&b[i]);
	}
	for(int i=2;i<=n;i++)
	{
		a[i-1]=b[i]-b[i-1]+88;
	}
	a[n]=0;
	calcsa(n,88*2);
	/*for(int i=1;i<=n-1;i++)
	{
		cout<<a[rk[i]]<<' ';
	}
	cout<<endl;
	for(int i=1;i<=n-1;i++)
	{
		cout<<' '<<h[i]<<' ';
	}
	cout<<endl;
	for(int i=1;i<=n-1;i++)
	{
		cout<<' '<<sa[i]<<' ';
	}
	cout<<endl;*/
	int lf=0,rg=n;
	while(lf<rg)
	{
		//cout<<lf<<endl;
		int mid=(lf+rg+1)/2;
		if(jiancha(mid))
		{
			//cout<<lf<<endl;
			lf=mid;
		}
		else
		{
			rg=mid-1;
		}
	}
	if(lf+1<5)
	{
		cout<<'0'<<endl;
	}
	else
	{
		cout<<lf+1<<endl;
	}
}
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值