【题解】Cover/绝世好题/挂饰

#4620. Cover

题目描述
你在一个坐标系中获得N个点,它们需要一个或多个矩形进行覆盖,以满足下面这些条件:

·每个矩形的边与坐标轴平行

·每个矩形的中心位于原点,即点(0,0)

·每个给定的点位于矩形的内部或边界上。

当然,你可以使用一个矩形覆盖所有的点,但这个矩形的面积可能非常大。我们的目标是找到所需的矩形,使这些矩形的面积总和最小。

题解:本题只需要对问题进行转化即可

  • 每个点都可以等效成一象限的点,因为每个矩形对每个象限是等效的
  • 考虑点a,b,若a.x<b.x&&a.y<b.y,则可以删除点a
  • 将所有点按x大小排序,再处理完无用的点,设 d p [ j ] dp[j] dp[j]表示前 j j j个点的花费
  • 方程: d p [ i ] = m a x ( b [ i ] . x ∗ b [ j ] . y + d p [ j − 1 ] ) , 1 < = j < = i dp[i]=max(b[i].x*b[j].y+dp[j-1]),1<=j<=i dp[i]=max(b[i].xb[j].y+dp[j1]),1<=j<=i
#include<cstdio> #include<algorithm> #include<cstring> #include<iostream> #define ll long long using namespace std; const int maxn=5005; struct node{ ll x,y; }a[maxn],b[maxn]; ll n,cnt,dp[maxn]; bool cmp(node p,node q) { return p.x<q.x; } int main() { scanf("%lld",&n); for(int i=1;i<=n;i++) { scanf("%lld%lld",&a[i].x,&a[i].y); if(a[i].x<0) a[i].x*=-1; if(a[i].y<0) a[i].y*=-1; } sort(a+1,a+1+n,cmp); for(int i=n;i>=1;i--) { if(a[i].y>b[cnt].y) b[++cnt].x=a[i].x,b[cnt].y=a[i].y; } reverse(b+1,b+1+cnt); for(int i=1;i<=cnt;i++) { for(int j=1;j<=i;j++) { ll t=b[i].x*b[j].y+dp[j-1]; if(dp[i]==0||t<dp[i]) dp[i]=t; } } printf("%lld",dp[cnt]*4); }

「JOISC 2014 Day4」挂饰

代码略显复杂

#include<cstdio> #include<algorithm> #include<cstring> using namespace std; const int maxn=2000*2000; struct node{ int v,w; }a[2005]; int n,ans,dp[maxn+5],sum=1; //dp[j]表示至少有j个挂坠的最大利润(j是背包容量) bool cmp(node x,node y) { return x.v>y.v; } int main() { for(int i=2;i<=maxn;i++) dp[i]=-0x3f3f3f3f; dp[0]=dp[1]=0;//本身有一个挂坠(记录可行性) scanf("%d",&n); for(int i=1;i<=n;i++) { scanf("%d%d",&a[i].v,&a[i].w); a[i].v--; if(a[i].v>0) sum+=a[i].v; } sort(a+1,a+1+n,cmp); for(int i=1;i<=n;i++) {//用可行状态去刷表(注意顺序) if(a[i].v>=0) { for(int j=sum;j>=0;j--) { dp[j+a[i].v]=max(dp[j+a[i].v],dp[j]+a[i].w); } } else { for(int j=-a[i].v;j<=sum;j++) { dp[j+a[i].v]=max(dp[j+a[i].v],dp[j]+a[i].w); } } } for(int i=0;i<=sum;i++) { ans=max(ans,dp[i]); } printf("%d",ans); }

「BZOJ4300」绝世好题

#include<cstdio> #include<cstring> #include<algorithm> #include<iostream> using namespace std; int dp[32],n,b,ans; //令dp[i]表示数列到目前为止**最后一项第i位为1**的最大子序列长度 int main() { scanf("%d",&n); for(int i=1;i<=n;i++) { scanf("%d",&b); int k=1;//max(dp[j]+1),(1<<j)&b==1 for(int j=0;j<=30;j++) { if((1<<j)&b) k=max(dp[j]+1,k); } for(int j=0;j<=30;j++) { if((1<<j)&b) dp[j]=max(dp[j],k); } ans=max(ans,k); } printf("%d",ans); }

__EOF__

本文作者仰望星空的蚂蚁
本文链接https://www.cnblogs.com/cqbzly/p/17530419.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   仰望星空的蚂蚁  阅读(3)  评论(0编辑  收藏  举报  
编辑推荐:
· 记一次.NET内存居高不下排查解决与启示
· 探究高空视频全景AR技术的实现原理
· 理解Rust引用及其生命周期标识(上)
· 浏览器原生「磁吸」效果!Anchor Positioning 锚点定位神器解析
· 没有源码,如何修改代码逻辑?
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」
点击右上角即可分享
微信分享提示