Monetone_Queue_2

已于 2024-04-02 19:25:29 修改
阅读量1k 收藏 18
点赞数 29
分类专栏： # 动态规划 文章标签： dp monotone_queue
于 2024-04-01 21:23:57 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/2301_78527293/article/details/137247218
版权

Monetone_Queue_2

Luogu P1725 琪露诺

f i f_i fi​ 表示在位置 i i i 取得得最大冰冻指数

f i = max ⁡ i − r i − l f ? + a i f_i=\max_{i-r}^{i-l}f_?+a_i fi​=maxi−ri−l​f?​+ai​

单调队列维护最大值即可，唯一需要注意不可达状态需要初始化为 − ∞ -\infty −∞

void solve(){
    cin >> n >> l >> r;
    for(int i = 0; i <= n; i ++){
        cin >> a[i];
    }
    for(int i = 0; i <= n + r; i ++){
        f[i] = -2e18;
    }
    f[0] = 0;
    for(int i = l; i <= n + r; i ++){
        while(h <= t && f[i - l] >= f[q[t]]) t -- ;
        q[++ t] = i - l;
        while(h <= t && q[h] < i - r) h ++;
        f[i] = f[q[h]] + a[i];
    }
    cout << *max_element(f + n, f + n + r + 1);
}

Luogu P3422 LOT-a jOURNEY TO Mars

破环成链，把原数组再复制一遍

1、顺时针情况

用 s i s_i si​ 表示从 1 到 i 每个结点的油量与距离之差的前缀和

则以 i 为起点的前缀和为 s j − s i − 1 s_j-s_{i-1} sj​−si−1​

从点 i 出发，对于任意 i ≤ j ≤ i + n − 1 i\leq j\leq i+n-1 i≤j≤i+n−1 , 都有 s j − s i − 1 ≥ 0 s_j-s_{i-1}\geq 0 sj​−si−1​≥0 , 则可绕行一周。

从大到小遍历，窗口向左滑动，用单调队列维护 s j s_j sj​ 的最小值

2、逆时针情况

用 s i s_i si​ 表示从 2n 到 i 的每个结点的油量与距离之差的后缀和

则以 i 为起点的后缀和为 s j − s i + 1 s_j-s_{i+1} sj​−si+1​

从点 i 出发，对于任意 i − n + 1 ≤ j ≤ i i-n+1\leq j\leq i i−n+1≤j≤i , 都有 s j − s i + 1 ≥ 0 s_j-s_{i+1}\geq 0 sj​−si+1​≥0 , 则可绕行一周。

从小到大遍历，窗口向右滑动，用单调队列维护 s j s_j sj​ 的最小值

void solve(){
	cin >> n;
	for(int i = 1; i <= n; i ++){
		cin >> p[i] >> d[i];
	}
	// 先处理正序行驶
	for(int i = 1; i <= n; i ++){
		s[i] = s[i + n] = p[i] - d[i];
	}
	for(int i = 1; i <= 2 * n; i ++){
		s[i] += s[i - 1];
	}
	// 窗口左滑
	h = 0, t = -1;
	for(int i = 2 * n; i >= 1; i --){
		while(h <= t && s[i] <= s[q[t]]) t --;
		q[++ t] = i;
		if(q[h] > i + n - 1) h ++;
		if(i <= n && s[q[h]] - s[i - 1] >= 0){
			fg[i] = true;
		}
	}

	// 再处理逆序行驶
	h = 0, t = -1;
	d[0] = d[n];
	for(int i = n; i >= 1; i --){
		s[i] = s[i + n] = p[i] - d[i - 1];
	}
	for(int i = 2 * n; i >= 1; i --){
		s[i] += s[i + 1];
	}
	// 窗口右滑
	for(int i = 1; i <= 2 * n; i ++){
		while(h <= t && s[i] <= s[q[t]]) t --;
		q[++ t] = i;
		while(h <= t && q[h] < i - n) h ++;
		if(i >= n + 1 && s[q[h]] - s[i + 1] >= 0){
			fg[i - n] = true;
		}
	}

	// 输出答案
	for(int i = 1; i <= n; i ++){
		cout << (fg[i] ? "TAK" : "NIE") << '\n';
	}
}

CF 372C Watching Fireworks Is Fun

单调队列优化DP, 设 f i , j f_{i,j} fi,j​ 为当前在燃放第 i i i 个烟花，位于位置 j j j 的最大幸福度

d i s = d × Δ t dis=d\times \Delta t dis=d×Δt

转移依据 ： 上一次位于什么位置

优化 : f i , j = max ⁡ j − d i s j + d i s f i − 1 , ? + b i − ∣ a i − j ∣ f_{i,j}=\max_{j-dis}^{j+dis}f_{i-1,?}+b_i-|a_i-j| fi,j​=maxj−disj+dis​fi−1,?​+bi​−∣ai​−j∣

用单调队列维护上一层窗口的最大值，注意本题需要优化空间

由于只用到上一维度，循环数组即可

#include<bits/stdc++.h>
using namespace std;
#define int long long
int n, m, d;
int f[2][150010];
int hd = 0, tl = -1, q[150010];
int a[310], b[310], t[310];
signed main(){
    // ios::sync_with_stdio(false); cin.tie(0), cout.tie(0);
    cin >> n >> m >> d;
    for(int i = 1; i <= m; i ++){
        cin >> a[i] >> b[i] >> t[i];
    }
    int now, last;
    for(int i = 1; i <= m; i ++){
        now = i & 1, last = (i - 1) & 1;
        int max_dis = d * (t[i] - t[i - 1]);
        hd = 0, tl = -1;
        for(int j = 1; j <= n; j ++){
            int ed = min(j + max_dis, n);
            if(j == 1){ // 初始化单调队列
                for(int k = 1; k <= ed; k ++){
                    while(hd <= tl && f[last][k] >= f[last][q[tl]]) tl --;
                    q[++ tl] = k;
                }
            }
            else{ // 一次加入一个
                while(hd <= tl && f[last][ed] >= f[last][q[tl]]) tl --;
                q[++ tl] = ed;
                while(hd <= tl && q[hd] < j - max_dis) hd ++;
            }
            f[now][j] = f[last][q[hd]] + b[i] - abs(a[i] - j);
        }
    }
    cout << *max_element(f[m & 1] + 1, f[m & 1] + n + 1) << '\n';
    return 0;
}

Luogu P2252 [NOI2005] 瑰丽华尔兹

先枚举 K K K 个区间，再枚举行或列，对于每行或者每列可以转化为序列的转移，当前位置 (x, y) 的滑行距离由定长区间内的最长距离转移，每个时间区间的长度即为窗口长度

状态 ： f [ x ] [ y ] f[x][y] f[x][y] 表示在位置 ( x , y ) (x,y) (x,y) 所获得的最长距离

转移 : f [ x ] [ y ] = m a x ( f [ x ′ ] [ y ′ ] + i − p o s ) f[x][y]=max(f[x^{'}][y^{'}]+i-pos) f[x][y]=max(f[x′][y′]+i−pos) , i − t i m ≤ j ≤ i i-tim\leq j\leq i i−tim≤j≤i

#include<bits/stdc++.h>
using namespace std;
#define int long long
int n, m, x, y, k, ans = -2e18;
int f[210][210];
char g[210][210];
int dx[5] = {0, -1, 1, 0, 0};
int dy[5] = {0, 0, 0, -1, 1};
int h, t;
struct node{int val, pos;} q[41000];
/* 出发点, 时间长度, 方向, 枚举距离*/
void dp(int x, int y, int tim, int d, int len){
	h = 0, t = -1;
	// 对于 i : [i - tim, i]
	for(int i = 1; i <= len; i ++){
		if(g[x][y] == 'x'){
			h = 0, t = -1;
		}
		else{
			while(h <= t && f[x][y] >= q[t].val + i - q[t].pos) t --;
			q[++ t] = node{f[x][y], i};
			if(q[h].pos < i - tim) h ++;
			f[x][y] = q[h].val + i - q[h].pos;
			ans = max(ans, f[x][y]);
		}
		x += dx[d], y += dy[d];
	}
}
// 
void solve(){
	cin >> n >> m >> x >> y >> k;
	for(int i = 1; i <= n; i ++){
		for(int j = 1; j <= m; j ++){
			cin >> g[i][j];
			f[i][j] = -2e18;
		}
	}
	f[x][y] = 0;
	for(int i = 1; i <= k; i ++){
		int s, t, d;
		cin >> s >> t >> d;
		int tim = t - s + 1;
		if(d == 1) for(int i = 1; i <= m; i ++) dp(n, i, tim, d, n);
		else if(d == 2) for(int i = 1; i <= m; i ++) dp(1, i, t - s + 1, d, n);
		else if(d == 3) for(int i = 1; i <= n; i ++) dp(i, m, t - s + 1, d, m);
		else for(int i = 1; i <= n; i ++) dp(i, 1, t - s + 1, d, m);
	}
	cout << ans;
}

signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0), cout.tie(0);

	int T = 1;
	while (T --){
		solve();
	}

	return 0;
}

Luogu P2569 股票交易

f[i][j] 表示第 i 天后，拥有 j 张股票的最大收益

  1. 第 i 天凭空买入，f[i][j]=-j*ap

  2. 第 i 天不买不卖, f[i][j]=f[i-1][j]

  3. 第 i 天在第 i-w-1 天的基础上买入 :

    f[i][j]=max{f[i-w-1][k]-ap*(j-k)}=max{f[i-w-1][k]+ap*k}-ap*j , j-as<=k<=j

  4. 第 i 天在第 i-w-1 的基础上卖出 :

    f[i][j]=max{f[i-w-1][k]+bp*(k-j)}=max{f[i-w-1][k]+bp*k}-bp*j , j<=k<=j+bs

#include<bits/stdc++.h>
using namespace std;
#define int long long
int n, m, w;
int f[2100][2100], q[2100];
void solve(){
    cin >> n >> m >> w;
    for(int i = 0; i <= n; i ++){
        for(int j = 0; j <= m; j ++){
            f[i][j] = -2e18;
        }
    }
    for(int i = 1; i <= n; i ++){
        int ap, bp, as, bs;
        cin >> ap >> bp >> as >> bs;
        // 1. 凭空买
        for(int j = 0; j <= as; j ++){
            f[i][j] = -j * ap;
        }
        // 2. 不买不卖
        for(int j = 0; j <= m; j ++){
            f[i][j] = max(f[i][j], f[i - 1][j]);
        }
        if(i <= w) continue; // 不能转移就跳过
        // 3. 买入 (窗口右滑)
        int h = 0, t = -1; 
        for(int j = 0; j <= m; j ++){
            while(h <= t && q[h] < j - as) h ++;
            while(h <= t && f[i - w - 1][j] + ap * j >= f[i - w - 1][q[t]] + ap * q[t]) t --;
            q[++ t] = j;
            f[i][j] = max(f[i][j], f[i - w - 1][q[h]] - ap * (j - q[h]));
        }
        // 4. 卖出 (窗口左滑)
        h = 0, t = -1;
        for(int j = m; j >= 0; j --){
            while(h <= t && q[h] > j + bs) h ++;
            while(h <= t && f[i - w - 1][j] + bp * j >= f[i - w - 1][q[t]] + bp * q[t]) t --;
            q[++ t] = j;
            f[i][j] = max(f[i][j], f[i - w - 1][q[h]] + bp * (q[h] - j));
        }
    }
    // 不留任何股票就是最优答案
    cout << f[n][0] << '\n';
}

signed main(){
    ios::sync_with_stdio(false);
    cin.tie(0), cout.tie(0);

    int T = 1;
    while (T --){
        solve();
    }

    return 0;
}
running_Turtle1
关注
专栏目录
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值