【SAM】BZOJ2882-工艺

【题目大意】

求一个循环数列的最小表示法。

【思路】

最小表示法的正解：★

SAM乱搞，和前面的POJ那道一样。然而MLE了，当作学习一下map的用法^ ^

map的使用方法（来源：☆）

一、map的说明  
  1   头文件 
  #include   <map> 
  
  2   定义 
  map<string,   int>   my_Map; 
  或者是typedef     map<string,   int>   MY_MAP; 
  MY_MAP   my_Map; 
  
  3   插入数据 
  (1)   my_Map["a"]   =   1; 
  (2)   my_Map.insert(map<string,   int>::value_type("b",2)); 
  (3)   my_Map.insert(pair<string,int>("c",3)); 
  (4)   my_Map.insert(make_pair<string,int>("d",4)); 
  
  4   查找数据和修改数据 
  (1)   int   i   =   my_Map["a"]; 
            my_Map["a"]   =   i; 
  (2)   MY_MAP::iterator   my_Itr; 
            my_Itr.find("b"); 
            int   j   =   my_Itr->second; 
            my_Itr->second   =   j; 
  不过注意，键本身是不能被修改的，除非删除。 
  
  5   删除数据 
  (1)   my_Map.erase(my_Itr); 
  (2)   my_Map.erase("c"); 
  还是注意，第一种情况在迭代期间是不能被删除的，道理和foreach时不能删除元素一样。 
  
  6   迭代数据 
  for   (my_Itr=my_Map.begin();   my_Itr!=my_Map.end();   ++my_Itr)   {} 
  
  7   其它方法 
  my_Map.size()               返回元素数目 
  my_Map.empty()       判断是否为空 
  my_Map.clear()           清空所有元素 
  可以直接进行赋值和比较：=,   >,   >=,   <,   <=,   !=   等等 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
const int MAXN=300002;
int n,k;
struct SAM
{
    int tot,last;
    int step[MAXN<<2],pre[MAXN<<2];
    map<int,int> next[MAXN<<2];
    inline int newNode(int cnt)
    {
        step[++tot]=cnt;
        pre[tot]=0;
        next[tot].clear();
        return tot;
    }
    
    inline void extend(int x)
    {
        int p=last;
        int np=newNode(step[last]+1);
        map<int,int>::iterator pos;
        while (p) 
        {
            pos=next[p].find(x);
            if (pos!=next[p].end()) break;
            next[p].insert(pair<int,int>(x,np));
            p=pre[p];
        }
        if (!p) pre[np]=1;
            else 
            {
                int q=pos->second;
                if (step[np]==step[p]+1) pre[np]=q;
                else
                {
                    int nq=newNode(step[p]+1);
                    next[nq]=next[q];
                    pre[nq]=pre[q];
                    pre[q]=pre[np]=nq;
                    while (p)
                    {
                        pos=next[p].find(x);
                        if (pos->second!=q) break;
                        pos->second=nq;
                        p=pre[p];
                    }
                }
            }
            last=np;
    }
    
    inline void clear()
    {
        tot=0;
        last=newNode(0);
    }
}suf;
int a[MAXN];

void init()
{
    scanf("%d",&n);
    suf.clear();
    for (int i=0;i<n;i++)
    {
        scanf("%d",&a[i]);
    }    
    for (int i=0;i<n;i++)
        suf.extend(a[i]);
    for (int i=0;i<n;i++)
        suf.extend(a[i]);
}

void solve()
{
    int j=1;
    for (int i=0;i<n;i++)
    {
        map<int,int>::iterator pos=suf.next[j].begin();
        printf("%d",pos->first);
        if (i!=n-1) printf(" ");
        j=pos->second;
    }
}

int main()
{
    init();
    solve();
    return 0;
}