1571 例题3 凸多边形的划分(LOJ10149) 区间动规思想40分 区间动归 dp 两定点一动点50分 int128使用100分 uint128使用100分

最新推荐文章于 2024-07-17 17:34:40 发布
最新推荐文章于 2024-07-17 17:34:40 发布
阅读量744 收藏
点赞数 1
分类专栏： 信息学奥赛一本通 提高篇 文章标签： 动态规划
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mrcrack/article/details/122331030
版权

总目录

在线测评地址(ybt)

在线测评地址(LOJ)

1.区间动规思想40分

ybt

未通过

测试点结果内存时间
测试点1答案错误620KB3MS
测试点2答案正确616KB4MS
测试点3答案正确608KB2MS
测试点4答案正确620KB2MS
测试点5答案正确616KB2MS
测试点6答案错误616KB2MS
测试点7答案错误612KB2MS
测试点8答案错误616KB2MS
测试点9答案错误616KB1MS
测试点10答案错误616KB2MS

LOJ

 问题比较陌生，样例还是要好好模拟

 该题的数据范围N≤50，每个点权值小于 10^9

数据极限(50-2)*10^9*10^9*10^9=4.8*10^28，unsigned long long(2^64=18446744073709551616)也必溢出，需使用高精度算法。

在上述模拟过程，感觉纯暴力即可AC.

编写，提交40分，代码如下：

#include <bits/stdc++.h>
#define ULL unsigned long long
#define maxn 110
using namespace std;
ULL a[maxn],b[maxn],mx;
int main(){
	int i,j,n;
	scanf("%d",&n);
	for(i=1;i<=n;i++)scanf("%llu",&a[i]),a[i+n]=a[i];
	for(i=1;i<=n;i++){//以i作为三角形顶点
		for(j=i+1;j<=i+n-2;j++)
			b[i]+=a[j]*a[j+1];
		b[i]*=a[i];
	}
	mx=(ULL)18e18;
	for(i=1;i<=n;i++)mx=min(mx,b[i]);
	printf("%llu\n",mx);
	return 0;
}

2.区间动归 dp 两定点一动点 50分

ybt

未通过

测试点结果内存时间
测试点1答案正确624KB2MS
测试点2答案正确636KB3MS
测试点3答案正确680KB2MS
测试点4答案正确688KB2MS
测试点5答案正确684KB2MS
测试点6答案错误644KB1MS
测试点7答案错误676KB2MS
测试点8答案错误688KB2MS
测试点9答案错误692KB2MS
测试点10答案错误696KB2MS

LOJ

 哪出错？估计还是在划分上，猜测应该在n比较大时，漏考虑了些情况。

反例很快找到，如下图，以六边形为例，以a为顶点的三角形划分，是之前的想法，漏考虑了，以a为顶点，以d为顶点的三角形划分。

 感觉情况挺多，复杂多变，如何思考？

难想之处在于，三角形，有两顶点相邻，但第三个顶点存在跨越，这个难办，以什么为束缚？

六边形的思考过程：

 区间动归 dp 两定点一动点 50分 代码如下：

#include <bits/stdc++.h>
#define ULL unsigned long long
#define maxn 110
#define INF 18e18
using namespace std;
ULL a[maxn],dp[maxn][maxn],mx;
int main(){
	int n,lt,rt,k,len,i,j;
	scanf("%d",&n);
	for(i=1;i<=n;i++)scanf("%llu",&a[i]),a[i+n]=a[i];
	for(i=1;i<=2*n;i++)
		for(j=i+2;j<=2*n;j++)
			dp[i][j]=INF;
	for(len=2;len<n;len++)
		for(lt=1;lt<=2*n;lt++){//三角形，三个顶点 
			rt=lt+len;//lt,rt是定点 
			if(rt>2*n)break;
			for(k=1;k<len;k++){//lt+k是动点,注意k取值从1开始 
				dp[lt][rt]=min(dp[lt][rt],dp[lt][lt+k]+dp[lt+k][rt]+a[lt]*a[lt+k]*a[rt]);
			}
		}
	mx=INF;
	for(i=1;i<=n;i++)mx=min(mx,dp[i][i+n-1]);
	printf("%llu\n",mx);
	return 0;
}

3.区间动归 int128使用(不使用高精度)100分

ybt

通过

测试点结果内存时间
测试点1答案正确652KB2MS
测试点2答案正确684KB2MS
测试点3答案正确748KB2MS
测试点4答案正确740KB2MS
测试点5答案正确772KB2MS
测试点6答案正确692KB2MS
测试点7答案正确744KB2MS
测试点8答案正确780KB2MS
测试点9答案正确788KB2MS
测试点10答案正确780KB2MS

LOJ

花了一晚上时间，翻了2021CSP复赛代码，基本确认，比赛中int128是可以使用的。

比赛建议，高精尽量不编，因没有一定的熟练程度，比赛中的成功率是比较低的。

翻看了int128的读写，之后，独立编码，区间动归 int128使用(不使用高精度)100分 代码如下：

#include <bits/stdc++.h>
#define maxn 110
#define INF 1e30
using namespace std;
__int128 a[maxn],dp[maxn][maxn],mx;
inline __int128 read(){
	__int128 x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){//读出负号，去除其他无用符号 
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while('0'<=ch&&ch<='9'){//读取数字 
		x=x*10+ch-'0';
		ch=getchar();
	}
	return x*f;
}
inline void print(__int128 x){
	if(x<0){
		putchar('-');//打印负号 
		x=-x;//变为正数 
	}
	if(x){//请注意此处是if不是while 
		print(x/10);//自高位到低位输出，采用递归形式 
		putchar(x%10+'0');//以字符形式打印 
	}
}
int main(){
	int n,lt,rt,k,len,i,j;
	scanf("%d",&n);
	for(i=1;i<=n;i++)a[i]=read(),a[i+n]=a[i];
	for(i=1;i<=2*n;i++)
		for(j=i+2;j<=2*n;j++)
			dp[i][j]=INF;
	for(len=2;len<n;len++)
		for(lt=1;lt<=2*n;lt++){//三角形，三个顶点 
			rt=lt+len;//lt,rt是定点 
			if(rt>2*n)break;
			for(k=1;k<len;k++){//lt+k是动点,注意k取值从1开始 
				dp[lt][rt]=min(dp[lt][rt],dp[lt][lt+k]+dp[lt+k][rt]+a[lt]*a[lt+k]*a[rt]);
			}
		}
	mx=INF;
	for(i=1;i<=n;i++)mx=min(mx,dp[i][i+n-1]);
	print(mx);
	return 0;
}

__int128能写成_int128吗，编码后，报错如下：

 __int128能写成int128吗，编码后，报错如下：

测试int128最大值代码如下：

#include <bits/stdc++.h>
using namespace std;
__int128 mx;
inline __int128 read(){
	__int128 x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){//读出负号，去除其他无用符号 
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while('0'<=ch&&ch<='9'){//读取数字 
		x=x*10+ch-'0';
		ch=getchar();
	}
	return x*f;
}
inline void print(__int128 x){
	if(x<0){
		putchar('-');//打印负号 
		x=-x;//变为正数 
	}
	if(x){//请注意此处是if不是while 
		print(x/10);//自高位到低位输出，采用递归形式 
		putchar(x%10+'0');//以字符形式打印 
	}
}
int main(){
	mx=1;
	printf("mx=");
	print((mx<<127)-1);
	printf("\n");
}

上述代码对应的输入输出数据如下：

mx=170141183460469231731687303715884105727

计算器中按出的2^127-1=1.7014118346046923173168730371588e+38基本吻合。

以下代码在处理数据溢出时有问题，亟待解决

#include <bits/stdc++.h>
using namespace std;
__int128 mx;
inline __int128 read(){
	__int128 x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){//读出负号，去除其他无用符号 
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while('0'<=ch&&ch<='9'){//读取数字 
		x=x*10+ch-'0';
		ch=getchar();
	}
	return x*f;
}
inline void print(__int128 x){
	if(x<0){
		putchar('-');//打印负号 
		x=-x;//变为正数 
	}
	if(x){//请注意此处是if不是while 
		print(x/10);//自高位到低位输出，采用递归形式 
		putchar(x%10+'0');//以字符形式打印 
	}
}
int main(){
	mx=1;
	printf("mx=");
	print(mx<<127);
	printf("\n");
}

请看对应的输入输出数据：

mx=--17014118346046923173168730371588410572(

开始有2个'-'号，结尾有一个'('号

上述代码证明，int128最大值是2^127-1=1.7*10^38

4.区间动归 uint128使用(不使用高精度)100分

ybt

通过

测试点结果内存时间
测试点1答案正确652KB4MS
测试点2答案正确672KB3MS
测试点3答案正确744KB2MS
测试点4答案正确736KB2MS
测试点5答案正确776KB2MS
测试点6答案正确680KB2MS
测试点7答案正确744KB2MS
测试点8答案正确776KB2MS
测试点9答案正确776KB2MS
测试点10答案正确780KB2MS

LOJ

再接再励，继续攻克uint128.

区间动归 uint128使用(不使用高精度)100分 代码如下：

#include <bits/stdc++.h>
#define maxn 110
#define INF 1e30
using namespace std;
__uint128_t a[maxn],dp[maxn][maxn],mx;
inline __uint128_t read(){
	__uint128_t x=0;
	char ch=getchar();
	while(ch<'0'||ch>'9'){//去除其他无用符号 
		ch=getchar();
	}
	while('0'<=ch&&ch<='9'){//读取数字 
		x=x*10+ch-'0';
		ch=getchar();
	}
	return x;
}
inline void print(__uint128_t x){
	if(x){//请注意此处是if不是while 
		print(x/10);//自高位到低位输出，采用递归形式 
		putchar(x%10+'0');//以字符形式打印 
	}
}
int main(){
	int n,lt,rt,k,len,i,j;
	scanf("%d",&n);
	for(i=1;i<=n;i++)a[i]=read(),a[i+n]=a[i];
	for(i=1;i<=2*n;i++)
		for(j=i+2;j<=2*n;j++)
			dp[i][j]=INF;
	for(len=2;len<n;len++)
		for(lt=1;lt<=2*n;lt++){//三角形，三个顶点 
			rt=lt+len;//lt,rt是定点 
			if(rt>2*n)break;
			for(k=1;k<len;k++){//lt+k是动点,注意k取值从1开始 
				dp[lt][rt]=min(dp[lt][rt],dp[lt][lt+k]+dp[lt+k][rt]+a[lt]*a[lt+k]*a[rt]);
			}
		}
	mx=INF;
	for(i=1;i<=n;i++)mx=min(mx,dp[i][i+n-1]);
	print(mx);
	return 0;
}

__uint128_t能写成__uint128吗，报错如下:

 __uint128_t能写成_uint128_t吗，报错如下:

uint128最大值时多少，测试代码如下：

#include <bits/stdc++.h>
#define maxn 110
#define INF 1e30
using namespace std;
__uint128_t mx;
inline __uint128_t read(){
	__uint128_t x=0;
	char ch=getchar();
	while(ch<'0'||ch>'9'){//去除其他无用符号 
		ch=getchar();
	}
	while('0'<=ch&&ch<='9'){//读取数字 
		x=x*10+ch-'0';
		ch=getchar();
	}
	return x;
}
inline void print(__uint128_t x){
	if(x){//请注意此处是if不是while 
		print(x/10);//自高位到低位输出，采用递归形式 
		putchar(x%10+'0');//以字符形式打印 
	}
}
int main(){
	mx=1;
	print((mx<<127)-1+(mx<<127));
	return 0;
}

上述代码对应的输入输出数据如下：

340282366920938463463374607431768211455

计算器中对应的2^128-1=3.4028236692093846346337460743177e+38与上述数据基本吻合。

目前在比赛中，没有见到uint128,稳妥起见，还是使用int128.

该题习得了什么？

陌生场景下的建模能力。

int128,uint128的读写，个人还是建议，比赛时要尽量避免高精度算法，尽量用int128(最大值在1.7*10^38),uint128(最大值在3.4*10^38).

目前在比赛中，没有见到uint128,稳妥起见，还是使用int128.

mrcrack
关注
专栏目录
凸多边形的划分
Averent的博客
07-06 3503
题目描述给定一具有N个顶点（从1到N编号）的凸多边形，每个顶点的权均已知。问如何把这个凸多边形划分成N-2个互不相交的三角形，使得这些三角形顶点的权的乘积之和最小？输入输出格式输入格式：第一行 顶点数N（N&lt;50）。 第二行 N个顶点（从1到N）的权值，权值为小于32768的整数。输出格式：一行，为各三角形顶点的权的乘积之和最小值。输入输出样例输入样例#1： 5 121 122 123 2...
LOJ10151. 「一本通 5.1 练习 2」分离与合体【区间DP+BFS】
dingxingdi的博客
07-08 525
【题目描述】 经过在机房里数日的切磋，LYD 从杜神牛那里学会了分离与合体，出关前，杜神牛给了他一个测试…… 杜神牛造了 nn 个区域，他们紧邻着排成一行，编号 1…n1…n。在每个区域里都放着一把 OI 界的金钥匙，每一把都有一定的价值，LYD 当然想得到他们了。然而杜神牛规定 LYD 不能一下子把他们全部拿走，而是每次只可以拿一把。为了尽快得到所有金钥匙，LYD 自然就用上了刚学的分离与合体特...
1571：【例 3】凸多边形的划分
weixin_45931113的博客
04-11 1200
时间限制: 1000 ms 内存限制: 524288 KB 提交数: 562 通过数: 150 【题目描述】 给定一个具有 N 个顶点的凸多边形，将顶点从 1 至 N 标号，每个顶点的权值都是一个正整数。将这个凸多边形划分成 N−2 个互不相交的三角形，试求这些三角形顶点的权值乘积和至少为多少。 【输入】 输入第一行为顶点数 N 第二行依次为顶点 1 至顶点 N 的权值。 ...
一本通OJ-P1571凸多边形的划分
weixin_38168139的博客
08-05 598
题目 题目链接 测试得分：　　30 主要算法：　　动态规划，区间DP，高精类 题干： 　　　　区间DP板子题 应试策略： 确定状态//f[i][j]表示的i-j这一凸边形划分的最小值 枚举区间长度FORa(l,1,n) 枚举区间开头for(int i=1,in=2*n-l+1;i<=in;i++) 状态转移f...
loj10149. 「一本通 5.1 例 3」凸多边形的划分
weixin_30315723的博客
08-16 244
思路： 　　区间dp，数值范围较大，需写高精度·。 　　以下是两种dp的方法，都可以AC #include<cstdio> #include<iostream> #include<cstdio> #include<cstring> using namespace std; const int maxn = 65; in...
区间dp(三)——#10149. 「一本通 5.1 例 3」凸多边形的划分
hhhcbw的博客，欢迎各位来访
09-18 460
题目链接：https://loj.ac/problem/10149 解题思路 凸多边形的划分也是区间dp的经典题目，dp[i][j]表示i~j进行凸多边形划分最少花费。dp[i][j]=max(dp[i][j],dp[i][k]+dp[k][j]+a[i]*a[k]*a[j])。a[]为每个点的权值。 我们选k，然后与i,j构成三角形，就可以把图分成三块，将问题分为紫色的两块子问题继续解决，且无后效性。 代码实现，肯定爆long long，我是用__int128 因为写起来简单QAQ ，当然也可以用高精度
LOJ 一本通提高篇5.1区间类动态规划 例题+练习
xxxLJYing的博客
08-03 656
是最讨厌的DP啊/哭/，但是总考不是/大哭/，我有什么办法/哭死了/。 好吧，DP还好几种：区间、树形、数位..@$#%4*（我才没有骂人）切腹自尽。 其实区间DP的用法还是比较单一的？（至少下面的题是?）我还不是一样做不出来。 目录 #10147. 「一本通 5.1 例 1」石子合并 #10148. 「一本通 5.1 例 2」能量项链 #10149. 「一本通 5.1 例 3」凸多边...
LOJ10151 分离与合体 [区间DP]
C2022YuRL的博客
07-20 231
LOJ 题意就是给你一个区间，每次进入一个区域 kkk （除区间末尾之外）。举个例子，一开始你进入区域 kkk ，就将 [1,n][1,n][1,n] 分成了 [1,k][1,k][1,k]，[k+1,n][k+1,n][k+1,n] 两个区间，还会获得（合并后所在区间左右端区域里金钥匙价值之和）∗*∗（之前分离的时候所在区域的金钥匙价值），在这里就是 （1号价值+n号价值）∗*∗ （k号价值）。然后在各自区间内再选任意一个区域发生分离，直到每个所在的区间只剩下了一个区域。 首先这道题有两个问，最大价值和方
[LOJ#2478][九省联考2018]林克卡特树（树形DP+带权二分）
女装的你，如此好看！
10-14 308
Address 洛谷P4383 BZOJ5252 LOJ#2478 Solution 简版题意：在一棵 nnn 个节点的边带权树上删掉 kkk 条边再加上 kkk 条权为 000 的边组成一棵新树，最大化新树的直径。 容易发现，如果删掉的 kkk 条边给定，那么答案就是这 k+1k+1k+1 个连通块的直径之和。 于是问题转化成在树上选出 k+1k+1k+1 条不相交的路径的最大权值和。点视为退化...
【代码超详解】LOJ #132. 树状数组 3 ：区间修改，区间查询
COFACTOR
02-07 617
一、题目描述 样例 样例输入 5 10 2 6 6 1 1 2 1 4 1 2 5 10 2 1 3 2 2 3 1 2 2 8 1 2 3 7 1 4 4 10 2 1 2 1 4 5 6 2 3 4 样例输出 15 34 32 33 50 二、算法分析说明与代码编写指导 如果不了解树状数组，建议先阅读之前的题解： https://blog.csdn.net/COFACTOR/arti...
1571：【例 3】凸多边形的划分——典型的区间dp问题，值得一看
最新发布
CylMK的博客
07-17 900
给定一个具有N个顶点的凸多边形，将顶点从1至N标号，每个顶点的权值都是一个正整数。将这个凸多边形划分成N−2个互不相交的三角形，试求这些三角形顶点的权值乘积和至少为多少。
LOJ #10149. 「一本通 5.1 例 3」凸多边形的划分 解题报告
日居月诸的博客
09-19 453
文章目录题目链接解题思路详细代码 题目链接 LOJ #10149. 「一本通 5.1 例 3」凸多边形的划分 解题思路 这道题显而易见是一道区间DP 我们将所有点逆时针编号 设f[l][r]f[l][r]f[l][r]表示从lll到逆时针rrr这一部分的图形（一个以lll、rrr和路径上其他点组成的多边形）所能达到的最大权值和。 易得： f[i][i+1]=0,else f[i][j]=...
「一本通 5.1 例 3」凸多边形的划分
bz02_2023f2的博客
06-30 357
给定一个具有 NN 个顶点的凸多边形，将顶点从 11 至 NN 标号，每个顶点的权值都是一个正整数。将这个凸多边形划分成 N-2N−2 个互不相交的三角形，试求这些三角形顶点的权值乘积和至少为多少。对于 100\%100% 的数据，有 N\le 50N≤50，每个点权值小于 10^9109。输出仅一行，为这些三角形顶点的权值乘积和的最小值。第二行依次为顶点 11 至顶点 NN 的权值。输入第一行为顶点数 NN。
LOJ P10149 凸多边形的划分 题解
weixin_30894583的博客
08-06 173
Analysis 区间dp+压位高精 dp五分钟，高精两小时 #include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #define maxn 110 #define INF 2147483647 using namespac...
背包问题求方案数量
weixin_53469462的博客
05-15 275
有 N 件物品和一个容量是 V 的背包。每件物品只能使用一次。 第 i 件物品的体积是 vi，价值是 wi。 求解将哪些物品装入背包，可使这些物品的总体积不超过背包容量，且总价值最大。 输出 最优选法的方案数。注意答案可能很大，请输出答案模 109+7的结果。 输入格式 第一行两个整数，N，V，用空格隔开，分别表示物品数量和背包容积。 接下来有 NN 行，每行两个整数 vi,wi，用空格隔开，分别表示第 i 件物品的体积和价值。 输出格式 输出一个整数，表示 方案数 模 1e9+7的结果。
凸多边形的划分（区间dp）
qq_52093121的博客
04-12 177
题目传送门 #include<bits/stdc++.h> using namespace std; typedef long long ll; const ll N=100,INF=1LL<<63,M=500; ll f[N][N]; int n; ll a[N]; bool cmp(vector<ll>& a,vector<ll>& b) { while(a.size()&&a.back()==0) a.pop
凸多边形划分
Rayhirox的博客
02-12 1219
又是卡特兰数。这道题需要注意的是一条边和两条边是不能构成多边形的，因此我们的计算应该从3开始，即将n减2，再算出相应的卡特兰数代码：#include&lt;iostream&gt; #include&lt;cstdio&gt; #include&lt;algorithm&gt; #include&lt;cstring&gt; #include&lt;string&gt; #include&lt;c...
LOJ-6279-数列分块入门3(分块, 二分)
05-24
题目描述 给定一个长度为n的序列a，有m次询问。每次询问给定l, r, k，要求区间[l, r]中第k小的数。 输入格式 第一行两个整数n,m。 接下来一行n个整数，表示序列a。 接下来m行，每行三个整数l,r,k。 输出格式 对于每个询问，输出一行一个数，表示区间[l, r]中第k小的数。 输入样例 5 3 1 3 2 4 5 1 5 2 2 4 1 3 5 3 输出样例 2 3 5 算法1 (分块) $O(n\sqrt n)$ 我们可以预处理出每一个块的信息，每个块又分为若干小块，每个小块中的数都小于等于块中的最大值，可以用二分查找确定小块中满足条件的数的个数，然后用根号分治处理。 时间复杂度 每次查询的时间复杂度是 $O(\sqrt n\log n)$ 的，总时间复杂度是 $O(n\sqrt n\log n)$ 的。 参考文献 算法竞赛进阶指南 C++ 代码 算法2 (二分答案) $O(n\log^2n)$ 我们二分答案，然后统计序列中小于等于这个答案的数的个数，统计方法可以用归并排序，在时间复杂度 $O(n\log n)$ 的时间内完成。 时间复杂度 每次查询的时间复杂度是 $O(n\log^2n)$ 的，总时间复杂度是 $O(mn\log^2n)$ 的。 参考文献 算法竞赛进阶指南 C++ 代码
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值