CSP202209-5 高维亚空间超频物质变压缩技术

已于 2023-04-17 17:06:28 修改
阅读量934 收藏 2
点赞数 2
文章标签： 算法 c++ 图论
于 2023-04-17 17:04:54 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_53899788/article/details/130203563
版权
该文探讨了在CSP202209-5问题中，如何对具有体积和质量属性的黄金进行分组压缩，以最小化总代价。提出了从暴力解法到利用斜率优化和单调队列实现的时间复杂度优化策略，最终达到高效解决大规模问题的目标。
摘要由CSDN通过智能技术生成

CSP202209-5高维亚空间超频物质变压缩技术

题意：

给定 n n n 块黄金，每个黄金有体积 v i v_i vi​。将黄金分组进行压缩，每一组内的黄金编号连续，压缩一组黄金的代价为 ( s − L ) 2 (s-L)^2 (s−L)2， s s s 为改组黄金的体积和。

每块黄金有质量 m i m_i mi​，分组必须满足 组之间编号最大的的黄金质量递增。

询问分段压缩的最小代价。

解析：

子任务1，2

令 f i f_i fi​ 为前 i i i 块黄金都压缩完成的最小代价。

则 f i = min ⁡ 0 < j < i { f j + ( s i − s j − L ) 2 } ( m i > m j ) f_i = \min\limits_{0<j<i}\Big\{ f_j + (s_i-s_j-L)^2\Big\}(m_i>m_j) fi​=0<j<imin​{fj​+(si​−sj​−L)2}(mi​>mj​)， s s s 为 v v v 的前缀和

初始条件 f 0 = 0 f_0 = 0 f0​=0，答案为 f n f_n fn​。

此时时间时间复杂度为 O ( n 2 ) O(n^2) O(n2)

可以通过前两个测试点

代码如下：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef double db;
#define fi first
#define se second
#define debug(x) cerr << #x << ": " << (x) << endl
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
const int maxn = 1e6+10;
const int maxm = 1e5+10;
const double eps = 1e-12;
const ll INF = 0x3f3f3f3f3f3f3f3f;
typedef pair<int, int> pii;
#define int ll



ll n, L;
ll v[maxn], m[maxn], s[maxn];
ll f[maxn];


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

	
	cin >> n >> L;
	for(int i = 1; i <= n; i++){
		cin >> s[i];
		s[i] += s[i-1];
	}
	for(int i = 1; i <= n; i++)
		cin >> m[i];
		
	memset(f, INF, sizeof(f));
	f[0] = 0;
	for(int i = 1; i <= n; i++){
		for(int j = 0; j < i; j++){ //上一组的终点
			if(m[i] <= m[j])
				continue; 
			f[i] = min(f[i], f[j]+(s[i]-s[j]-L) * (s[i]-s[j]-L));
		}
	}
	cout << f[n] << endl;
	return 0;

}


子任务3

因为 m i = i m_i = i mi​=i，所以对于 i i i，可以从任意 j ( j < i ) j(j<i) j(j<i) 进行转移。

状态转移方程： f i = min ⁡ 0 < j < i { f j + ( s i − s j − L ) 2 } f_i = \min\limits_{0<j<i}\Big\{ f_j + (s_i-s_j-L)^2\Big\} fi​=0<j<imin​{fj​+(si​−sj​−L)2}因为有 s i × s j s_i\times s_j si​×sj​ 项，可以考虑斜率优化。

整理一下： f j + s j 2 + 2 L s j = 2 s i × s j + f i − s i 2 + 2 L s i f_j+s_j^2+2Ls_j = 2s_i\times s_j+f_i-s_i^2+2Ls_i fj​+sj2​+2Lsj​=2si​×sj​+fi​−si2​+2Lsi​确实是可以斜率优化的。 y = f j + s j 2 + 2 L s j , k = 2 s i , x = s j y = f_j+s_j^2+2Ls_j, k = 2s_i, x = s_j y=fj​+sj2​+2Lsj​,k=2si​,x=sj​

而且 k , x k,x k,x 都是单调的，所以可以用单调队列维护下凸壳。

代码如下：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef double db;
#define fi first
#define se second
#define debug(x) cerr << #x << ": " << (x) << endl
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
const int maxn = 1e6+10;
const int maxm = 1e5+10;
const double eps = 1e-12;
const ll INF = 0x3f3f3f3f3f3f3f3f;
typedef pair<int, int> pii;
#define int ll

#define x(a) (s[a])
#define y(a) (f[a]+s[a]*s[a]+2*L*s[a]) 
#define k(a) (2*s[a])

ll n, L;
ll v[maxn], m[maxn], s[maxn], p[maxn];
int q[maxn], hh = 1, tt = 0;
ll f[maxn];

long double slope(int a, int b){
	long double x1 = (long double)x(a);
	long double y1 = (long double)y(a);
	
	long double x2 = (long double)x(b);
	long double y2 = (long double)y(b);
	
	if(fabs(x2-x1) < eps)
		return y2 > y1 ? 1e18 : -1e18;
	else
		return (y2-y1)/(x2-x1);
}

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

	
	cin >> n >> L;
	for(int i = 1; i <= n; i++){
		cin >> s[i];
		s[i] += s[i-1];
	}
	for(int i = 1; i <= n; i++)
		cin >> m[i];
	
	if(n <= 2000){	
		memset(f, INF, sizeof(f));
		f[0] = 0;
		for(int i = 1; i <= n; i++){
			for(int j = 0; j < i; j++){
				if(m[i] <= m[j])
					continue; 
				f[i] = min(f[i], f[j]+(s[i]-s[j]-L) * (s[i]-s[j]-L));
			}
		}
		cout << f[n] << endl;
		return 0;
	}
	else{
		q[++tt] = 0;
		f[0] = 0;	
		for(int i = 1; i <= n; i++){
     	   while(hh < tt && slope(q[hh], q[hh+1]) <= k(i))
     	   	hh++;     
    	    int j = q[hh];
    	    f[i] = f[j] + (s[i]-s[j]-L) * (s[i]-s[j]-L);
     	   while(hh < tt && slope(q[tt-1], q[tt]) >= slope(q[tt], i))
     	   	tt--;				
     	   q[++tt] = i;
		}	
		cout << f[n] << endl;
		return 0;
	}
	return 0;	
}


子任务4

i i i 的决策点 j j j 必须满足 m j < m i m_j < m_i mj​<mi​，增加一维偏序，可以CDQ分治。

根据黄金质量分为左右两部分。当左侧的状态值计算完成，用左侧的状态更新右侧的状态，必然满足质量的限制。

用左侧值更新右侧值时，左侧的 x x x 单调，右侧处理成 k k k 单调的情况，可以继续用单调队列维护下凸壳。

时间复杂度为 O ( n log ⁡ 2 n ) O(n\log^2 n) O(nlog2n)

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef double db;
#define fi first
#define se second
#define debug(x) cerr << #x << ": " << (x) << endl
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
const int maxn = 1e6+10;
const int maxm = 1e5+10;
const double eps = 1e-12;
const ll INF = 0x3f3f3f3f3f3f3f3f;
typedef pair<int, int> pii;
#define int ll

#define x(a) (s[a])
#define y(a) (f[a]+s[a]*s[a]+2*L*s[a]) 
#define k(a) (2*s[a])

ll n, L;
ll m[maxn], s[maxn];
int q[maxn], hh = 1, tt = 0;
ll f[maxn];
struct node{
	int m, id;
	ll k, x, y;
}a[maxn], tmp[maxn];
bool cmpm(node a, node b){
	if(a.m == b.m)
		return a.id < b.id;
	else
		return a.m < b.m;
}
bool cmpx(node a, node b){
	return a.x < b.x;
}
long double slope(int i, int j){
	long double x1 = (long double)a[i].x;
	long double y1 = (long double)a[i].y;
	
	long double x2 = (long double)a[j].x;
	long double y2 = (long double)a[j].y;
	
	if(fabs(x2-x1) < eps)
		return y2 > y1 ? 1e18 : -1e18;
	else
		return (y2-y1)/(x2-x1);
}

void cdq(int ml, int mr, int pl, int pr){
	if(pl > pr)
		return;
	if(ml == mr){
		for(int i = pl; i <= pr; i++){
			int j = a[i].id;
			a[i].y = f[j] + s[j]*s[j] + 2 * L * s[j];
		}		
		return; 
	}
	
	int mmid = (ml+mr) >> 1;
	int pmid = upper_bound(a+pl, a+1+pr, (node){mmid, n+1, 0, 0}, cmpm)-a-1;
	cdq(ml, mmid, pl, pmid);	
	int hh = 1, tt = 0;
	for(int i = pl; i <= pmid; i++){
		while(hh < tt && slope(q[tt-1], q[tt]) >= slope(q[tt], i))
			tt--;
		q[++tt] = i;
	}
	
	if(pr > pmid){
		sort(a+pmid+1, a+pr+1, cmpx);
		for(int i = pmid+1; i <= pr; i++){
			while(hh < tt && slope(q[hh], q[hh+1]) <= a[i].k)
				hh++;
			if(hh <= tt){
				int id = a[i].id;
				int j = a[q[hh]].id;
				if(id < j)
					continue;
				f[id] = min(f[id], f[j]+(s[id]-s[j]-L)*(s[id]-s[j]-L));
			}
		}
		sort(a+pmid+1, a+pr+1, cmpm);
	}
	
	cdq(mmid+1, mr, pmid+1, pr);	
	sort(a+pl, a+pr+1, cmpx);
}
signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);

	
	cin >> n >> L;
	for(int i = 1; i <= n; i++){
		cin >> s[i];
		s[i] += s[i-1];
	}
	for(int i = 1; i <= n; i++)
		cin >> m[i];
		
	memset(f, INF, sizeof(f));
	f[1] = (s[1]-L) * (s[1]-L);
	for(int i = 1; i <= n; i++){
		a[i].id = i;
		a[i].m = m[i];
		a[i].k = 2 * s[i];
		a[i].x = s[i];
		f[i] = (s[i]-L) * (s[i]-L);
	}
	sort(a+1, a+1+n, cmpm); //按质量排序 
	
	cdq(1, n, 1, n);
	cout << f[n] << endl;
	return 0;

}
hydqiln
关注
csp 2022 - 09 压轴题 高维亚空间超频物质变压缩技术
Tommyjia0405的博客
01-22 1148
斜率优化 + 分治降维 解决9月份最后一题，想了好久
CCF-CSP 202209-2-何以包邮？ (Python,C++) 穷举 一维二维 动态规划 满分实现
m0_50078259的博客
03-06 1723
CCF-CSP 202209-2-何以包邮？ (Python ，C++实现) 穷举，动态规划，满分实现
CSP202209 03.18 JAVA 如此编码 何以包邮？ 防疫大数据 吉祥物投票 高维亚空间超频物质变压缩技术 持续更新欢迎交流..
weixin_45322373的博客
03-18 689
既然不能通过枚举所有的情况来找到答案 ，那肯定是有一种途径来找到最优解，所以自然的联想到了dp（动态规划）；可以转化为先求书的总价sum，再用sum减去x得y，问题就变成了寻找不超过y（达到包邮条件）的最大数字组合，再用总价减去这个数字组合，即可得到最终答案。寻找不超过y的最大数字组合，这种类型就是典型的背包问题。递推关系式：v[i]：各商品价值 w[i]：各商品重量 （这里书的重量就是它的价值）子问题：当前背包容量j，前i个物品最佳组合对应的价值。处理数据时按照风险名单判定来处理。
CSP202209-5---＞高维亚空间超频物质变压缩技术(50分代码)
weixin_51954774的博客
03-05 430
CSP刷题
CSP 202209题解：如此编码，何以包邮，防疫大数据，吉祥物投票，高维亚空间超频物质变压缩技术
fmy_xfk的博客
09-18 3514
CSP202209全部满分题解：如此编码，何以包邮，防疫大数据，吉祥物投票，高维亚空间超频物质变压缩技术
ccf-csp历年满分题解 + 考点汇总 + 经验总结 + 常用代码模板（已完成63题，持续更新中...）
一只可爱的小猴子的博客
03-19 7767
3. 当题目中没有过多的说明某个操作或数据可能出现的情况，但确实可能会出现递归处理等情况，可以先不考虑递归的情况，后续实在没有其他错误的时候，再来修改这个操作的处理情况（比如说，这道题目中，cancel这个记录有没有可能操作的是一条cancel记录，出现递归的情况，但是。//右对齐，字符宽度为m，如果要输出的字符宽度>=m，则全部输出，如果字符宽度
CCF-CSP历年真题大全附题解(202312已更)
热门推荐
weixin_53919192的博客
05-22 15万+
各位朋友，历年的题目你们要是有不同的解法想和大家进行分享的，可以私聊我发我题目编号和代码，我也可以更新到文章中，给需要的朋友多点参考~~ CCF-CSP真题拿来练手，如果对自己没有拿高分的期望的话，可以就琢磨着写一二三题就好，三题可以就拿最基础的情况，再慢慢添东西更新代码，持续更新，博主的算力有限，努力提升进步中，写一道发一道，目前先从简单慢慢做，有想法的可以一起讨论呀~~
CSP-防疫大数据202209-3
m0_50435987的博客
10-01 1632
CSP
2022 CSP-J1 CSP-S1 初赛 资料集（2022.07.08).pdf
07-08
知识点 5：CSP-J1 CSP-S1 初赛复习方法 本资源提供了 CSP-J1 CSP-S1 初赛的复习方法，旨在帮助选手更好地备战考试。 知识点 6：CSP-J1 CSP-S1 初赛知识汇总 本资源提供了 CSP-J1 CSP-S1 初赛的知识汇总，旨在帮助...
NOIP CSP-J CSP-S 初赛 第1轮 学习资料集（J）-2023.08.22.pdf
08-22
【CSP-J CSP-S 初赛】是中国计算机学会（CCF）组织的“全国青少年信息学奥林匹克联赛”（NOIP）的一部分，旨在检验初中生和高中生的信息技术能力，特别是编程和算法解决能力。赛事分为两个级别，CSP-J（入门级，...
2023 CSP-J2 CSP-S2 复赛 第2轮 真题讲解.pdf
10-24
【CSP-J2 CSP-S2】是中国计算机学会（CCF）组织的计算机编程能力认证，主要面向中学生，旨在提升他们的编程技能和信息学奥林匹克竞赛水平。CSP-J2是入门级，适合初级选手；CSP-S2是提高级，针对有一定基础的参赛者。...
dp算法练习【7】
MADAO的博客
09-08 1362
给你一个按排序的整数数组nums，返回组成的新数组，要求也按排序。平方后，数组变为 [16,1,0,9,100]排序后，数组变为 [0,1,9,16,100]非常简单，思路平方后排序。优化：插入法。
【数据结构】排序算法系列——希尔排序（附源码+图解）
Skrrapper的博客
09-07 1055
希尔排序实际上是个相当复杂的排序算法，这主要是跟它的步长序列gap到底该如何取、后续应该减小有关。时，数组已经接近有序的了，就整体而言，最后一次整体的插入排序就可以大大提高效率——我们从插入排序的时间复杂度分析也可以看出，可以看到每次减小gap的规律是将原先的gap/2，但事实上这只是其中一种处理方法，并不说明这是最优解。，使得平均情况下的时间复杂度能够在O(n^1.3)左右，而不好的选择则可能导致接近最坏情况的性能。中，我们会将整体数据一分为多份，进行散布式的插入排序，这时候每一个子序列之间的间隙就是。
【Python 千题 —— 算法篇】寻找两个正序数组的中位数
繁依Fanyi的博客
09-07 1250
合并排序法：易于理解，但时间复杂度较高，适合小规模数据。双指针法：通过双指针合并数组，时间复杂度 O(m + n)，适合中等规模数据。二分查找法：通过二分查找，在较短数组上进行划分，时间复杂度为 O(log(min(m, n)))，是处理大规模数据的最佳选择。对于这种涉及高效搜索和排序的题目，掌握二分查找的应用至关重要。通过本题目，我们不仅能够提升对数组和中位数的理解，还能够学到如何优化时间复杂度，使算法在处理大规模数据时更加高效。
洛谷刷题之B2089 数组逆序重存放
resco的博客
09-07 1180
洛谷刷题之B2089 数组逆序重存放
算法练习题22——解码
hello77的blogggg 祝你得偿所愿
09-11 296
小明有一串很长的英文字母，其中可能包含大写和小写字母。在这串字母中，有很多字母是连续且重复的。为了压缩这串字母的表达形式，小明决定将连续相同的字母用“字母 + 出现次数”的形式表示。
二分系列(二分查找)9/12
最新发布
2301_78191305的博客
09-12 333
排序数组->使用二分。第一个位置:>=；第二个位置:<=又因为<=target的话可以先看作>target-1 再看作>=(target+1)-1；所以直接写一个>=的函数。
算法 | 基础 | 出现奇数次的数字
蘑菇蘑菇不会开花的专栏
09-07 716
本篇是关于异或（）运算的运用。后期看算法过程中如果再碰到异或的都会收录到本篇中。
C语言程序设计-练习篇
2402_85207276的博客
09-09 535
·打印整数二进制的奇数位和偶数位 ·用循环打印"*"组成的X形图案 ·获得月份天数
CSP 202209-2
08-31
CSP 202209-2是一道题目，题目的名称是"何以包邮？"。该题的时间限制为1.0秒，内存限制为512.0MB。该题目的描述可以在CCF-CSP真题附题解大全中找到。同时，还提供了一种解法，该解法是直接无脑解的方法，可以得到70分。这个解法是通过输入购物车中的图书数量n和包邮最低值x，然后依次输入每本书的价格，对所有书本的价格进行排序，遍历每一个价格，判断是否小于最小的包邮值，如果小于，则将该价格存入集合temp中，如果不小于，则更新最小消费金额，并结束循环。最后打印出最小消费金额min_money。需要注意的是，该解法虽然能得到70分，但在写的过程中没有对每本书做标记，可能导致一些书重复购买了两次，不符合题意。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* *3* [CCF-CSP真题《202209-2 何以包邮?》思路+python满分题解](https://blog.csdn.net/weixin_53919192/article/details/127173638)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值