【数学建模】【APIO2015】Palembang Bridges
Description
一条东西走向的穆西河将巴邻旁市一分为二,分割成了区域 A 和区域 B。
Input
输入的第一行包含两个正整数 K 和 N,分别表示桥的上限数量和居民的数量。
Output
输出仅为一行,包含一个整数,表示 D1+D2+⋯+DN 的最小值。
Sample Input
B 0 A 4
B 1 B 3
A 5 B 7
B 2 A 6
B 1 A 7
Sample Output
HINT
K=1或K=2
题解
/*自己搞出来的一道题,真是好感动*/
为了方便,在同一侧的直接处理掉。
K=1时,显然就是取中位数。
K=2时,考虑每一对(S,T),(绿蓝代表两座桥),那么肯定是要选两条红线短的那一边。不过这样还不够好考虑,我们不如让着两条直线伸长相等的距离,延伸到中点。
于是对于每一对(S,T),直接看Mid到两墙的距离就可以选择了。
进一步得出,如果我们把它们按中点排序,那么最后的局面肯定是,左部分都去桥1,右部分都去桥2,有一个分割点。
那么我们可以枚举分割点,分别计算两部分的最优桥,很神奇的把K=2分成了两个K=1,得到n^2的算法。
显然最优桥我们是可以动态维护的。
考虑点i(一对ST),把它从右部分加入左部分。
那么它对于最优桥的影响是可以讨论的,过程很傻逼可以自己试一试。
我得到的结果是,考虑左右部分肯定都是偶数,那么另左部分中位数取(len/2),右部分中位数取(len/2+1),这么做比较方便。
然后加点到左部分的时候,因为我们已经按中点排序,所以只可能是一左一右或者两右(左右是相对当前最优桥而言的)。
那么如果是一左一右,最优桥不发生变化,对ans新的贡献也就是这个点的贡献。
如果是两右,会使最优桥右移一位,但对于之前已加入的点是没有影响的,ans的改变还是只有这个点。
于是用两个树状数组模拟两边,每次求最优桥(这个也用树状数组求第k小数),然后维护ans也就是计算当前点的影响就行了。
复杂度O(nlogn)。说不清还是看代码吧。
主要考察数学建模、对数据结构的应用,感觉这题还是蛮好的。
代码
代码能力各种逗啊QwQ 调了好久 但调出来后真是好久没觉得这么爽了
#include<cstdio> #include<algorithm> #define ll long long using namespace std; const int maxn=1e5+5;int abs(int x){return x>0?x:-x;}struct point{int x,y,idx,idy,mid;bool operator <(const point&aa)const {return mid<aa.mid;} }a[maxn];int b[maxn*2],l,len,k,n; ll ans;void prepare(){for(int i=1;i<=l;i++)b[++len]=a[i].x,b[++len]=a[i].y;sort(b+1,b+len+1);for(int i=1;i<=l;i++){if(a[i].x>a[i].y) swap(a[i].x,a[i].y);a[i].idx=lower_bound(b+1,b+len+1,a[i].x)-b;a[i].idy=lower_bound(b+1,b+len+1,a[i].y)-b;} }int S1[maxn*2],S2[maxn*2];int lowbit(int o){return o&-o;} int add(int o,int k,int* C){while(o<=len){C[o]+=k;o+=lowbit(o);} }int find(int k,int *C){int ret=0,cnt=0;for(int i=17;i>=0;i--){ret+=(1<<i);if(ret>len||cnt+C[ret]>=k) ret-=(1<<i);else{cnt+=C[ret];}}return ret+1; }ll solve(){ll ans1=0,ans2=0,ansx=0;prepare();int A=0,B=b[len/2+1];for(int i=1;i<=len;i++)ans2+=abs(B-b[i]);for(int i=1;i<=l;i++)add(a[i].idx,1,S2),add(a[i].idy,1,S2);ansx=ans2;for(int i=1;i<=l;i++){int l1=i*2,l2=len-i*2;add(a[i].idx,1,S1);add(a[i].idy,1,S1);add(a[i].idx,-1,S2);add(a[i].idy,-1,S2);ans2-=abs(a[i].x-B)+abs(a[i].y-B);A=b[find(l1/2,S1)],B=b[find(l2/2+1,S2)];ans1+=abs(a[i].x-A)+abs(a[i].y-A);if(ans1+ans2<ansx) ansx=ans1+ans2;}return ansx; }int main(){char p,q;int u,v;scanf("%d%d",&k,&n);if(k==1){ans=n;for(int i=1;i<=n;i++){scanf("\n%c %d %c %d",&p,&u,&q,&v);if(p==q) ans+=abs(u-v),ans--;else b[++l]=u,b[++l]=v;}sort(b+1,b+l+1);int A=b[l/2];for(int i=1;i<=l;i++)ans+=abs(A-b[i]);printf("%lld\n",ans);}else{ans=n;for(int i=1;i<=n;i++){l++;scanf("\n%c %d %c %d",&p,&a[l].x,&q,&a[l].y);if(p==q) ans+=abs(a[l].x-a[l].y)-1,l--;else a[l].mid=a[l].x+a[l].y;}sort(a+1,a+l+1);printf("%lld",ans+solve());}return 0; }
转载于:https://www.cnblogs.com/xkui/p/4536421.html
【数学建模】【APIO2015】Palembang Bridges相关推荐
- 2021 年高教社杯全国大学生数学建模竞赛A题分析
2021 年高教社杯全国大学生数学建模竞赛A题分析 题目 赛题分析 前言 问题一分析 问题二分析 问题三分析 题目 A 题 "FAST"主动反射面的形状调节 中国天眼--500 米 ...
- 2021 年第十一届 MathorCup 高校数学建模挑战赛A题分析
2021 年第十一届 MathorCup 高校数学建模挑战赛A题分析 题目:自动驾驶中的车辆调头问题 问题分析 前言 综述 问题一分析 问题二分析 问题三分析 问题四分析 问题五分析 问题六分析 题目 ...
- 全国计算机建模三等奖,2009年全国数学建模真题(论文获国家 三等奖).doc
2009年全国数学建模真题(论文获国家 三等奖) 制动器试验台的控制方法分析 摘要 制动器的作用在于使行驶减速或者停止,其设计是车辆设计中最重要的环节之一,直接影响着人身和车辆的安全.本文分析了制动器 ...
- 《数学建模:基于R》一一2.2 方差分析
本节书摘来自华章计算机<数学建模:基于R>一书中的第2章,第2.2节,作者:薛 毅 更多章节内容可以访问云栖社区"华章计算机"公众号查看. 2.2 方差分析 方差分析是 ...
- 数学建模优化模型简单例题_数学建模之优化模型:存储模型
点击上方「蓝字」关注我们 最近,为申报市级精品课程,我为我校"数学建模与科学计算"课程录制了讲课视频,下面是3.1节优化模型的第一个例子:存储模型.敬请大家批评指正! 优化模型是数 ...
- 资料分享:数学建模资料分享 -- 神经网络部分
背景 周日的时候,为数学建模俱乐部的同学们进行了一场有关人工神经网络方面的分享.虽然在这个方面有一些积累,但过于零散,所以拿了一堆文件拼凑成整体的内容. 幸亏自己有分享的习惯,学会一些知识就写下来,这 ...
- 资料分享:数学建模资料分享 -- 图论部分
背景 今天上午,在教六第一阶梯教室为数学建模俱乐部的同学们分享了有关图论的基本知识和应用. 课后,为同学们留了一个算法实现的小练习,大家可以先做一下.在本图文的末尾处,我把上课的资料以及代码分享出来, ...
- 如何准备数学建模竞赛!
昨天早晨,我到教十一实验室的时候遇到史会峰老师,他说正准备给学生们进行数学建模的培训.今天早晨,我又遇到了孔令才老师,他同样也说准备给学生们进行数学建模的培训.看到这么多同事在做这个事情,想想自己也应 ...
- 全国大学生数学建模竞赛中,哈工大被禁用MATLAB
整理 | 屠敏 出品 | AI科技大本营(ID:rgznai100) AK47 VS "小米加步枪"同台竞技,最终会有什么样的结果? 目前我们尚未可知,但是这样的"竞赛& ...
最新文章
- 【学术相关】Github大盘点!2021年最值得推荐的38篇AI论文
- 分享一个HTML5画布实现的超酷文字弹跳球效果
- CORS(跨域资源共享)
- (一)带有图像到图像转换的移动风格迁移
- RN和React路由详解及对比
- sql as关键字_SQL AS关键字概述和示例
- Ubuntu Server 安装桌面untiy
- 《一秒学会C++》异步回调函数(C++11)
- mysql group by rowid_MySQL里面的group by问题浅析
- py+webdriver学习之旅6:cookie
- 创建一个带目录的Word模板
- PyQt5-使用QDialog显示通用对话框
- cad插入块_【CAD教程】创建块和块应用步骤
- 逻辑、认识论和本体论“三统一”
- 小马哥---高仿苹果6 R7226刷机拆机主板图与开机识别图 修改WiFi信号增强
- 微信扫码支付:此公众号并没有这些scope的权限,错误码:10005
- IIS服务器安全配置[摘]
- 好记性不如烂笔头(2)
- 【渗透测试自学系列】—黑客是如何利用常用的搜索引擎的?
- 如何区分实用新型专利和发明专利
热门文章
- mvc在页面上显示PDF
- 31全志r58平台Android4.4.2下打开USB摄像头
- bs4 CSS选择器
- 鸟哥Linux私房菜知识点总结3到5章
- 引用与指针的异同-基础篇
- hdu 1116 欧拉回路 并查集 一组字符串能否首尾相连成一个字符串
- 【Notes】《Thinking in Java》【Chapter 11】Part II
- UVA11825 黑客的攻击 Hackers' Crackdown 状压DP,二进制,子集枚举
- [转]bootstrap table 动态列数
- 利用 :placeholder-shown 选择器实现 label 浮动效果