【hdu 4333】Revolving Digits
【链接】http://acm.hdu.edu.cn/showproblem.php?pid=4333
【题意】
就是给你一个数字,然后把最后一个数字放到最前面去,经过几次变换后又回到原数字,问在这些数字中,比原数字小的,相等的,大的分别有多少个。比如341-->134-->413-->341,所以和原数字相比,比原数字小的有一个,相等的有一个,大的有一个,要去重
【题解】
把原数字在右边复制一份,合成一个长度为2len的字符串T,然后做一下扩展KMP.算出T的每个位置的extend[i]值;
在extend[1..len]中,如果extend[i]>=len,就表明从这个位置开始的长度为len的字符串和原串是一样的,否则,比较
T[extend[i]]和T[i+extend[i]]的大小就能比较出大小了.
(也即从高位开始比较);
这样的复制字符的方式能够得到所有的len个数字.
但是这样还没有去掉重复的情况;
比如
1212
复制成了12121212
取的时候会取出以下4个字符串
1212
2121 1212 2121
也即每个数字会重复出现两次;
这个重复出现的次数,和原字符的最小循环节出现的次数对应;
也即len/(len-Next[len]) 当len%(len-Next[len])==0时
否则每个数字只出现一次。
(也即没有循环子串)
最小循环节的讲解在这KMP求最小循环节
扩展KMP算法的Next数组和原KMP算法的Next数组是不一样的。
【错的次数】
0
【反思】
扩展KMP的应用
【代码】
#include <bits/stdc++.h> using namespace std; #define lson l,m,rt<<1 #define rson m+1,r,rt<<1|1 #define LL long long #define rep1(i,a,b) for (int i = a;i <= b;i++) #define rep2(i,a,b) for (int i = a;i >= b;i--) #define mp make_pair #define pb push_back #define fi first #define se second #define ms(x,y) memset(x,y,sizeof x) #define ri(x) scanf("%d",&x) #define rl(x) scanf("%lld",&x) #define rs(x) scanf("%s",x) #define oi(x) printf("%d",x) #define ol(x) printf("%lld",x) #define oc putchar(' ') #define os(x) printf(x) #define all(x) x.begin(),x.end() #define Open() freopen("F:\\rush.txt","r",stdin) #define Close() ios::sync_with_stdio(0) typedef pair<int,int> pii; typedef pair<LL,LL> pll; const int dx[9] = {0,1,-1,0,0,-1,-1,1,1}; const int dy[9] = {0,0,0,-1,1,-1,1,-1,1}; const double pi = acos(-1.0); const int N=20e4; int Next[N+10],extend[N+10],lens,lent,f[N]; char S[N+10],T[N+10]; void makenext(int m){ int a = 0; Next[0] = lens; while(a < lens - 1 && S[a] == S[a + 1]) a++; Next[1] = a; a = 1; for(int k = 2; k < lens; k ++) { int p = a + Next[a] - 1,L = Next[k - a]; if( (k - 1) + L >= p) { int j = (p - k + 1) > 0 ? (p - k + 1) : 0; while(k + j < lens && S[k + j] == S[j]) j++; Next[k] = j; a = k; } else Next[k] = L; } } void GetNext(const char *T){ int a=0; int MinLen = lens < lent ? lens : lent; while(a < MinLen && S[a] == T[a] ) a++; extend[0]=a; a=0; for(int k=1;k < lent;k++){ int p=a+extend[a]-1,L = Next[k-a]; if((k-1)+L>=p){ int j=(p-k+1)>0? (p-k+1):0; while(k + j < lent && T[k+j] == S[j]) j++; extend[k]=j; a=k; } else extend[k]=L; } } void getf(const char *s){ f[0] = 0,f[1] = 0; int j = 0; rep1(i,1,lens-1){ while (j > 0 && s[i]!=s[j]) j = f[j]; f[i+1] = (s[i]==s[j]?++j:0); } } int main(){ //Open(); int TT,kk = 0; ri(TT); while (TT--){ rs(S); lens = strlen(S),lent = 2*lens; makenext(lens); getf(S); int ju = lens-f[lens],tt; if (lens%ju==0) tt = lens/ju; else tt = 1; rep1(i,0,lens-1) T[i] = S[i]; rep1(i,lens,2*lens-1) T[i] = T[i-lens]; GetNext(T); int a = 0,b = 0,c = 0; rep1(i,0,lens-1) if (extend[i]>=lens) b++; else if (S[extend[i]] < T[i+extend[i]]) a++; else if (S[extend[i]] > T[i+extend[i]]) c++; os("Case ");oi(++kk);os(": ");oi(c/tt);oc;oi(b/tt);oc;oi(a/tt); puts(""); } return 0; }