[SDOI2016] 生成魔咒 - 后缀数组,平衡树,STL,时间倒流
[SDOI2016] 生成魔咒
Description
初态串为空,每次在末尾追加一个字符,动态维护本质不同的子串数。
Solution
考虑时间倒流,并将串反转,则变为每次从开头删掉一个字符,即每次从后缀集合中删掉一个后缀。
预处理出后缀数组和高度数组后,用平衡树维护所有后缀集合(按照后缀排序),要删除一个后缀 \(S[sa[p],n]\) 时,找到它在平衡树上的前驱 \(u\) 和后继 \(v\) ,如果都存在,那么这一步的贡献就是
\]
约定 \(h[p]\) 表示 \(S[sa[p],n]\) 与 \(S[sa[p-1],n]\) 的 LCP 长度。
如果 \(u\) 或 \(v\) 不存在,则当作 \(LCP\) 为零处理,仍然成立。
求 \(LCP\) 可以暴力用 ST 表维护。但考虑到这里每次删除操作最多只会再影响一个元素,我们可以顺便记录一下,即当我们删除 \(p\) 的时候令 \(h[v] = Min(h[p],h[v])\) 即可。
#include <bits/stdc++.h>
using namespace std;
#define int long long
int n,m,sa[1000005],y[1000005],u[1000005],v[1000005],o[1000005],r[1000005],h[1000005],T;
int str[1000005];
long long ans;
map <int,int> mp;
set <int> s;
vector <int> an;
signed main()
{
scanf("%lld",&n);
m=n;
for(int i=1; i<=n; i++)
scanf("%lld",&str[i]);
reverse(str+1,str+n+1);
for(int i=1; i<=n; i++)
mp[str[i]]++;
int ind=0;
for(map<int,int>::iterator it=mp.begin(); it!=mp.end(); ++it)
it->second=++ind;
for(int i=1; i<=n; i++)
str[i]=mp[str[i]];
for(int i=1; i<=n; i++)
u[str[i]]++;
for(int i=1; i<=m; i++)
u[i]+=u[i-1];
for(int i=n; i>=1; i--)
sa[u[str[i]]--]=i;
r[sa[1]]=1;
for(int i=2; i<=n; i++)
r[sa[i]]=r[sa[i-1]]+(str[sa[i]]!=str[sa[i-1]]);
for(int l=1; r[sa[n]]<n; l<<=1)
{
memset(u,0,sizeof u);
memset(v,0,sizeof v);
memcpy(o,r,sizeof r);
for(int i=1; i<=n; i++)
u[r[i]]++, v[r[i+l]]++;
for(int i=1; i<=n; i++)
u[i]+=u[i-1], v[i]+=v[i-1];
for(int i=n; i>=1; i--)
y[v[r[i+l]]--]=i;
for(int i=n; i>=1; i--)
sa[u[r[y[i]]]--]=y[i];
r[sa[1]]=1;
for(int i=2; i<=n; i++)
r[sa[i]]=r[sa[i-1]]+((o[sa[i]]!=o[sa[i-1]])||(o[sa[i]+l]!=o[sa[i-1]+l]));
}
{
int i,j,k=0;
for(int i=1; i<=n; h[r[i++]]=k)
for(k?k--:0,j=sa[r[i]-1]; str[i+k]==str[j+k]; k++);
}
ans=(long long)n*(long long)(n+1)/(long long)2;
for(int i=1; i<=n; i++)
ans-=(long long)h[i];
an.push_back(ans);
for(int i=1; i<=n; i++)
s.insert(i);
for(int i=1; i<=n; i++)
{
int p=r[i],u=0,v=0;
set<int>::iterator it,it1,it2;
it=s.find(p);
it1=it;
it2=it;
if(it1!=s.begin())
{
--it1;
u=*it1;
}
if(it2!=s.end())
{
++it2;
if(it2!=s.end())
v=*it2;
}
int tmp=max(h[p],h[v]);
ans -= n-i+1 - tmp;
h[v]=min(h[p],h[v]);
s.erase(it);
an.push_back(ans);
}
for(int i=n-1; i>=0; --i)
printf("%lld\n",an[i]);
}
[SDOI2016] 生成魔咒 - 后缀数组,平衡树,STL,时间倒流的更多相关文章
- 【bzoj4516】[Sdoi2016]生成魔咒 后缀数组+倍增RMQ+STL-set
题目描述 魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1.2 拼凑起来形成一个魔咒串 [1,2].一个魔咒串 S 的非空字串被称为魔咒串 S 的生成魔咒. 例如 S=[1,2 ...
- BZOJ 4516: [Sdoi2016]生成魔咒——后缀数组、并查集
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4516 题意 一开始串为空,每次往串后面加一个字符,求本质不同的子串的个数,可以离线.即长度为 ...
- BZOJ.4516.[SDOI2016]生成魔咒(后缀数组 RMQ)
题目链接 后缀自动机做法见这(超好写啊). 后缀数组是可以做的: 本质不同的字符串的个数为 \(子串个数-\sum_{ht[i]}\),即 \(\frac{n(n+1)}{2}-\sum_{ht[i] ...
- BZOJ 4516: [Sdoi2016]生成魔咒(后缀数组)
传送门 解题思路 题目其实就是动态维护本质不同的串的个数.考虑到只有加数字的操作,所以可以用后缀数组.题目是每次往后加数字,这样不好处理,因为每次加数字之后所有的后缀都会改变.所以要转化一下思路,就是 ...
- BZOJ4516: [Sdoi2016]生成魔咒(后缀数组 set RMQ)
题意 题目链接 Sol 毒瘤SDOI 终于有一道我会做的题啦qwq 首先,本质不同的子串的个数 $ = \frac{n(n + 1)}{2} - \sum height[i]$ 把原串翻转过来,每次就 ...
- BZOJ4516: [Sdoi2016]生成魔咒 后缀自动机
#include<iostream> #include<cstdio> #include<cstring> #include<queue> #inclu ...
- BZOJ 4516: [Sdoi2016]生成魔咒 [后缀自动机]
4516: [Sdoi2016]生成魔咒 题意:询问一个字符串每个前缀有多少不同的子串 做了一下SDOI2016R1D2,题好水啊随便AK 强行开map上SAM 每个状态的贡献就是\(Max(s)-M ...
- liberOJ #2033. 「SDOI2016」生成魔咒 后缀数组
#2033. 「SDOI2016」生成魔咒 题目描述 魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1 11.2 22 拼凑起来形成一个魔咒串 [1,2] [1, 2] ...
- [bzoj4516][Sdoi2016]生成魔咒——后缀自动机
Brief Description 魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1.2 拼凑起来形成一个魔咒串 [1,2]. 一个魔咒串 S 的非空字串被称为魔咒串 S 的生 ...
随机推荐
- sql关系型运算符优先级高到低为:not >and> or
今天在做项目的时候发现一个查询的结果不太对. 随后拿出sql仔细端详一番,where条件中发现一个条件本应该是 …… xx in (‘13’,‘14’)……,却写成了…… xx = ‘13’ or x ...
- Qt实践-基于C++的UI设计基础
学习网址:http://c.biancheng.net/view/1824.html 我的代码: 头文件qwdlgmanual.h: #ifndef QWDLGMANUAL_H #define QWD ...
- PP: Extracting statisticla graph features for accurate and efficient time series classification
Problem: TSC, time series classification; Traditional TSC: find global similarities or local pattern ...
- echarts 【图表的基本使用】
一.柱状图 <!DOCTYPE html> <html> <head> <meta charset="utf-8"> <tit ...
- 一个margin就可以让块状元素响应居中,很实用
之前总结过水平居中的很多方法,但今天在<css世界>这本书里看到margin的一个特性,一行代码就搞定很实用,分享一下 margin: auto能在块级元素设定宽高之后自动填充剩余宽高.m ...
- mysql数据库环境配置及部分问题
亲身经历了MySQL初学者的痛苦,把主要问题和解决方法整理一下. 一.解压版环境配置 1.把压缩包解压到某盘符下. 解压后在类似“E:\mysql-5.7.22-winx64”这个文件夹中可以看到以上 ...
- C语言 typedef struct _STUDENT {}STUDENT,*PSTUDENT;
#include <stdio.h> #include <stdlib.h> #include <string.h> //给stuct _STUDENT 起一个别名 ...
- win10系统家庭版升级到专业版
win10家庭版升级专业版密钥:VK7JG-NPHTM-C97JM-9MPGT-3V66T4N7JM-CV98F-WY9XX-9D8CF-369TT FMPND-XFTD4-67FJC-HDR8C-3 ...
- JavaScript 13 Ajax技术(未完)
<body> <!-- 添加文档主体内容 --> <header> <nav>JavaScript - Ajax - 读取XML文件</nav&g ...
- js监听textarea 内容的变化,并计算内容的长度c
html代码如下: <div class="customer-message"> <label for="customerMessage&quo ...