[贪心+后缀数组] poj 3623 Best Cow Line, Gold

该博客探讨了如何利用贪心算法和后缀数组在O(1)时间内解决POJ 3623问题。问题要求从给定的字符序列中,每次选择首字符或尾字符,找到能构造出的最小字典序字符串。通过比较由首字符和尾字符构建的字符串字典序,结合后缀数组的快速查询特性,可以高效地得出解决方案。
摘要由CSDN通过智能技术生成

题意:

给N个字符,每次只能取第一个或者最后一个,问构成最小字典序的串是什么。

思路:

贪心,每次取字典序最小的串,的首字符。

其实就是l和r比较一下构成的串那个字典序小。

这里运用后缀数组要实现O(1)的查询。

将原串反拼进串内做da。

然后根据ra判断。

代码:

#include"cstdlib"
#include"cstdio"
#include"cstring"
#include"cmath"
#include"queue"
#include"algorithm"
#include"iostream"
#include"map"
#define N 70003
using namespace std;
char ans[N];
int v[N],sa[N],ra[N],height[N];
int Log[N];
int wa[N],wb[N],wv[N],wws[N];
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void da(int n,int m)
{
    int i,j,p,*x=wa,*y=wb;
    for(i=0; i<m; i++) wws[i]=0;
    for(i=0; i<n; i++) wws[x[i]=v[i]]++;
    for(i=1; i<m; i++) wws[i]+=wws[i-1];
    for(i=n-1; i>=0; i--) sa[--wws[x[i]]]=i;
    for(j=1,p=1; p<n; j*=2,m=p)
    {
        for(i=n-j,p=0; 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++) wws[i]=0;
        for(i=0; i<n; i++) wws[wv[i]]++;
        for(i=1; i<m; i++) wws[i]+=wws[i-1];
        for(i=n-1; i>=0; i--) sa[--wws[wv[i]]]=y[i];
        for(swap(x,y),p=1,i=1,x[sa[0]]=0; i<n; i++) x[sa[i]]=cmp(y,sa[i],sa[i-1],j)?p-1:p++;
    }
    return ;
}
void gethei(int n)
{
    int i,j,k=0;
    for(i=1; i<=n; i++) ra[sa[i]]=i;
    for(i=0; i<n; i++)
    {
        if(k) k--;
        j=sa[ra[i]-1];
        while(v[i+k]==v[j+k]) k++;
        height[ra[i]]=k;
    }
    return ;
}
int main()
{
    int n;
    while(scanf("%d",&n)!=-1)
    {
        int len=n;
        for(int i=0;i<n;i++)
        {
            char x[12];
            scanf("%s",x);
            v[i]=x[0]-'A'+3;
        }
        for(int i=n-2;i>=0;i--) v[len++]=v[i];
        v[len]=0;
        da(len+1,40);
        gethei(len);
        int cnt=0;
        int l=0,r=n-1;
        while(cnt!=n)
        {
            if(ra[l]>ra[(n-1)*2-r])
            {
                ans[cnt++]=v[r]+'A'-3;
                r--;
            }
            else
            {
                ans[cnt++]=v[l]+'A'-3;
                l++;
            }
        }
        for(int i=0;i<cnt;i++)
        {
            printf("%c",ans[i]);
            if((i+1)%80==0) puts("");
        }
        if(cnt%80!=0) puts("");
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值