1. 程式人生 > >jzoj 3467. 【NOIP2013模擬聯考7】最長上升子序列(lis) dfs+lis+手工棧

jzoj 3467. 【NOIP2013模擬聯考7】最長上升子序列(lis) dfs+lis+手工棧

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; }