string str = "bacbababadababacambabacaddababacasdsd";
string ptr = "ababaca";

/*目标字符串是ababaca*/
NEXT[0] = -1;
NEXT[1] = 0;"a"
NEXT[2] = 0;"ab"
NEXT[3] = 1;"aba"
NEXT[4] = 2;"abab"
NEXT[5] = 3;"ababa"
NEXT[6] = 0;:ababac"NEXT[7] = 1;""ababaca"

#include<cstdio>
#include<iostream>
#include<string>
using namespace std;
const int maxn = 1000002;
string S;//主串
string T;//模式串int NEXT[maxn];//当前位置如果不匹配则指针回溯的位置即最长前后缀相同的长度
int slen, tlen;void getnext()//生成next数组
{//a b a c a b a d a b a c a b aint j, k;//j是模式串的当前角标,k是在模式串当前角标之前已经匹配的前缀后缀的长度
    j = 0; k = -1; NEXT[0] = -1;/*已知的next[j]和k一位置T[0]T[1]...T[k - 1]=T[j - k]T[j - k + 1]....T[j - 1];推算next[j+1]则将T[j]与T[k]比较（即判断该位置是否能接着前后缀相同)若T[k] == T[j] 则T[0] ..... T[k] = T[j - k] ... T[j-1]即next[++j] = ++k;意味着已匹配的前后缀长度加1之后，T[0..j]拥有k+1的相同前后缀也就是在j+1位与其不匹配时，就返回第k+1位，即相同的前缀的后一位若T[k] != T[j] 即不匹配时，不可能一棍子打死从0开始，退一步找一个相对短的匹配前缀k代表的是目前为止相同前后缀的长度，若当前j不匹配，则将k跳回next[k]处可以理解为T数组自己对自己用next已经记录的部分匹配,既然k处不匹配了相当于前缀在k处与自己失配，将k前移回上一个能够匹配的前缀长度，判断是否能适配现在的T[j]就像已经适配的前缀是ababac 后缀是ababab，发现b和c不相同，回溯一次前缀为ababa最后一位a仍然与当前T[j] = b，不相同，再回溯一次为b与当前相同，就继续开始判定了这里说一下默认的是next[k] = k-1;这一点通过相同时可以很容易看出来，但是不够全面，所以使用k = next[k];而不是k--*/
    while (j < tlen){if (k == -1 || T[j] == T[k])NEXT[++j] = ++k;elsek = NEXT[k];  }
}
int KMP_index()//读取出现的位置（第一个）
{int i = 0, j = 0;//i指向主串，j指向模式串，i不会回退
    getnext();while (i < slen && j < tlen){if (j == -1 || S[i] == T[j]){i++;j++;}else{j = NEXT[j];}}if (j == tlen)return i - tlen;else return -1;
}
int KMP_Count()//出现的次数
{int ans = 0;int i, j = 0;if (slen == 1 && tlen == 1){if (S[0] == T[0])return 1;else return 0;}getnext();for (i = 0; i < slen; i++){while (j > 0 && S[i] != T[j])j = NEXT[j];if (S[i] == T[j])j++;if (j == tlen){ans++;j = NEXT[j];}}return ans;
}
int main()
{int TT;int i, cc;cin >> TT;while (TT--){cin >> S >> T;slen = S.length();tlen = T.length();cout << "模式串T在主串S中首次出现的位置是: " << KMP_index() << endl;cout << "模式串T在主串S中出现的次数为: " << KMP_Count() << endl;for (int i = 0; i < tlen; i++)cout << NEXT[i] << " ";}return 0;}

string T = "aaaabaaaa";

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
//next[i] : x[i .... m - 1]与x[0 .... m-1]的最长公共前缀
//extend[i]: y[i ... n - 1] 与x[0 .... m - 1]的最长公共前缀
void get_next(char *x, int m, int *next)//m为子串长度，x为目标子串
{next[0] = m;int j = 0;//初始化next[1]while (j + 1 < m && x[j] == x[j + 1]) j++;next[1] = j;//推导之后的值int k = 1;for (int i = 2; i < m; i++){int p = next[k] + k - 1;int L = next[i - k];if (i + L < p + 1) next[i] = L;else{j = max(0, p - i + 1);while (i + j < m && x[i + j] == x[j]) j++;next[i] = j;k = i;}}}
void ExKMP(char *x, int m, char *y, int n, int *next, int *extend)
{get_next(x, m, next);int j = 0;while (j < n && j < m && x[j] == y[j]) j++;extend[0] = j;int k = 0;for (int i = 1; i < n; i++){int p = extend[k] + k - 1;int L = next[i - k];if (i + L < p + 1) extend[i] = L;else{j = max(0, p - i + 1);while (i + j < n && j < m && y[i + j] == x[j]) j++;extend[i] = j;k = i;}}
}

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;

const int maxn = 110010;

//转换原始串int INIT(char *st, char *tmp){    int i, len = strlen(st);    tmp[0] = '@';//字符串开头增加一个特殊字符，防止越界    for (i = 1; i <= 2 * len; i += 2)    {        tmp[i] = '#';        tmp[i + 1] = st[i / 2];    }    tmp[2 * len + 1] = '#';    tmp[2 * len + 2] = '$';//字符串结尾加一个字符，防止越界    tmp[2 * len + 3] = 0;    return 2 * len + 1;//返回转换字符串的长度}//Manacher算法计算过程int MANACHER(char *st, int len, int *Len){    int mx = 0, ans = 0, po = 0;//mx即为当前计算回文串最右边字符的最大值    for (int i = 1; i <= len; i++)    {        if (mx>i)            Len[i] = min(mx - i, Len[2 * po - i]);//在Len[j]和mx-i中取个小        else            Len[i] = 1;//如果i>=mx，要从头开始匹配        while (st[i - Len[i]] == st[i + Len[i]])            Len[i]++;        if (Len[i] + i>mx)//若新计算的回文串右端点位置大于mx，要更新po和mx的值        {            mx = Len[i] + i;            po = i;        }        ans = max(ans, Len[i]);    }    return ans - 1;//返回Len[i]中的最大值-1即为原串的最长回文子串额长度 }int main(){    char str[maxn];//原字符串    while (cin >> str)    {        char tmp[maxn << 1];//转换后的字符        int Len[maxn << 1];        int len = INIT(str, tmp);

        cout << MANACHER(tmp, len, Len) << endl;    }}