Description

Input

一行,一个字符串S

Output

一行,一个整数,表示所求值

Sample Input

cacao

Sample Output

54

HINT

2<=N<=500000,S由小写英文字母组成

Source

好吧,看来品酒大会是蒯了这个题的套路。。。今天做了三个套路一模一样的题目。。。

看到这个题直接枚举lcp的长度L,然后就又是品酒大会并查集从大到小合并的套路了。。。

每个height相当于是要维护lcp为L的对数,以及len的和,合并的时候瞎jb乘法原理一下。。。

哎,思维越来越僵化了。。。话说没想到nlogn^2能跑过50w。。。

// MADE BY QT666
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=600050;
int gi(){int x=0,flag=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-') flag=-1;ch=getchar();}while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();return x*flag;
}
int sa[N],len,y[N],rk,rnk[N],height[N],fa[N];
char a[N];
ll ansz[N],ansl[N],Ans,sz[N],le[N];
struct data{int fir,sec,id;
}x[N];
bool cmp(const data &a,const data &b){if(a.fir==b.fir) return a.sec<b.sec;else return a.fir<b.fir;
}
void work2(){rk=1;y[x[1].id]=rk;for(int i=2;i<=len;i++){if(x[i-1].fir!=x[i].fir||x[i-1].sec!=x[i].sec) rk++;y[x[i].id]=rk;}
}
void work(){sort(x+1,x+1+len,cmp);work2();for(int i=1;i<=len;i<<=1){for(int j=1;j+i<=len;j++) x[j].fir=y[j],x[j].sec=y[j+i],x[j].id=j;for(int j=len-i+1;j<=len;j++) x[j].fir=y[j],x[j].sec=0,x[j].id=j;sort(x+1,x+1+len,cmp);work2();if(rk==len) break;}
}
void get_height(){int kk=0;for(int i=1;i<=len;i++) rnk[sa[i]]=i;for(int i=1;i<=len;i++){if(kk) kk--;int j=sa[rnk[i]-1];while(a[i+kk]==a[j+kk]) kk++;height[rnk[i]]=kk;}
}
struct Data{int l,r,hi;
}g[N];
bool cmp2(const Data &a,const Data &b){return a.hi>b.hi;
}
int find(int x) {if(x!=fa[x]) fa[x]=find(fa[x]);return fa[x];
}
void merge(int x,int y,int hi){ansz[hi]+=sz[x]*sz[y];ansl[hi]+=le[x]*sz[y]+le[y]*sz[x];fa[x]=y;sz[y]+=sz[x];le[y]+=le[x];
}
int main(){scanf("%s",a+1);len=strlen(a+1);for(int i=1;i<=len;i++) x[i].id=i,x[i].fir=x[i].sec=a[i]-'a'+1;work();for(int i=1;i<=len;i++) sa[y[i]]=i;for(int i=1;i<=len;i++) sz[i]=1,fa[i]=i,le[i]=len-i+1;get_height();for(int i=2;i<=len;i++) g[i-1]=(Data){sa[i-1],sa[i],height[i]};sort(g+1,g+len,cmp2);for(int i=1;i<=len;i++){int x=find(g[i].l),y=find(g[i].r);if(x!=y) merge(x,y,g[i].hi);}for(int i=0;i<=len;i++) Ans+=ansl[i]-2*ansz[i]*i;printf("%lld\n",Ans);return 0;
}

  

转载于:https://www.cnblogs.com/qt666/p/7105983.html

bzoj 3238: [Ahoi2013]差异相关推荐

  1. BZOJ 3238: [Ahoi2013]差异 [后缀数组 单调栈]

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 2326  Solved: 1054 [Submit][Stat ...

  2. bzoj 3238: [Ahoi2013]差异(后缀数组+单调栈)

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 3443  Solved: 1562 [Submit][Stat ...

  3. @bzoj - 3238@ [Ahoi2013]差异

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 给定一个长度为 n 的字符串 S,令 Ti 表示它从第 i 个字 ...

  4. BZOJ 3238 [Ahoi2013]差异

    后缀自动机+后缀树 有一种不难想到的后缀数组的做法,详见黄学长题解 然而我做这题的目的主要是熟悉一下后缀自动机和后缀树- - 一个有趣的性质:反串的后缀自动机的parent树就是正串的后缀树,其树上的 ...

  5. 3238: [Ahoi2013]差异

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 2106  Solved: 953 [Submit][Statu ...

  6. 【BZOJ】 3238: [Ahoi2013]差异

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3238 求:$${\sum _{i=1}^{n-1}\sum _{j=i+1}^{n}len ...

  7. 【BZOJ】3238: [Ahoi2013]差异

    [题意]给定长度为n的小写字母字符串,令Ti表示以i开头的后缀,求Σ[Ti+Tj-2*lcp(Ti,Tj)],1<=i<j<=n. [算法]后缀自动机 [题解]Σ(Ti+Tj)只与n ...

  8. 【bzoj 3238】差异

    传送门~ 解题思路 求后缀数组以后搞个单调栈. 因为height数组的意义让这个单调栈很烦,注意要理清思路再写. 代码: #include<algorithm> #include<i ...

  9. BZOJ3238: [Ahoi2013]差异

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 4840  Solved: 2298 [Submit][Stat ...

最新文章

  1. JQuery插件:ScrollTo平滑滚动到页面指定位置
  2. udp重发机制_UDP 协议
  3. codevs 1507 酒厂选址
  4. Python如何读取文件
  5. 深入对比数据科学工具箱:Python和R的异常处理机制
  6. 【Bug】MQ消息与事务提交
  7. 添加服务oracle,oracle 11g(四)给oracle添加为系统服务(脚本)
  8. LeetCode(8):字符串转整数(atoi)
  9. idea springboot启动报SLF4J:Failed to load class “org.slf4j.impl.StaticLoggerBinder”
  10. linux蓝牙打印机驱动安装失败怎么办,打印机驱动安装失败怎么办?
  11. Q版京剧脸谱来喽——黑脸
  12. flutter图片聊天泡泡_Flutter极致的业务封装——各类聊天气泡(一)
  13. 矩阵的乘法口诀(一)
  14. 美股 基础知识2-常见问答
  15. 软件需求分析方法收集
  16. Linux man的使用以及汉化
  17. NVivo最新v1.0重磅来袭!不一样的用户界面,不一样的体验!
  18. Java对图片Base64转码--HTML对Base64解码 [ Java加强版 ]
  19. GBASE数据库助力城市轨道交通数据信息国产化
  20. c语言vs2017入门自学,写给自己的用VS2017学C语言[3]

热门文章

  1. RA维持期达标式减停生物制剂(DRESS研究)的18个月数据
  2. [Android]使用Gradle提交自己开源Android库到Maven中心库
  3. android常用网址
  4. 安全应对MySQL攻击
  5. 《花开在眼前》听得泪流满面
  6. Ghost in IE6.web标准网页IE6中的幽灵。
  7. Mike Krueger 加入Mono团队
  8. python 日志输出变量_python中记录打印的log模块logging的用法实例
  9. 浏览器的cookie和服务器存储的session有什么关联,下面妙文为你解答
  10. cron每月1号_微信服务号按粉丝标签分组群发消息怎样实现?