洛谷 P1284 三角形牧场

题目描述

和所有人一样，奶牛喜欢变化。它们正在设想新造型的牧场。奶牛建筑师Hei想建造围有漂亮白色栅栏的三角形牧场。她拥有N(3≤N≤40)块木板，每块的长度Li(1≤Li≤40)都是整数，她想用所有的木板围成一个三角形使得牧场面积最大。

请帮助Hei小姐构造这样的牧场，并计算出这个最大牧场的面积。

输入输出格式

输入格式：

第1行：一个整数N

第2..N+1行：每行包含一个整数，即是木板长度。

输出格式：

仅一个整数：最大牧场面积乘以100然后舍尾的结果。如果无法构建，输出-1。

输入输出样例

输入样例#1：

5
1
1
3
3
4

输出样例#1：

692

说明

样例解释：692=舍尾后的（100×三角形面积），此三角形为等边三角形，边长为4。

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

DP+海伦公式~

暴搜可以过4个点，加读入优化5个。

第一次见到这样的DP，居然被分到普及+……我还是好好学动规吧……

用f[i][j]表示是否可以组成两条边，长度分别是i和j。然后类似背包DP，海伦公式计算面积，好像并没有卡精度。

（刚开始以为WA了6个点是因为没开long long，事实证明加上以后还是WA……这道题用int就可以了，请自动忽视long long~读入优化也可以不要~）

#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;
#define ll long long

int n,a[41],tot,f[1601][1601],kkz;
ll ans;

ll cal(double u,double v,double k)
{
    double kkz=(u+v+k)/2;
    return sqrt(kkz*(kkz-u)*(kkz-v)*(kkz-k))*100;
}

int pan(int u,int v,int k)
{
	int totnum=u+v+k,maxx=max(u,max(v,k));
	if(maxx*2<totnum) return 1;
	return 0;
}

int read()
{
    int f=1,totnum=0;char s=getchar();
    while(s<'0' || s>'9') {if(s=='-') f=-1;s=getchar();}
    while(s>='0' && s<='9') {totnum=totnum*10+s-'0';s=getchar();}
    return f*totnum;
}

int main()
{
    n=read();
    for(int i=1;i<=n;i++)
    {
        a[i]=read();tot+=a[i];
    }
    kkz=tot/2;f[0][0]=1;ans=-1;
    for(int z=1;z<=n;z++)
      for(int i=kkz;i>=0;i--)
        for(int j=kkz;j>=0;j--)
        {
        	if(i>=a[z]) f[i][j]|=f[i-a[z]][j];
        	if(j>=a[z]) f[i][j]|=f[i][j-a[z]];
        	if(f[i][j] && pan(i,j,tot-i-j)) ans=max(ans,cal(i,j,tot-i-j));
		}
	printf("%lld\n",ans);
    return 0;
}