Triangular Pastures POJ - 1948

Triangular Pastures POJ - 1948

sum表示木条的总长。a[i]表示第i根木条长度。ans[i][j][k]表示用前i条木条，摆成两条长度分别为j和k的边是否可能。

那么ans[i][j][k]=ans[i-1][j-a[i]][k] || ans[i-1][j][k-a[i]]

可以用滚动数组优化。

最后在ans[n]中枚举i和j，如果ans[n][i][j]为true，再算出sum-i-j的长度，判断是否存在分别以i,j,sum-i-j为三边长的三角形。如果存在，再用海伦公式算出三角形面积，用面积去更新最大面积。

错误记录：
没有在dp的时候判边界，也就是没有判j>=a[i]和k>=a[i]。

 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 typedef long long LL;
 bool ans[][];
 LL a[],n,sum;
 double p,anss;
 int main()
 {
     LL i,j,k,x,y,z;
     scanf("%lld",&n);
     for(i=;i<=n;i++)
     {
         scanf("%lld",&a[i]);
         sum+=a[i];
     }
     ans[][]=;
     for(i=;i<=n;i++)
         for(j=sum;j>=;j--)
             for(k=sum;k>=;k--)
             {
                 if(j>=a[i])
                     ans[j][k]|=ans[j-a[i]][k];
                 if(k>=a[i])
                     ans[j][k]|=ans[j][k-a[i]];
             }
     for(x=;x<=sum;x++)
         for(y=x;y<=sum-x-;y++)
         {
             if(!ans[x][y])    continue;
             z=sum-x-y;
             if(y>z)    break;
             if(x+y<=z) continue;
             p=(double)(x+y+z)/2.0;
             anss=max(anss,sqrt(p*(p-x)*(p-y)*(p-z)));
         }
     if(anss!=0.0)
         printf("%lld",(long long)(anss*(double)));
     else
         printf("-1");
     return ;
 }