poj 1948 dp（组成面积最大的三角形）

题意：给定n(n<=40)个木棍的长度（长度<=40），现在要求将所有的木棍都用上，求能围成的最大的三角形的面积。

思路：首先有了三角形的三边长，求面积显然用海伦公式。其次要考虑怎么组成三角形，想了半天觉得复杂度都会超。后来看题解，发现做法基本都是O(n*sum*sum)即O(40*800*800)的，按说一分钟的时限应该是会超的。反正如果是这样就好写了。dp[i][j]表示用边长i和j（其余一个的边长也就固定了）能否组成三角形。dp只维护左下三角形。最后把能够组成三角形的方案扫一遍求面积即可。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <queue>
using namespace std;
#define INF 0x3fffffff
#define clr(s,t) memset(s,t,sizeof(s))
#define N 45
int s[N],n,p,dp[805][805],sum=0;
int res=0;
int check(int x,int y){
    return x>=y && x<=sum;
}
double area(int x,int y){
    int z = p-x-y;
    double tmp = p/2.;
    if(p<=2*x || p<=2*y || p<=2*z)
        return 0;
    return sqrt(tmp*(tmp-x)*(tmp-y)*(tmp-z));
}
int main(){
    int i,j,k;
    scanf("%d",&n);
    for(i = 1;i<=n;i++){
        scanf("%d",&s[i]);
        sum += s[i];
    }
    p = sum;
    sum >>= 1;
    clr(dp, 0);
    dp[0][0] = 1;
    for(k = 1;k<=n;k++)
        for(i = sum;i>=0;i--)
            for(j = i;j>=0;j--)
                if(dp[i][j]){
                    if(check(i+s[k],j))
                        dp[i+s[k]][j] = 1;
                    if(check(i,j+s[k]))
                        dp[i][j+s[k]] = 1;
                    if(check(j+s[k],i))
                        dp[j+s[k]][i] = 1;
                }
    for(i = 1;i<=sum;i++)
        for(j = 1;j<=i;j++)
            if(dp[i][j])
                res = max(res,(int)(100*area(i,j)));
    printf("%d\n",res?res:-1);
    return 0;
}