bzoj3173: [Tjoi2013]最长上升子序列(fhqtreap)
这题用fhqtreap可以在线。
fhqtreap上维护以i结尾的最长上升子序列,数字按从小到大加入, 因为前面的数与新加入的数无关, 后面的数比新加入的数小, 所以新加入的数对原序列其他数的值没有影响。
然后就可以直接fhqtreap模拟啦
忘记所有答案取max WA了半天T T
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#define lt tree[x].ls
#define rt tree[x].rs
using namespace std;
const int maxn=, seed=, inf=1e9+;
struct treap{int sum, mx, rnd, size, ls, rs;}tree[maxn];
int n, x, tott, root, ans;
void read(int &k)
{
int f=; k=; char c=getchar();
while(c<'' || c>'') c=='-' && (f=-), c=getchar();
while(c<='' && c>='') k=k*+c-'', c=getchar();
k*=f;
}
inline void build(int &x)
{
tree[x=++tott].rnd=rand()<<|rand();
tree[x].sum=tree[x].mx=tree[x].size=;
}
inline int max(int a, int b){return a>b?a:b;}
inline void up(int x)
{
if(!x) return;
tree[x].size=tree[lt].size+tree[rt].size+;
tree[x].mx=max(tree[x].sum, max(tree[lt].mx, tree[rt].mx));
}
void split(int x, int &l, int &r, int k)
{
if(!k) l=, r=x;
else if(k==tree[x].size) l=x, r=;
else if(k<=tree[lt].size) r=x, split(lt, l, lt, k), up(x);
else l=x, split(rt, rt, r, k-tree[lt].size-), up(x);
}
void merge(int &x, int l, int r)
{
if(!l || !r) x=l+r;
else if(tree[l].rnd<tree[r].rnd) x=l, merge(rt, rt, r), up(x);
else x=r, merge(lt, l, lt), up(x);
}
inline int solve(int pos)
{
int x, y, tmp;
build(tmp); split(root, x, y, pos);
tree[tmp].sum=tree[tmp].mx=max(, tree[x].mx)+;
merge(x, x, tmp); merge(root, x, y);
return ans=max(ans, tree[tmp].sum);
}
int main()
{
srand(seed); read(n); tree[].mx=-inf; tree[].rnd=inf;
for(int i=;i<=n;i++) read(x), printf("%d\n", solve(x));
}
貌似还可以用BIT, 过会再补
先把所有数设成1,1表示<=x的数,0表示>x的数,然后从n到1扫一遍,对于每个要求的位置pos,就是找到前缀为pos的位置插入,至于怎么找到前缀为pos的位置,就用BIT求第k大的方法就好了。
求出原序列之后, 求每一个数结尾的LIS就好了
注意:BIT求第k大的时候add是单点修改
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
using namespace std;
const int maxn=;
int n, mx, ans, x;
int tree1[maxn], tree2[maxn], pos[maxn], p[maxn], f[maxn];
void read(int &k)
{
int f=; k=; char c=getchar();
while(c<'' || c>'') c=='-' && (f=-), c=getchar();
while(c<='' && c>='') k=k*+c-'', c=getchar();
k*=f;
}
inline int max(int a, int b){return a>b?a:b;}
inline void add1(int x, int delta){for(;x<=mx;x+=x&-x) tree1[x]+=delta;}
inline void add2(int x, int delta){for(;x<=n;x+=x&-x) tree2[x]=max(tree2[x], delta);}
inline int query2(int x){int sum=; for(;x;x-=x&-x) sum=max(sum, tree2[x]); return sum;}
int main()
{
read(n);
for(mx=;mx<n;mx<<=);
for(int i=;i<=mx;i++)
tree1[i]+=(i<=n), i+(i&-i)<=mx&&(tree1[i+(i&-i)]+=tree1[i]);
for(int i=;i<=n;i++) read(p[i]), p[i]++;
for(int i=n;i;i--)
{
int l=, r=mx, k=;
while(l<r)
{
int mid=(l+r)>>;
k+=tree1[mid];
if(k<p[i]) l=mid+;
else k-=tree1[mid], r=mid;
}
pos[l]=i;
add1(l, -);
}
for(int i=, tmp;i<=n;i++)
f[pos[i]]=query2(pos[i]-)+, add2(pos[i], f[pos[i]]);
for(int i=;i<=n;i++) ans=max(ans, f[i]), printf("%d\n", ans);
}
bzoj3173: [Tjoi2013]最长上升子序列(fhqtreap)的更多相关文章
- [BZOJ3173][Tjoi2013]最长上升子序列
[BZOJ3173][Tjoi2013]最长上升子序列 试题描述 给定一个序列,初始为空.现在我们将1到N的数字插入到序列中,每次将一个数字插入到一个特定的位置.每插入一个数字,我们都想知道此时最长上 ...
- bzoj3173[Tjoi2013]最长上升子序列 平衡树+lis
3173: [Tjoi2013]最长上升子序列 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2253 Solved: 1136[Submit][S ...
- bzoj3173: [Tjoi2013]最长上升子序列(树状数组+二分倒推)
3173: [Tjoi2013]最长上升子序列 题目:传送门 题解: 好题! 怎么说吧...是应该扇死自己...看错了两次题: 每次加一个数的时候,如果当前位置有数了,是要加到那个数的前面,而不是直 ...
- BZOJ3173 TJOI2013最长上升子序列(Treap+ZKW线段树)
传送门 Description 给定一个序列,初始为空.现在我们将1到N的数字插入到序列中,每次将一个数字插入到一个特定的位置.每插入一个数字,我们都想知道此时最长上升子序列长度是多少? Input ...
- bzoj千题计划316:bzoj3173: [Tjoi2013]最长上升子序列(二分+树状数组)
https://www.lydsy.com/JudgeOnline/problem.php?id=3173 插入的数是以递增的顺序插入的 这说明如果倒过来考虑,那么从最后一个插入的开始删除,不会对以某 ...
- BZOJ3173:[TJOI2013]最长上升子序列(Splay)
Description 给定一个序列,初始为空.现在我们将1到N的数字插入到序列中,每次将一个数字插入到一个特定的位置.每插入一个数字,我们都想知道此时最长上升子序列长度是多少? Input 第一行一 ...
- BZOJ3173 TJOI2013最长上升子序列(splay)
容易发现如果求出最后的序列,只要算一下LIS就好了.序列用平衡树随便搞一下,这里种一棵splay. #include<iostream> #include<cstdio> #i ...
- 【LG4309】【BZOJ3173】[TJOI2013]最长上升子序列
[LG4309][BZOJ3173][TJOI2013]最长上升子序列 题面 洛谷 BZOJ 题解 插入操作显然用平衡树就行了 然后因为后面的插入对前面的操作无影响 就直接在插入完的序列上用树状数组求 ...
- BZOJ 3173: [Tjoi2013]最长上升子序列
3173: [Tjoi2013]最长上升子序列 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1524 Solved: 797[Submit][St ...
随机推荐
- 「日常训练」 Fire!(UVA-11624)
与其说是训练不如说是重温.重新写了Java版本的代码. import java.util.*; import java.math.*; import java.io.BufferedInputStre ...
- Python爬虫与反爬虫(7)
[Python基础知识]Python爬虫与反爬虫(7) 很久没有补爬虫了,相信在白蚁二周年庆的活动大厅比赛中遇到了关于反爬虫的问题吧 这节我会做个基本分享. 从功能上来讲,爬虫一般分为数据采集,处理, ...
- eclipse查看源代码问题
最近分析源代码时,eclipse总是出错,显示org.eclipse.core.runtime.CoreException,解决方法:在builderpath点击 add external jars, ...
- Codeforces Round #504 (rated, Div. 1 + Div. 2, based on VK Cup 2018 Final) E. Down or Right
从(1,1,n,n)每次只变一个坐标,进行询问. 如果问到对角线有距离限制, 再从(1,1,n/2,n/2)询问到(n/2,n/2,n,n) 记住前半部分贪心忘上走,后本部分贪心往右走 因为最后的路线 ...
- mysql数据导到本地
需求: 把mysql查询结果导出到txt(其他格式亦可),放在本地,供下一步使用 首先网上查了下,select * from driver into outfile 'a.txt'; 前面是你的sql ...
- 使用Firebug或chrome-devToolBar深入学习javascript语言核心
使用Firebug和chrome-devToolBar调试页面样式或脚本是前端开发每天必做之事.这个开发神器到底能给我们带来哪些更神奇的帮助呢?这几天看的一些资料中给了我启发,能不通过Firebug和 ...
- MySQL双主复制
原文发表于cu:2017-06-12 本文简单介绍MySQL双主复制原理及1个简单是双主复制验证. 一.MySQL双主复制原理 1. 双主复制原理 master-master复制的两台服务器,既是ma ...
- libCurl 初步认识 - cur easy
cur easy接口简洁明了,主接口4个,辅接口5个. 主接口 初始化 + 配参数 + 执行 + 销毁 初始化 CURL* curl_easy_init() 获得CURL句柄,返回值需要判空. 配参数 ...
- day04 list tuple (补)
今日内容: 1. 列表 2. 列表的增删改查 3. 列表的嵌套 4. 元组和元组嵌套 5. range 列表 列表: 能装对象的对象. 有顺序的(按照我们添加的顺序保存) 在代码中使用[]表示列表. ...
- ES6的新特性(1)——ES6 的概述
ES6 的概述 首先,感谢马伦老师的ES6新特性的教程. ECMAScript 和 JavaScript 的关系是 ECMAScript 和 JavaScript 的关系是,前者是后者的规格,后者是前 ...