POJ 题目3581 Sequence（后缀数组+离散化）

最新推荐文章于 2020-10-10 18:46:59 发布

Jogging_Clown

最新推荐文章于 2020-10-10 18:46:59 发布

阅读量623

点赞数 1

分类专栏：后缀数组

本文链接：https://blog.csdn.net/yu_ch_sh/article/details/48056369

版权

后缀数组专栏收录该内容

25 篇文章 0 订阅

订阅专栏

Sequence

Time Limit: 5000MS		Memory Limit: 65536K
Total Submissions: 5738		Accepted: 1227
Case Time Limit: 2000MS

Description

Given a sequence, {A₁, A₂, ..., A_n} which is guaranteedA₁> A₂, ..., A_n, you are to cut it into three sub-sequences and reverse them separately to form a new one which is the smallest possible sequence in alphabet order.

The alphabet order is defined as follows: for two sequence {A₁,A₂, ..., A_n} and {B₁,B₂, ..., B_n}, we say {A₁,A₂, ..., A_n} is smaller than {B₁,B₂, ..., B_n} if and only if there exists suchi ( 1 ≤ i ≤ n) so that we have A_i <B_{i and}A_j = B_j for each j <i.

Input

The first line contains n. (n ≤ 200000)

The following n lines contain the sequence.

Output

output n lines which is the smallest possible sequence obtained.

Sample Input

Sample Output

Hint

{10, 1, 2, 3, 4} -> {10, 1 | 2 | 3, 4} -> {1, 10, 2, 4, 3}

Source

POJ Founder Monthly Contest – 2008.04.13, Yao Jinyu

解析转自：http://blog.csdn.net/whyorwhnt/article/details/9840641

题意：给定一个数列，将其分为非空的三段。然后分别将三段的数字翻转连接在一起组成新的数列。输出处理后字典序最小的结果。

思路：首先注意必须写成单Case的形式！！！while(scanf("%d",&n)!=EOF)会RE或者WA。。。

因为要翻转，所以读入时直接反着读。由于没有告诉数值的上限，所以接下来离散化。

对于第一段：第一段翻转之后必须最小，因此对翻转过的数列求sa，则排名在前面且sa[i]>1(因为另两段不能为空)的第一个i为第一段。

接下来处理后两段：显然要让第二段的开头尽可能小，但是切去第二段后剩下的部分会直接接在第二段后面，所以不能直接求sa。为保证第三段接在第二段后面时整体最小，所以需要把剩余的串复制到剩余串的后面，再求sa.注意数组要开2倍。

以下举例转自 poj3581 Sequence(后缀数组) - Damonbaby的日志 - 网易博客

刚开始的时候觉得没有必要将剩下的数字复制一片贴在剩下的数字后面，所以wa了很多次。。。
最后在discuss里面看见一组数据：
9
8 4 -1 5 0 5 0 2 3
第一步：
3 2 0 5 0 5 -1 4 8 对应输出 -1 4 8
第二步
3 2 0 5 0 5（开始的时候我并没有复制一遍）对应输出：0 5
第三步
3 2 0 5 对应输出： 3 2 0 5
可以看见这样做是不对的。。
必须要将剩下的字符串复制一遍贴在后面，然后再来求后缀数组。。。
正解：
第一步：
3 2 0 5 0 5 -1 4 8 对应输出 -1 4 8
第二步
3 2 0 5 0 5 3 2 0 5 0 5 对应输出： 0 5 0 5；
第三步
3 2 对应输出：3 2；

ac代码

Problem: 3581  User: kxh1995 
Memory: 8020K  Time: 2032MS 
Language: C++  Result: Accepted

#include<stdio.h>                 
#include<string.h>                 
#include<algorithm> 
#include<stdlib.h>                
#include<iostream>                
#define min(a,b) (a>b?b:a)             
#define max(a,b) (a>b?a:b)     
#define INF 0xfffffff             
using namespace std;                     
int sa[400010],Rank[400010],rank2[400010],c[400010],*x,*y,s[400010],k,temp[400010],num[400010];     
void cmp(int n,int sz)          
{          
    int i;          
    memset(c,0,sizeof(c));          
    for(i=0;i<n;i++)          
        c[x[y[i]]]++;          
    for(i=1;i<sz;i++)          
        c[i]+=c[i-1];          
    for(i=n-1;i>=0;i--)          
        sa[--c[x[y[i]]]]=y[i];          
}          
void build_sa(int *s,int n,int sz)          
{          
    x=Rank,y=rank2;          
    int i,j;          
    for(i=0;i<n;i++)          
        x[i]=s[i],y[i]=i;          
    cmp(n,sz);          
    int len;          
    for(len=1;len<n;len<<=1)          
    {          
        int yid=0;          
        for(i=n-len;i<n;i++)          
        {          
            y[yid++]=i;          
        }          
        for(i=0;i<n;i++)          
            if(sa[i]>=len)          
                y[yid++]=sa[i]-len;          
            cmp(n,sz);          
			swap(x,y);          
			x[sa[0]]=yid=0;          
			for(i=1;i<n;i++)          
			{          
				if(y[sa[i-1]]==y[sa[i]]&&sa[i-1]+len<n&&sa[i]+len<n&&y[sa[i-1]+len]==y[sa[i]+len])          
					x[sa[i]]=yid;          
				else         
					x[sa[i]]=++yid;          
			}          
			sz=yid+1;          
			if(sz>=n)          
				break;          
    }          
    for(i=0;i<n;i++)          
        Rank[i]=x[i];          
}          
struct s
{
	int val,id;
}a[200020];
int cmp(const void *a,const void *b)
{
	return (*(struct s *)a).val-(*(struct s *)b).val;
}
int main()
{
	int n;
	scanf("%d",&n);
		int i;
		for(i=n-1;i>=0;i--)
		{
			scanf("%d",&num[i]);
			a[i].val=num[i];
			a[i].id=i;
		}
		qsort(a,n,sizeof(a[0]),cmp);
		temp[0]=1;
		for(i=1;i<n;i++)
		{
			temp[i]=(a[i].val==a[i-1].val?temp[i-1]:i+1);
		}
		for(i=0;i<n;i++)
		{
			s[a[i].id]=temp[i];
		}
		s[n]=0;
		build_sa(s,n+1,n+5);
		int minn=INF;
		for(i=1;i<=n;i++)
		{
			if(sa[i]>1)
				break;
		}
		int m=sa[i];
		for(i=m;i<n;i++)
		{
			printf("%d\n",num[i]);
		}
		for(i=0;i<m;i++)
		{
			s[m+i]=s[i];
		}
		
		int nn=m*2;
		s[nn]=0;
		//for(i=0;i<nn;i++)
		//	{
		//		printf("****%d\n",s[i]);
		//	}
		build_sa(s,nn+1,n+5);
		for(i=1;i<=nn;i++)
		{
			if(sa[i]>0&&sa[i]<m)
				break;
		}
		int j;
		for(j=sa[i];j<m;j++)
		{
			printf("%d\n",num[j]);
		}
		for(j=0;j<sa[i];j++)
		{
			printf("%d\n",num[j]);
		}
	//	while(1);
}