hdu5135 Little Zu Chongzhi's Triangles 状压dp

前几天突然得知要去打一站区域赛，时间紧迫，只能拿前几年区域赛的真题来练习了。

此题是2014年广州站第二水题，开始看错了题意，以为两个三角形可以有公共边，后来看了题解发现是两个三角形不可以有公共边

状压dp，数据范围小，直接用1<<12位做状压，用第几根木棍就把第几位设为1，总状态也不多，保存合法状态，直接dp就好了，dp的时候注意公共边的要舍去

AC代码

<span style="font-family:KaiTi_GB2312;font-size:18px;">#include <iostream>
#include <queue>
#include <stack>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
double dp[1 << 12];//只有12根木棍，所以只需要1左移12位做状态
int n,a[12];
double cal(int a, int b, int c)
{
    if(a + b <= c) return 0.0;
    double p;
    p = (a + b + c) * 0.5;
    double ans = sqrt(p*(p-a)*(p-b)*(p-c));
    return ans;
}
vector<int> v;    //用来存放合法状态
int main()
{
    while(scanf("%d", &n) && n)
    {
        memset(dp, 0, sizeof(dp));
        for(int i = 0; i < n; i++)
            scanf("%d", &a[i]);
        v.clear();
        for(int i = 0; i < n; i++)
            for(int j = i + 1; j < n; j++)
            {
                for(int  k= j + 1; k < n; k++)
                {
                    int st = (1<<i)|(1<<j)|(1<<k);    //获取当前状态的位置
                    dp[st] = cal(a[i],a[j],a[k]);
                    if(a[i] + a[j] > a[k]) v.push_back(st);
                }
            }
        for(int i = 0; i < (1 << n); i++)     //枚举所有状态
            for(int j = 0; j < v.size(); j++)
            {
                if(i & v[j]) continue;        //因为一根木棍不能属于两个三角形
                dp[i|v[j]] = max(dp[i|v[j]], dp[i] + dp[v[j]]);
            }
        printf("%.2lf\n",dp[(1<<n)-1]);     //最终状态，用上了所有木棍
    }
}
</span>