CodeForces - 960F[动态开点线段树优化dp]详解

题意：给一张有向图，每条边有边权与编号，求一条最长的路径，这条路径的边权与编号都是递增的。（编号指输入顺序）

首先我们回忆一下普通得LIS得做法：就是dp[i]以第i个结尾得最长上升子序列的长度，那么同理我们就可以推出这道题：(假设边是从u -> v)dp[v][w]就是以结点v结尾并且权值不大于w的LIS的最大长度
那么状态转移方程就是dp[v][w] = max(dp[u][1,2,3,4…w-1])很明显如果暴力去循环的化时间复杂度就会爆炸那么我们想一下如何去优化

---

优化：首先我们应如何在对于dp的第二维我们应该想如何快速的找到(dp[u][1,2,3,4…w-1])之中的最大值，由于是选边的编号必须是递增的所以我们dp[][1,2,3,4…w-1]这个数组应支持修改，假设dp[v][w]找到了答案那么我们还要去更新dp[v][>=w]的所有数值因为每一次只能接在前面的点后面

---

1.对于第二维找动态区间最大值很明显就是线段树
2.对于第一维我们就可以开n个线段树，由于空间开销大我们直接动态开点去做就好了
下面代码：

#include <iostream>
#include <cstdio>
#include <stack>
#include <sstream>
#include <vector>
#include <map>
#include <cstring>
#include <deque>
#include <cmath>
#include <deque>
#include <iomanip>
#include <queue>
#include <algorithm>
#include <set>
#define mid ((l + r) >> 1) 
#define Lson rt << 1, l , mid
#define Rson rt << 1|1, mid + 1, r
#define ms(a,al) memset(a,al,sizeof(a))
#define _for(i,a,b) for( int i = (a); i < (b); ++i)
#define _rep(i,a,b) for( int i = (a); i <= (b); ++i)
#define for_(i,a,b) for( int i = (a); i >= (b); -- i)
#define rep_(i,a,b) for( int i = (a); i > (b); -- i)
#define lowbit(x) ((-x) & x)
#define IOS std::ios::sync_with_stdio(0); cin.tie(0); cout.tie(0)
#define INF 0x3f3f3f3f
#define hash Hash
#define all(x) x.begin(),x.end()
#define next Next
#define f first
#define s second
using namespace std;
const int N = 2e5 + 10;
const double eps = 1e-9;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
inline LL read() {
   LL s=0,w=1;
   char ch=getchar();
   while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
   while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
   return s*w;
}
struct node 
{
    int val;
    int l, r;
}tr[N * 40];
int idx, root[N];
inline void update(int &now, int l, int r, int pos, int val)
{
    if(!now) now = ++ idx;
    if(l == r)
    {
        tr[now].val = max(tr[now].val,val);
        return;
    }
    if(pos <= mid)
        update(tr[now].l,l,mid,pos,val);
    else 
        update(tr[now].r,mid+1,r,pos,val);
    tr[now].val = max(tr[tr[now].l].val, tr[tr[now].r].val);
}

inline int query(int now, int l, int r, int posl, int posr)
{
    if(!now) return 0;
    if(posl <= l && posr >= r) return tr[now].val;
    int res = 0;
    if(posl <= mid)
       res = max(res,query(tr[now].l,l,mid,posl,posr));
    if(posr > mid)
      res = max(res,query(tr[now].r,mid+1,r,posl,posr));
     return res;
}

int main()
{
    IOS;
    int n, m;
    int res = 0;
    cin >> n >> m;
    while(m --)
    {
        int l, r, w;
        cin >> l >> r >> w;
        // if(!w) continue;
        int t = query(root[l],1,1e5+10,1,w) + 1;
        // cout << t << endl;
        res = max(t,res);
        update(root[r],1,1e5+10,w+1,t);
    }
    cout << res ;
    return 0;
}