noi.ac NA531 【神树和物品】

今日成就:本来以为过了这题,然后被mcfx发现写假并针对地造了一组hack数据之后FST了。

复杂度什么的咱也不会证,咱也不会卡,被hack之后只能FST。

是个决策单调性sb题,但是由于太菜不怎么会写决策单调性于是胡搞了一波,然后还是被mcfx卡了\(\mathcal{O}(n^2)\),成功\(100\rightarrow44\)

单调队列写丑了我也没什么办法...sb挂sb题...

没啥说的了,直接扔代码:

#include<bits/stdc++.h>

using namespace std;

const int N=3e5+5;

int n,a[N],x[N],y[N];

void qread(int &xx){
	xx=0;int ch=getchar();
	while(ch<'0'||ch>'9'){
		ch=getchar();
	}
	while(ch>='0'&&ch<='9'){
		xx=xx*10+ch-'0';
		ch=getchar();
	}
}

long long f[N];

long long calc(int l,int r){
	int u=abs(a[r]-x[l]);
	return 1LL*u*u*u+1LL*y[l]*y[l]*y[l];
}

int dq[N],dqfr=1,dqen;

int bs(int l,int r,int pos1,int pos2){
	int mid;
	while(l<r){
		mid=(l+r)>>1;
		f[pos1]+calc(pos1+1,mid)>=f[pos2]+calc(pos2+1,mid)?r=mid:l=mid+1;
	}
	return l;
}

int main(){
	qread(n);
	for(int i=1;i<=n;i++){
		qread(a[i]);
	}
	for(int i=1;i<=n;i++){
		qread(x[i]);
	}
	for(int i=1;i<=n;i++){
		qread(y[i]);
	}
	dq[++dqen]=0;
	for(int i=1;i<=n;i++){
		while(dqfr<dqen&&f[dq[dqfr]]+calc(dq[dqfr]+1,i)>=f[dq[dqfr+1]]+calc(dq[dqfr+1]+1,i)){
			++dqfr;
		}
		f[i]=f[dq[dqfr]]+calc(dq[dqfr]+1,i);
		while(dqen>dqfr&&bs(1,n+1,dq[dqen],i)<=bs(1,n+1,dq[dqen-1],dq[dqen])){
			--dqen;
		}
		dq[++dqen]=i;
	}
	printf("%lld\n",f[n]);
	return 0;
}
posted @ 2019-07-05 15:10  __BLUESKY007  阅读(328)  评论(1编辑  收藏  举报