[bzoj4881][DP]线段游戏

Description

quailty和tangjz正在玩一个关于线段的游戏。在平面上有n条线段，编号依次为1到n。其中第i条线段的两端点坐
标分别为(0,i)和(1,p_i)，其中p_1,p_2,…,p_n构成了1到n的一个排列。quailty先手，他可以选择一些互不相交
的线段，将它们拿走，当然他也可以一条线段也不选。然后tangjz必须拿走所有剩下的线段，若有两条线段相交，
那么他就输了，否则他就赢了。注意若quailty拿走了全部线段，那么tangjz也会胜利。quailty深深喜欢着tangjz
，所以他不希望tangjz输掉游戏，请计算他有多少种选择线段的方式，使得tangjz可以赢得游戏。

Input

第一行包含一个正整数n(1<=n<=100000)，表示线段的个数。
第二行包含n个正整数p_1,p_2,…,p_n(1<=p_i<=n)，含义如题面所述。

Output

输出一行一个整数，即tangjz胜利的方案数，因为答案很大，请对998244353取模输出。

Sample Input

5

1 2 4 5 3

Sample Output

8

题解

我的做法怎么跟别人的都不一样啊…
感觉人家的就是神仙我就是弟弟呀…
首先容易发现，选出的子序列没有交点等价于这个子序列中没有逆序对
实际上问题就变为了要你把这个序列分成两个上升子序列
似乎只有我是傻乎乎地直接做这个问题的…
考虑dp
设$g[i]$表示选了$i$且合法，并且没有被选的子序列的最大值小于$a[i+1]$的方案数
转移可以逆向一下，不枚举选了的子序列的前一个，我们枚举没有被选的子序列的最后一个
不妨设枚举的这个位置是$u$，显然在$[u+1,i]$这一段区间要求单调递增
合法的时候直接累加$g[u]$就行了
然后又发现合法的$u$是连续一段的，维护一个$L$表示$[L,i]$这一段单调递增，显然$u$的取值范围在$[L-1,i-1]$这一段
树状数组做一下就好了…
最后统计答案的时候反向扫一遍
我真是个弟弟

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
#include<map>
#include<bitset>
#include<set>
#define LL long long
#define mp(x,y) make_pair(x,y)
#define pll pair<long long,long long>
#define pii pair<int,int>
using namespace std;
inline int read()
{
	int f=1,x=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int stack[20];
inline void write(LL x)
{
	if(x<0){putchar('-');x=-x;}
    if(!x){putchar('0');return;}
    int top=0;
    while(x)stack[++top]=x%10,x/=10;
    while(top)putchar(stack[top--]+'0');
}
inline void pr1(int x){write(x);putchar(' ');}
inline void pr2(LL x){write(x);putchar('\n');}
const int MAXN=100005;
const int mod=998244353;

int s[MAXN],a[MAXN],g[MAXN],n;
int lowbit(int x){return x&-x;}
void change(int x,int c){if(!x)return ;for(;x<=n+1;x+=lowbit(x))s[x]=(s[x]+c)%mod;}
int qry(int x){int ret=0;for(;x>=1;x-=lowbit(x))ret=(ret+s[x])%mod;return ret;}

int main()
{
	n=read();
	for(int i=1;i<=n;i++)a[i]=read();
	a[n+1]=n+1;a[n+2]=n+2;
	int L;g[1]=1;
	change(a[1],g[1]);
	if(a[1]<a[2])L=1;else L=2;
	bool tf=false;
	for(int i=2;i<=n;i++)
	{
		if(a[i]<a[i-1])tf=true;
		g[i]=qry(a[i+1]);if(!tf)g[i]=(g[i]+1)%mod;
		if(a[i]<a[i+1])
		{
			if(a[i]>a[i-1])change(a[i],g[i]);
			else 
			{
				for(int j=L-1;j<i;j++)change(a[j],-g[j]);
				change(a[i],g[i]);L=i;change(a[i-1],g[i-1]);
			}
		}
		else
		{
			for(int j=L-1;j<i;j++)change(a[j],-g[j]);
			change(a[i],g[i]);L=i+1;
		}
	}
	int ans=0;int u;
	for(int i=n;i>=1;i--)
	{
		ans=(ans+g[i])%mod;
		if(a[i]>a[i+1])break;
		u=i-1;
	}
	if(!u){ans=(ans+1)%mod;}
	pr2(ans);
	return 0;
}