JZOJ 3648. 【GDOI2014】beyond
生活随笔
收集整理的這篇文章主要介紹了
JZOJ 3648. 【GDOI2014】beyond
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
Description
Input
第一行:包含一個整數N。
第二行:包含一個長度為N的字符串,字符串中只包含小寫字母。
第三行:包含一個長度為N的字符串,字符串中只包含小寫字母。
Output
輸出答案只包含一個數字L,表示圓環最大可能有的格子數。
Sample Input
輸入1:
5
abcdx
cdabz
輸入2:
4
abcd
cdab
Sample Output
輸出1:
4
輸出2:
4
Data Constraint
對于20% 的數據,1 <= N <= 5,000
對于50% 的數據,1 <= N <= 600,000
對于100% 的數據,1 <= N <= 2,000,000
Solution
簡化題意:求兩個串前綴的最大循環同構子串的長度。
于是對于兩個串互相各做一遍擴展KMP,求出 Extend 數組。
掃一遍,如果滿足 ext1[ext2[i]+1]≥i?1 即說明滿足條件(畫圖理解一下)。
時間復雜度 O(N) 。
Code
#include<cstdio> using namespace std; const int N=2e6+5; int n,ans,mx,id; int next[N],ext1[N],ext2[N]; char s1[N],s2[N]; inline int max(int x,int y) {return x>y?x:y; } inline void exkmp(char *s,char *t,int *ext) {next[1]=n;mx=id=0;for(int i=2;i<=n;i++){int j=next[i-id+1];if(i+j>mx){j=max(mx-i+1,0);while(i+j<=n && t[i+j]==t[j+1]) j++;}next[i]=j;if(i+j-1>mx) mx=i+j-1,id=i;}mx=id=0;for(int i=1;i<=n;i++){int j=next[i-id+1];if(i+j>mx){j=max(mx-i+1,0);while(i+j<=n && s[i+j]==t[j+1]) j++;}ext[i]=j;if(i+j-1>mx) mx=i+j-1,id=i;} } int main() {scanf("%d",&n);scanf("%s",s1+1),scanf("%s",s2+1);exkmp(s1,s2,ext1),exkmp(s2,s1,ext2);for(int i=1;i<=n;i++){int j=ext1[i];if(ext2[j+1]>=i-1) ans=max(ans,i+j-1);}printf("%d",ans);return 0; }總結
以上是生活随笔為你收集整理的JZOJ 3648. 【GDOI2014】beyond的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 2682. 【WC2012选拔12.17
- 下一篇: JZOJ 3600. 【CQOI2014