POJ 1509 Glass Beads---最小表示法

题意：

　　T组数据，每组数据给出一个字符串，求这个字符串的最小表示发(只要求输出起始位置坐标)

---

　　SAM入门题(检测板子是否正确)。

　　将字符串S加倍丢进SAM中,然后走字符串长度次，每次贪心的沿最小的边走，然后答案就是:sam.e[po].len-len+1

---

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<cstring>
using namespace std;
#define maxn 10010
#define llg long long
#define SIZE 26 
#define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
llg n,m,T;

char s[maxn];

struct SAM
{
    struct 
    {
        llg len,f/*parent边*/,ch[SIZE];
        void init()
        {
            len=0,f=-1;
            memset(ch,0xff,sizeof(ch));
        }
    }e[maxn<<1];
    llg idx,last;
    
    void init() {idx=last=0; e[idx++].init();}
    
    int newnode() {e[idx].init(); return idx++;}
    
    void add(llg c)
    {
        int end=newnode(),tmp=last;
        e[end].len=e[last].len+1;
        for (;tmp!=-1 && e[tmp].ch[c]==-1;tmp=e[tmp].f){e[tmp].ch[c]=end;}//跳parent tree的边，找到第一个具有c出边的点并停下，对于没有的点连出一条边权是c的边指向当前点（end）。
        if (tmp==-1) e[end].f=0;//如果没有任何一个点有权值为c的出边，则说明了相应字符c是第一次出现在自动机中。
        else
        {
            llg nxt=e[tmp].ch[c];
            if (e[tmp].len+1==e[nxt].len) e[end].f=nxt;//如果找到的第一个具有出边c的点的出边c左指向点的len=lastlen+1，则直接把新建点的parent边连向nxt点。
            else
            {
                llg np=newnode();
                e[np]=e[nxt];
                e[np].len=e[tmp].len+1;//新建点np
                e[nxt].f=e[end].f=np;
                for (;tmp!=-1 && e[tmp].ch[c]==nxt;tmp=e[tmp].f) {e[tmp].ch[c]=np;}//沿parent边往祖先走把所有原本连向nxt的带边权c的边改为连向np
            }//如果不满足MAX(s)+1=MAX(s[tmp][c])则新建点
        }
        last=end;
    }
}sam;

int main()
{
    yyj("SAM");
    cin>>T;
    while (T--)
    {
        sam.init();
        scanf("%s",s);
        llg len=strlen(s);
        for (llg i=0;i<len;i++) sam.add(s[i]-'a');
        for (llg i=0;i<len;i++) sam.add(s[i]-'a');
        llg po=0;
        for (llg i=0;i<len;i++)
            for (llg j=0;j<26;j++)
                if (sam.e[po].ch[j]!=-1)
                {
                    po=sam.e[po].ch[j];
                    break;
                }
        printf("%lld
",sam.e[po].len-len+1);
    }
    return 0;
}