五月月赛 寻宝 exkmp + 主席树
: 寻宝
时间限制: Sec 内存限制: MB
提交: 解决:
[提交] [状态] [讨论版] [命题人:admin]
题目描述
采蘑菇的小西佬找到了一张上古年间的藏宝图,上面画着m座连绵不断的山,他决定去地图上记载的地点探险,可当他到达时,他发现当地其实有n座山,并且由于年代久远有些山也改变了形状,他不能准确的找到地图上藏宝的地点,于是他决定从一些连续的山中随意选择一座山作为地图上的第一个点,当他选定了起点以后,若当前的山与地图上记载山的高度相同他就会继续前往下一座山,否则就停下来,直到他到达地图上的最后一座山位置。我们假设如果他经过了大于k座山他就能找到宝藏(如果第一座山高度就不相同,视为经过0座山),现在的问题就是他有多大的可能性找到宝藏。
题目共有q次询问,每次询问有两个整数l,r代表采蘑菇的小西佬选择的区间,并且会给出第一次询问的k,之后的每个询问的k按照下面的式子算出
k=(ans*10086+666)%x; ans代表上一次询问的答案,x=min(m,r-l+1); 输入
第一行输入三个数n,m,k(<=n,m,k<=1e5); 第二行输入n个数代表n座山的高度(<=a[i]<=1e9) 第二行输入m个数代表m座山的高度(<=b[i]<=1e9) 第四行输入一个数q;(<=q<=1e5) 接下来每行输入两个整数l,r;(<=l<=r<=n) 输出
对于每个询问请输出对应的概率,以逆元的形式输出(mod=1e9+7);
题面
用exkmp预处理出每个点可以向后面跑的长度,由于k是变化的,所以需要主席树来回答。
#include <bits/stdc++.h> using namespace std;
#define pb push_back
#define fi first
#define se second
#define debug(x) cerr<<#x << " := " << x << endl;
#define bug cerr<<"-----------------------"<<endl;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
//const int inf = 0x3f3f3f3f;
const int mod = 1e9+;
/**********showtime************/ const int N = 1e5 + ;
int nxt[N], ex[N];
int S[N], P[N];
int a[N];
int n,m,k;
void GETNEXT() {
int i = , j, po, len=m;
nxt[] = len;
while(P[i] == P[i+] && i+ < len) i++;
nxt[] = i;
po = ;
for(i = ; i < len; i++) {
if(nxt[i-po] + i < nxt[po] + po)
nxt[i] = nxt[i-po];
else {
j=nxt[po] + po - i;
if(j < ) j = ;
while(i + j < len && P[j] == P[j+i])
j++;
nxt[i] = j;
po = i;
}
}
}
void EXKMP()
{
int i = , j, po, len = n, l2 = m;
while(S[i] == P[i] && i < l2 && i < len) i++;
ex[] = i;
po = ;
for(i = ; i < len; i++)
{
if(nxt[i-po] + i < ex[po] + po) ex[i]=nxt[i-po];
else {
j = ex[po] + po - i;
if(j < ) j = ;
while(i + j < len && j < l2 && S[j+i] == P[j]) j++;
ex[i] = j;
po = i;
}
}
} struct node{
int le, ri;
int cnt;
} tree[N * ];
int tot = ;
int root[N];
int build(int le, int ri) {
int p = ++tot;
if(le == ri) {
tree[p].cnt = ;
return p;
}
int mid = (le + ri) >> ;
tree[p].le = build(le, mid);
tree[p].ri = build(mid+, ri);
tree[p].cnt = tree[tree[p].le].cnt + tree[tree[p].ri].cnt;
return p;
}
int ins(int now, int le, int ri, int x, int val) {
int p = ++tot;
tree[p] = tree[now];
if(le == ri) {
tree[p].cnt += val;
return p;
}
int mid = (le + ri) >> ;
if(x <= mid) {
tree[p].le = ins(tree[now].le, le, mid, x, val);
}
else {
tree[p].ri = ins(tree[now].ri, mid+, ri, x, val);
}
tree[p].cnt = tree[tree[p].le].cnt + tree[tree[p].ri].cnt;
return p;
}
int query(int L, int R, int le, int ri, int rt1, int rt2) {
if(L > R) return ;
if(le >= L && ri <= R) {
return tree[rt2].cnt - tree[rt1].cnt;
} int res = ;
int mid = (le + ri) >> ;
if(mid >= L) res += query(L, R, le, mid, tree[rt1].le, tree[rt2].le);
if(mid < R) res += query(L, R, mid+, ri, tree[rt1].ri, tree[rt2].ri);
return res;
}
ll ksm(ll a, ll n) {
ll res = ;
while(n > ) {
if(n & ) res = res * a % mod;
a = a * a % mod;
n = n>>;
}
return res;
}
int main(){
scanf("%d%d%d", &n, &m, &k);
for(int i=; i<n; i++) scanf("%d", &S[i]);
for(int i=; i<m; i++) scanf("%d", &P[i]); GETNEXT();
EXKMP();
root[] = build(, n); for(int i=; i < n; i++){
root[i+] = ins(root[i], , n, ex[i]+, );
} int q; scanf("%d", &q);
int ans = ;
for(int i=; i<=q; i++) { int l, r;
scanf("%d%d", &l, &r);
if(i > ) {
int x = min(m, r - l + );
k = (1ll*ans * + ) % x;
}
int c = query(k+, n, , n, root[l-], root[r]);
// cout<<c << " " << k << endl;
ans = 1ll*c * ksm(r - l + , mod-) % mod;
printf("%d\n", ans);
} return ;
}
五月月赛 寻宝 exkmp + 主席树的更多相关文章
- BZOJ5361[Lydsy1805月赛]对称数——主席树+随机化
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=5361 好神的一道题啊! 容易看出来是要用维护权值的数据结构,因此树链剖分首先pass掉. ...
- bzoj3207--Hash+主席树
题目大意: 给定一个n个数的序列和m个询问(n,m<=100000)和k,每个询问包含k+2个数字:l,r,b[1],b[2]...b[k],要求输出b[1]~b[k]在[l,r]中是否出现. ...
- bzoj1901--树状数组套主席树
树状数组套主席树模板题... 题目大意: 给定一个含有n个数的序列a[1],a[2],a[3]--a[n],程序必须回答这样的询问:对于给定的i,j,k,在a[i],a[i+1],a[i+2]--a[ ...
- BZOJ 3626: [LNOI2014]LCA [树链剖分 离线|主席树]
3626: [LNOI2014]LCA Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2050 Solved: 817[Submit][Status ...
- BZOJ 1146: [CTSC2008]网络管理Network [树上带修改主席树]
1146: [CTSC2008]网络管理Network Time Limit: 50 Sec Memory Limit: 162 MBSubmit: 3522 Solved: 1041[Submi ...
- BZOJ 2588: Spoj 10628. Count on a tree [树上主席树]
2588: Spoj 10628. Count on a tree Time Limit: 12 Sec Memory Limit: 128 MBSubmit: 5217 Solved: 1233 ...
- BZOJ 1901: Zju2112 Dynamic Rankings[带修改的主席树]【学习笔记】
1901: Zju2112 Dynamic Rankings Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 7143 Solved: 2968[Su ...
- [bzoj3932][CQOI2015][任务查询系统] (主席树)
Description 最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分.超级计算机中的 任务用三元组(Si,Ei,Pi)描述,(Si,Ei,Pi)表示任务从第Si ...
- [bzoj2588][count on a tree] (主席树+lca)
Description 给定一棵N个节点的树,每个点有一个权值,对于M个询问(u,v,k),你需要回答u xor lastans和v这两个节点间第K小的点权.其中lastans是上一个询问的答案,初始 ...
随机推荐
- Windows 使用 helm3 和 kubectl
简介: 主要原因是,我不会 vim ,在 linux 上修改 charts 的很蹩脚,所以就想着能不能再 windows 上执行 helm 命令,将 charts install linux 上搭建的 ...
- case和decode的用法(行转列)
创建了一张成绩表,如下图所示: 在oracle中,这两个函数我们都可以使用,代码及结果如下: decode用法: select Name,decode(Subject,'语文',1,'数学',2,'英 ...
- Could not launch "APP_NAME" process launch failed: 4294967295
真机调试忽然遇到这个问题, Could not launch "APP_NAME" process launch failed: 如图所示: 模拟器上能正常调试………… 这个问题还 ...
- Linux 文件系统相关的基本概念
本文介绍 Linux 文件系统相关的基本概念. 硬盘的物理组成 盘片硬盘其实是由单个或多个圆形的盘片组成的,按照盘片能够容纳的数据量,分为单盘(一个硬盘里面只有一个盘片)或多盘(一个硬盘里面有多个盘片 ...
- 实用小工具推荐 OpenWrite
[实用小工具推荐]给技术同学们推荐一款比较好用的工具,可以实现一稿多发,主流的技术渠道基本涵盖了:https://www.openwrite.cn/ 因为工作的关系,认识了很多做技术公众号的小伙伴,同 ...
- bit、byte、kb、mb、g的区别
1Byte=8bit1KB=1024Byte(字节)=8*1024bit1MB=1024KB1GB=1024MB1TB=1024GB bit是计算机数据的最小单元.要么是0,要么是1. byte 关键 ...
- HashMap、Hash Table、ConcurrentHashMap
这个这个...本王最近由于开始找实习工作了,所以就在牛客网上刷一些公司的面试题,大多都是一些java,前端HTML,js,jquery,以及一些好久没有碰的算法题,说实话,有点难受,其实在我不知道的很 ...
- 【JDK】JDK源码分析-ReentrantLock
概述 在 JDK 1.5 以前,锁的实现只能用 synchronized 关键字:1.5 开始提供了 ReentrantLock,它是 API 层面的锁.先看下 ReentrantLock 的类签名以 ...
- Netty学习(十)-Netty文件上传
今天我们来完成一个使用netty进行文件传输的任务.在实际项目中,文件传输通常采用FTP或者HTTP附件的方式.事实上通过TCP Socket+File的方式进行文件传输也有一定的应用场景,尽管不是主 ...
- 大话 Spring Session 共享
javaweb中我们项目稍微正规点,都会用到单点登录这个技术.实现它的方法各家有各界的看法.这几天由于公司项目需求正在研究.下面整理一下最近整理的心得. 简介 在分布式项目中另我们头疼的是多项目之间的 ...