【牛客NOIP模拟】牛半仙的妹子序列【DP】【Segment Tree Beats】

题意：排列的极长上升子序列个数

$n\le 2\times 10^5$

显然有个 dp

$$f_n=\sum _{i<n,p_i<p_n,\nexists i<j<n\text{ s.t. }{p}_i<p_j<p_n}{f}_i$$

相当于要对 $i<n,p_i<p_n$ 维护一个单调栈。

这样并不好搞，但我们换个角度，每个 $i$ 产生贡献的 $n$ 也是关于 $p_i$ 的单调栈的形式。

（仔细想想好像不是单调栈，感性理解好了）

考虑怎么维护这东西。用一棵权值线段树来维护当前每个位置的栈顶的值。

添加一个数 $v$ 时，视为在所有 $p_i<v$ 的栈 $i$ 中压入了 $v$，相当于是对 $[1,v]$ 取 $\min$ 。

然后栈顶值为 $v$ 的位置可以产生贡献，用最大值判一下就可以了。

最后这个位置新开一个栈，即把 $v$ 改成 $+\infty$，并把贡献放到树上。

用 SGB 维护即可。

因为只有单点修改，所以复杂度为 $O(n\log n)$

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#define MAXN 200005
using namespace std;
const int MOD=998244353,INF=0x3f3f3f3f;
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
inline int read()
{
	int ans=0;
	char c=getchar();
	while (!isdigit(c)) c=getchar();
	while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
	return ans;
}
int f[MAXN];
#define lc p<<1
#define rc p<<1|1
struct node{int mx,se,sum;}t[MAXN<<2];
int lzy[MAXN<<2];
inline node merge(node a,node b)
{
	if (a.mx==b.mx)	a.se=max(a.se,b.se),a.sum=add(a.sum,b.sum);
	else
	{
		if (a.mx<b.mx) swap(a,b);
		a.se=max(a.se,b.mx);
	}
	return a;
}
inline void pushlzy(int p,int v){t[p].mx=min(t[p].mx,v),lzy[p]=min(lzy[p],v);}
inline void pushdown(int p)
{
	if (lzy[p]<INF)
	{
		pushlzy(lc,lzy[p]),pushlzy(rc,lzy[p]);
		lzy[p]=INF;
	}
}
void modify(int p,int l,int r,int k,int v)
{
	if (l==r) return (void)(t[p].mx=INF,t[p].sum=v);
	int mid=(l+r)>>1;
	pushdown(p);
	if (k<=mid) modify(lc,l,mid,k,v);
	else modify(rc,mid+1,r,k,v);
	t[p]=merge(t[lc],t[rc]);
}
void modify(int p,int l,int r,int ql,int qr,int v)
{
	if (ql<=l&&r<=qr&&v>t[p].se) return pushlzy(p,v);
	if (qr<l||r<ql) return;
	int mid=(l+r)>>1;
	pushdown(p);
	modify(lc,l,mid,ql,qr,v),modify(rc,mid+1,r,ql,qr,v);
	t[p]=merge(t[lc],t[rc]);
}
node ans;
void query(int p,int l,int r,int ql,int qr)
{
	if (ql<=l&&r<=qr)
	{
		if (ql==l) ans=t[p];
		else ans=merge(ans,t[p]);
		return;
	}
	if (qr<l||r<ql) return;
	pushdown(p);
	int mid=(l+r)>>1;
	query(lc,l,mid,ql,qr),query(rc,mid+1,r,ql,qr);
}
int main()
{
	memset(lzy,0x3f,sizeof(lzy));
	int n=read();
	for (int i=1;i<=n;i++)
	{
		int v=read();
		modify(1,1,n,1,v,v),query(1,1,n,1,v);
		f[i]=(ans.mx==v? ans.sum:0);
		if (!f[i]) f[i]=1;
		modify(1,1,n,v,f[i]);
	}
	printf("%d\n",t[1].sum);
	return 0;
}