bzoj4850[Jsoi2016]灯塔

省选题可以到网上找一下题意

我们可以先假设i<j,再假设i>j

就能神奇的推出一些式子

就能做出来了

大概就是用单调栈+二分

用一个结构体(p,l,r)表示在(l,r)区间中p最优

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
inline int max(int x,int y)
{
	return x>y?x:y;
}
inline int abs(int x)
{
	return max(x,-x);
}
inline void swap(int &x,int &y)
{
	x^=y^=x^=y;
}
inline int read()
{
	int x=0,f=1,ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-'){f=-1;}ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int a[100005];
inline double calc(int i,int j)
{
	return a[i]+sqrt(abs(i-j))-a[j];
}
struct stu
{
	int p,l,r;
	stu(){}
	stu(int p,int l,int r):p(p),l(l),r(r){}
};
inline int find(stu t,int x)
{
	int l=t.l,r=t.r;
	while(l<r)
	{
		int mid=(l+r)>>1;
		if(calc(t.p,mid)<=calc(x,mid)) r=mid;
		else l=mid+1;
	}
	return l;
}
int n;
stu q[100005];
inline void solve(double dp[])
{
	int head=1,tail=0;
	for(int i=1;i<=n;i++)
	{
		q[head].l++;
		if(head<=tail&&q[head].r<q[head].l) head++;
		if(head>tail||(calc(q[tail].p,n)<calc(i,n)))
		{
			while(head<=tail&&calc(q[tail].p,q[tail].l)<calc(i,q[tail].l)) tail--;
			if(head>tail) q[++tail]=stu(i,i,n);
			else
			{
				int t=find(q[tail],i);
				q[tail].r=t-1;
				q[++tail]=stu(i,t,n);
			}
		}
		dp[i]=calc(q[head].p,i);
	}
}
double f[100005],g[100005];
int main()
{
	n=read();
	for(int i=1;i<=n;i++) a[i]=read();
	solve(f);
	for(int i=1;i<=n/2;i++) swap(a[i],a[n-i+1]);
	solve(g);
	for(int i=1;i<=n;i++) printf("%d\n",max(0,ceil(max(f[i],g[n-i+1]))));
	return 0;
}

 

转载于:https://www.cnblogs.com/gcyyzf/p/6772608.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值