题目链接:https://atcoder.jp/contests/abc127/tasks/abc127_f

题目大意

  初始状态下$f(x) = 0$,现在有 2 种模式的询问,第一种以“1 a b”的形式,需要进行操作$f(x) = f(x) + |x - a| + b$;第二种以“2”的形式,求使得 f(x) 取得最小值的 x 取值和 f(x) 值,如果有多个 x,输出任意一个即可。

分析

  考虑第一种询问已经出现了 k 次,现在遇到第二种询问。此时$f(x) = \sum_{i = 1}^k (|x - a_i| + b_i)$,可以发现$b_i$完全是独立于 x 的,所以只要讨论$\sum_{i = 1}^k |x - a_i|$的最小值即可。
  于是问题就退化成在数轴上找一个点,使得它到给定的 k 个点的距离之和最短。
  当点个数为奇数,最优解是中位数。 
  当点个数为偶数,最优解是位于中间的两个数中任取一个。
  所以这个问题的本质就是让你动态找中位数,维护中位数可以用双堆。

代码如下

  1. #include <bits/stdc++.h>
  2. using namespace std;
  3.  
  4. #define INIT() ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
  5. #define Rep(i,n) for (int i = 0; i < (n); ++i)
  6. #define For(i,s,t) for (int i = (s); i <= (t); ++i)
  7. #define rFor(i,t,s) for (int i = (t); i >= (s); --i)
  8. #define ForLL(i, s, t) for (LL i = LL(s); i <= LL(t); ++i)
  9. #define rForLL(i, t, s) for (LL i = LL(t); i >= LL(s); --i)
  10. #define foreach(i,c) for (__typeof(c.begin()) i = c.begin(); i != c.end(); ++i)
  11. #define rforeach(i,c) for (__typeof(c.rbegin()) i = c.rbegin(); i != c.rend(); ++i)
  12.  
  13. #define pr(x) cout << #x << " = " << x << " "
  14. #define prln(x) cout << #x << " = " << x << endl
  15.  
  16. #define LOWBIT(x) ((x)&(-x))
  17.  
  18. #define ALL(x) x.begin(),x.end()
  19. #define INS(x) inserter(x,x.begin())
  20.  
  21. #define ms0(a) memset(a,0,sizeof(a))
  22. #define msI(a) memset(a,inf,sizeof(a))
  23. #define msM(a) memset(a,-1,sizeof(a))
  24.  
  25. #define MP make_pair
  26. #define PB push_back
  27. #define ft first
  28. #define sd second
  29.  
  30. template<typename T1, typename T2>
  31. istream &operator>>(istream &in, pair<T1, T2> &p) {
  32. in >> p.first >> p.second;
  33. return in;
  34. }
  35.  
  36. template<typename T>
  37. istream &operator>>(istream &in, vector<T> &v) {
  38. for (auto &x: v)
  39. in >> x;
  40. return in;
  41. }
  42.  
  43. template<typename T1, typename T2>
  44. ostream &operator<<(ostream &out, const std::pair<T1, T2> &p) {
  45. out << "[" << p.first << ", " << p.second << "]" << "\n";
  46. return out;
  47. }
  48.  
  49. inline int gc(){
  50. static const int BUF = 1e7;
  51. static char buf[BUF], *bg = buf + BUF, *ed = bg;
  52.  
  53. if(bg == ed) fread(bg = buf, , BUF, stdin);
  54. return *bg++;
  55. }
  56.  
  57. inline int ri(){
  58. int x = , f = , c = gc();
  59. for(; c<||c>; f = c=='-'?-:f, c=gc());
  60. for(; c>&&c<; x = x* + c - , c=gc());
  61. return x*f;
  62. }
  63.  
  64. typedef long long LL;
  65. typedef unsigned long long uLL;
  66. typedef pair< double, double > PDD;
  67. typedef pair< int, int > PII;
  68. typedef pair< string, int > PSI;
  69. typedef set< int > SI;
  70. typedef vector< int > VI;
  71. typedef vector< PII > VPII;
  72. typedef map< int, int > MII;
  73. typedef pair< LL, LL > PLL;
  74. typedef vector< LL > VL;
  75. typedef vector< VL > VVL;
  76. typedef priority_queue< int > PQIMax;
  77. typedef priority_queue< int, VI, greater< int > > PQIMin;
  78. const double EPS = 1e-;
  79. const LL inf = 0x7fffffff;
  80. const LL infLL = 0x7fffffffffffffffLL;
  81. const LL mod = 1e9 + ;
  82. const int maxN = 2e5 + ;
  83. const LL ONE = ;
  84. const LL evenBits = 0xaaaaaaaaaaaaaaaa;
  85. const LL oddBits = 0x5555555555555555;
  86.  
  87. int Q;
  88. PQIMax maxH;
  89. PQIMin minH;
  90. LL sumB, sumLA, sumRA, ans;
  91.  
  92. int main(){
  93. INIT();
  94. cin >> Q;
  95. while(Q--) {
  96. int x, a, b;
  97. cin >> x;
  98. if(x == ) {
  99. cin >> a >> b;
  100. sumB += b;
  101. if(maxH.size() == minH.size()) {
  102. if(maxH.empty()) {
  103. maxH.push(a);
  104. sumLA += a;
  105. }
  106. else if(a >= maxH.top()) {
  107. minH.push(a);
  108. sumRA += a;
  109. }
  110. else {
  111. maxH.push(a);
  112. sumLA += a;
  113. }
  114. }
  115. else if(maxH.size() > minH.size()) {
  116. if(a >= maxH.top()) {
  117. minH.push(a);
  118. sumRA += a;
  119. }
  120. else {
  121. maxH.push(a);
  122. sumLA += a;
  123. sumLA -= maxH.top();
  124. minH.push(maxH.top());
  125. sumRA += maxH.top();
  126. maxH.pop();
  127. }
  128. }
  129. else {
  130. if(a < minH.top()) {
  131. maxH.push(a);
  132. sumLA += a;
  133. }
  134. else {
  135. minH.push(a);
  136. sumRA += a;
  137. sumRA -= minH.top();
  138. maxH.push(minH.top());
  139. sumLA += minH.top();
  140. minH.pop();
  141. }
  142. }
  143. }
  144. else {
  145. int x = maxH.top();
  146. ans = sumRA - sumLA + sumB;
  147. if(maxH.size() > minH.size()) ans += maxH.top();
  148. else if(maxH.size() < minH.size()) {
  149. ans -= minH.top();
  150. x = minH.top();
  151. }
  152. cout << x << " " << ans << endl;
  153. }
  154. }
  155. return ;
  156. }

AtCoder ABC 127F Absolute Minima的更多相关文章

  1. ATCODER ABC 099

    ATCODER ABC 099 记录一下自己第一场AK的比赛吧...虽然还是被各种踩... 只能说ABC确实是比较容易. A 题目大意 给你一个数(1~1999),让你判断它是不是大于999. Sol ...

  2. Atcoder ABC 141

    Atcoder ABC 141 A - Weather Prediction SB题啊,不讲. #include<iostream> #include<cstdio> #inc ...

  3. Atcoder ABC 139E

    Atcoder ABC 139E 题意: n支球队大循环赛,每支队伍一天只能打一场,求最少几天能打完. 解法: 考虑抽象图论模型,既然一天只能打一场,那么就把每一支球队和它需要交手的球队连边. 求出拓 ...

  4. Atcoder ABC 139D

    Atcoder ABC 139D 解法: 等差数列求和公式,记得开 $ long long $ CODE: #include<iostream> #include<cstdio> ...

  5. Atcoder ABC 139C

    Atcoder ABC 139C 题意: 有 $ n $ 个正方形,选择一个起始位置,使得从这个位置向右的小于等于这个正方形的高度的数量最多. 解法: 简单递推. CODE: #include< ...

  6. Atcoder ABC 139B

    Atcoder ABC 139B 题意: 一开始有1个插口,你的插排有 $ a $ 个插口,你需要 $ b $ 个插口,问你最少需要多少个插排. 解法: 暴力模拟. CODE: #include< ...

  7. Atcoder ABC 139A

    Atcoder ABC 139A 题意: 给你两个字符串,记录对应位置字符相同的个数 $ (n=3) $ 解法: 暴力枚举. CODE: #include<iostream> #inclu ...

  8. atcoder abc 244

    atcoder abc 244 D - swap hats 给定两个 R,G,B 的排列 进行刚好 \(10^{18}\) 次操作,每一次选择两个交换 问最后能否相同 刚好 \(10^{18}\) 次 ...

  9. AtCoder ABC 250 总结

    AtCoder ABC 250 总结 总体 连续若干次一样的结果:30min 切前 4 题,剩下卡在 T5 这几次卡在 T5 都是一次比一次接近, 什么 dp 前缀和打挂,精度被卡,能水过的题连水法都 ...

随机推荐

  1. Delphi动态添加控件

    {动态添加导航} var Panl:Tpanel; MainPage,Subpage:TPageControl; TabSheet1: TTabSheet; ToolBar2: TToolBar; S ...

  2. UVA - 10347 - Medians(由三中线求三角形面积)

    AC代码: #include<cstdio> #include<cmath> #include<algorithm> #include<iostream> ...

  3. Gerrit(1): Manage Projects

    1) Register an openid account https://login.ubuntu.com/+login 2) Custom settings set SSH pubkey set ...

  4. Centos 能ping通域名和公网ip但是网站不能够打开,服务器拒绝了请求。打开80端口解决。

    博客搬迁,给你带来的不便,敬请谅解! http://www.suanliutudousi.com/2017/10/29/centos-%E8%83%BDping%E9%80%9A%E5%9F%9F%E ...

  5. Jmeter+ant

    1.下载 ant,解压到非中文目录,并配置环境变量,不会的自行 google 2.将 jmeter 中 extras 子目录里的 ant-jmeter-1.1.1.jar 复制到 ant 中的 lib ...

  6. ArcGis相接面补节点c#

    相接(Touch)面执行切割后 新面与原相接面会缺少公共节点. private void AddPointToTouchesPolygon(IFeatureCursor newFeatureCurso ...

  7. D3.js绘制 颜色:RGB、HSL和插值 (V3版本)

    颜色和插值   计算机中的颜色,常用的标准有RGB和HSL.   RGB:色彩模式是通过对红(Red).绿(Green).蓝(Blue)三个颜色通道相互叠加来得到额各式各样的颜色.三个通道的值得范围都 ...

  8. 如何恢复误删的OneNote页面

    今天不小心把半个月的日记删掉了!(为了减少页面数量,每个月的日记会放在同一个页面上). 幸运的是OneNote有自动备份功能,喜极而泣. 操作方法来自微软支持 打开丢失了最近笔记的笔记本. 单击“文件 ...

  9. linux安装lolcat实现彩色文字输出信息

    [root@localhost ~]# mount /dev/sr0 /media/[root@localhost ~]# rpm -ivh epel-release-latest-7.noarch. ...

  10. Eureka注册中心是什么?

    Eureka注册中心是什么? Eureka注册中心是什么? Eureka是Netflix开发的服务发现组件,本身是一个基于REST的服务.Spring Cloud将它集成在其子项目spring-clo ...