hdu5442(2015长春网络赛F题)

最新推荐文章于 2018-08-19 11:17:45 发布
最新推荐文章于 2018-08-19 11:17:45 发布
阅读量435 收藏
点赞数
分类专栏： # hdu代码 # 模板 # 字符串 # 后缀数组 文章标签： hdu 网络赛 长春 后缀数组 字符串
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ms961516792/article/details/48524167
版权
hdu代码 同时被 3 个专栏收录
35 篇文章 0 订阅
订阅专栏
模板
13 篇文章 0 订阅
订阅专栏
字符串
7 篇文章 0 订阅
订阅专栏

题意：

给出一个字符串，只由‘a’~‘z'组成，字符串是一个首尾相接的串。我们要找到一个起点，顺时针或者逆时针的读这个串，找到字典序最大的读法，如果有多种，输出起点坐标小的那个，如果起点坐标一样，输出顺时针的那个。


思路：

把串变成倍长的，反着再处理一次。

坑死我了。。。用罗穗骞的da算法超时，用dc3RE，用了个网上找的，过了。


代码：

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <set>
#include <algorithm>
#include <math.h>
using namespace std;
typedef long long ll;
const int N=1e5;
int sa[N],height[N],Rank[N],tmp[N],top[N];
char s[N];
void dc3(int n){
    int i,j,len,na;
    na=(n<256?256:n);
    memset(top,0,na*sizeof(int));
    for(i=0;i<n;i++)top[Rank[i]=s[i]&0xff]++;
    for(i=1;i<na;i++)top[i]+=top[i-1];
    for(i=0;i<n;i++)sa[--top[Rank[i]]]=i;
    for(len=1;len<n;len<<=1){
        for(i=0;i<n;i++){
            j=sa[i]-len;if(j<0)j+=n;
            tmp[top[Rank[j]]++]=j;
        }
        sa[tmp[top[0]=0]]=j=0;
        for(i=1;i<n;i++){
            if(Rank[tmp[i]]!=Rank[tmp[i-1]]||
                Rank[tmp[i]+len]!=Rank[tmp[i-1]+len])
                top[++j]=i;
            sa[tmp[i]]=j;
        }
        memcpy(Rank,sa,n*sizeof(int));
        memcpy(sa,tmp,n*sizeof(int));
        if(j>=n-1)break;
    }
}

char rt[N];

void _print(int ans1,int ans2,int n)
{
    int t1=ans1;
    int t2=ans2;
    t2+=n;
    for(int i=0;i<n;i++){
        if(rt[t1]>rt[t2])
        {
            printf("%d 0\n",ans1+1);
            return ;
        }
        if(rt[t1]<rt[t2])    
        {
            printf("%d 1\n",ans2+1);
            return ;
        }
        t1++;
        if(t1>n+n)
            break;
        t2--;
        if(t2<0)
            break;
    }
    if(ans1<ans2)
        printf("%d 0\n",ans1+1);
    else if(ans2<ans1)
        printf("%d 1\n",ans2+1);
    else
        printf("%d 0\n",ans1+1);
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&n);
        scanf("%s",s);
        for(int i=0;i<n;i++)
        {
            s[i+n]=s[i];
            rt[i]=rt[i+n]=s[i];
        }
        s[n+n]=0;
        dc3(n+n+1);
//        for(int i=0;i<n+n;i++)
//            printf("1--%d\n",Rank[i]);
        int ans1=0;
        for(int i=1;i<n;i++)
        {
            if(Rank[i]>Rank[ans1])
                ans1=i;
        }
        for(int i=0;i<n;i++)
        {
            s[i]=rt[n-i-1];
            s[i+n]=rt[n-i-1];
        }
        s[n+n]='z'+1;
        dc3(n+n+1);
        int ans2=0;
        for(int i=1;i<n;i++)
        {
            if(Rank[i]>Rank[ans2])
                ans2=i;
        }
        ans2=n-1-ans2;
        _print(ans1,ans2,n);
//        printf("%d %d\n",ans1,ans2);
//        for(int i=0;i<n+n;i++)
//            printf("2--%d\n",Rank[i]);
    }
    return 0;
}


附上TLE的da算法代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

#define maxn 20005
int wa[maxn],wb[maxn],wv[maxn],_ws[maxn];
int cmp(int *r,int a,int b,int l)
{return r[a]==r[b]&&r[a+l]==r[b+l];}
void da(char *r,int *sa,int n,int m)
{
     int i,j,p,*x=wa,*y=wb,*t;
     for(i=0;i<m;i++) _ws[i]=0;
     for(i=0;i<n;i++) _ws[x[i]=r[i]]++;
     for(i=1;i<m;i++) _ws[i]+=_ws[i-1];
     for(i=n-1;i>=0;i--) sa[--_ws[x[i]]]=i;
     for(j=1,p=1;p<n;j*=2,m=p)
     {
       for(p=0,i=n-j;i<n;i++) y[p++]=i;
       for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
       for(i=0;i<n;i++) wv[i]=x[y[i]];
       for(i=0;i<m;i++) _ws[i]=0;
       for(i=0;i<n;i++) _ws[wv[i]]++;
       for(i=1;i<m;i++) _ws[i]+=_ws[i-1];
       for(i=n-1;i>=0;i--) sa[--_ws[wv[i]]]=y[i];
       for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
       x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
     }
     return;
}
int _rank[maxn],_height[maxn];
void cal_height(char *r,int *sa,int n)
{
     int i,j,k=0;
     for(i=1;i<=n;i++) _rank[sa[i]]=i;
     for(i=0;i<n;_height[_rank[i++]]=k)
     for(k?k--:0,j=sa[_rank[i]-1];r[i+k]==r[j+k];k++);
     return;
}

char r[maxn],rt[maxn];
int sa[maxn];

void _print(int ans1,int ans2,int n)
{
    int t1=ans1;
    int t2=ans2;
    t2+=n;
    for(int i=0;i<n;i++){
        if(rt[t1]>rt[t2])
        {
            printf("%d 0\n",ans1+1);
            return ;
        }
        if(rt[t1]<rt[t2])    
        {
            printf("%d 1\n",ans2+1);
            return ;
        }
        t1++;
        t2--;
    }
    if(ans1<ans2)
        printf("%d 0\n",ans1+1);
    else if(ans2<ans1)
        printf("%d 1\n",ans2+1);
    else
        printf("%d 0\n",ans1+1);
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&n);
        scanf("%s",r);
        for(int i=0;i<n;i++)
        {
            r[i+n]=r[i];
            rt[i]=rt[i+n]=r[i];
        }
        r[n+n]=0;
        da(r,sa,n+n+1,200);
        cal_height(r,sa,n+n);
//        for(int i=0;i<n+n;i++)
//            printf("1--%d\n",_rank[i]);
        int ans1=0;
        for(int i=1;i<n;i++)
        {
            if(_rank[i]>_rank[ans1])
                ans1=i;
        }
        for(int i=0;i<n;i++)
        {
            r[i]=rt[n-i-1];
            r[i+n]=rt[n-i-1];
        }
        r[n+n]='z'+1;
        da(r,sa,n+n+1,200);
        cal_height(r,sa,n+n);
        int ans2=0;
        for(int i=1;i<n;i++)
        {
            if(_rank[i]>_rank[ans2])
                ans2=i;
        }
        ans2=n-1-ans2;
        _print(ans1,ans2,n);
//        printf("%d %d\n",ans1,ans2);
//        for(int i=0;i<n+n;i++)
//            printf("2--%d\n",_rank[i]);
    }
    return 0;
}


附上RE的dc3算法代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

#define maxn 1000003

#define F(x) ((x)/3+((x)%3==1?0:tb))
#define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)
int wa[maxn],wb[maxn],wv[maxn],_ws[maxn];
int c0(char *r,int a,int b)
{
    return r[a]==r[b]&&r[a+1]==r[b+1]&&r[a+2]==r[b+2];
}
int c12(int k,char *r,int a,int b)
{
    if(k==2) 
        return r[a]<r[b] || r[a]==r[b] && c12(1,r,a+1,b+1);
    else 
        return r[a]<r[b] || r[a]==r[b] && wv[a+1]<wv[b+1];
}

void sort(char *r,int *a,int *b,int n,int m)
{
     int i;
     for(i=0;i<n;i++) wv[i]=r[a[i]];
     for(i=0;i<m;i++) _ws[i]=0;
     for(i=0;i<n;i++) _ws[wv[i]]++;
     for(i=1;i<m;i++) _ws[i]+=_ws[i-1];
     for(i=n-1;i>=0;i--) b[--_ws[wv[i]]]=a[i];
     return;
}
void dc3(char *r,int *sa,int n,int m)
{
     int i,j,*san=sa+n,ta=0,tb=(n+1)/3,tbc=0,p;
     char *rn=r+n;
     r[n]=r[n+1]=0;
     for(i=0;i<n;i++) if(i%3!=0) wa[tbc++]=i;
     sort(r+2,wa,wb,tbc,m);
     sort(r+1,wb,wa,tbc,m);
     sort(r,wa,wb,tbc,m);
     for(p=1,rn[F(wb[0])]=0,i=1;i<tbc;i++)
     rn[F(wb[i])]=c0(r,wb[i-1],wb[i])?p-1:p++;
     if(p<tbc) dc3(rn,san,tbc,p);
     else for(i=0;i<tbc;i++) san[rn[i]]=i;
     for(i=0;i<tbc;i++) if(san[i]<tb) wb[ta++]=san[i]*3;
     if(n%3==1) wb[ta++]=n-1;
     sort(r,wb,wa,ta,m);
     for(i=0;i<tbc;i++) wv[wb[i]=G(san[i])]=i;
     for(i=0,j=0,p=0;i<ta && j<tbc;p++)
     sa[p]=c12(wb[j]%3,r,wa[i],wb[j])?wa[i++]:wb[j++];
     for(;i<ta;p++) sa[p]=wa[i++];
     for(;j<tbc;p++) sa[p]=wb[j++];
     return;
}

int _rank[maxn],_height[maxn];
void cal_height(char *r,int *sa,int n)
{
     int i;
     for(i=1;i<=n;i++) _rank[sa[i]]=i;
     return;
}

char r[maxn],rt[maxn];
int sa[maxn];

void _print(int ans1,int ans2,int n)
{
    int t1=ans1;
    int t2=ans2;
    t2+=n;
    for(int i=0;i<n;i++){
        if(rt[t1]>rt[t2])
        {
            printf("%d 0\n",ans1+1);
            return ;
        }
        if(rt[t1]<rt[t2])    
        {
            printf("%d 1\n",ans2+1);
            return ;
        }
        t1++;
        t2--;
    }
    if(ans1<ans2)
        printf("%d 0\n",ans1+1);
    else if(ans2<ans1)
        printf("%d 1\n",ans2+1);
    else
        printf("%d 0\n",ans1+1);
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&n);
        scanf("%s",r);
        for(int i=0;i<n;i++)
        {
            r[i+n]=r[i];
            rt[i]=rt[i+n]=r[i];
        }
        r[n+n]=0;
        dc3(r,sa,n+n+1,200);
        cal_height(r,sa,n+n);
//        for(int i=0;i<n+n;i++)
//            printf("1--%d\n",_rank[i]);
        int ans1=0;
        for(int i=1;i<n;i++)
        {
            if(_rank[i]>_rank[ans1])
                ans1=i;
        }
        for(int i=0;i<n;i++)
        {
            r[i]=rt[n-i-1];
            r[i+n]=rt[n-i-1];
        }
        r[n+n]='z'+1;
        dc3(r,sa,n+n+1,200);
        cal_height(r,sa,n+n);
        int ans2=0;
        for(int i=1;i<n;i++)
        {
            if(_rank[i]>_rank[ans2])
                ans2=i;
        }
        ans2=n-1-ans2;
        _print(ans1,ans2,n);
//        printf("%d %d\n",ans1,ans2);
//        for(int i=0;i<n+n;i++)
//            printf("2--%d\n",_rank[i]);
    }
    return 0;
}





NeverMoreH
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
hdu杭电网络编程结课报告 聊天室
02-16
hdu杭电网络编程结课报告 聊天室
HDU 递归题详解大全（含代码）
02-11
蟠桃记 1 折线分割平面 2 不容易系列之一 2 骨牌铺方格 3 不容易系列之(3)—— LELE的RPG难题 3 Children’s Queue 3 献给杭电五十周年校庆的礼物 3 钥匙计数之二 3 钥匙计数之一 3 母牛的故事 3 ...
hdu5442
bigsungod的专栏
09-18 534
题意： 有一个len长度的环，求某一位置pos，以该点为头，其字典序最大（正逆顺序都行）。 如果有且只有一个 ，输出其 开头的下标， 下标从1 开始， 再输出0 表示顺时针  1表示逆时针 如果有多个（包括在正序与逆序），输出 开头下标最小的那个。 另外，就是顺时针 逆时针 一样且下标相等的情况  算成顺时针的。 思路：字符串的最大表示法 + kmp，kmp针对逆序的 因为跑完最大表
HDU5442（字符串的最大表示法或者后缀数组）
DS00HY的小博客
08-19 258
题目网址 https://cn.vjudge.net/problem/HDU-5442  字符串的最大表示法 https://wenku.baidu.com/view/b0ef1be7a6c30c2258019ede.html ，这个博客里面代码https://blog.csdn.net/zy691357966/article/details/39854359 字符串的最大表示： 题解方法转...
hdu 5442 Favorite Donut（最小表示法+kmp）
不慌不忙、不急不躁
09-13 869
题目链接：hdu 5442 Favorite Donut 正序逆序分别做一遍最小表示法，但是因为逆序求出最小表示法下标对应的是为最大值，所以我们用kmp求出循环节长度，然后用下标取模即可。 #include #include #include #include #include #include #include using namespace std; con
hdu5442Favorite Donut(好题)
_zidaoziyan的博客
09-22 559
Favorite Donut Time Limit: 1500/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others) Total Submission(s): 1961    Accepted Submission(s): 482 Problem Description Lulu has a s
hdu5442 后缀数组
zzucaicai的专栏
09-18 343
第一道AC的后缀数组，总算也是深入理解了各个数组的意义 height数组连续>＝length，就说明两个顺序相同，保证字母序最大的前提下，先开始位置小 再顺时针优先于逆时针 #include #include #include using namespace std; const int N = int(2e5)+10; inline int max(int a ,int b){
HDU杭电 计算机网络实验报告
03-11
压缩包包含十份报告，已经通过验收，实验内容：交换机、生成树、静态路由、NAT等完全根据教材实验要求
杭电HDU2050题的ac程序
04-10
一个十分简单的程序，能够ac杭电hdu的第2050题，无注释，简单明了
【HDU5442】 Favorite Donut（后缀数组）
u011332631的专栏
09-15 688
题意，给你一个长度为n的字符串，让你选择从某一个地方断开（可以正着取也可以反着取）。问你最大的字典序是从哪里断开。 首先先将原串复制一遍，然后用后缀数组求出字典序最大的位置，如果字典序相同则求出最靠前的位置。 然后再将该串翻转，得到反向取字典序最大的串，再直接比较就好了。 code： #include #include #include using namespace std; const
hdu 5442 (后缀数组)
diaocuiguo2493的博客
09-20 86
稍微学习了下第一次用后缀数组- - ， 强行凑出答案 ， 感觉现在最大的问题是很多算法都不知道 ，导致有的题一点头绪都没有（就像本题）。 /*推荐 《后缀数组——处理字符串的有力工具》——罗穗骞 后缀数组sa， sa[ i ] = a表示字符串从第a个开始到结尾的字典序排序为i 本题是绕成了一个环，所以我将字符串重复一遍再用后缀数组模板跑一次。 顺序的话，取最后一个即可...
HDU 5442 (串的最大表示+KMP）
84970000的专栏
09-13 548
题意：一个环形的字符串长度为n，求出长度n的字典序最大的子串，如果多个符合的子串，取顺时针起始字符下标最小的子串，如果从该起始字符的顺时针和逆时针得到结果相同，取顺时针的那个。 做法：这道题据说可以用后缀数组、后缀自动机之类的重型兵器解决，也可以用字符串的最大表示法加KMP这种相对轻量级的工具来解决。字符串的最大表示可以参考http://blog.csdn.net/zy691357966/art
hdu 5442 Favorite Donut 最小表示法+KMP 2015长春网络赛
万事屋
09-13 1672
题意： 有一个len长度的环，问有没有 最大的 长度为 len 的 串在 这个环里。 如果有的话，且只有一个 ，输出其 开头的下标， 下标从1 开始， 再输出0 表示 顺时针 1 表示逆时针 如果多个，输出 开头下标最小的那个。 如果 还有 ，就是顺时针 逆时针 一样的情况 输出 下标 0。 做法： 很容易想到最小表示法 然后 正序直接用最大表示法 ，可以得到 最大字典序，最小下标的 开头位置。 逆序 就把串反下，用最大表示法，可以得到 最大字典序， 也是下标最小的 开头位置， 但是因
2013长春现场赛
菜鸟成长ing
12-07 769
2013长春现场赛水题 虽然我想不出思路看别人题解懂的  #include #include #include #include using namespace std; const int maxn = 105; const int inf = 1<<30; int n,m,pos; __int64 s[maxn],p[maxn],price[maxn]; int find( __i
HDU 5442 Favorite Donut（最大表示法）
aozil_yang的博客
08-03 477
Favorite Donut Time Limit: 1500/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others) Total Submission(s): 2746    Accepted Submission(s): 667 Problem Description Lulu has a swe
hdu 5442
u013570474的专栏
09-28 360
最小表示法的推广 -》 最大表示 分别求 正序跟逆序的，但是由于这道题要起点最小，那么逆序的时候不能找到第一个就返回。要继续往后找#include <iostream> #include <string.h> #include <stdlib.h> #include <stdio.h> #include <map>using namespace std; #define MAXN 20010 ch
hdu5442 Favorite Donut（最大表示法+KMP）
Code of Duty
09-14 582
头一次比赛中遇到最大(最小)表示法的题目 因为这个算法可以把字典序最大(最小)的串的最小的起始位置求出来，所以以为正着跑一遍  反着跑一遍就OK了。 其实反着跑的时候我是把字符串翻转的，所以求出来的其实是最大的起始位置 -  - WA了N次  还以为是题意理解错了   -  - #include using namespace std; int n; int f(char wo
HDU5442 最小（大）表示法
Since_natural_ran
12-02 316
最小（大）表示法​ 一个字符串扩增到两倍的时候，在0~n-1的范围中都可以通过循环找到与原字符串同构的字符串，但是其中有两个特殊的字符串，字典序最大和字典序最小（数字也可以当做字典序来处理），直接暴力找到最小或者最大的算法复杂度是O(n^2) ，有一种最小（大）表示法可以解决此问题。百度论文：https://wenku.baidu.com/view/c6c5e7335a8102d276a22fa
HDU5442 Favorite Donut（KMP+最大表示法）
hcx11333的博客
08-17 191
题意：给一个环，要找出里面的最大的串，输出串最大时起点的下标（从1开始）和顺时针（0）还是逆时针（1）。如果有多种情况就输出最小的下标，如果顺时针和逆时针相同就选择顺时针。 思路：正逆各用一次最大表示法即可。对于逆时针，因为要考虑字符串的循环对下标的影响，因此在与顺时针方向比较之前先做一次KMP，求出循环节长度。 #include #include #include #incl
hdu-1569(网络流)-最大点权独立集
最新发布
05-30
这是一道经典的网络流问题，可以使用最小割模型进行建模。 首先，将所有点拆成两个点，分别表示该点被选中和不被选中两种状态。对于每个点，如果它被选中，则将它的点权加入最终结果中。接下来，对于每条边 $(u,v)$...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值