POJ 1948 Triangular Pastures (二维01背包)

题目链接：http://poj.org/problem?id=1948

题目大意：给最多40根木棍，每根长度不超过40，要用完所有的木棍构成面积最大的三角形，求出最大的面积。

思路：设f[i][j]为第一条边长度为i，第二条边长度为j是否可行

状态转移方程：

if(j>=a[i])

f[j][k]=f[j][k]||f[j-a[i]][k]

if(k>=a[i])

f[j][k]=f[j][k]||f[j][k-a[i]]

#include
#include
#include
using namespace std;

int f[810][810];
int a[50];

void swap(int &a, int &b){
int tmp = a;
a = b;
b = tmp;
}

int san(int a,int b,int c)
{
if(a > b)   swap(a,b);
if(a > c)   swap(a,c);
if(b > c)   swap(b,c);

if(a + b < c)   return -1;
double p = (a + b + c) * 1.0 / 2;
double ret;
ret = sqrt(p * (p-a) * (p-b) * (p-c)) * 100;

return (int)ret;
}

int main(){
int n;
scanf("%d",&n);
int sum = 0;
for (int i = 0; i < n; i ++){
scanf("%d",&a[i]);
sum += a[i];
}
memset(f, 0, sizeof(f));
f[0][0] = 1;
int ban = sum / 2;
for (int k = 0; k < n; k ++){
for (int i = ban; i >= 0; i --){
for (int j = ban; j >= 0; j --){
if (i - a[k] >= 0){
f[i][j] = f[i][j] | f[i - a[k]][j];
}
if (j - a[k] >= 0){
f[i][j] = f[i][j] | f[i][j - a[k]];
}
}
}
}
int max = -1;
for (int i = 1; i <= ban; i ++){
for (int j = 1; j <= ban; j ++){
if (f[i][j]){
if (san(i, j, sum - i - j) > max){
max = san(i, j, sum - i - j);
}
}
}
}
printf("%d\n", max);
return 0;
}