Codeforces Round #219 (Div. 1) C. Watching Fireworks is Fun(单调队列优化dp)

本文介绍了如何运用动态规划和单调队列优化来解决一个关于烟花燃放时的最大收益问题。给定一定数量的点和烟花,每个烟花在特定时间和地点燃放,目标是在烟花燃放时最大化收益,而每秒最多能移动一定距离。通过建立状态转移方程并利用单调队列维护最大值,实现了高效求解。
摘要由CSDN通过智能技术生成

题目链接
题意:

给出一条线段,有 n n n 个点,编号为 1 1 1 n n n ,有 m m m 个烟花,第 i i i 个烟花会在 t i t_i ti 时刻在 a i a_i ai处燃放,假如此时你的位置为 x x x,那么收益就是 b i − ∣ a i − x ∣ b_i-|a_i-x| biaix 。每秒最多可以移动 d d d单位,求为能获得的最大收益。

题解:

先考虑 d p dp dp ,设 d p [ i ] [ j ] dp[i][j] dp[i][j]表示第 i i i 个烟花燃放时处在 j j j 位置的最大收益,那么可以写出转移方程:

d p [ i ] [ j ] = max ⁡ ( d p [ i − 1 ] [ k ] ) + b i − ∣ a i − j ∣ dp[i][j]= \max (dp[i-1][k]) +b_i-|a_i-j| dp[i][j]=max(dp[i1][k])+biaij

j − ( t [ i ] − t [ i − 1 ] ) ∗ d ≤ k ≤ j + ( t [ i ] − t [ i − 1 ] ) ∗ d j-(t[i]-t[i-1])*d \leq k \leq j+(t[i]-t[i-1])*d j(t[i]t[i1])dkj+(t[i]t[i1])d

很明显,这肯定不能暴力求解 d p [ i − 1 ] [ k ] dp[i-1][k] dp[i1][k] ,既然是要找最大值,那么可以利用单调队列维护最大值即可。

代码:

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<stack>
#include<set>
#include<ctime>
#define iss ios::sync_with_stdio(false)
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
typedef pair<int,int>pii;
const int MAXN=2e5+5;
const int mod=1e9+7;
const int inf=0x3f3f3f3f;
int n,m,d;
struct node
{
    ll a,b,t;
}f[MAXN];
ll dp[2][MAXN];
ll q[MAXN];
bool cmp(node x,node y)
{
    return x.t<y.t;
}
int main()
{
    memset(dp,-inf,sizeof dp);
    cin>>n>>m>>d;
    for(int i=1;i<=m;i++)
    {
        cin>>f[i].a>>f[i].b>>f[i].t;
    }
    sort(f+1,f+1+m,cmp);
    int fl=1;
    for(int i=1;i<=n;i++) dp[0][i]=0;
    for(int i=1;i<=m;i++)
    {
        int k=1;
        int l=1,r=0;
        for(int j=1;j<=n;j++)
        {
            ll maxx=j+(f[i].t-f[i-1].t)*d;
            if(maxx>n) maxx=n;
            ll minn=j-(f[i].t-f[i-1].t)*d;
            if(minn<1) minn=1;
            while(k<=maxx)
            {
                while(l<=r&&dp[fl^1][k]>dp[fl^1][q[r]]) r--;
                q[++r]=k;
                k++;
            }
            while(l<=r&&q[l]<minn) l++;
            dp[fl][j]=dp[fl^1][q[l]]+f[i].b-abs(f[i].a-j);
        }
        fl=(fl^1);
    }
    ll ans=-1e18;
    for(int i=1;i<=n;i++)
    {
        ans=max(ans,dp[fl^1][i]);
    }
    printf("%lld\n",ans);
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值