2021牛客暑期多校训练营2 - G - League of Legends（ dp + 单调队列 ）

题目链接：点击进入

题目

题意

n 个区间，分成 k 组，每组取相交的公共长度，使所有组的相交长度和最大

思路

将区间分成两部分：一部分是包含其它区间的大区间，对于这类区间，他既可以单独算作一组，其长度计入答案；又可以并入它所包含的那组，不影响答案。
另一部分就是一些没有包含关系的区间。我们可以对这些区间进行排序 ( 左端点排序 ) ，然后进行 dp 。
规定 dp [ i ] [ j ] 代表前 i 个区间分成 j 组的交和的最大长度。假设区间形式 ( $l_i$ , $r_i$ )
dp [ i ] [ j ] = max ( dp [ k - 1 ] [ j - 1 ] + r [ k ] - l [ i ] ) 且 r [ k ] > l [ i ] ;
max 里的意思就是，前 k - 1 个区间构成 j - 1 组，后面第 k , k + 1 … i 这 i - k + 1 个区间构成一组，同时要求这些区间要有交，否则这个式子就错了，就成了加上负数了。( 因为是左端点排序，同时也没有包含关系，那么求这 i - k + 1 个区间的交集，我只需要看最左边的区间 ( k ) 跟最右边的区间 ( i ) 交集即可，其它区间的交最终都会归到这一交集 )。对于 max 这一部分，我们可以用单调队列来维护。队尾更新 dp [ k - 1 ] [ j - 1 ] + r [ k ] 的最大值，同时从队首把 r [ k ] <= l [ i ] 的 k 都弹掉，最后队首维护的就是最大值。
最后还需要枚举剩下的大区间，对于这些区间，排序后求前缀和，更新不足 k 组的 dp情况 ，取最大值为最终答案。

代码

//#pragma GCC optimize(3)//O3
//#pragma GCC optimize(2)//O2
#include<iostream>
#include<string>
#include<map>
#include<set>
//#include<unordered_map>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<stack>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<fstream>
#define X first
#define Y second
#define best 131
#define INF 0x3f3f3f3f3f3f3f3f
#define pii pair<int,int>
#define lowbit(x) x & -x
#define inf 0x3f3f3f3f
//#define int long long
//#define double long double
//#define rep(i,x,y) for(register int i = x; i <= y;++i)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const int maxn=5e3+10;
const int mod=998244353;
const double eps=1e-9;
const int N=5e3+10;
 
inline int read()
{
    int k = 0, f = 1 ;
    char c = getchar() ;
    while(!isdigit(c)){if(c == '-') f = -1 ;c = getchar() ;}
    while(isdigit(c)) k = (k << 1) + (k << 3) + c - 48 ,c = getchar() ;
    return k * f ;
}
 
struct node
{
    int l,r;
    bool operator<(const node s)const
    {
        if(l==s.l)
            return r>s.r;
        else
            return l<s.l;
    }
}p[maxn],a[maxn];
int n,k,seg[maxn],cnt,dp[maxn][maxn];
int q[maxn],tot; 
  
int main()
{
//  ios::sync_with_stdio(false);
//  cin.tie(0);cout.tie(0);
	cin>>n>>k;
	for(int i=1;i<=n;i++) cin>>p[i].l>>p[i].r;
	sort(p+1,p+n+1);
	int maxx=inf,ans=0;
	for(int i=n;i>=1;i--)
	{
		if(p[i].r<maxx)
		{
			maxx=p[i].r;
			a[++cnt]=p[i]; 
		}
		else
			seg[++tot]=p[i].r-p[i].l;
	}
	sort(a+1,a+cnt+1);
	memset(dp,-0x3f,sizeof(dp));
	dp[0][0]=0;
	for(int j=1;j<=k;j++)
	{
		int l=1,r=0;
		for(int i=1;i<=cnt;i++)
		{
			while(l<=r&&dp[q[r]-1][j-1]+a[q[r]].r<=dp[i-1][j-1]+a[i].r) r--;
			q[++r]=i;
			while(l<=r&&a[q[l]].r<=a[i].l) l++; 
			if(l<=r) dp[i][j]=dp[q[l]-1][j-1]+a[q[l]].r-a[i].l;
		}
	}
	sort(seg+1,seg+tot+1,greater<int>());
	for(int i=1;i<=tot;i++) seg[i]+=seg[i-1];
	for(int i=0;i<=tot&&k-i>=1;i++)
		ans=max(ans,dp[cnt][k-i]+seg[i]);
    printf("%d",ans);
    return 0;
}