YbtOJ 冲刺 NOIP2020 模拟赛 Day10 D.弱者对决【区间dp】

滑稽大佬

于 2020-12-04 09:38:17 发布

阅读量399

点赞数

分类专栏：区间dp

本文链接：https://blog.csdn.net/qq_35786326/article/details/110622307

版权

区间dp 专栏收录该内容

5 篇文章

订阅专栏

本文介绍了一种通过区间贡献优化来解决特定类型问题的算法。该算法通过对询问区间进行离散化处理，并使用动态规划的方式，记录每个询问区间的贡献值，以达到最大化最终答案的目的。文章详细介绍了状态定义、状态转移方程以及优化后的复杂度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目：

传送门

题意：

有 $m$ 个形如 ${a_i,b_i,c_i\}$ 的三元组，表示在 $a_i\sim b_i$ 的区间内， $x_{min}$ 如果 $c_i$ ，则此询问有 $x_{min}$ 的贡献
问如何安排 $x_i$ 的大小能使得最后答案最大

分析：

首先我们可以发现， $x_i$ 一定是某个 $c_i$ ，因为如果 $c_{i1}<x_i<c_{i2}$ ，那么对于 $i_1$ 并不会产生更多的贡献，而 $i_2$ 虽然产生了，但显然比 $x_i=c_{i2}$ 的贡献小
据此，我们阔以考虑将 $c$ 离散化处理
我们需要记录的量有区间 $(l, r)$ 、 $x_{min}$ 、 $x_{min}$ 的位置
那么设 $f_{l,r,w}$ 表示对于所有 $l<=a_i\&b_i<=r$ 的询问，当 $x_{min}=w$ 时取得最大贡献， $k$ 为 $l\sim r$ 间的某个位置，表示 $w$ 的位置
同时统计下 $w$ 会对哪些询问产生贡献，显然需要满足 $l<=a_i<=k<=b_i<=r\&c_i>=w$ ，设这样的询问有 $cnt_w$ 个
有 $f_{l,r,w}=w*cnt_w+max\{f_{l,k-1,h}\}+max\{f_{k+1,r,h}\}\ |\ h>=w$
直接转移复杂度为 $O(n^3m)$
但注意到是取最大值，我们可以对每个 $f_{l,r}$ 搞一个后缀最大值，转移的时候直接调用
此时复杂度来到 $O(n^2m)$ ，可以解决此题
对于每个 $x_i$ 的值，我们只需要记录下每个最大值是从 $w$ 为哪个值转移来的

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>  
#include<string>
#include<algorithm>
#include<queue>
#define LL long long
using namespace std;
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
int a[4005],b[4005],c[4005],cnt[4005],bc[4005];
int n,m,t;
void S(int l,int r,int k)
{
	memset(cnt,0,sizeof(cnt));
	for(int i=1;i<=m;i++)
	  if(l<=a[i]&&a[i]<=k&&k<=b[i]&&b[i]<=r)
	    cnt[c[i]]++;
	for(int i=t;i;i--) cnt[i]+=cnt[i+1];
	return;
}
int f[55][55][4005],kf[55][55][4005],g[55][55][4005],kg[55][55][4005];
void get(int l,int r,int w)
{
	if(!l||!r||l>r) return;
	int k=kf[l][r][w];
	get(l,k-1,kg[l][k-1][w]);
	printf("%d ",bc[w]);
	get(k+1,r,kg[k+1][r][w]); 
	return;
}
int main()
{
	n=read();m=read();
	for(int i=1;i<=m;i++) a[i]=read(),b[i]=read(),bc[i]=c[i]=read();
	sort(bc+1,bc+1+m);
	t=unique(bc+1,bc+1+m)-bc-1;
	for(int i=1;i<=m;i++) c[i]=lower_bound(bc+1,bc+1+t,c[i])-bc;
	for(int len=1;len<=n;len++)
	  for(int l=1;l<=n;l++)
	  {
	  	int r=l+len-1;
	  	if(r>n) continue;
	  	for(int k=l;k<=r;k++)
	  	{
	  		S(l,r,k);
	  		for(int w=1;w<=t;w++)
	  		  if(f[l][r][w]<g[l][k-1][w]+g[k+1][r][w]+bc[w]*cnt[w]||!f[l][r][w])
	  		    f[l][r][w]=g[l][k-1][w]+g[k+1][r][w]+bc[w]*cnt[w],kf[l][r][w]=k;
		}
		for(int w=t;w;w--)
		  if(g[l][r][w+1]>f[l][r][w]) g[l][r][w]=g[l][r][w+1],kg[l][r][w]=kg[l][r][w+1];
	  	  else g[l][r][w]=f[l][r][w],kg[l][r][w]=w;
	  }
	printf("%d\n",g[1][n][1]);
	get(1,n,kg[1][n][1]);
	return 0;
}