(http://acm.hust.edu.cn/vjudge/contest/view.action?cid=105904#problem/A)
题意:给定字符串s,问s循环后的最小串。即s为“BBACD”,答案为“ACDBB”。
解法:最小表示法直接就能解,为了写板子,又用后缀自动机写了一遍。
最小表示法:
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <complex>
//#include <tr1/unordered_set>
//#include <tr1/unordered_map>
#include <bitset>
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define inf (1e18)
#define debug(a) cout << #a" = " << (a) << endl;
#define debugarry(a, n) for (int i = 0; i < (n); i++) { cout << #a"[" << i << "] = " << (a)[i] << endl; }
#define clr(x, y) memset(x, y, sizeof x)
#define LL long long
#define ll long long
#define uLL unsigned LL
using namespace std;
const int maxn = 10000 * 2 + 30;
int Find(char *s){
int n=strlen(s),p=0,q=1;
while(p<n&&q<n){
int i;
for(i=0;i<n;i++)
if( s[(p+i)%n] != s[(q+i)%n] )
break;
if( s[(p+i)%n] > s[(q+i)%n] )
p=p+i+1;
else
q=q+i+1;
if(p==q) q++;
}
return min(p,q);
}
char s[maxn];
int main(){
// freopen("input.txt","r",stdin);
int T;
scanf("%d",&T);
while(T--){
scanf("%s",s);
printf("%d\n",Find(s)+1);
}
return 0;
}
后缀自动机:
其实就是把s复制一遍,然后直接从头往后遍历字典序小最小的那个就好啦。
有一个问题是,因为题目要求输出为之最靠前的答案,我的做法是取出这个串之后,在做一遍kmp,麻烦一点。有的人直接用找的的那个位置的step值作为答案,没理解为什么这样是对的。。。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <queue>
//#include <tr1/unordered_set>
//#include <tr1/unordered_map>
#include <bitset>
//#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define inf 1e9
#define debug(a) cout << #a" = " << (a) << endl;
#define debugarry(a, n) for (int i = 0; i < (n); i++) { cout << #a"[" << i << "] = " << (a)[i] << endl; }
#define clr(x, y) memset(x, y, sizeof x)
#define ll long long
#define ull unsigned long long
#define FOR(i,a,b) \
for(i=a;a<b?i<=b:i>=b;a<b?i++:i--)
const int maxn = 10000*5+100;
struct suffix_automaton{
int son[maxn][27],pre[maxn],step[maxn];
int last,si;
void pushback(int v)
{
step[++si]=v;
pre[si]=-1;
clr(son[si],-1);
}
void Extend(int ch)
{
pushback(step[last]+1);
int p=last,np=si;
while( p!=-1 && son[p][ch]==-1 )
son[p][ch]=np,
p=pre[p];
if(p==-1) pre[np]=0;
else{
int q=son[p][ch];
if(step[q]==step[p]+1) pre[np]=q;
else{
pushback(step[p]+1);
int nq=si;
memcpy(son[nq],son[q],sizeof son[q]);
pre[nq]=pre[q];
pre[q]=pre[np]=nq;
while(p!=-1&&son[p][ch]==q)
son[p][ch]=nq,
p=pre[p];
}
}
last=np;
}
void init()
{
si=last=0;
clr(son[0],-1);
pre[0]=-1;
step[0]=0;
}
int le[maxn];
void cca()
{
for(int i=si;i>=0;i--)
{
le[i]=0;
for(int j=0;j<26;j++)
if(son[i][j]!=-1)
le[i] = max(le[i],le[son[i][j]]+1);
}
}
int Find(int p,int len,int cnt)
{
if( cnt == len ) return step[p];
for(int i=0;i<26;i++)
if( son[p][i]!=-1&&le[son[p][i]]+cnt>=len-1 )
return Find(son[p][i],len,cnt+1);
}
}su;
void getFail(char *P,int *f)
{
int m=strlen(P);
f[0]=0;
f[1]=0;
for(int i=1;i<m;i++)
{
int j=f[i];
while(j&&P[i]!=P[j]) j=f[j];
f[i+1]=P[i]==P[j]?j+1:0;
}
}
int Find(char *T,char *P,int *f)
{
int n=strlen(T),m=strlen(P);
getFail(P,f);
int j=0;
for(int i=0;i<n;i++){
while(j&&P[j]!=T[i]) j=f[j];
if(P[j]==T[i]) j++;
if(j==m) return i-m+1;
}
return -1;
}
char s[maxn],ss[maxn];
int f[maxn];
int main()
{
//freopen("input.txt","r",stdin);
int T;
scanf("%d",&T);
while(T--)
{
scanf("%s",s);
int len=strlen(s);
for(int i=0;i<len;i++)
s[i+len]=s[i];
s[2*len]=0;
su.init();
//debug(s);
for(int i=0;i<len*2;i++)
su.Extend(s[i]-'a');
su.cca();
int ans=su.Find(0,len,0)-len;
printf("%d\n",ans+1);
}
return 0;
}