jzoj 3467. 【NOIP2013模擬聯考7】最長上升子序列(lis) dfs+lis+手工棧
阿新 • • 發佈:2019-02-11
Input
輸入檔案lis.in的第一行有一個正整數n,表示操作個數。接下來n行每行有兩個整數op,x。如果op為0,則表示新增x這個數字;如果op為1,則表示回到第x次操作之後。
Output
對於每次操作,在輸出檔案lis.out中輸出一個答案,表示當前最長上升子序列的長度
Data Constraint
30%的資料 n<=1000
另外20%的資料沒有第二個操作
80%的資料 n<=200000
100%的資料 n<=500000且所有輸入的數字都是長整型範圍內的非負整數
Solution
部分分給的很良心啊,感覺這才是模擬賽應該有的樣子
可以發現每次字首都是不變的,把它看成一顆樹,也就是說0操作等同於拓展一個結點,1操作等同於在x處拓展一個新的兒子。那麼我們把操作離線下來dfs一遍順便統計lis就可以了。
注意直接遞迴會爆棧要手動模擬,還有就是INF的大小問題
手工棧賊!難!打!
Code
#include <stdio.h>
#include <string.h>
#include <stack>
#define rep(i,st,ed) for (int i=st;i<=ed;i++)
#define fill(x,t) memset(x,t,sizeof(x))
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define INF 2147483647
#define N 1000015
#define E 1000015
struct edge{int y,next;}e[E];
struct data{int x,now,tmp,rec;};
std::stack<int>stack;
int ans[N],a[N],t[N],q[N],f[N],s[N],mx;
int cur[N],ls[N],edgeCnt=0;
int tmp[N],rec[N];
bool vis[N];
void addEdge(int x,int y) {
e[++edgeCnt]=(edge){y,ls[x]}; cur[x]=ls[x]=edgeCnt;
}
int find(int x) {
int l=0,r=mx;
while (l<=r) {
int mid=(l+r)>>1;
if (s[mid]>=a[x]) r=mid-1;
else l=mid+1;
}
return l;
}
/*
void dfs(int x) {
t[++t[0]]=a[x];
int tmp=s[f[t[0]]=find(x)];
int rec=mx;
mx=max(f[t[0]],mx);
s[f[t[0]]]=min(s[f[t[0]]],a[x]);
ans[x]=mx;
for (int i=ls[x];i;i=e[i].next) {
dfs(e[i].y);
}
s[f[t[0]]]=tmp;
mx=rec;
t[t[0]--]=0;
}
*/
void DFS() {
t[0]=0;
fill(vis,0);
fill(f,0);
stack.push(0);
mx=0;
while (!stack.empty()) {
int x=stack.top();
if (!vis[x]) {
t[++t[0]]=a[x];
tmp[stack.size()]=s[f[t[0]]=find(x)];
rec[stack.size()]=mx;
mx=max(mx,f[t[0]]);
s[f[t[0]]]=min(s[f[t[0]]],a[x]);
// printf("%d\n", s[mx]);
ans[x]=mx;
vis[x]=1;
}
if (!cur[x]) {
s[f[t[0]]]=tmp[stack.size()];
mx=rec[stack.size()];
t[t[0]--]=0;
stack.pop();
} else {
stack.push(e[cur[x]].y);
cur[x]=e[cur[x]].next;
}
}
}
int main(void) {
freopen("lis.in","r",stdin);
freopen("lis.out","w",stdout);
int n,now=0,tail=0;
scanf("%d",&n);
s[0]=0;
rep(i,1,n) {
s[i]=INF;
int opt,x;
scanf("%d%d",&opt,&x);
if (opt==0) addEdge(now,++tail),a[now=tail]=x;
else now=q[x];
q[i]=now;
}
rep(i,n+1,N-1) s[i]=INF;
// dfs(0);
DFS();
rep(i, 1, n) {
printf("%d\n", ans[q[i]]-1);
}
return 0;
}