BZOJ 3203 Sdoi2013 保护出题人 凸包+三分

12 篇文章 0 订阅
2 篇文章 0 订阅

题目大意:太长自己看
sumi 表示第 i 个僵尸以及之前的僵尸的体力总和,disi表示第 i 个僵尸与房屋的初始距离
我们发现我们能消灭一个僵尸当且仅当y>=sumidisi
那么我们要求的显然就是 max{sumidisi}
我们将一个僵尸抽象成一个点 sumidisi ,那么我们发现每个回合僵尸之间的相对位置是不变的
因此我们可以维护一个凸包,三分即可

#include <cstdio>
#include <cstring>
#include <iomanip>
#include <iostream>
#include <algorithm>
#define M 100100
#define EPS 1e-7
using namespace std;
typedef long double ld;
struct Point{
    ld x,y;
    Point() {}
    Point(ld _,ld __):
        x(_),y(__) {}
    friend ld Get_Slope(const Point &p1,const Point &p2)
    {
        return (p1.y-p2.y)/(p1.x-p2.x);
    }
}O,stack[M];int top;
int n;
long long d,a[M],sum[M];
ld ans;
void Insert(const Point &p)
{
    while( top>=2 && Get_Slope(stack[top-1],stack[top]) >= Get_Slope(stack[top],p) - EPS )
        top--;
    stack[++top]=p;
}
ld Trisection()
{
    int i,l=1,r=top;
    while(r-l>=3)
    {
        int l_mid=(l+l+r)/3,r_mid=(l+r+r)/3;
        ld _l=Get_Slope(stack[l_mid],O);
        ld _r=Get_Slope(stack[r_mid],O);
        if(_l>_r)
            r=r_mid;
        else
            l=l_mid; 
    }
    ld re=0;
    for(i=l;i<=r;i++)
        re=max(re,Get_Slope(stack[i],O));
    return re;
}
int main()
{
    int i;
    long long x;
    cin>>n>>d;
    for(i=1;i<=n;i++)
    {
        #ifdef PoPoQQQ
            scanf("%I64d%I64d",&a[i],&x);
        #else
            scanf("%lld%lld",&a[i],&x);
        #endif
        sum[i]=sum[i-1]+a[i];
        Insert(Point(i*d,sum[i-1]));
        O=Point(i*d+x,sum[i]);
        ans+=Trisection();
    }
    cout<<fixed<<setprecision(0)<<ans<<endl;
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值