CF39C-Moon Craters【dp】

1|0正题

题目链接:https://www.luogu.com.cn/problem/CF39C


1|1题目大意

坐标轴上有n个圆,给出每个圆的位置ci和半径ri

要求选出最多的圆使得他们不相交,求方案。

1n2000


1|2解题思路

转换为选出最多的不交区间,然后离散一下,按照区间长度排序。

fi,j表示区间i内最右边的圆边界是r时的最多数量。

然后转移就可以了。

时间复杂度O(n2)


1|3code

#include<cstdio> #include<cstring> #include<algorithm> #include<vector> using namespace std; const int N=2100,inf=2147483647; int n,l[N],r[N],f[N][N<<1],fr[N][N<<1]; int cnt,p[N],b[N<<1]; vector<int> q[N<<1]; bool cmp(int x,int y) {return (long long)r[x]-l[x]<=(long long)r[y]-l[y];} void print(int x,int r){ if(r<l[x])return; if(fr[x][r]){ printf("%d ",fr[x][r]); print(fr[x][r],r); print(x,l[fr[x][r]]); } else print(x,r-1); } int main() { scanf("%d",&n); for(int i=1;i<=n;i++){ int c,R; scanf("%d%d",&c,&R); l[i]=c-R;r[i]=c+R;p[i]=i; b[++cnt]=c-R;b[++cnt]=c+R; } l[++n]=-inf;r[n]=inf;p[n]=n; b[++cnt]=-inf;b[++cnt]=inf; sort(b+1,b+1+cnt); cnt=unique(b+1,b+1+cnt)-b-1; sort(p+1,p+1+n,cmp); for(int xx=1;xx<=n;xx++){ int x=p[xx]; l[x]=lower_bound(b+1,b+1+cnt,l[x])-b; r[x]=lower_bound(b+1,b+1+cnt,r[x])-b; f[x][l[x]]=1; for(int i=l[x]+1;i<=r[x];i++){ f[x][i]=f[x][i-1]; for(int j=0;j<q[i].size();j++){ int y=q[i][j]; if(l[y]<l[x])continue; if(f[x][l[y]]+f[y][r[y]]>f[x][i]) f[x][i]=f[x][l[y]]+f[y][r[y]],fr[x][i]=y; } } q[r[x]].push_back(x); } printf("%d\n",f[n][r[n]]-1); print(n,r[n]); return 0; }

__EOF__

本文作者QuantAsk
本文链接https://www.cnblogs.com/QuantAsk/p/15181789.html
关于博主:退役OIer,GD划水选手
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   QuantAsk  阅读(73)  评论(0编辑  收藏  举报
编辑推荐:
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
· Linux系列:如何用heaptrack跟踪.NET程序的非托管内存泄露
阅读排行:
· TypeScript + Deepseek 打造卜卦网站:技术与玄学的结合
· 阿里巴巴 QwQ-32B真的超越了 DeepSeek R-1吗?
· 【译】Visual Studio 中新的强大生产力特性
· 张高兴的大模型开发实战:(一)使用 Selenium 进行网页爬虫
· 【设计模式】告别冗长if-else语句:使用策略模式优化代码结构
点击右上角即可分享
微信分享提示