zoukankan      html  css  js  c++  java
  • hdu5442(2015长春赛区网络赛1006)后缀数组+KMP /最小表示法?

    题意:给定一个由小写字母组成的长度为 n 的字符串,首尾相连,可以从任意一个字符开始,顺时针或逆时针取这个串(长度为 n),求一个字典序最大的字符串的开始字符位置和顺时针或逆时针。如果有多个字典序最大的字符串,优先选择开始位置靠前的,如果开始位置相同,优先选择顺时针。

    这种字符串的问题,第一反应是后缀数组,为了达到首尾相连的目的,所以先复制了一个两倍长的该字符串,然后再将串倒置,也弄成两倍长度,得到顺逆时针的两倍长的串,并对这两个字符串分别做后缀数组,得到 sa 数组(该串字典序第几小的后缀的开始字符是第几个)由于字符串的后缀长度必然不等,因此所有后缀都有确定的名次。

    由于我们需要的是该串的循环n个串,所以对于每个 2*n 的串,它的后缀中我们需要的只有前 n 个,而每个串我们需要的也只有它的前 n 个字符,那么就出现了有可能有多个后缀的前 n 个字符是相同的,但是由于后缀的存在,使这些串的名次就有了差别。因此 sa 数组中我们得到的字典序最大的一个串,它的前 n 个字符只能确定是字典序最大的其中之一,因此我们就需要在 2*n 长度的串中找到第一个符合我们要求的串。所以我从得到的字典序最大的串中取出前 n 个字符,用这个串作为模式串在 2*n 串中 KMP 。

    在顺时针的串中匹配到开始位置最靠前的第一个相同的串,而因为逆时针的串是将原串的所有字符导致,所以在逆时针串中开始位置越靠前,则它在原串中的开始位置就越靠后,因此我们需要KMP匹配找出的就是在 2*n 串中开始位置最靠后并且在 n 前的一个相同串。

    最后再将得到的两个串比较他们的各个情况,输出最符合要求的一个。

    恩,rank竟然变成了关键字,导致我当时CE了一发。

    后来在被问到的时候我仔细再想了一下,觉得顺时针貌似并不需要用KMP,因为多个前 n 个字符相同的后缀,由于开始位置靠前的后缀长度会更长,所以字典序就会更大,所以我们取一个字典序最大的后缀,它在前 n 个字符相同的后缀串中开始位置也应该是更靠前的,所以我们用后缀数组得到的字典序最大的串就应该是我们需要的一个串,暂时我还没有想到反例。

    但是逆时针却必须要匹配一下,因为字典序最大的,在 2*n 的串中开始位置更靠前,但在原串中开始位置就会越靠后,所以就要匹配得到在 2*n 串中最靠后的一个。例子就是 abcabc ,逆时针取的时候字典序最大的后缀开始位置是原串中的最后一个 c ,因此我们需要找出第三个字符 c 。所以我用了KMP。

    但是据说有个做法是最小表示法,然而我并不知道那是什么,所以我还是用后缀数组+KMP过的Orz……

      1 #include <iostream>
      2 #include <string.h>
      3 #include <algorithm>
      4 #include <stdio.h>
      5 using namespace std;
      6 const int MAXN=40010;
      7 
      8 int sa[MAXN];//SA数组,表示将S的n个后缀从小到大排序后把排好序的
      9 //的后缀的开头位置顺次放入SA中
     10 int t1[MAXN],t2[MAXN],c[MAXN];//求SA数组需要的中间变量,不需要赋值
     11 int rk[MAXN],height[MAXN];
     12 void build_sa(int s[],int n,int m)
     13 {
     14     int i,j,p,*x=t1,*y=t2;
     15     for(i=0;i<m;i++)c[i]=0;
     16     for(i=0;i<n;i++)c[x[i]=s[i]]++;
     17     for(i=1;i<m;i++)c[i]+=c[i-1];
     18     for(i=n-1;i>=0;i--)sa[--c[x[i]]]=i;
     19     for(j=1;j<=n;j<<=1)
     20     {
     21         p=0;
     22         for(i=n-j;i<n;i++)y[p++]=i;//后面的j个数第二关键字为空的最小
     23         for(i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
     24         for(i=0;i<m;i++)c[i]=0;
     25         for(i=0;i<n;i++)c[x[y[i]]]++;
     26         for(i=1;i<m;i++)c[i]+=c[i-1];
     27         for(i=n-1;i>=0;i--)sa[--c[x[y[i]]]]=y[i];
     28         swap(x,y);
     29         p=1;x[sa[0]]=0;
     30         for(i=1;i<n;i++)
     31             x[sa[i]]=y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+j]==y[sa[i]+j]?p-1:p++;
     32         if(p>=n)break;
     33         m=p;//下次基数排序的最大值
     34     }
     35 }
     36 
     37 void getHeight(int s[],int n)
     38 {
     39     int i,j,k=0;
     40     for(i=0;i<=n;i++)rk[sa[i]]=i;
     41     for(i=0;i<n;i++)
     42     {
     43         if(k)k--;
     44         j=sa[rk[i]-1];
     45         while(s[i+k]==s[j+k])k++;
     46         height[rk[i]]=k;
     47     }
     48 }
     49 
     50 char str[MAXN];
     51 int s0[MAXN];
     52 int s1[MAXN];
     53 int ans0[MAXN];
     54 int ans1[MAXN];
     55 int pp1[MAXN],pp0[MAXN];
     56 
     57 int check(int n){
     58     for(int i=0;i<n;++i){
     59         if(ans0[i]<ans1[i])return 1;
     60         else if(ans0[i]>ans1[i])return 0;
     61     }
     62     return -1;
     63 }
     64 
     65 int main(){
     66     int T;
     67     scanf("%d",&T);
     68     while(T--){
     69         int n;
     70         scanf("%d",&n);
     71         scanf("%s",str);
     72         for(int i=0;i<=n;i++)s0[i]=str[i]-'a'+1;
     73         for(int i=0;i<=n;i++)s0[n+i]=str[i]-'a'+1;
     74         s0[2*n]=0;
     75 /*        for(int i=0;i<=2*n;++i)printf("%d ",s0[i]);
     76         printf("
    ");*/
     77         build_sa(s0,2*n+1,28);
     78 /*        for(int i=0;i<=2*n;++i)printf("%d ",sa[i]);
     79         printf("
    ");*/
     80         int p0,p1;
     81         for(int i=2*n;i>=1;--i){
     82             if(sa[i]<n){p0=sa[i];break;}
     83         }
     84         for(int i=0;i<n;++i){
     85             ans0[i]=s0[p0+i];
     86 //            printf("%d ",ans0[i]);
     87         }
     88 /*        printf("
    ");
     89         printf("
    ");*/
     90 
     91         for(int i=0,j=n-1;i<n;i++,j--)s1[j]=str[i]-'a'+1;
     92         for(int i=0,j=n-1;i<n;i++,j--)s1[n+i]=str[j]-'a'+1;
     93         s1[2*n]=0;
     94 /*        for(int i=0;i<=2*n;++i)printf("%d ",s1[i]);
     95         printf("
    ");*/
     96         build_sa(s1,2*n+1,28);
     97 /*        for(int i=0;i<=2*n;++i)printf("%d ",sa[i]);
     98         printf("
    ");*/
     99         for(int i=2*n;i>=1;--i){
    100             if(sa[i]<n){p1=sa[i];break;}
    101         }
    102         for(int i=0;i<n;++i){
    103             ans1[i]=s1[p1+i];
    104 //            printf("%d ",ans1[i]);
    105         }
    106 
    107 /*        printf("
    ");
    108         printf("
    ");
    109 */
    110         int cc=check(n);
    111 //        printf("c%d
    ",cc);
    112         if(cc==-1){
    113             int i,j;
    114             pp0[0]=pp0[1]=0;
    115             for(i=1;i<n;++i){
    116                 j=pp0[i];
    117                 while(j&&ans0[i]!=ans0[j])j=pp0[j];
    118                 pp0[i+1]=ans0[i]==ans0[j]?j+1:0;
    119             }
    120             j=0;
    121             for(int i=0;i<2*n;++i){
    122                 while(j&&s0[i]!=ans0[j])j=pp0[j];
    123                 if(s0[i]==ans0[j])j++;
    124                 if(j==n){
    125                     p0=i-n+1;
    126                     break;
    127                 }
    128             }
    129             pp1[0]=pp1[1]=0;
    130             for(i=1;i<n;++i){
    131                 j=pp1[i];
    132                 while(j&&ans1[i]!=ans1[j])j=pp1[j];
    133                 pp1[i+1]=ans1[i]==ans1[j]?j+1:0;
    134             }
    135             j=0;
    136             for(int i=0;i<2*n;++i){
    137                 while(j&&s1[i]!=ans1[j])j=pp1[j];
    138                 if(s1[i]==ans1[j])j++;
    139                 if(j==n){
    140                     if(i-n+1<n)p1=i-n+1;
    141                 }
    142             }
    143             p1=n-(p1+1);
    144         //    printf("p0:%d p1:%d
    ",p0,p1);
    145             if(p0<p1){
    146                 printf("%d 0
    ",p0+1);
    147             }
    148             else if(p0>p1){
    149                 printf("%d 1
    ",p1+1);
    150             }
    151             else printf("%d 0
    ",p0+1);
    152 
    153         }
    154         else if(cc==0){
    155             int i,j;
    156             pp0[0]=pp0[1]=0;
    157             for(i=1;i<n;++i){
    158                 j=pp0[i];
    159                 while(j&&ans0[i]!=ans0[j])j=pp0[j];
    160                 pp0[i+1]=ans0[i]==ans0[j]?j+1:0;
    161             }
    162             j=0;
    163             for(int i=0;i<2*n;++i){
    164                 while(j&&s0[i]!=ans0[j])j=pp0[j];
    165                 if(s0[i]==ans0[j])j++;
    166                 if(j==n){
    167                     p0=i-n+1;
    168                     break;
    169                 }
    170             }
    171             printf("%d 0
    ",p0+1);
    172         }
    173         else if(cc==1){
    174             int i,j;
    175             pp1[0]=pp1[1]=0;
    176             for(i=1;i<n;++i){
    177                 j=pp1[i];
    178                 while(j&&ans1[i]!=ans1[j])j=pp1[j];
    179                 pp1[i+1]=ans1[i]==ans1[j]?j+1:0;
    180             }
    181             j=0;
    182             for(int i=0;i<2*n;++i){
    183                 while(j&&s1[i]!=ans1[j])j=pp1[j];
    184                 if(s1[i]==ans1[j])j++;
    185                 if(j==n){
    186                     if(i-n+1<n)p1=i-n+1;
    187             //        p0=i-n+1;
    188                 }
    189             }
    190             p1=n-(p1+1);
    191             printf("%d 1
    ",p1+1);
    192         }
    193     }
    194     return 0;
    195 }
    View Code
  • 相关阅读:
    接口:
    抽象类:
    构造方法与重载
    get方法和set方法
    类的成员变量与方法、构造方法
    构造方法
    Java面向对象学习-----类的成员变量2
    Java面向对象学习-----类的成员变量
    java字符串分割
    android 使用讯飞人脸识别api报错:java.lang.UnsatisfiedLinkError
  • 原文地址:https://www.cnblogs.com/cenariusxz/p/4808079.html
Copyright © 2011-2022 走看看