【jzoj4920】【降雷皇】【动态规划】【二分答案】【数据结构】

题目大意

求一个序列的最长上升子序列及方案数。

解题思路

首先维护一个数组，t[i]表示长度为i的最长上升子序列末尾最小是多少，这样就可以求出第一个答案。每个长度再用权值线段树维护每个末尾的方案数，转移方案数时求一下比它小的末尾方案数之和。

code

#include<cmath>#include<cstdio>#include<algorithm>#define LL long long#define min(a,b) ((a<b)?a:b)#define max(a,b) ((a>b)?a:b)#define fo(i,j,k) for(int i=j;i<=k;i++)#define fd(i,j,k) for(int i=j;i>=k;i--)using namespace std;int const maxn=1e5,inf=1e9;int n,type,t[maxn+10],lson[maxn*16+10],rson[maxn*16+10],cntpon=maxn;LL f[maxn*16+10],mo=123456789;LL qury(int now,int l,int r,int pos){    int mid=(l+r)>>1;    if(l==r)return f[now];    if(pos<=mid)return (lson[now])?qury(lson[now],l,mid,pos):0;    else return (f[lson[now]]+((rson[now])?qury(rson[now],mid+1,r,pos):0))%mo;}void add(int now,int l,int r,int pos,LL val){    int mid=(l+r)>>1;    f[now]=(f[now]+val)%mo;    if(l==r)return;    if(pos<=mid){        if(!lson[now])lson[now]=++cntpon;        add(lson[now],l,mid,pos,val);    }else{        if(!rson[now])rson[now]=++cntpon;        add(rson[now],mid+1,r,pos,val);    }}int main(){    //freopen("hamon.in","r",stdin);    //freopen("hamon.out","w",stdout);    freopen("d.in","r",stdin);    freopen("d.out","w",stdout);    scanf("%d%d",&n,&type);    int a,l,r,mid,mx=0;    fo(i,mx+1,n)t[i]=inf;    fo(i,1,n){        scanf("%d",&a);        l=0,r=mx;        while(l!=r){            mid=(l+r+1)>>1;            if(t[mid]>=a)r=mid-1;            else l=mid;        }        if(a<t[l+1])t[l+1]=a;        if(l+1>mx)mx=l+1;        add(l+1,1,maxn,a,max(1,qury(l,1,maxn,a-1)));    }    printf("%d\n",mx);    if(type)printf("%lld\n",qury(mx,1,maxn,maxn));    return 0;}