bzoj1132 [POI2008]Tro 计算几何

题意:有n个点,求任意三个点组成的三角形面积之和。

一眼叉积,然后没有把式子写下来所以没有发现前缀和性质= =。。。
然后又被斜率卡了,用double直接除就会炸,int强制转double就能对是什么鬼。。
具体的话就是用叉积计算面积,对于一个点k来说,他的贡献等于(sum.x-k.y)*(sum.y-k.x),那么枚举第二个点,加进来按照斜率排个序以后直接做。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e5+5;
typedef long long ll;
struct node
{
    int x,y;
    double k;
    node(){}
    node(int a,int b):x(a),y(b)
    {
        k=x?(double)y/x:1e10;
    }
    bool operator < (const node &Y)const
    {
        if (x!=Y.x)return x<Y.x;
        return y<Y.y;
    }
}a[N],p[N];
int n,tot;
ll ans;
bool cmp(const node &x,const node &y)
{
    return x.k<y.k;
}
int main()
{
    scanf("%d",&n);
    fo(i,1,n)scanf("%d%d",&a[i].x,&a[i].y);
    sort(a+1,a+1+n);
    fo(i,1,n)
    {
        ll sumx=0,sumy=0;
        tot=0;
        fo(j,i+1,n)p[++tot]=node(a[j].x-a[i].x,a[j].y-a[i].y);
        sort(p+1,p+1+tot,cmp);
        fo(j,1,tot)
        {
            ans+=sumx*p[j].y-sumy*p[j].x;
            sumx+=p[j].x;
            sumy+=p[j].y; 
        }
    }
    printf("%lld.%d\n",ans>>1,ans&1?5:0);  
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值