ACM模版

描述


题解

计算几何问题……不难，就是麻烦，精度问题也需要着重注意，注意人家输入精确到 10−12 ，而不是拼接时精确到 10−12 ！

我是通过判断面积是否可以构成正方形（与最大边符合），三角形是否有五个，四边形是否有两个，七个多边形一共23条边排序后是否符合七巧板的规格等等，当然，我这个也不是最简单的，完全可以通过判断多边形之间各边关系来判定是否可行，具体的需要自己去发掘了，七巧板中边与边的关系实在是丰富，O(∩_∩)O哈哈~

记得上一次听到七巧板这个名词大概是在小学三年级了吧，好怀念啊！

代码

#include <iostream>
#include <algorithm>
#include <cmath>

using namespace std;

const int MAXN = 7;
const int MAXA = 23;
const double ESP = 1e-7;

double res = 0;

struct Lpoint
{
    double x, y;
};  // 点

/* * 要求按照逆时针方向输入多边形顶点 * 可以是凸多边形或凹多边形 */
double area_of_polygon(int vcount, Lpoint plg[])
{
    int i;
    double s;
    if (vcount < 3)
    {
        return 0;
    }
    s = plg[0].y * (plg[vcount - 1].x - plg[1].x);
    for (i = 1; i < vcount; i++)
    {
        s += plg[i].y * (plg[(i - 1)].x - plg[(i + 1) % vcount].x);
    }
    return s / 2;
}

Lpoint plg[MAXN];
double area[MAXA];

int main(int argc, const char * argv[])
{
    int vcount;
    int cnt_3 = 0;
    int cnt_4 = 0;
    int cnt_area = 0;
    for (int i = 0; i < MAXN; i++)
    {
        cin >> vcount;
        if (vcount == 3)
        {
            cnt_3++;
        }
        else
        {
            cnt_4++;
        }
        for (int j = 0; j < vcount; j++)
        {
            cin >> plg[j].x >> plg[j].y;
        }
        for (int j = 0; j < vcount; j++)
        {
            area[cnt_area++] = pow(plg[j].x - plg[(j + 1) % vcount].x, 2)
                             + pow(plg[j].y - plg[(j + 1) % vcount].y, 2);
        }
        res += area_of_polygon(vcount, plg);
    }

    if (cnt_3 != 5 || cnt_4 != 2)
    {
        cout << "NO\n";
        return 0;
    }

    sort(area, area + cnt_area);

    if (fabs(area[9] - area[0]) > ESP || fabs(area[10] - area[9]) < ESP)
    {
        cout << "NO\n";
        return 0;
    }
    if (fabs(area[15] - area[10]) > ESP || fabs(area[16] - area[15]) < ESP)
    {
        cout << "NO\n";
        return 0;
    }
    if (fabs(area[20] - area[16]) > ESP || fabs(area[21] - area[20] < ESP))
    {
        cout << "NO\n";
        return 0;
    }
    if (fabs(area[22] - area[21]) > ESP)
    {
        cout << "NO\n";
        return 0;
    }

    puts(fabs(area[MAXA - 1] - res) < ESP ? "YES" : "NO");

    return 0;
}

参考

《Liuctic计算几何库》