hdu 4573 Throw the Stones

本文介绍了一种计算三维空间中凸包体积的方法。通过定义结构体来表示顶点和面,并利用向量运算来判断点的位置关系,最终实现对三维凸包体积的有效计算。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

求三维凸包,其中求体积的部分值得借鉴

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define MAXN 10005
const double eps = 1e-9;
int n, fc_cnt, leg[MAXN][MAXN];
double mx_ff, dis_ff;

inline int sig(double x)
{
    return (x>eps) - (x<-eps);
}
struct node
{
    double x, y, z;
    node(double a = 0.0, double b = 0.0, double c = 0.0):x(a),y(b),z(c) {}
    
    node operator - (node& a) const
    {
        return node(x-a.x, y-a.y, z-a.z);
    }
    node operator + (node& a) const
    {
        return node(x+a.x, y+a.y, z+a.z);
    }
    node operator ^ (node a) const // this * a * sin($)
    {
        return node(y*a.z-a.y*z, z*a.x-a.z*x, x*a.y-a.x*y);
    }
    double operator * (node a) const // this * a * cin($)
    {
        return x*a.x + y*a.y + z*a.z;
    }
    node cross(node a, node b)
    {
        return a - *this^b - *this;
    }
    double L()
    {
        return sqrt(x*x + y*y + z*z);
    }
    double vv6(node a, node b)
    {
        return (a^b)*(*this);
    }
} P[MAXN];
struct FACE
{
    int a, b, c;
    int isOK;
    FACE(int x = 0, int y = 0, int z = 0, int ik = 0):a(x),b(y),c(z), isOK(ik){}
    int visit(node mm)
    {
        return sig(P[a].cross(P[b], P[c])*(mm - P[a])) > 0;
    }
}facc[MAXN*8];
void deal(int o, int a, int b);
void dfs(int f, int o)
{
    facc[f].isOK = 0;
    dis_ff -= P[facc[f].a].vv6(P[facc[f].b], P[facc[f].c]);
    deal(o, facc[f].b, facc[f].a);
    deal(o, facc[f].a, facc[f].c);
    deal(o, facc[f].c, facc[f].b);
}
void deal(int o, int a, int b)
{
    int i = leg[a][b];
    if (facc[i].isOK)
    {
        if (facc[i].visit(P[o]))
            dfs(i, o);
        else
        {
            dis_ff += P[b].vv6(P[a], P[o]);
            FACE add(b, a, o, 1);
            leg[b][a] = leg[a][o] = leg[o][b] = fc_cnt;
            facc[fc_cnt++] = add;
        }
    }
}
void _update(int idx)
{
    for (int i = 0; i< fc_cnt; i++)
    {
        if (facc[i].isOK && facc[i].visit(P[idx]))
        {
            dfs(i, idx);
            break;
        }
    }
}
double getV()
{
    double v = 0;
    for (int i = 0; i<fc_cnt; i++)
    {
        if (facc[i].isOK)
        {
            v += P[facc[i].a].vv6(P[facc[i].b], P[facc[i].c]);
        }
    }
    return v/6;
}
void solve()
{
    mx_ff = dis_ff = 0;
    int kg = 1, flg = 1;
    int i, j;
    for (i = 1; i< n; i++)
    {
        if (flg == 1)
        {
            flg += sig((P[0]-P[i]).L());
            if (flg > 1)
                swap(P[1], P[i]);
        }
        else if (flg == 2)
        {
            flg += sig((P[0].cross(P[1], P[i])).L());
            if (flg > 2)
                swap(P[i], P[2]);
        }
        else if (flg == 3)
        {
            flg += sig(P[0].cross(P[1], P[2]) * (P[i]-P[0])) != 0;
            if (flg > 3)
            {
                swap(P[3], P[i]);
                fc_cnt = 0;
                for (j = 0; j< 4; j++)
                {
                    FACE add((j+1)%4, (j+2)%4, (j+3)%4, 1);
                    if (add.visit(P[j]))
                        swap(add.b, add.c);
                    leg[add.a][add.b] = leg[add.b][add.c] = leg[add.c][add.a] = fc_cnt;
                    facc[fc_cnt++] = add;
                }
                for (j = 4; j<= i; j++)
                {
                    _update(j);
                }
                kg = i+1;
                mx_ff = getV();
            }
        }
        else
        {
            dis_ff = 0;
            _update(i);
            dis_ff /= 6;
            if (sig(dis_ff - mx_ff) > 0)
            {
                mx_ff = dis_ff;
                kg = i+1;
            }
        }
    }
    printf("%d %.2lf\n", kg, mx_ff);
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
#endif
    int cs = 0;
    while (scanf("%d", &n) != EOF)
    {
        for (int i = 0; i< n; ++i)
        {
            scanf("%lf%lf%lf", &P[i].x, &P[i].y, &P[i].z);
        }
        printf("Case #%d:\n", ++cs);
        solve();
    }
    return 0;
}


资源下载链接为: https://pan.quark.cn/s/22ca96b7bd39 在 IT 领域,文档格式转换是常见需求,尤其在处理多种文件类型时。本文将聚焦于利用 Java 技术栈,尤其是 Apache POI 和 iTextPDF 库,实现 doc、xls(涵盖 Excel 2003 及 Excel 2007+)以及 txt、图片等格式文件向 PDF 的转换,并实现在线浏览功能。 先从 Apache POI 说起,它是一个强大的 Java 库,专注于处理 Microsoft Office 格式文件,比如 doc 和 xls。Apache POI 提供了 HSSF 和 XSSF 两个 API,其中 HSSF 用于读写老版本的 BIFF8 格式(Excel 97-2003),XSSF 则针对新的 XML 格式(Excel 2007+)。这两个 API 均具备读取和写入工作表、单元格、公式、样式等功能。读取 Excel 文件时,可通过创建 HSSFWorkbook 或 XSSFWorkbook 对象来打开相应格式的文件,进而遍历工作簿中的每个 Sheet,获取行和列数据。写入 Excel 文件时,创建新的 Workbook 对象,添加 Sheet、Row 和 Cell,即可构建新 Excel 文件。 再看 iTextPDF,它是一个用于生成和修改 PDF 文档的 Java 库,拥有丰富的 API。创建 PDF 文档时,借助 Document 对象,可定义页面尺寸、边距等属性来定制 PDF 外观。添加内容方面,可使用 Paragraph、List、Table 等元素将文本、列表和表格加入 PDF,图片可通过 Image 类加载插入。iTextPDF 支持多种字体和样式,可设置文本颜色、大小、样式等。此外,iTextPDF 的 TextRenderer 类能将 HTML、
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值