Codeforces Round #661 (Div. 3) F. Yet Another Segments Subset

dp[i][j]表示左边界是i右边界是j时最多包含的线段个数。

传统的区间dp是n^3的，考虑如何优化。

如果存在一条线段正好左端点是i右端点是j，我们可以直接对这个线段计数，因为这个线段不会影响其他线段。

第一个转移是dp[i+1][j]，左端点直接向右移动一格。

第二个转移是切割[i,j]这个区间，但是你不需要O(n)地切割，使用以i为左端点的线段的右端点切割即可。

考虑复杂度的方式：每条线段至多被枚举n次。因此复杂度是O(n^2)。

#include<bits/stdc++.h>
#define pb push_back
#define fi first
#define se second
#define sz(x)  (int)x.size()
#define cl(x)  x.clear()
#define all(x)  x.begin() , x.end()
#define rep(i , x , n)  for(int i = x ; i <= n ; i ++)
#define per(i , n , x)  for(int i = n ; i >= x ; i --)
#define mem0(x)  memset(x , 0 , sizeof(x))
#define mem_1(x)  memset(x , -1 , sizeof(x))
#define mem_inf(x)  memset(x , 0x3f , sizeof(x))
#define debug(x)  cerr << #x << " = " << x << '\n'
#define ddebug(x , y)  cerr << #x << " = " << x << "   " << #y << " = " << y << '\n'
#define ios std::ios::sync_with_stdio(false) , cin.tie(0)
using namespace std ;
typedef long long ll ;
typedef long double ld ;
typedef pair<int , int> pii ;
typedef pair<ll , ll> pll ;
typedef double db ;
const int mod = 998244353 ;
const int maxn = 6000 + 10 ;
const int inf = 0x3f3f3f3f ;
const double eps = 1e-6 ; 
int n ;
int l[maxn] , r[maxn] ;
int q[maxn] ;
int cnt ;
vector<int> pos[maxn] ;
int dp[maxn][maxn] ;
//int num[maxn][maxn] ;
int num(int l , int r)
{
    return upper_bound(all(pos[l]) , r) - lower_bound(all(pos[l]) , r) ;
}
int dfs(int l , int r)
{
    if(dp[l][r] != -1)  return dp[l][r] ;
    int c = num(l , r) ;
    if(l > r)  return 0 ;
    if(l == r)  
    {
        dp[l][r] = c ;
        return c ;
    } 
    int res = dfs(l + 1 , r) ;
    for(auto u : pos[l])
    {
        if(u >= r)  break ;
        res = max(res , dfs(l , u) + dfs(u + 1 , r)) ;
    }
    dp[l][r] = c + res ;
    return dp[l][r] ;
}
int main()
{
    ios ;
    int T ;
    cin >> T ;
    while(T --)
    {
        cin >> n ;
        int now = 0 ;
        rep(i , 1 , n)  cin >> l[i] >> r[i] , q[++ now] = l[i] , q[++ now] = r[i] ;
        sort(q + 1 , q + now + 1) ;
        cnt = unique(q + 1 , q + now + 1) - q - 1 ;
        int mx = 0 ;
        rep(i , 1 , n)  l[i] = lower_bound(q + 1 , q + cnt + 1 , l[i]) - q , r[i] = lower_bound(q + 1 , q + cnt + 1 , r[i]) - q ;
        rep(i , 1 , n)  mx = max({mx , l[i] , r[i]}) ;
        rep(i , 1 , mx)  cl(pos[i]) ;
        rep(i , 1 , n)  pos[l[i]].pb(r[i]) ;
        rep(i , 1 , mx)  sort(all(pos[i])) ;
        //rep(i , 1 , n)  num[l[i]][r[i]] ++ ;
        rep(i , 1 , mx)  rep(j , 1 , mx)  dp[i][j] = -1 ;
        dfs(1 , mx) ;
        //rep(i , 1 , n)  num[l[i]][r[i]] -- ;
        cout << dp[1][mx] << '\n' ;
    }
    return 0 ;
}