poj 3294 Girls' research (manacher,回文串)

Note: This article is available in Chinese only. 本文暂无英文版本。 View original
poj 3294 题意：先做个简单替换，然后求替换后的字符串的最长回文串，以及这个最长回文串的开始和结束位置。思路：manacher。需要注意的是，返回下标的时候如果字符串长度为偶数，那么中间是没有字符的，需要特判一下（我的做法是 left+(ans%2==0)）。
代码实现
  1/* ***********************************************
  2Author :111qqz
  3Created Time :2016年04月18日 星期一 19时40分06秒
  4File Name :code/hdu/3294.cpp
  5************************************************ */
  6
  7#include <cstdio>
  8#include <cstring>
  9#include <iostream>
 10#include <algorithm>
 11#include <vector>
 12#include <queue>
 13#include <set>
 14#include <map>
 15#include <string>
 16#include <cmath>
 17#include <cstdlib>
 18#include <ctime>
 19#define fst first
 20#define sec second
 21#define lson l,m,rt<<1
 22#define rson m+1,r,rt<<1|1
 23#define ms(a,x) memset(a,x,sizeof(a))
 24typedef long long LL;
 25#define pi pair < int ,int >
 26#define MP make_pair
 27
 28using namespace std;
 29const double eps = 1E-8;
 30const int dx4[4]={1,0,0,-1};
 31const int dy4[4]={0,-1,1,0};
 32const int inf = 0x3f3f3f3f;
 33const int N=4E5+7;
 34char st[N],st2[N];
 35int p[N];
 36char ch;
 37
 38void manacher(char *s)
 39{
 40    int len = strlen(s);
 41
 42    for ( int i = len ;  i>= 0 ; i--)
 43    {
 44    s[i*2+2] = s[i];
 45    s[i*2+1] = '#';
 46    }
 47
 48    int id = 0 ;
 49    int maxlen = 0 ;
 50    int center = 0 ;
 51    s[0]='%';
 52    for ( int i = 2 ; i < 2*len+1 ; i++)
 53    {
 54    if (id+p[id]>i) p[i] = min(p[2*id-i],p[id]+id-i);
 55    else p[i] = 1;
 56
 57    while (s[p[i]+i]==s[i-p[i]]) p[i]++;
 58
 59    if (id+p[id]<i+p[i]) id = i;
 60    if (p[i]>maxlen)
 61    {
 62        maxlen = p[i];
 63        center = i;
 64    }
 65    }
 66    int ans = (maxlen-1);
 67   // cout<<"ans:"<<ans<<" center:"<<center<<endl;
 68    if (ans<2)
 69    {
 70    puts("No solution!");
 71    }
 72    else
 73    {
 74    int left = center/2-1-ans/2+(ans%2==0);
 75    int right = center/2-1+ans/2;
 76    printf("%d %d\n",left,right);
 77    for ( int  i = left ; i <= right ; i++) printf("%c",st2[i]);
 78    printf("\n");
 79    }
 80
 81}
 82int main()
 83{
 84    #ifndef  ONLINE_JUDGE
 85    freopen("code/in.txt","r",stdin);
 86  #endif
 87
 88    while (scanf("%c %s\n",&ch,st)!=EOF)
 89    {
 90        int k = ch-'a';
 91        int len = strlen(st);
 92        for ( int i = 0 ; i < len ; i++)
 93        {
 94        int tmp = (st[i]-'a'-k+26);
 95        st[i] = char(tmp+'a');
 96        st2[i] = st[i];
 97        }
 98    //    cout<<"st:"<<st<<endl;
 99        manacher(st);
100    }
101
102  #ifndef ONLINE_JUDGE
103  fclose(stdin);
104  #endif
105    return 0;
106}
107
108
109
110
111
112
113
114
115#include <bits/stdc++.h>
116#define fst first
117#define sec second
118#define lson l,m,rt<<1
119#define rson m+1,r,rt<<1|1
120#define ms(a,x) memset(a,x,sizeof(a))
121typedef long long LL;
122#define pi pair < int ,int >
123#define MP make_pair
124
125using namespace std;
126const double eps = 1E-8;
127const int dx4[4]={1,0,0,-1};
128const int dy4[4]={0,-1,1,0};
129const int inf = 0x3f3f3f3f;
130const int N=5E5+79;
131int t[2*N];
132char st[N];
133int p[2*N];  //p[i]为以第i个位置的字符为中心的回文串的半径长，默认为1.
134vector<int>vec[N];
135int a[2*N];
136char Ma[2*N];
137int Mp[2*N];
138int Manacher( char *s)
139{
140    int len = strlen(s);
141    int id = 0; //id表示之前位置的能延伸到最远位置的字符串的中心位置。
142    int maxlen = 0 ;  //maxlen是为了更新答案而用。。。就是记录了一个最大值。。。
143    int mx = 0 ;//mx为当前位置之前的回文串能延伸到的最远位置 即：max(p[j]+j) (j<i)
144                //如果知道之前能延伸到最远位置的字符串的中心位置的下标为id，那么就是p[id]+id;
145    for ( int i = len ; i >= 0 ; i--)  //插入'#'是为了将字符串长度奇偶性不同的统一考虑。
146    {
147        Ma[i+i+2] = s[i];
148        Ma[i+i+1] = '#';
149    }
150    Ma[0]='*'; //s[0] ='*'以及用字符数组最后默认的s[len+len+2]='\0'是为了下面while 循环的时候不溢出。。
151                //两边的字符一定要不一样。。用string的话记得两边都加字符。。。
152    for ( int i = 2 ; i < 2*len+1 ; i++)
153    {
154        if (Mp[id]+id>i) Mp[i] = min(Mp[2*id-i],Mp[id]+id-i);
155        else Mp[i] = 1;
156        while (Ma[i-Mp[i]]==Ma[i+Mp[i]]) Mp[i]++;
157        if (id+Mp[id]<i+Mp[i]) id = i;
158        if (maxlen<Mp[i]) maxlen = p[i];
159    }
160    return 2*len+2;
161    //这道题是问最长回文串的长度。。。如果是问回文串是什么的话。。。根据id和maxlen也可以构造出来。。。
162}
163
164int lowbit( int x) { return x&(-x);}
165void update ( int x,int delta)
166{
167    for ( int i = x; i < 2*N ; i+=lowbit(i))
168    t[i]+=delta;
169}
170int Sum( int x)
171{
172    int res = 0 ;
173    for ( int i = x ; i >= 1 ; i-=lowbit(i))
174    res = res + t[i];
175    return res;
176}
177int main()
178{
179        #ifndef  ONLINE_JUDGE
180       freopen("./in.txt","r",stdin);
181  #endif
182    int T;
183    cin>>T;
184    while (T--)
185    {
186        ms(t,0);
187        for ( int i = 0 ; i < N ; i++) vec[i].clear();
188        scanf("%s",st);
189        int n = strlen(st);
190        int len = Manacher(st);
191        int k = 1;
192        for ( int i = 2 ; i < len ; i+=2)
193        {
194        a[k] = Mp[i]/2-1;
195
196        printf("k:%d a[k]：%d p[i]:%d\n",k,a[k],Mp[i]);
197        vec[k-a[k]].push_back(k);
198        k++;
199        }
200        LL ans = 0 ;
201        for ( int i = 1 ; i <= n  ; i++)
202        {
203//      cout<<"a[i]:"<<a[i]<<endl;
204        for ( int j = 0 ; j < vec[i].size() ; j++)
205            update(vec[i][j],1);
206        ans += Sum(min(i+a[i],n)) -Sum(i);
207        }
208//      cout<<Sum(2*N)<<endl;
209        cout<<ans<<endl;
210    }
211  #ifndef ONLINE_JUDGE
212  fclose(stdin);
213  #endif
214    return 0;
215}
Related
