HDU 5785 Interesting manacher + 延迟标记
题意:给你一个串,若里面有两个相邻的没有交集的回文串的话,设为S[i...j] 和 S[j+1...k],对答案的贡献是i*k,就是左端点的值乘上右端点的值。
首先,如果s[x1....j]、s[x2....j]、s[x3....j]....s[xn....j]、是回文串,而且s[j+1...y1]、s[j+1...y2]、s[j+1...y3]、也是回文串,那么这些串对答案的贡献一共就是(x1+x2+...+xn)*(y1+y2+....+yn)
所以想到了用cntL[i]表示:以i这个字符为结尾的回文串的左端点之和,cntR[i]表示:以i这个字符为开始的回文串的右端点之和。
所以答案就是for i=1 to lenstr-1 ans += cntL[i]*cntR[i+1];
通俗点说吧,根据manacher,p[i]就表示以i这个字符为中心的回文串的半径长度,所以有s[i-p[i]...i+p[i]]是一个回文串
所以cntL[i+p[i]]就要加上i-p[i],然后因为s[i-p[i]+1...i+p[i]-1]也是一个回文串,所以cntL[i+p[i]-1]就要加上i-p[i]+1等等之类的。所以就是在[i..i+p[i]]这个区间上,依次加上 i、i-1、i-2、.....i-p[i],这样一个公差为-1的等差数列。这个可以用延迟标记来做
先看简单的,在[L,R]上加上一个数val,我们只需要cnt[L] += val; cnt[R+1] -= val; 然后O(n)扫描,即可得到这个数组的所有值。但是现在是等差数列呢?其实是差不多的,只不过要开多一个数组而已。
首先,如果是在[begin,end]上加上一个以val为首相的等差数列,可以这样考虑。首先cnt[begin] += val;然后可以算出结束的时候的值是多少把?cnt[end+1] -= calc_end_val;即可。然后中间的公差,因为可能有叠加现象,我们要开多个数组subL[i]表示这个位置有多少个-1,这个subL[i]就像一开始的最简单的这样传递下去即可。就是subL[begin] += 1; subL[end+1] -= 1;
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
#include <cstdio> #include <cstdlib> #include <cstring> #include <cmath> #include <algorithm> using namespace std; #define inf (0x3f3f3f3f) typedef long long int LL; #include <iostream> #include <sstream> #include <vector> #include <set> #include <map> #include <queue> #include<stack> #include <string> const int maxn = 1000000+20; const int MOD = 1000000007; int cntL[maxn],cntR[maxn]; int subL[maxn],subR[maxn]; char str[maxn<<1]; int p[maxn<<1]; void addL (int begin,int end,int val) { if (end<begin) return ; cntL[begin] += val; cntL[end+1] -= begin-end+val; subL[begin+1]++; subL[end+1]--; cntL[begin] %= MOD; cntL[end+1] = (cntL[end+1]+MOD)%MOD; subL[begin+1]%=MOD; subL[end+1]=(subL[end+1]+MOD)%MOD; return ; } void addR (int begin,int end,int val) { if (end<begin) return ; cntR[begin] += val; cntR[end+1] -= begin-end+val; subR[begin+1]++;//记录有多少个数列覆盖的,记录公差 subR[end+1]--; cntR[begin] %= MOD; cntR[end+1] = (cntR[end+1]+MOD)%MOD; subR[begin+1]%=MOD; subR[end+1]=(subR[end+1]+MOD)%MOD; return ; } int manacher (char str[],int lenstr) { str[0]='*'; for (int i=lenstr;i>=0;i--) { str[i+i+2]=str[i+1]; str[i+i+1]='#'; } int id=0,maxlen=0; for (int i=2;i<=2*lenstr+1;i++) { if (p[id]+id>i) { p[i]=min(p[id]+id-i,p[2*id-i]); } else p[i]=1; while (str[i+p[i]] == str[i-p[i]]) ++p[i]; if (p[id]+id<p[i]+i) id=i; maxlen=max(maxlen,p[i]); } return maxlen-1; } void init_arr(int gg) { for (int i=1;i<=gg;++i) { subL[i] += subL[i-1]; cntL[i] += cntL[i-1]-subL[i]; subR[i] += subR[i-1]; cntR[i] += cntR[i-1]-subR[i]; subL[i] %= MOD; subR[i] %= MOD; cntL[i]=(cntL[i]+MOD)%MOD; cntR[i]=(cntR[i]+MOD)%MOD; } return ; } void init () { memset(cntL,0,sizeof cntL); memset(cntR,0,sizeof cntR); memset(subL,0,sizeof subL); memset(subR,0,sizeof subR); return ; } void work () { init(); int lenstr = strlen(str+1); manacher(str,lenstr); for (int i=2;i<=2*lenstr+1;++i) { int len = p[i]-1; if (len==0) continue; if (i&1) //'#'所在位置 { int pos = i/2;//截断 //len必须是偶数 int t = len/2; int begin = pos-t+1; //这些细节要慢慢算啦,慢慢推公式 int end = t+pos; addL(pos+1,end,pos); addR(begin,pos,end); } else //字母所在位置 { int pos = i/2; len--; len/=2; int begin = pos-len; int end = pos+len; addL(pos,end,pos); addR(begin,pos,end); } } init_arr(lenstr); LL ans = 0; for (int i=1;i<=lenstr-1;++i) { ans += 1LL*cntL[i]*cntR[i+1]; ans %= MOD; } printf ("%I64d\n",ans); return ; } int main () { #ifdef local freopen("data.txt","r",stdin); #endif while(scanf("%s",str+1)!=EOF) work(); return 0; }
posted on 2016-08-17 22:29 stupid_one 阅读(182) 评论(0) 编辑 收藏 举报