NC 54585 DP + BIT

159 篇文章 1 订阅
题意

传送门 NC 54585

题解

d p [ i ] [ j ] dp[i][j] dp[i][j] 代表以 a [ i ] a[i] a[i] 结尾长度为 j j j 的子序列,则有递推式
d p [ i ] [ j ] = ∑ i ′ < i 且 a [ i ′ ] < a [ i ] d p [ i ′ ] [ j − 1 ] dp[i][j] = \sum\limits_{i'<i 且 a[i']<a[i]}dp[i'][j-1] dp[i][j]=i<ia[i]<a[i]dp[i][j1] B I T [ k ] BIT[k] BIT[k] a [ i ] a[i] a[i] 为索引,维护长度为 k k k 的子序列的数量和。

#include <bits/stdc++.h>
using namespace std;
#define maxn 500005
#define maxk 12
typedef long long ll;
const int mod = 998244353;
int n, k, a[maxn], dp[maxk], bit[maxn][maxk];

void compress(int *x, int n)
{
    vector<int> xs(n);
    for (int i = 0; i < n; i++)
    {
        xs[i] = x[i];
    }
    sort(xs.begin(), xs.end());
    xs.erase(unique(xs.begin(), xs.end()), xs.end());
    for (int i = 0; i < n; i++)
    {
        x[i] = 1 + lower_bound(xs.begin(), xs.end(), x[i]) - xs.begin();
    }
}

void add(int i, int j, int x)
{
    while (i <= n)
    {
        bit[i][j] = (bit[i][j] + x) % mod;
        i += i & -i;
    }
}

int sum(int i, int j)
{
    int s = 0;
    while (i > 0)
    {
        s = (s + bit[i][j]) % mod;
        i -= i & -i;
    }
    return s;
}

int main()
{
    scanf("%d%d", &n, &k);
    for (int i = 0; i < n; i++)
    {
        scanf("%d", a + i);
    }
    compress(a, n);
    int res = 0;
    for (int i = 0; i < n; i++)
    {
        dp[1] = 1;
        for (int j = 2; j <= k; j++)
        {
            dp[j] = sum(a[i] - 1, j - 1);
        }
        for (int j = 1; j < k; j++)
        {
            add(a[i], j, dp[j]);
        }
        res = (res + dp[k]) % mod;
    }
    printf("%d\n", res);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值