[TJOI2014]上升子序列(树状数组)

2986. 上升子序列 - AcWing题库

P3970 [TJOI2014]上升子序列 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

#include <bits/stdc++.h>
#include <iostream>
#include<unordered_map>
#define x first
#define y second

#define  ios ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);

using namespace std;

typedef  long  long LL ;
typedef  unsigned long  long ULL ;
typedef pair<int,int> PII ;
typedef pair<LL,int> PLI ;
const int N = 1e6+10,mod = 1e9 + 7 ;
const LL INF2 = 1e17;


map<LL,LL> ma;
int n,m;
LL a[N],b[N];
bool f[N];// 记录一个数 是否是第一次出现
LL lastans[N],tr[N];
// lastans 这个数上一次出现时对 答案的 影响     

int lowbit(int x)
{
    return x & -x;
}


void add(int x,LL v)
{
    for(int i=x;i <=n ;i+= lowbit(i))
    tr[i] = (tr[i] + v ) % mod;
}

LL query(LL x)
{
    LL ans =0;
    for(int i=x;i;i-= lowbit(i))
    ans = (ans + tr[i]) % mod ;
    return ans ;
}


inline void solve()
{
    cin >> n;
    for(int i= 1; i<= n; i++ ) cin >> a[i],b[i] = a[i];
    sort(b + 1, b + 1 + n);
    int tot = unique(b + 1, b + 1 +n ) -(b + 1 );

    for(int i=1;i<= tot; i ++ ) ma[b[i]] = i;
    // ma[i]  i这个数在原序列的 从小到大的排序的排名
    LL cnt =0;
    for(int i=1;i <= n;i ++ )
    {
        int j= ma[a[i]];
        if(!f[j])  // 当前这个数是第一次出现
        {
            LL x = query(j-1);// 比当前这个数 小的数的个数x
            cnt = (cnt +x ) % mod;
            add(j,x+ 1);// 在当前这个数的 位置 插入 x + 1

            lastans[j] =x;
            f[j] = 1;
            continue;
        }
        LL x = query(j-1);
        cnt = (cnt + x - lastans[j] + mod ) % mod;
        //之前 cnt 已经加过了 j 这个排名 第一次 出现的影响 ,当前j 的影响 包含之前的影响要减去
        add(j,(((x - lastans[j]) % mod ) + mod ) % mod );
        lastans[j] = x;
    }

    cout << cnt << endl;
}
signed main()
{
    ios
    int T=1;
//    cin>>T;
    while(T -- )
    {
        solve();

    }


    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
树状数组可以用来解决最长公共子序列问题。下面是使用树状数组优化的最长公共子序列求解算法。 首先,我们需要将两个序列分别离散化,将每个数映射到一个连续的整数区间内,然后将它们分别存储在两个数组中。 接着,我们定义一个二维数组`dp`,其中`dp[i][j]`表示序列1中前i个数和序列2中前j个数的最长公共子序列长度。则有以下状态转移方程: ```c if (a[i] == b[j]) dp[i][j] = dp[i-1][j-1] + 1; else dp[i][j] = max(dp[i-1][j], dp[i][j-1]); ``` 其中,`a`和`b`分别是两个离散化后的序列。 时间复杂度为O(n^2)。 然后,我们可以使用树状数组来优化这个算法,将时间复杂度降至O(nlogn)。 我们需要定义一个树状数组`c`,其中`c[i]`表示序列1中前i个数中最后一个数在序列2中出现的位置。然后,我们可以用二分查找来找到序列1中第i个数在序列2中出现的最晚位置,即`c[i]`。 接着,我们可以通过遍历序列1中的每个数,用树状数组更新`c`数组,并根据`c`数组和状态转移方程来更新`dp`数组。 具体来说,对于序列1中的第i个数,我们可以用二分查找在序列2中找到它出现的最晚位置`pos`,然后用树状数组将`pos`更新为i。接着,我们可以遍历序列2中的每个数,如果它在序列1中出现过,则可以根据状态转移方程来更新`dp`数组。 时间复杂度为O(nlogn)。 以下是使用树状数组优化的最长公共子序列求解算法的完整代码实现: ```c #include <stdio.h> #include <stdlib.h> #include <string.h> #define MAX_N 100000 int a[MAX_N + 10], b[MAX_N + 10]; int c[MAX_N + 10]; int dp[MAX_N + 10][2]; int n, m; int lowbit(int x) { return x & (-x); } void update(int x, int val) { while (x <= n) { c[x] = max(c[x], val); x += lowbit(x); } } int query(int x) { int res = 0; while (x) { res = max(res, c[x]); x -= lowbit(x); } return res; } int main() { scanf("%d %d", &n, &m); for (int i = 1; i <= n; i++) { scanf("%d", &a[i]); } for (int i = 1; i <= m; i++) { scanf("%d", &b[i]); } // 离散化 int k = 1; for (int i = 1; i <= n; i++) { for (int j = k; j <= m; j++) { if (a[i] == b[j]) { a[i] = j; k = j + 1; break; } } } // 初始化 memset(c, 0, sizeof(c)); memset(dp, 0, sizeof(dp)); // 动态规划求解 for (int i = 1; i <= n; i++) { int pos = query(a[i]); dp[i][0] = dp[i-1][1]; dp[i][1] = dp[i-1][1]; if (pos > 0) { dp[i][1] = max(dp[i][1], dp[pos-1][0] + i - pos + 1); } update(a[i], i); } printf("%d\n", dp[n][1]); return 0; } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值