[Manacher 线段树 离线] 2015 计蒜之道 初赛 第三场 商品推荐走马灯
来源:互联网 发布:js confirm 样式 编辑:程序博客网 时间:2024/05/21 10:40
题目大意是求区间内所有回文子串的权值和
具体做法是离线处理询问 一个回文中心会对包含它的询问产生影响 而在询问左半边的中心最先被左边界限制 右半边也是这样
那么把询问拆成左右两部分
这里我还有个去重的过程 因为若一个回文中心是个字符而非间隙 那么我会把这个中心的权值算两边 那就减去好了
官方题解:http://blog.jisuanke.com/?p=56
题目询问的是一个区间内所有的回文串的权值和,因为考虑的对象是回文串,所以我们可以从回文中心下手。首先用Manacher 算法处理出每一个回文中心能够向两边延伸的最长回文串长度;又题目是支持离线的,所以我们可以把询问离线,然后考虑每一个回文中心对询问的影响。
对于一个询问,我们要考虑的显然只有那些在区间内部的回文中心,并且在区间左半部分的回文中心在向两边延伸时先碰到区间左端点;在区间右半部分的回文中心在向两边延伸时先碰到区间右端点,所以可以把每一个询问区间从中间拆开成两个,分别统计答案就行了。
接下来以统计左半部分区间的答案为例,右半部分类似。按区间的右端点排序,把位于当前区间右端点左边的回文中心用一棵线段树维护起来以方便统计答案。具体的维护方法是考虑这个回文中心到线段树上每一个点构成的回文串的权值和,这个值是对应原串的一个区间和,若处理出前缀和数组,就可以方便地在线段树上进行区间更新了,同时询问答案就是在线段树上询问一段区间的和。
用同样的方法处理一遍左半区间和右半区间,就能得出完整的答案。时间复杂度 O( mlogn ) 。
#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; } return *p1++;}inline void read(int &x){ char c=nc(),b=1; for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1; for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;}const int N=200005;int n,a[N],f[N];int m,ql[N],qr[N];ll sum[N];struct SEG{ int TH,M; ll T[N<<2],S[N<<2],F[N<<2]; int C[N<<2],L[N<<2]; void Build(int n,ll *sum){ for (M=1,TH=0;M<n+2;M<<=1,TH++); for (int i=1;i<=n;i++) T[M+i]=F[M+i]=0,S[M+i]=sum[i],C[M+i]=0,L[M+i]=1; for (int i=M-1;i;i--){ T[i]=0; F[i]=0; C[i]=0; S[i]=S[i<<1]+S[i<<1|1]; L[i]=L[i<<1]+L[i<<1|1]; } } void Pushdown(int rt){ int p; for (int i=TH;i;i--) if (C[p=rt>>i]){T[p<<1]+=F[p]*L[p<<1]-(ll)C[p]*S[p<<1]; C[p<<1]+=C[p]; F[p<<1]+=F[p];T[p<<1|1]+=F[p]*L[p<<1|1]-(ll)C[p]*S[p<<1|1]; C[p<<1|1]+=C[p]; F[p<<1|1]+=F[p];C[p]=F[p]=0; } } void Add(int s,int t,ll r){ for (Pushdown(s+=M-1),Pushdown(t+=M+1);s^t^1;){ if (~s&1) T[s^1]+=r*L[s^1]-S[s^1],F[s^1]+=r,C[s^1]++; if ( t&1) T[t^1]+=r*L[t^1]-S[t^1],F[t^1]+=r,C[t^1]++; T[s>>=1]=T[s<<1]+T[s<<1|1]; T[t>>=1]=T[t<<1]+T[t<<1|1]; } while (s>>=1) T[s]=T[s<<1]+T[s<<1|1]; } ll Sum(int s,int t){ if (s>t) return 0; ll ret=0; for (Pushdown(s+=M-1),Pushdown(t+=M+1);s^t^1;s>>=1,t>>=1){ if (~s&1) ret+=T[s^1]; if ( t&1) ret+=T[t^1]; } return ret; }}Seg;struct event{ int l,r,idx; event(int l=0,int r=0,int idx=0):l(l),r(r),idx(idx) { } bool operator < (const event &B) const{ return r<B.r; }}eve[N];int tot;ll Ans[N];namespace BIT{ ll c1[N],c2[N]; int maxn; inline void init(int n){ maxn=n; for (int i=1;i<=n;i++) c1[i]=c2[i]=0; } inline void add(int x,ll r){ for (int i=x;i<=maxn;i+=i&-i) c1[i]+=r,c2[i]+=r*x; } inline void add(int s,int t,ll r){ add(s,r); add(t+1,-r); } inline ll sum(int x){ ll ret1=0,ret2=0; for (int i=x;i;i-=i&-i) ret1+=c1[i],ret2+=c2[i]; return ret1*(x+1)-ret2; } ll sum(int s,int t){ if (s>t) return 0; return sum(t)-sum(s-1); }}int main(){ int t=0,pnt; freopen("t.in","r",stdin); freopen("t.out","w",stdout); read(n); read(m); a[0]=-10000; a[++t]=-1000; for (int i=1;i<=n;i++) read(a[++t]),a[++t]=-1000; n=t; a[n+1]=-10000; int id=1; f[1]=1; for (int i=2;i<=n;i++){ if (i<id+f[id]) f[i]=min(f[2*id-i],id+f[id]-i); else f[i]=1; while (a[i-f[i]]==a[i+f[i]]) f[i]++; if (id+f[id]<i+f[i]) id=i; } for (int i=1;i<=m;i++) read(ql[i]),read(qr[i]),ql[i]=2*ql[i]-1,qr[i]=2*qr[i]+1; tot=0; for (int i=1;i<=m;i++) eve[++tot]=event(ql[i],((ql[i]+qr[i])>>1)-1,i); sort(eve+1,eve+tot+1); for (int i=1;i<=n;i++) sum[i]=sum[i-1]+(a[i]==-1000?0:a[i]); Seg.Build(n,sum-1); pnt=1; BIT::init(n>>1); for (int i=1;i<=n;i++){ Seg.Add(i-f[i]+1,i,sum[i]); if (~i&1) BIT::add((i>>1)-(f[i]>>1)+1,i>>1,a[i]); while (pnt<=tot && eve[pnt].r<i) pnt++; while (pnt<=tot && eve[pnt].r==i){ Ans[eve[pnt].idx]+=Seg.Sum(eve[pnt].l,eve[pnt].r); Ans[eve[pnt].idx]-=BIT::sum(((eve[pnt].l-1)>>1)+1,eve[pnt].r>>1); pnt++; } } reverse(a+1,a+n+1); reverse(f+1,f+n+1); tot=0; for (int i=1;i<=m;i++) eve[++tot]=event(n-qr[i]+1,n-((ql[i]+qr[i])>>1)+1,i); sort(eve+1,eve+tot+1); for (int i=1;i<=n;i++) sum[i]=sum[i-1]+(a[i]==-1000?0:a[i]); Seg.Build(n,sum-1); pnt=1; BIT::init(n>>1); for (int i=1;i<=n;i++){ Seg.Add(i-f[i]+1,i,sum[i]); if (~i&1) BIT::add((i>>1)-(f[i]>>1)+1,i>>1,a[i]); while (pnt<=tot && eve[pnt].r<i) pnt++; while (pnt<=tot && eve[pnt].r==i){ Ans[eve[pnt].idx]+=Seg.Sum(eve[pnt].l,eve[pnt].r); Ans[eve[pnt].idx]-=BIT::sum(((eve[pnt].l-1)>>1)+1,eve[pnt].r>>1); pnt++; } } for (int i=1;i<=m;i++) printf("%lld\n",Ans[i]); return 0;}
0 0
- [Manacher 线段树 离线] 2015 计蒜之道 初赛 第三场 商品推荐走马灯
- 2015 计蒜之道 初赛 商品推荐走马灯(dp)
- 计蒜之道 初赛 第三场 题解 Manacher o(n)求最长公共回文串 线段树
- 2017 计蒜之道 初赛 第三场
- 2017 计蒜之道 初赛 第三场 A+B
- 2017 计蒜之道 初赛 第三场
- 2017计蒜之道第六场初赛 微软大楼设计方案(中等)(线段树)
- [非旋转Treap 权值线段树] 2016 计蒜之道 初赛 第五场 腾讯的新游戏
- 计蒜之道 初赛 第三场--腾讯手机地图 题解
- 2016计蒜之道初赛第三场----百度的科学计算器【模拟】
- 2016 计蒜之道 初赛 第三场 百度帐号的选取方案(中等) KMP
- 2017计蒜之道程序设计大赛初赛第三场题解
- 2017 计蒜之道 初赛 第三场-A. 腾讯课堂的物理实验
- 2017 计蒜之道初赛第三场 A.腾讯课堂的物理实验
- [最大密度子图] 2017 计蒜之道 初赛 第三场 腾讯狼人杀
- 2017 计蒜之道 初赛 第三场 A. 腾讯课堂的物理实验
- 2017 计蒜之道 初赛 第三场 腾讯狼人杀(简单)(暴搜)
- [最大密度子图] 2017 计蒜之道 初赛 第三场. 腾讯狼人杀
- JSQLParser源码分析之CreateTable
- DB 查询分析器7.01 新增的保存执行结果到多个文件功能
- inline内联的用法与作用
- instanceof判定类的实例
- Unity3d Android SDK接入解析(一)Unity3d 与 Android之间的互相调用
- [Manacher 线段树 离线] 2015 计蒜之道 初赛 第三场 商品推荐走马灯
- Android RecyclerView
- 利用scp 远程上传下载文件/文件夹,ssh上传下载文件
- RelativeLayout相对布局
- ecshop模板动态引入lbi
- android下获取无线wif信号、ssid、MAC等操作类
- 票据:汇票业务解析
- Flex布局——色子的小案例
- 根文件系统相关操作