uvalive4513

https://vjudge.net/problem/UVALive-4513

终于做出来了。。。。。。。。。

各种sb错误,最后对拍出来了,还没改对。。。。。。。。。。。。。。。。。。。。。

快半天的时间。。。。。。。。。。

看上去第一眼就猜出了一个结论:这些相同的字符串肯定是按后缀排序的顺序存在的。

于是就想出了一个方法:求出sa和lcp,然后扫描每个长为m的区间,rmq其中最小的lcp和最大的sa作为答案(扫描是指枚举左端点,右端点是左端点+m-2)

这样复杂度nlogn,可以接受,用单调队列?(不清楚,没学好),于是我用了线段树。

(请无视)然后各种wa,各种改,搞得我怀疑自己的方法是错的,然后看了看题解,好像都是二分判定。再想了想,貌似是对的,但是不确定,就翻了翻别人的代码,看见和自己的想法差不多,就确定是对的了,然后就各种调试,看见一个又一个sb错误:update1和update搞反找了半天。。。然后又对拍拍了半天拍出来了组错误数据,发现sa[0]=0这样的串会错,于是就改掉,对了。。。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 80010
struct data
{
    int val,pos;
};
int n,m,k;
char s[N];
int Rank[N],sa[N],temp[N],lcp[N],tree[N<<2],tree1[N<<2];
bool cp(int i,int j)
{
    if(Rank[i]!=Rank[j]) return Rank[i]<Rank[j];
    int ri=i+k<n?Rank[i+k]:-1;
    int rj=j+k<n?Rank[j+k]:-1;
    return ri<rj;
}
void update(int l,int r,int x,int pos,int num)
{
    if(l==r)
    {
        tree[x]=num;
        return;
    }
    if(pos>(l+r)/2) update((l+r)/2+1,r,x*2+1,pos,num);
    else update(l,(l+r)/2,x*2,pos,num);
    tree[x]=min(tree[x*2],tree[x*2+1]);
}
int query(int l,int r,int x,int a,int b)
{
    if(l>b||r<a) return N;
    if(l>=a&&r<=b) return tree[x];
    int i=query(l,(l+r)/2,x*2,a,b);
    int j=query((l+r)/2+1,r,x*2+1,a,b);
    return min(i,j);
}
void update1(int l,int r,int x,int pos,int Pos)
{
    if(l==r)
    {
        tree1[x]=Pos;
        return;
    }
    if(pos>(l+r)/2) update1((l+r)/2+1,r,x*2+1,pos,Pos);
    else update1(l,(l+r)/2,x*2,pos,Pos);
    tree1[x]=max(tree1[x*2],tree1[x*2+1]);
}
int query1(int l,int r,int x,int a,int b)
{
    if(l>b||r<a) return -N;
    if(l>=a&&r<=b) return tree1[x];
    int i=query1(l,(l+r)/2,x*2,a,b);
    int j=query1((l+r)/2+1,r,x*2+1,a,b);
    return max(i,j);
}
void Sa()
{
    n=strlen(s);
    for(int i=0;i<n;i++)
    {
        Rank[i]=s[i]; sa[i]=i;
    }
    for(k=1;k<=n;k*=2)
    {
        sort(sa,sa+n,cp);
        temp[sa[0]]=1;
        for(int i=1;i<n;i++) temp[sa[i]]=temp[sa[i-1]]+(cp(sa[i-1],sa[i]));
         for(int i=0;i<n;i++) Rank[i]=temp[i];
    }
}
void Lcp()
{
    n=strlen(s);
    for(int i=0;i<n;i++) Rank[sa[i]]=i;
    int h=0;
    for(int i=0;i<n;i++)
    {
        int j=sa[Rank[i]-1];
        if(Rank[i]<=0) continue;
        if(h>0) h--;
        for(;((i+h<n)&&(j+h<n));h++)
        {
            if(s[i+h]!=s[j+h]) break;
        }
        lcp[Rank[i]-1]=h;
    }
    update1(1,n,1,1,sa[0]);
    for(int i=0;i<n-1;i++)
    {
        update(1,n,1,i+1,lcp[i]);
        update1(1,n,1,i+2,sa[i+1]);
    }
}
void solve()
{
    memset(sa,0,sizeof(sa));
    memset(lcp,0,sizeof(lcp));
    memset(tree,0x3f3f,sizeof(tree));
    memset(tree1,0,sizeof(tree1));
    memset(Rank,0,sizeof(Rank));
    Sa(); Lcp();
    int l=1;
    data ans; ans.val=0; ans.pos=0;
    for(int i=1;i+m-2<n;i++)
    {
        data x;
        x.val=query(1,n,1,i,i+m-2);
        x.pos=query1(1,n,1,i,i+m-1);
        if(x.val>ans.val||(x.val==ans.val&&x.pos>ans.pos))
        {
            ans.val=x.val;
            ans.pos=x.pos;
        }
    }
    if(ans.val==0) printf("none\n");
    else printf("%d %d\n",ans.val,ans.pos);
}
int main()
{
    while(1)
    {
        scanf("%d",&m);
        if(!m) break;
        scanf("%s",s);
        if(m==1)
        {
            printf("%d %d\n",strlen(s),0);
            continue;
        }
        n=strlen(s);
        solve();
    }
    return 0;
}

时间: 2024-10-12 20:36:31

uvalive4513的相关文章

UVALive4513 Stammering Aliens(哈希法 | 后缀数组)

题目链接:http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=12580 [思路] 求出现次数不小于k次的最长可重复子串和最后的出现位置. 法一: 后缀数组,二分长度,划分height.时间复杂度为O(nlogn) 法二: Hash法.构造字符串的hash函数,二分长度,求出hash(i,L)后排序,判断是否存在超过k个相同hash 值得块即可.时间为O(nlog2n). 注意划分height一定要精确且如果m=1需要特判(对于

9月——都已经9月了还不好好刷题?。。

UVALive4513 (字符串hash 题目:在一个串中找到至少出现m次的最长的串. 思路:直接hash,然后二分答案. /* * @author: Cwind * http://www.cnblogs.com/Cw-trip/ */ //#pragma comment(linker, "/STACK:102400000,102400000") #include <iostream> #include <map> #include <algorithm&