题目描述

传送门

题目大意:两个长度为r的子串相同称为r相似。两个r相似的子串的价值是两个子串开头位置的价值乘积。
求r相似的子串数量和最大的代价(r=0…n-1)

题解

后缀自动机parent树上两个点的lca的len是两个字符串的最长公共后缀,但是这个要求的是开头位置的乘积,所以我就考虑将串反置,建立后缀自动机并构建parent树。
然后在parent树上进行树形DP ,我们知道在parent树上的每个节点代表的子串都是长度区间为(len(fa),len(i)]的一段区间,但是我们其实在做的时候只需要考虑len(i)就可以了,最后做完了,在将ans[i]+=ans[i+1].这是为什么呢?你想啊能匹配到i+1那么长那么他一定有i那么长并且满足条件的子串啊。
那么怎么DP?其实也好做,因为我们要求的是最长公共后缀严格是len的也就是不允许两个节点出自该点的同一个子树,那么我们每加入一个子树就将子树的size给他,然后算下一个子树的时候,用累加的值计算答案即可。最大代价的计算方法也是同样的道理。
需要注意的是子树上的某些点不是真是存在的,即没有对应子串的节点,这些点在统计的时候需要特别注意。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 600003
#define inf 1000000000
#define LL long long
using namespace std;
int point[N],nxt[N],v[N],ch[N][30],fa[N],l[N];
int n,root,np,cnt,last,p,nq,q,tot,mp[N],mark[N];
LL a[N],mx[N],mn[N],ans[N],cnt1[N],size[N];
char s[N];
void add(int x,int y)
{tot++; nxt[tot]=point[x]; point[x]=tot; v[tot]=y;//cout<<x<<" "<<y<<" "<<mp[y]<<endl;
}
void extend(int x)
{int c=s[x]-'a'+1;p=last; np=++cnt; last=np; mark[np]=1; mp[np]=x;l[np]=l[p]+1;for (;p&&!ch[p][c];p=fa[p]) ch[p][c]=np;if (!p) fa[np]=root; else {q=ch[p][c];if (l[q]==l[p]+1) fa[np]=q;else {nq=++cnt; l[nq]=l[p]+1;memcpy(ch[nq],ch[q],sizeof(ch[nq]));fa[nq]=fa[q];fa[q]=fa[np]=nq;for (;ch[p][c]==q;p=fa[p]) ch[p][c]=nq;}}
}
void dfs(int x)
{if (!mark[x]) mx[x]=mn[x]=inf,size[x]=0;else size[x]=1,mx[x]=mn[x]=a[mp[x]];for (int i=point[x];i;i=nxt[i]) {dfs(v[i]);cnt1[l[x]]+=size[x]*size[v[i]];if (mx[x]!=inf)ans[l[x]]=max(ans[l[x]],max(mx[v[i]]*mx[x],mn[v[i]]*mn[x])),mx[x]=max(mx[x],mx[v[i]]),mn[x]=min(mn[x],mn[v[i]]);else mx[x]=mx[v[i]],mn[x]=mn[v[i]];size[x]+=size[v[i]];    }
}
int main()
{scanf("%d",&n);scanf("%s",s+1); reverse(s+1,s+n+1);last=root=++cnt;for (int i=1;i<=n;i++) extend(i);for (int i=1;i<=cnt;i++) add(fa[i],i);for (int i=1;i<=n;i++) scanf("%lld",&a[n-i+1]);memset(ans,128,sizeof(ans));dfs(1);for (int i=n-1;i>=0;i--) cnt1[i]+=cnt1[i+1],ans[i]=max(ans[i],ans[i+1]);for (int i=0;i<=n-1;i++) if (cnt1[i]) printf("%lld %lld\n",cnt1[i],ans[i]);else printf("0 0\n");
} 

bzoj 4199: [Noi2015]品酒大会 (后缀自动机+DP)相关推荐

  1. BZOJ.4199.[NOI2015]品酒大会(后缀自动机 树形DP)

    BZOJ 洛谷 后缀数组做法. 洛谷上SAM比SA慢...BZOJ SAM却能快近一倍... 只考虑求极长相同子串,即所有后缀之间的LCP. 而后缀的LCP在后缀树的LCA处.同差异这道题,在每个点处 ...

  2. BZOJ 4199 [Noi2015]品酒大会(后缀自动机 + parent树上统计)

    4199: [Noi2015]品酒大会 Time Limit: 10 Sec  Memory Limit: 512 MB Submit: 1598  Solved: 905 [Submit][Stat ...

  3. [BZOJ]4199: [Noi2015]品酒大会(后缀数组+笛卡尔树)

    Time Limit: 10 Sec  Memory Limit: 512 MB Description Input Output Sample Input 10 ponoiiipoi 2 1 4 7 ...

  4. UOJ #131 BZOJ 4199 luogu P2178【NOI2015】品酒大会 (后缀自动机、树形DP)

    UOJ #131 BZOJ 4199 luogu P2178[NOI2015]品酒大会 (后缀自动机.树形DP) 水是水,但是写出了不少问题,因此写一发博客. https://www.luogu.or ...

  5. [BZOJ4199][NOI2015]品酒大会-后缀数组

    品酒大会 题目描述 一年一度的"幻影阁夏日品酒大会"隆重开幕了.大会包含品尝和趣味挑战 两个环节,分别向优胜者颁发"首席品酒家"和"首席猎手" ...

  6. 洛谷P2178 [NOI2015]品酒大会 后缀数组+单调栈

    P2178 [NOI2015]品酒大会 题目链接 https://www.luogu.org/problemnew/show/P2178 题目描述 一年一度的"幻影阁夏日品酒大会" ...

  7. 洛谷 P2178 [NOI2015]品酒大会 解题报告

    P2178 [NOI2015]品酒大会 题目描述 一年一度的"幻影阁夏日品酒大会"隆重开幕了.大会包含品尝和趣味挑战 两个环节,分别向优胜者颁发"首席品酒家"和 ...

  8. [NOI2015]品酒大会

    题目描述 一年一度的"幻影阁夏日品酒大会"隆重开幕了.大会包含品尝和趣味挑战 两个环节,分别向优胜者颁发"首席品酒家"和"首席猎手"两个奖项 ...

  9. hdu5343 后缀自动机+dp

    给定两个串,分别截取字串X和Y,连接组成X+Y,求不同的X+Y的方案数. 对于X+Y,如果重复的部分其实就是从同一个X+Y的某个地方断开弄成不同的X和Y,那么只要使得X和X+Y匹配得最长就行了. 因此 ...

  10. BZOJ2806(后缀自动机+DP)

    题目:http://61.187.179.132/JudgeOnline/problem.php?id=2806 题意:给定一个由M个01串组成的字典.依据这个字典和一个阀值L,可以断言一个01串是否 ...

最新文章

  1. string替换_GEE数据类型—String,Number
  2. laravel log 对象_swoole运行模式加速laravel应用的详细介绍(life)
  3. c 开发服务器的性能,高性能网络编程(三):下一个10年,是时候考虑C10M并发问题了...
  4. C++带参数控制台(实例)
  5. Python+socket实现TCP套接字服务端自由限速
  6. 「一道面试题」ArrayList和LinkedList有什么区别及使用场景
  7. 《善用佳软:高效能人士的软件应用之道》一2.3 思维导图:思想脉络的逻辑结构...
  8. 易我数据恢复向导 v2.0 官网
  9. 什么是台式机主板的“双通道”认知误区?
  10. 编译原理(八)消除空产生式
  11. 前端js中文转拼音(例:张三转为ZhangSan)
  12. hdu 3018 图 欧拉回路 并查集
  13. 设置屏幕常亮 switch
  14. 万字长文读透 Redis
  15. 奖励稀疏_好奇心解决稀疏奖励任务
  16. 评分卡:WOE、IV、PSI计算及ROC和KS曲线
  17. 活动预告|CodeWisdom软件供应链系列学术报告:第5期(鲍凌峰 浙江大学)
  18. css超级无敌不会,一大堆疑问大全
  19. 想升职?不妨和老板争一争!
  20. 《SEM运营36计》第10计:笑里藏刀

热门文章

  1. SAP内部培训效果考察表存档
  2. “当前页面的脚本发生错误”如何解决
  3. 数据结构 课堂测验错题_200911
  4. matlab速成学习
  5. Word转成PDF格式会变吗?教你3个免费方法
  6. Python 绘制正多边形
  7. ue4 中KismetProceduralMeshLibrary SliceProceduralMesh的原理
  8. 如何设置qq说说展示时间_QQ空间说说新增定时发表功能
  9. linux V4L2子系统——v4l2架构(5)之v4l2_device与v4l2_subdev异步机制
  10. 高数笔记(六):拉格朗日中值定理,柯西中值定理,洛必达法则,泰勒定理