【BZOJ2251】外星联络 后缀数组

原创 2016年06月01日 19:59:14

题目大意是找出一个字符串中所有出现次数大于1的字串并且按字典序输出出现次数

height数组上乱搞。。。

为了防止统计重复,我们统计时可以统计从sa[i]开始的字串,也就是长度Len>=height[i]+1的字串,这样可以保证不重复。

/**************************************************************
    Problem: 2251
    User: cqyzhb
    Language: C++
    Result: Accepted
    Time:120 ms
    Memory:1444 kb
****************************************************************/
 
#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
#define MAXN 6005
int N,A[MAXN],M;
char s[MAXN];
void Init(){    scanf("%d",&N); scanf("%s",s+1); for(int i=1;i<=N;i++)A[i]=s[i]-'0'+1; M=2;return ;  }
int wa[MAXN],wb[MAXN],wc[MAXN],sa[MAXN],h[MAXN],rank[MAXN];
void Debug()
{
    for(int i=1;i<=N;printf("%d ",A[i]),i++);
    putchar('\n');
    for(int i=1;i<=N;printf("%d ",sa[i]),i++);
}
void make_sa()
{
    int *t,*x=wa,*y=wb,i,k,cnt;
    for(i=1;i<=M;i++)wc[i]=0;
    for(i=1;i<=N;i++)wc[x[i]=A[i]]++;
    for(i=1;i<=M;i++)wc[i]+=wc[i-1];
    for(i=N;i>=1;i--)sa[wc[x[i]]--]=i;
    for(k=1;k<=N;k=k*2)
    {
        cnt=0;
        for(i=N-k+1;i<=N;i++)y[++cnt]=i;
        for(i=1;i<=N;i++)if(sa[i]>k)y[++cnt]=sa[i]-k;
        for(i=1;i<=M;i++)wc[i]=0;
        for(i=1;i<=N;i++)wc[x[y[i]]]++;
        for(i=1;i<=M;i++)wc[i]+=wc[i-1];
        for(i=N;i>=1;i--)sa[wc[x[y[i]]]--]=y[i];
        for(t=x,x=y,y=t,cnt=1,x[sa[1]]=1,i=2;i<=N;i++)
        x[sa[i]]= (y[sa[i-1]]==y[sa[i]])&&(y[sa[i-1]+k]==y[sa[i]+k]) ? cnt : ++cnt ;
        M=cnt; if(cnt==N)break;
    }
    for(i=1;i<=N;i++)rank[sa[i]]=i;
    k=0;
//  Debug();
    for(i=1;i<=N;h[rank[i]]=k,i++)
    {
        if(k)k--;
        for(int j=sa[rank[i]-1];A[i+k]==A[j+k];k++);
    }
    return ;
}
int main()
{
    //freopen("in.txt","r",stdin);
    Init();
    make_sa();
    for(int i=1;i<=N;i++)
    {
        for(int j=h[i]+1;sa[i]+j-1<=N;j++)
        {
            int l=i,r;
        //  for(l=i;h[l]>=j && l>=1;l--);
            for(r=i+1;h[r]>=j && r<=N;r++);
            if(r-l>1)printf("%d\n",r-l);
        }
    }
    return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

BZOJ 2251: [2010Beijing Wc]外星联络|后缀数组

后缀数组求出height数组 然后扫一边 暴力统计出解就可以了 #include #include #include #include #include #include #include #i...
  • ws_yzy
  • ws_yzy
  • 2016年01月12日 15:30
  • 440

BZOJ 2251 Beijing WC 2010 外星联络 后缀数组

给出一个01串,求:出现次数大于1 的子串所出现的次数。 不过数据范围有点迷。。 后缀数组模板题。 构建出后缀数组以后,每个后缀都会产生子串。 比如ababa:a aba ababa ba b...

BZOJ2251: [2010Beijing Wc]外星联络 后缀数组

2251: [2010Beijing Wc]外星联络 Time Limit: 30 Sec  Memory Limit: 256 MB Submit: 791  Solved: 475 [Subm...
  • Oakley_
  • Oakley_
  • 2017年01月12日 20:20
  • 307

BZOJ2251 外星联络 (Trie树)

题目大意给出一个01串,要求按照字典序输出在串中出现两次以上的子串的出现次数。题解我是在hzwer学长的后缀数组专题中看到这道题的,但是我并没有很理解黄学长的暴力是怎么思考的…于是我在另一个地方得到了...

bzoj2251 外星联络

trie树的dfs
  • sdfzyhx
  • sdfzyhx
  • 2016年06月19日 23:07
  • 527

[2010Beijing Wc]外星联络 后缀数组

题目要求求出所有的子串,考虑到每个子串都必定是一个后缀的前缀,可以使用后缀数组。         对于后缀数组sa[],发现按顺序i=0 to n-1;j=sa[] to n-1枚举,其中每个子串必...

【bzoj 1923】[Sdoi2010]外星千足虫(高斯消元)

执笔生宣,浅墨入画,绘一幅丹青风华

[题解]bzoj1923(SDOI2010)外星千足虫

DescriptionInput第一行是两个正整数 N, M。 接下来 M行,按顺序给出 Charles 这M次使用“点足机”的统计结果。每行 包含一个“01”串和一个数字,用一个空格隔开。“01”串...

[BZOJ1923][Sdoi2010]外星千足虫(高斯消元)

题目描述传送门题解高斯消元求解异或方程组 和普通消元的思路差不多,同样是找到一个可用的交换,然后把原来的加减运算全部改为异或 如果a(i,i)无法交换成0的话,说明这里有可能出现无解或多解,那么回...

【bzoj1923】[Sdoi2010]外星千足虫

高斯消元
  • KikiDMW
  • KikiDMW
  • 2017年03月19日 20:05
  • 246
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【BZOJ2251】外星联络 后缀数组
举报原因:
原因补充:

(最多只允许输入30个字)