【bzoj1535】[POI2005]Sza-Template【乱搞】

这道题我用一个思维难度和代码难度都比较低的算法过了,感觉性价比还是蛮高的qwq
我们首先设 f [ i ] s [ i . . n ] s [ 1.. n ] 匹配的最大长度。
对于怎么求 f ,我是直接二分+哈希乱搞求的。
考虑怎样的一个前缀的长度 l ,可以作为模板。
对于所有 f [ i ] l i ,也就是从i开始可以和前缀完全匹配,我们把这些 i 从小到大排好。显然要满足任何两个相邻的 i ,它们的差必须小于等于 l
现在就好办了。我们按照每个位置的 f 值从小到大加进来。然后维护两个 s e t ,一个存位置,一个存相邻两个的距离,随便搞搞就好。
时间复杂度 O ( n l o g 2 n ) 。差一点就超时了= =

#include<bits/stdc++.h>
using namespace std;
typedef unsigned long long ull;
const int N=500005;
int n,ans;
ull base[N],h[N];
char s[N];
set<int> s1;
multiset<int> s2;
ull gethash(int l,int r){
    return h[r]-h[l-1]*base[r-l+1];
}
struct data{
    int x,v;
}a[N];
bool cmp(data a,data b){
    return a.v>b.v;
}
void insert(int x){
    s1.insert(x);
    set<int>::iterator it1=s1.lower_bound(x);
    it1--;
    set<int>::iterator it2=s1.upper_bound(x);
    s2.erase(s2.find(*it2-*it1));
    s2.insert(x-*it1);
    s2.insert(*it2-x);
}
int main(){
    scanf("%s",s+1);
    n=strlen(s+1);
    base[0]=1;
    for(int i=1;i<=n;i++){
        base[i]=base[i-1]*48271;
        h[i]=h[i-1]*48271+s[i]-'a';
    }
    for(int i=1;i<=n;i++){
        int l=0,r=n-i+1,mid,ans;
        while(l<=r){
            mid=(l+r)/2;
            if(gethash(1,mid)==gethash(i,i+mid-1)){
                ans=mid;
                l=mid+1;
            }else{
                r=mid-1;
            }
        }
        a[i].x=i;
        a[i].v=ans;
    }
    sort(a+1,a+n+1,cmp);
    s1.insert(0);
    s1.insert(n+1);
    s2.insert(n+1);
    for(int i=n,j=1;i>=1;i--){
        while(a[j].v>=i){
            insert(a[j].x);
            j++;
        }
        multiset<int>::iterator it=s2.end();
        it--;
        if(*it<=i){
            ans=i;
        }
    }
    printf("%d\n",ans);
    return 0;
}

猜你喜欢

转载自blog.csdn.net/ez_2016gdgzoi471/article/details/81810320