转博客大法好

自己画一画看一看,就会体会到这个设置关键点的强大之处了.

CODE(sa)

O(nlogn)→1436msO(nlogn)\to 1436msO(nlogn)→1436ms

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
template<class T>inline void read(T &num) {
register char ch; register int flg = 1;
while(!isdigit(ch=getchar()))if(ch=='-')flg=-flg;
for(num=0; isdigit(ch); num=num*10+ch-'0', ch=getchar());
num *= flg;
}
const int MAXN = 5e4+5; int x[MAXN], y[MAXN], c[MAXN], Log[MAXN]; int n, s[MAXN], b[MAXN], tot;
struct SA {
int sa[MAXN], rk[MAXN], h[MAXN], f[MAXN][16];
inline void Get_Sa(int m) {
for(int i = 1; i <= m; ++i) c[i] = 0;
for(int i = 1; i <= n; ++i) ++c[x[i]=s[i]];
for(int i = 2; i <= m; ++i) c[i] += c[i-1];
for(int i = n; i >= 1; --i) sa[c[x[i]]--] = i;
for(int k = 1; k <= n; k<<=1) {
int p = 0;
for(int i = n-k+1; i <= n; ++i) y[++p] = i;
for(int i = 1; i <= n; ++i) if(sa[i]>k) y[++p] = sa[i]-k;
for(int i = 1; i <= m; ++i) c[i] = 0;
for(int i = 1; i <= n; ++i) ++c[x[i]];
for(int i = 2; i <= m; ++i) c[i] += c[i-1];
for(int i = n; i >= 1; --i) sa[c[x[y[i]]]--] = y[i];
swap(x, y);
x[sa[1]] = 1; p = 1;
for(int i = 2; i <= n; ++i)
x[sa[i]] = (y[sa[i]] == y[sa[i-1]] && y[sa[i]+k] ==y[sa[i-1]+k]) ? p : ++p;
if((m=p) == n) break;
}
}
inline void Get_Height() {
int k = 0;
for(int i = 1; i <= n; ++i) rk[sa[i]] = i;
for(int i = 1; i <= n; ++i) {
if(rk[i] == 1) continue;
if(k) --k;
int j = sa[rk[i]-1];
while(i+k <= n && j+k <= n && s[i+k] == s[j+k]) ++k;
h[rk[i]] = k;
}
}
inline void rev() {
reverse(rk + 1, rk + n + 1);
for(int i = 1; i <= n; ++i) sa[i] = n-i+1;
}
inline void init() {
Get_Sa(n);
Get_Height();
}
inline void initST() {
for(int i = 1; i <= n; ++i) f[i][0] = h[i];
for(int j = 1; j <= Log[n]; ++j)
for(int i = 1; i <= n-(1<<j)+1; ++i)
f[i][j] = min(f[i][j-1], f[i+(1<<j-1)][j-1]);
}
inline int query(int l, int r) {
l = rk[l], r = rk[r];
if(l > r) swap(l, r);
++l;
int k = Log[r-l+1];
return min(f[l][k], f[r-(1<<k)+1][k]);
}
}sa1, sa2;
int m;
int main() {
read(n), read(m);
for(int i = 1; i <= n; ++i) read(s[i]);
--n;
for(int i = 1; i <= n; ++i) b[++tot] = (s[i]=s[i+1]-s[i]);
s[n+1] = 0;
sort(b + 1, b + tot + 1); tot = unique(b + 1, b + tot + 1) - b - 1;
for(int i = 1; i <= n; ++i) s[i] = lower_bound(b + 1, b + tot + 1, s[i]) - b; for(int i = 2; i <= n; ++i) Log[i] = Log[i>>1] + 1;
sa1.init(), sa1.initST();
reverse(s + 1, s + n + 1);
sa2.init(), sa2.rev(), sa2.initST();
LL ans = 0;
for(int len = 1; (len<<1)+m <= n; ++len)
for(int i = 1, j; (j=i+len+m) <= n; i += len) {
int l = min(sa2.query(i, j), len), r = min(sa1.query(i, j), len);
if(l + r > len) ans += l + r - len;
}
printf("%lld\n", ans);
}

CODE(hash)

O(nlog2n)→688msO(nlog^2n)\to 688msO(nlog2n)→688ms

实测hash速度是sa的2.5倍…

常数的幽怨…

市面上貌似找不到hash的代码的说…

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
template<class T>inline void read(T &num) {
register char ch; register int flg = 1;
while(!isdigit(ch=getchar()))if(ch=='-')flg=-flg;
for(num=0; isdigit(ch); num=num*10+ch-'0', ch=getchar());
num *= flg;
}
const int MAXN = 5e4+5;
const int p = 137;
int n, m, s[MAXN], h[MAXN], mul[MAXN];
inline int hsh(int l, int r) {
return h[r] - h[l-1]*mul[r-l+1];
}
inline int querypre(int i, int j, int up) {
int l = 0, r = min(min(i, j), up), mid;
while(l < r) {
mid = (l + r + 1) >> 1;
if(hsh(i-mid+1, i) == hsh(j-mid+1, j)) l = mid;
else r = mid-1;
}
return l;
}
inline int querysuf(int i, int j, int up) {
int l = 0, r = min(min(n-i+1, n-j+1), up), mid;
while(l < r) {
mid = (l + r + 1) >> 1;
if(hsh(i, i+mid-1) == hsh(j, j+mid-1)) l = mid;
else r = mid-1;
}
return l;
}
int main() {
read(n), read(m); mul[0] = 1;
for(int i = 1; i <= n; ++i) read(s[i]);
for(int i = 1; i < n; ++i) {
h[i] = h[i-1] * p + s[i+1]-s[i];
mul[i] = mul[i-1] * p;
}
s[n] = 0; --n;
LL ans = 0;
for(int len = 1; (len<<1)+m <= n; ++len)
for(int i = 1, j; (j=i+len+m) <= n; i += len) {
int l = querypre(i, j, len), r = querysuf(i, j, len);
if(l + r > len) ans += l + r - len;
}
printf("%lld\n", ans);
}

BZOJ 2119: 股市的预测 (Hash / 后缀数组 + st表)的更多相关文章

  1. BZOJ 2119 股市的预测(后缀数组)

    首先要差分+离散化. 然后就是求形如ABA的串有多少,其中B的长度确定为k. 我们用到了设置关键点的思想.我们枚举A的长度L.然后在\(1,1+L,1+L*2,1+L*3...\)设置关键点.然后我们 ...

  2. bzoj千题计划312:bzoj2119: 股市的预测(后缀数组+st表)

    https://www.lydsy.com/JudgeOnline/problem.php?id=2119 题意:将给定数组差分后,求ABA形式的字串个数,要求|B|=m,|A|>0 1.后缀数 ...

  3. 【BZOJ 2119】 2119: 股市的预测 (后缀数组+分块+RMQ)

    2119: 股市的预测 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 404  Solved: 188 Description 墨墨的妈妈热爱炒股,她 ...

  4. SPOJ 687 Repeats(后缀数组+ST表)

    [题目链接] http://www.spoj.com/problems/REPEATS/en/ [题目大意] 求重复次数最多的连续重复子串的长度. [题解] 考虑错位匹配,设重复部分长度为l,记s[i ...

  5. POJ 3693 Maximum repetition substring(后缀数组+ST表)

    [题目链接] poj.org/problem?id=3693 [题目大意] 求一个串重复次数最多的连续重复子串并输出,要求字典序最小. [题解] 考虑错位匹配,设重复部分长度为l,记s[i]和s[i+ ...

  6. BZOJ_4516_[Sdoi2016]生成魔咒_后缀数组+ST表+splay

    BZOJ_4516_[Sdoi2016]生成魔咒_后缀数组+ST表+splay Description 魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1.2 拼凑起来形成一个魔 ...

  7. UVA10829 L-Gap Substrings(后缀数组+ST表)

    后缀数组+ST表. 代填的坑. \(Code\ Below:\) #include <bits/stdc++.h> #define ll long long using namespace ...

  8. BZOJ 2119: 股市的预测 [后缀数组 ST表]

    2119: 股市的预测 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 331  Solved: 153[Submit][Status][Discuss ...

  9. BZOJ 2119: 股市的预测 SA

    2119: 股市的预测 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 434  Solved: 200[Submit][Status][Discuss ...

随机推荐

  1. hdoj4812 D Tree(点分治)

    题目链接:https://vjudge.net/problem/HDU-4812 题意:给定一颗带点权的树,求是否存在一条路经的上点的权值积取模后等于k,如果存在多组点对,输出字典序最小的. 思路: ...

  2. [转帖]中芯国际14nm秋季量产 7nm工艺或在2020年底问世

    中芯国际14nm秋季量产 7nm工艺或在2020年底问世 https://news.mydrivers.com/1/641/641087.htm 中芯正在发力.. 今年秋天 14nm两场 明年底 7n ...

  3. Oracle - 函数及多表关联

    函数一般是在数据上执行的,它给数据的转换和处理提供了方便.只是将取出的数据进行处理,不会改变数据库中的值.函数根据处理的数据分为单行函数和聚合函数(组函数),组函数又被称作聚合函数,用于对多行数据进行 ...

  4. winform messageBox.Show()

    MessageBox.Show(" 5 个参数...... ",     " 亮仔提示",     MessageBoxButtons.OKCancel,    ...

  5. 修改git默认的编辑器nano为vim的方法

    git默认的编辑器是nano,使用起来不易操作,下面介绍两种方法将git默认的编辑器修改为vim.  git config --global core.editor vim .git/config文件 ...

  6. ns nat rule

    ns nat rule NAT实现方式: NAT的实现方式有三种,即静态转换(Static Nat).动态转换(Dynamic Nat) 和 端口多路复用(OverLoad). 静态转换是指将内部网络 ...

  7. 点击登录页面成功后,后端返回数据需要保存,在另外一个页面,发送ajax请求的时候需要登录返回数据的其中的一部分当做参数然后拿到新的数据

    对于这个怎么操作首先我们要在登录的ajax请求中把后端的数据保存到sessionstorage中,代码如下 登录ajax $.ajax({ type:'post', url:xxxxxxxxx, da ...

  8. wabacus JaveEE开发框架

    http://www.wabacus.org/ css学习网站:http://www.divcss5.com/rumen/r422.shtml

  9. netty的断线重连问题

    手里的这个项目需要作为客户端,不断的接收服务端发来的数据,用的netty框架,但是一直存在一个问题,就是断线重连问题. 什么是断线重连呢? 就是我们这个客户端要保证一直与服务端保持连接,这样客户端才能 ...

  10. java 读取property

    Properties prop = new Properties(); String path = AlarmController.class.getResource("/").g ...