zoukankan      html  css  js  c++  java
  • poj 3294 Life Forms

    后缀数组的题目,把后缀连接起来,这个还是先二分答案,然后选取一段连续的height值,判断这些height代表的后缀有没有覆盖一半以上的字符串。

    得出答案的长度之后还要在枚举连续的heigh,判断有没有答案,有的话标记其中一个。

    最后再按照sa输出答案。这样就可以保证字典序。

    #include <iostream>
    #include <cstdio>
    #include <cstring>
    using namespace std;
    const int maxn=1e5+10000;
    int r[maxn];
    char a[1111];
    int *rank,height[maxn],sa[maxn];
    int wx[maxn],wy[maxn],cnt[maxn];
    int n,m;
    inline bool cmp(int *r,int a,int b,int l)
    {
        return r[a]==r[b]&&r[a+l]==r[b+l];
    }
    void da(int *r,int n,int m)
    {
        r[n+1]=0;
        int i,l,p,*x=wx,*y=wy,*t;
        memset(cnt,0,sizeof(int)*(m+1));
        for(int i=1;i<=n;i++) cnt[x[i]=r[i]]++;
        for(int i=1;i<=m;i++) cnt[i]+=cnt[i-1];
        for(int i=n;i>=1;i--) sa[cnt[x[i]]--]=i;
        for(l=1,p=1;p<n;l<<=1,m=p)
        {
            for(p=1,i=n-l+1;i<=n;i++) y[p++]=i;
            for(i=1;i<=n;i++)
            if(sa[i]>l)
            y[p++]=sa[i]-l;
            memset(cnt,0,sizeof(int)*(m+1));
            for(i=1;i<=n;i++) cnt[x[i]]++;
            for(i=1;i<=m;i++) cnt[i]+=cnt[i-1];
            for(i=n;i>=1;i--) sa[cnt[x[y[i]]]--]=y[i];
            for(t=x,x=y,y=t,p=1,x[sa[1]]=1,i=2;i<=n;i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],l)?p:++p;
        }
        rank=x;
        int j,k=0;
    	for(i=1;i<=n;height[rank[i++]]=k)
    	for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
    	return;
    }
    
    bool fla[1111];
    int c[maxn];
    bool ans[maxn];
    bool chk(int ret)
    {
        for(int i=2;i<=n;i++)
        if(height[i]>=ret)
        {
            int j=i+1;
            while(j<=n&&height[j]>=ret) j++;
            j--;
            memset(fla,0,sizeof(fla));
            for(int k=i-1;k<=j;k++)
            fla[c[sa[k]]]=1;
            int cnt=0;
            for(int k=1;k<=m;k++)
            cnt+=fla[k];
            if(cnt>m/2) return true;
            i=j;
        }
        return false;
    }
    
    int main()
    {
    //    freopen("in.txt","r",stdin);
        int ret;
        int tt=0;
        while(scanf("%d",&m),m)
        {
            if(tt++) printf("
    ");
            n=0;
            for(int i=1;i<=m;i++)
            {
                scanf("%s",a+1);
                ret=strlen(a+1);
                for(int j=1;j<=ret;j++)
                r[n+j]=a[j];
                for(int j=n+1;j<=n+ret;j++)
                c[j]=i;
                n+=ret;
                r[++n]=570+i;
            }
            da(r,n,680);
            int st=0,ed=1111,mid;
            while(st<ed)
            {
                mid=st+ed+1>>1;
                if(chk(mid)) st=mid;
                else ed=mid-1;
            }
            if(m==1)
            {
                for(int i=1;i<n;i++)
                printf("%c",r[i]);
                printf("
    ");
            }
            else if(st==0)
            printf("?
    ");
            else
            {
                memset(ans,0,sizeof(ans));
                for(int i=1;i<=n;i++)
                if(height[i]>=st)
                {
                    int j=i+1;
                    while(height[j]>=st) j++;
                    j--;
                    memset(fla,0,sizeof(fla));
                    for(int k=i-1;k<=j;k++)
                    fla[c[sa[k]]]=1;
                    int cnt=0;
                    for(int k=1;k<=m;k++)
                    cnt+=fla[k];
                    if(cnt>m/2)
                    {
                        ans[sa[i-1]]=1;
                    }
                    i=j;
                }
                for(int i=1;i<=n;i++)
                if(ans[sa[i]])
                {
                    for(int j=1;j<=st;j++)
                    printf("%c",r[j+sa[i]-1]);
                    printf("
    ");
                }
            }
        }
        return 0;
    }
    


  • 相关阅读:
    HybridAPP开发框架Ionic+AngularJS+Cordova搭建
    MySQL数据类型DECIMAL用法
    RabbitMQ用户管理
    Linux中Mysql root用户看不到mysql库问题解决方式
    RabbitMQ环境安装
    RabbitMQ简介
    Mysql根据一个基库生成其他库与其不同的库升级脚本
    关于中文乱码的解决方法(URL方式)
    java反射机制
    the jar file rt.jar has no source attachment
  • 原文地址:https://www.cnblogs.com/pangblog/p/3253861.html
Copyright © 2011-2022 走看看