标签：ch return int tag P2900 Land inline calc Acquisition

P2900 [USACO08MAR]Land Acquisition G

斜率优化 dp 。

首先，我们可以发现这个题可以贪心：对于一个长度和宽度都小于等于另外一个地的地，我们可以直接划分到另外那一个去，直接把当前这个忽略即可。

然后我们可以发现这个求一个最大值可以直接写成一个 dp 。

设按\(h\)降序排序后： \(dp[i]\) 表示把前 \(i\) 个点划分完毕的代价。

很容易发现此时 \(h\) 单调递减，\(w\) 单调递增。

方程就是 \(dp[i]=\min{dp[j]+(h[j+1]*w[i])}\) 。

可以直接李超树优化即可。

这里写了两个代码：

Code1

#include<bits/stdc++.h>
using namespace std;
template <typename T>
inline void read(T &x){
	x=0;bool f=false;char ch=getchar();
	while(!isdigit(ch)){f=(ch=='-');ch=getchar();}
	while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
	return ;
}
template <typename T>
inline void write(T x){
	if(x<0) x=-x,putchar('-');
	if(x>9) write(x/10);
	putchar(x%10^48);
	return ;
}
const int N=1e6+5,M=1e6;
#define int long long
struct node{
	int h,w;
	inline bool operator < (const node &B){return h==B.h?w>B.w:h>B.h;}
}a[N];
int n,m,top;
int tag[N<<2],k[N],b[N];
#define calc(i,x) (k[i]*x+b[i])
void Modify(int x,int l,int r,int d){
	if(l==r){
		if(calc(tag[x],l)>calc(d,l)) tag[x]=d;
		return ;
	}
	if(!tag[x]) return tag[x]=d,void();
	int mid=l+r>>1;
	if(calc(tag[x],mid)>calc(d,mid)){
		if(k[tag[x]]>k[d]) Modify(x<<1,l,mid,tag[x]),tag[x]=d;
		else Modify(x<<1|1,mid+1,r,tag[x]),tag[x]=d;
	}
	else{
		if(k[tag[x]]>k[d]) Modify(x<<1|1,mid+1,r,d);
		else Modify(x<<1,l,mid,d);
	}
	return ;
}
int Query(int x,int l,int r,int v){
	if(l==r) return calc(tag[x],v);
	int mid=l+r>>1,res=calc(tag[x],v);
	if(v<=mid) res=min(res,Query(x<<1,l,mid,v));
	else res=min(res,Query(x<<1|1,mid+1,r,v));
	return res;
}
//dp[i]=dp[j]+(a[i].w*a[j].h)
int dp[N];
signed main(){
	read(n);
	for(int i=1;i<=n;i++) read(a[i].h),read(a[i].w);
	sort(a+1,a+n+1);top=1;
	for(int i=1;i<=n;i++) if(a[top].w<a[i].w) a[++top]=a[i];
	n=top;k[0]=b[0]=M;
	k[n+1]=a[1].h,b[n+1]=0;
	Modify(1,1,M,n+1);
	for(int i=1;i<=n;i++){
		dp[i]=Query(1,1,M,a[i].w);
		k[i]=a[i+1].h,b[i]=dp[i];
		Modify(1,1,M,i);
	}
	write(dp[n]);
	return 0;
}

Code2

#include<bits/stdc++.h>
using namespace std;
template <typename T>
inline void read(T &x){
	x=0;bool f=false;char ch=getchar();
	while(!isdigit(ch)){f=(ch=='-');ch=getchar();}
	while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
	return ;
}
template <typename T>
inline void write(T x){
	if(x<0) x=-x,putchar('-');
	if(x>9) write(x/10);
	putchar(x%10^48);
	return ;
}
const int N=1e6+5,M=1e6;
#define int long long
struct node{
	int h,w;
	inline bool operator < (const node &B){return h==B.h?w>B.w:h>B.h;}
}a[N];
int n,m,top;
int tag[N<<2],k[N],b[N];
#define calc(i,x) (k[i]*x+b[i])
void Modify(int x,int l,int r,int d){
	if(l==r){if(calc(d,l)<calc(tag[x],l)) tag[x]=d;return ;}
	if(!tag[x]){tag[x]=d;return ;}
	int mid=l+r>>1;
	int Y1=calc(tag[x],mid),Y2=calc(d,mid);
	if(k[tag[x]]<k[d]){
		if(Y1<=Y2) Modify(x<<1,l,mid,d);
		else Modify(x<<1|1,mid+1,r,tag[x]),tag[x]=d;
	}
	else if(k[tag[x]]>k[d]){
		if(Y1<=Y2) Modify(x<<1|1,mid+1,r,d);
		else Modify(x<<1,l,mid,tag[x]),tag[x]=d;
	}
	else if(b[tag[x]]<b[d]) tag[x]=d;
	return ;
}
int Query(int x,int l,int r,int d){
	if(l==r) return calc(tag[x],d);
	int mid=l+r>>1;int res=calc(tag[x],d);
	if(d<=mid) res=min(res,Query(x<<1,l,mid,d));
	else res=min(res,Query(x<<1|1,mid+1,r,d));
	return res;
}
//dp[i]=dp[j]+(a[i].w*a[j+1].h)
int dp[N];
signed main(){
	read(n);
	for(int i=1;i<=n;i++) read(a[i].h),read(a[i].w);
	sort(a+1,a+n+1);top=1;
	for(int i=1;i<=n;i++) if(a[top].w<a[i].w) a[++top]=a[i];
	n=top;k[0]=b[0]=M;
	k[n+1]=a[1].h,b[n+1]=0;
	Modify(1,1,M,n+1);
	for(int i=1;i<=n;i++){
		dp[i]=Query(1,1,M,a[i].w);
		k[i]=a[i+1].h,b[i]=dp[i];
		Modify(1,1,M,i);
	}
	write(dp[n]);
	return 0;
}

标签：ch,return,int,tag,P2900,Land,inline,calc,Acquisition
来源： https://www.cnblogs.com/Akmaey/p/14677025.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。