最長(zhǎng)上升子序列

最長(zhǎng)上升子序列(Longest Increasing Subsequence)

int len[__];

int LIS(int a[],int n)
{
    int lis=0;
    for(int i=1;i<=n;++i)
    {
        int x=lower_bound(len+1,len+lis+1,a[i])-len;
        len[x]=a[i],lis=max(x,lis);
    }
    return lis;
}

最長(zhǎng)上升子序列方案數(shù)

struct node
{
    int x,p,ans;

    bool operator<(const node &b)const
    {
        if(x!=b.x)return x<b.x;
        return p<b.p;
    }
}dp[__];

int n,a[__],b[__],len[__];

/*
樹(shù)狀數(shù)組:
單點(diǎn)修改void add(int x,int val)
查詢(xún)前綴和int sum(int x)
*/

int main()
{
    sf("%d",&n);
    for(int i=1;i<=n;++i)
        sf("%d",&a[i]),b[i]=a[i];

  //離散化數(shù)組a[]
    sort(b+1,b+1+n);
    int m=unique(b+1,b+1+n)-b-1;
    fup(i,1,n)a[i]=lower_bound(b+1,b+1+m,a[i])-b;

  //求最長(zhǎng)上升子序列及其每個(gè)數(shù)的深度
    int lis=0;
    for(int i=1;i<=n;++i)
    {
        int x=lower_bound(len+1,len+lis+1,a[i])-len;
        len[x]=a[i];
        dp[i]={x,i,1};
        lis=max(x,lis);
    }

    sort(dp+1,dp+n+1);
  //x->第i層, y->第i+1層
    int x=1,y=0,ans=0;

  //初始化第一層, 使y指向第二層
    while(++y<=n && dp[y].x==1)
        if(lis==1)++ans;

    for(int i=1;i<lis;++i)
    {
        int z=y;
      //雙指針
        for(;y<=n && dp[y].x==i+1;++y)
        {
            for(;x<=n && dp[x].x==i && dp[x].p<dp[y].p;++x)
                T.add(a[dp[x].p],dp[x].ans);
            dp[y].ans=1ll*dp[y].ans*T.sum(a[dp[y].p]-1)%mod;
            if(i+1==lis && (ans+=dp[y].ans)>=mod)
                ans-=mod;
        }
      //清空樹(shù)狀數(shù)組
        for(--x;dp[x].x==i;--x)
            T.clear(a[dp[x].p]);
        x=z;
    }
    printf("%d\n",ans);
    return 0;
}
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀(guān)點(diǎn),簡(jiǎn)書(shū)系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

友情鏈接更多精彩內(nèi)容