【清华集训2017模拟12.09】塔
2017-12-14 16:52
344 查看
题目
有一个塔,他的名字叫做粽粑,粽粑的每一层都有一个颜色 .粽粑非常厉害,它在吸收天地精华之后会长高.粽粑的长高方式有两种:
1.在塔顶长出一层.
2.在塔底长出一层,即原来的第一层变成第二层,第二层变成第三层,以此类推,新长出来的是第一层.
粽粑有可能在某个时刻不是很开心,这个时候它会撤销它的前若干次长高.
你现在想知道粽粑长高的奥秘,于是找到了粽粑,发现它的入口上写着这么一句话:要进入粽粑,请找出一段最长的塔的区间,满足翻转后颜色不变.
粽粑会不断的长高(或撤销),所以它每次长高后你都要回答.为了你的方便,粽粑一开始的高度为0.
分析
我发现我的字符串学傻了⊙﹏⊙‖∣回文树模板题,但是卡不过。
假设上一次的答案为ans,
对于头尾,我们分别维护长度为ans的字符串的两个哈希值(这个长度为ans的字符串的哈希值高位分别在前或后),
看哈希值是否一样就能判断是否为回文串。每次长高只会让ans+1或+2,判断一下。
#pragma GCC optimize(3) #include <cmath> #include <iostream> #include <cstdio> #include <cstdlib> #include <cstring> #include <algorithm> #include <queue> #include <map> const int maxlongint=2147483647; const long long mo=1e9+7; const int N=10000005; using namespace std; int n,m,ans ,top; long long sum [4]; char t ; short int num[N*3]; long long mi ,ny,val; long long pow1(long long x,int y) { long long sum=1; for(;y;x=x*x%mo,y>>=1) if(y&1) sum=sum*x%mo; return sum; } int main() { scanf("%d\n",&n); ny=pow1(101,mo-2); int l=N,r=l-1; mi[0]=1; for(int i=1;i<=n;i++) mi[i]=mi[i-1]*101%mo; for(int i=1;i<=n;i++) { char c,c1,c2; c=getchar(),c1=getchar(),c2=getchar(); int x=((c1-'0')*10+c2-'0'+(ans[top]!=0?ans[top]:0))%100; if(c=='1') { num[++r]=x,t[++top]='1'; ans[top]=ans[top-1]; sum[top][0]=sum[top-1][0],sum[top][1]=sum[top-1][1],sum[top][2]=sum[top-1][2],sum[top][3]=sum[top-1][3]; long long v1=sum[top][0],v2=sum[top][1],len=ans[top]; if(len) v1=((v1-mi[len-1]*num[r-len]%mo+mo)*101%mo+num[r])%mo; if(len) v2=((v2-num[r-len]+mo)%mo*ny%mo+mi[len-1]*num[r])%mo; sum[top][0]=v1,sum[top][1]=v2; v1=(v1+mi[len]*num[r-len])%mo; v2=(v2*101+num[r-len])%mo; if(v1==v2) { ans[top]=len+1; sum[top][0]=v1,sum[top][1]=v2; sum[top][2]=(sum[top][2]+mi[len]*num[l+len])%mo; sum[top][3]=(sum[top][3]*101+num[l+len])%mo; } if(r-l+1>=len+2) { v1=(v1+mi[len+1]*num[r-len-1])%mo; v2=(v2*101+num[r-len-1])%mo; if(v1==v2) { ans[top]=len+2; sum[top][0]=v1,sum[top][1]=v2; sum[top][2]=(sum[top][2]+mi[len]*num[l+len])%mo; sum[top][3]=(sum[top][3]*101+num[l+len])%mo; sum[top][2]=(sum[top][2]+mi[len+1]*num[l+len+1]%mo)%mo; sum[top][3]=(sum[top][3]*101+num[l+len+1])%mo; } } } else if(c=='2') { num[--l]=x,t[++top]='0'; ans[top]=ans[top-1]; sum[top][0]=sum[top-1][0],sum[top][1]=sum[top-1][1],sum[top][2]=sum[top-1][2],sum[top][3]=sum[top-1][3]; long long v1=sum[top][2],v2=sum[top][3],len=ans[top]; if(len) v1=((v1-mi[len-1]*num[l+len]%mo+mo)*101%mo+num[l])%mo; if(len) v2=((v2-num[l+len]+mo)%mo*ny%mo+mi[len-1]*num[l])%mo; sum[top][2]=v1,sum[top][3]=v2; v1=(v1+num[l+len]*mi[len])%mo; v2=(v2*101+num[l+len])%mo; if(v1==v2) { ans[top]=len+1; sum[top][2]=v1,sum[top][3]=v2; sum[top][0]=(sum[top][0]+mi[len]*num[r-len])%mo; sum[top][1]=(sum[top][1]*101+num[r-len])%mo; } if(r-l+1>=len+2) { v1=(v1+num[l+len+1]*mi[len+1]%mo)%mo; v2=(v2*101+num[l+len+1])%mo; if(v1==v2) { ans[top]=len+2; sum[top][2]=v1,sum[top][3]=v2; sum[top][0]=(sum[top][0]+mi[len]*num[r-len])%mo; sum[top][1]=(sum[top][1]*101+num[r-len])%mo; sum[top][0]=(sum[top][0]+mi[len+1]*num[r-len-1])%mo; sum[top][1]=(sum[top][1]*101+num[r-len-1])%mo; } } } else { for(;x--;top--) { if(t[top]=='1') num[r--]=0; else num[l++]=0; sum[top][0]=sum[top][1]=sum[top][2]=sum[top][3]=0; } } val+=1ll*ans[top]; } printf("%lld\n",val); }
相关文章推荐
- 【清华集训2017模拟12.09】Tree
- [JZOJ5495]【清华集训2017模拟12.09】MiniumCut
- JZOJ 5496. 【清华集训2017模拟12.09】Tree
- 【JZOJ 5495】【清华集训2017模拟12.09】MiniumCut(最小割树)
- 5483. 【清华集训2017模拟11.26】简单路径
- [JZOJ5498]【清华集训2017模拟12.10】大佬的难题
- jzoj5316 【清华集训2017模拟8.19】merge
- [JZOJ5485]【清华集训2017模拟11.26】字符串
- 【JZOJ5317】【清华集训2017模拟8.19】func
- [jzoj]5484. 【清华集训2017模拟11.26】快乐树(树形DP)
- 【清华集训2017模拟12.10】回文串(回文树+树链剖分)
- 【JZOJ 5282】【清华集训2017模拟】等差划分数
- 【清华集训2017模拟】Create
- 【JZOJ 5276】【清华集训2017模拟】神奇的玩具
- [JZOJ5500]【清华集训2017模拟12.10】营养餐
- 【JZOJ 5284】【清华集训2017模拟】超级翻转
- 【JZOJ 5296】【清华集训2017模拟】Sequence
- JZOJ 5485. 【清华集训2017模拟11.26】字符串
- 【JZOJ 5295】【清华集训2017模拟】Create
- JZOJ5489. 【清华集训2017模拟11.28】海明距离