凸多边形三角划分(HNOI’97)

一、试题描述

给定一个具有N(N<50)个顶点(从1到N编号)的凸多边形,每个顶点的权均已知。问如何把这个凸多边形划分成N-2个互不相交的三角形,使得这些三角形顶点的权的乘积之和最小?

输入文件:第一行 顶点数N

          第二行 N个顶点(从1到N)的权值

输出格式:最小的和的值

          各三角形组成的方式

输入示例:5

121 122 123 245 231

输出示例:The minimum is :12214884

          The formation of 3 triangle:

          3 4 5, 1 5 3, 1 2 3

二、试题分析

这是一道很典型的动态规划问题。设F[I,J]I<J)表示从顶点I到顶点J的凸多边形三角剖分后所得到的最大乘积,我们可以得到下面的动态转移方程:

F[I,J]=Min{F[I,K]+F[K,J]+S[I]*S[J]*S[K]}     (I<K<J)

目标状态为:F[1,N]

但我们可以发现,由于这里为乘积之和,在输入数据较大时有可能超过长整形甚至实形的范围,所以我们还需用高精度计算,但这是大家的基本功,程序中就没有写了,请读者自行完成。

三、参考程序

Var S          :Array[1..50] Of Integer;

    F          :Array[1..50,1..50] Of Comp; 

    D           :Array[1..50,1..50] Of Byte;

    N            :Integer;

Procedure Init;    (输入数据)

Var I          :Integer;

Begin

Readln(N);

For I:=1 To N Do Read(S[I]);

End;

Procedure Dynamic;   (动态规划)

Var I,J,K      :Integer;

Begin

For I:=1 To N Do

    For J:=I+1 To N Do F[I,J]:=Maxlongint; (赋初始值)

For I:=N-2 Downto 1 Do

    For J:=I+2 To N Do

      For K:=I+1 To J-1 Do 

        If (F[I,J]>F[I,K]+F[K,J]+S[I]*S[J]*S[K]) Then

Begin

            F[I,J]:=F[I,K]+F[K,J]+S[I]*S[J]*S[K];

            D[I,J]:=K;     (记录父节点)

          End;          

End;

Procedure Print(I,J:Integer); (输出每个三角形)

Begin

If J=I+1 Then Exit;

Write(',',I,' ',J,' ',D[I,J]);

Out(I,D[I,J]);

Out(D[I,J],J);

End;

Procedure Out; (输出信息)

Begin

Assign(Output,'Output.Txt'); Rewrite(Output);

Writeln('The minimum is :',F[1,N]:0:0);

Writeln('The formation of ',N-2,' triangle:');

Write(1,' ',N,' 'D[1,N]);

Out(1,D[1,N]);

Out(D[1,N],N);

Close(Output);

End;

Begin (主程序)

Init;

Dynamic;

Out;

End.

 

根据引用\[1\]和引用\[2\]的描述,题目中的影魔拥有n个灵魂,每个灵魂有一个战斗力ki。对于任意一对灵魂对i,j (i<j),如果不存在ks (i<s<j)大于ki或者kj,则会为影魔提供p1的攻击力。另一种情况是,如果存在一个位置k,满足ki<c<kj或者kj<c<ki,则会为影魔提供p2的攻击力。其他情况下的灵魂对不会为影魔提供攻击力。 根据引用\[3\]的描述,我们可以从左到右进行枚举。对于情况1,当扫到r\[i\]时,更新l\[i\]的贡献。对于情况2.1,当扫到l\[i\]时,更新区间\[i+1,r\[i\]-1\]的贡献。对于情况2.2,当扫到r\[i\]时,更新区间\[l\[i\]+1,i-1\]的贡献。 因此,对于给定的区间\[l,r\],我们可以根据上述方法计算出区间内所有下标二元组i,j (l<=i<j<=r)的贡献之和。 #### 引用[.reference_title] - *1* *3* [P3722 [AH2017/HNOI2017]影魔(树状数组)](https://blog.csdn.net/li_wen_zhuo/article/details/115446022)[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^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [洛谷3722 AH2017/HNOI2017 影魔 线段树 单调栈](https://blog.csdn.net/forever_shi/article/details/119649910)[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^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值