[Swust OJ 1094]--中位数(巧用set,堆排序)
题目链接:http://acm.swust.edu.cn/problem/1094/
中位数(又称中值,英语:Median),统计学中的专有名词,代表一个样本、种群或概率分布中的一个数值,其可将数值集合划分为相等的上下两部分。对于有限的数集,可以通过把所有观察值高低排序后找出正中间的一个作为中位数。如果观察值有偶数个,则中位数不唯一,通常取最中间的两个数值的平均数作为中位数。
多组输入
第一行:一个正整数N (0<N<1000000)
第二行:N个正整数。(0=<A[i]<2^30)
每组数据先输出”Case X:”,X表示测试数据的编号,从1开始。
第二行输出N个数,第i个数对应数组前i个值的中位数。(精确到小数点后一位)
1
2
3
4
5
|
5
1 2 3 4 5
6
2 5 4 8 7 4
|
1
2
3
4
5
|
Case 1:
1.0 1.5 2.0 2.5 3.0
Case 2:
2.0 3.5 4.0 4.5 5.0 4.5
|
#include <iostream>
#include <set>
#include <cstdio>
using namespace std; #define rep(i,a,b) for(int i=a;i<=b;i++)
int t, n;
multiset<double>mpt; void mid_num(){
double x;
printf("Case %d:\r\n", ++t);
mpt.clear();
scanf("%lf", &x);
mpt.insert(x);
multiset <double> ::iterator it = mpt.begin();
printf("%.1lf", x);
rep(i, , n){
scanf("%lf", &x);
mpt.insert(x);
if (x < *it){
if (i & ) printf(" %.1lf", *(--it));
else {
multiset <double> ::iterator it1 = it;
printf(" %.1lf", (*it + *(--it1)) / 2.0);
}
}
else {
if (i & ) printf(" %.1lf", *it);
else {
++it;
multiset <double> ::iterator it1 = it;
printf(" %.1lf", (*it + *(--it1)) / 2.0); }
}
}
printf("\r\n");
} int main(){
while (~scanf("%d", &n))
mid_num();
return ;
}
堆排序的代码也贴出来吧(原谅我太懒没有优化这代码~~~)
#include <stdio.h>
#include <string.h>
int lstack[], ltot, rstack[], rtot, mid;
int Max(int a, int b){
return a > b ? a : b;
}
int Min(int a, int b){
return a<b ? a : b;
}
void lup(int step){
while (step != ){
if (lstack[step]>lstack[step / ])lstack[step] ^= lstack[step / ] ^= lstack[step] ^= lstack[step / ];
else break;
step = step / ;
}
if (step == && lstack[] > mid) lstack[] ^= mid ^= lstack[] ^= mid;
}
void ldown(){
if (ltot > && mid < lstack[]) mid ^= lstack[] ^= mid ^= lstack[];
else return;
int step = ;
while (step * < ltot){
if (step * + >= ltot){
if (lstack[step] < lstack[step * ])
lstack[step] ^= lstack[step * ] ^= lstack[step] ^= lstack[step * ], step = step * ;
else return;
}
else{
if (lstack[step] <= lstack[step * ] && lstack[step * + ] <= lstack[step * ])
lstack[step] ^= lstack[step * ] ^= lstack[step] ^= lstack[step * ], step = step * ;
else if (lstack[step] <= lstack[step * + ] && lstack[step * ] <= lstack[step * + ])
lstack[step] ^= lstack[step * + ] ^= lstack[step] ^= lstack[step * + ], step = step * + ;
else return;
}
}
}
void rup(int step){
while (step != ){
if (rstack[step]<rstack[step / ])rstack[step] ^= rstack[step / ] ^= rstack[step] ^= rstack[step / ];
else break;
step = step / ;
}
if (step == && rstack[]<mid) rstack[] ^= mid ^= rstack[] ^= mid;
}
void rdown(){
if (rtot> && mid>rstack[]) mid ^= rstack[] ^= mid ^= rstack[];
else return;
int step = ;
while (step * < rtot){
if (step * + >= rtot){
if (rstack[step] > rstack[step * ])
rstack[step] ^= rstack[step * ] ^= rstack[step] ^= rstack[step * ], step = step * ;
else return;
}
else{
if (rstack[step] >= rstack[step * ] && rstack[step * + ] >= rstack[step * ])
rstack[step] ^= rstack[step * ] ^= rstack[step] ^= rstack[step * ], step = step * ;
else if (rstack[step] >= rstack[step * + ] && rstack[step * ] >= rstack[step * + ])
rstack[step] ^= rstack[step * + ] ^= rstack[step] ^= rstack[step * + ], step = step * + ;
else return;
}
}
}
int main()
{
int t = , n, i;
while (scanf("%d", &n) != EOF){
printf("Case %d:\r\n", t++);
scanf("%d", &mid);
ltot = rtot = ;
printf("%.1lf", (double)mid);
for (i = ; i < n; i++){
if (i % ){
scanf("%d", &lstack[ltot]);
lup(ltot);
rdown();
ltot++;
printf(" %.1lf", (mid + lstack[]) / 2.0);
}
else{
scanf("%d", &rstack[rtot]);
rup(rtot);
ldown();
rtot++;
printf(" %.1lf", (double)mid);
}
}
printf("\r\n");
}
return ;
}
[Swust OJ 1094]--中位数(巧用set,堆排序)的更多相关文章
- [Swust OJ 404]--最小代价树(动态规划)
题目链接:http://acm.swust.edu.cn/problem/code/745255/ Time limit(ms): 1000 Memory limit(kb): 65535 Des ...
- [Swust OJ 649]--NBA Finals(dp,后台略(hen)坑)
题目链接:http://acm.swust.edu.cn/problem/649/ Time limit(ms): 1000 Memory limit(kb): 65535 Consider two ...
- SWUST OJ NBA Finals(0649)
NBA Finals(0649) Time limit(ms): 1000 Memory limit(kb): 65535 Submission: 404 Accepted: 128 Descri ...
- [Swust OJ 746]--点在线上(线段树解法及巧解)
题目链接:http://acm.swust.edu.cn/problem/746/ Time limit(ms): 1000 Memory limit(kb): 65535 fate是一个数学大牛 ...
- [Swust OJ 1026]--Egg pain's hzf
题目链接:http://acm.swust.edu.cn/problem/1026/ Time limit(ms): 3000 Memory limit(kb): 65535 hzf ...
- [swustoj 1094] 中位数
中位数(1094) 问题描述 中位数(又称中值,英语:Median),统计学中的专有名词,代表一个样本.种群或概率分布中的一个数值,其可将数值集合划分为相等的上下两部分.对于有限的数集,可以通过把所有 ...
- 九度OJ 1371 最小的K个数 -- 堆排序
题目地址:http://ac.jobdu.com/problem.php?pid=1371 题目描述: 输入n个整数,找出其中最小的K个数.例如输入4,5,1,6,2,7,3,8这8个数字,则最小的4 ...
- [Swust OJ 1023]--Escape(带点其他状态的BFS)
解题思路:http://acm.swust.edu.cn/problem/1023/ Time limit(ms): 5000 Memory limit(kb): 65535 Descript ...
- [Swust OJ 1125]--又见GCD(数论,素数表存贮因子)
题目链接:http://acm.swust.edu.cn/problem/1125/ Time limit(ms): 1000 Memory limit(kb): 65535 Descriptio ...
随机推荐
- 经典阅读-《Effective C++》Item1:视C++为一个联邦语言
C++已经是个多重范型编程语言(multiparadigm programming language),一个同时支持过程形式(procedural).面向对象形式(object-oriented).泛 ...
- Android 学习开发笔记《Service 与 Thread 的区别 》
很多时候,你可能会问,为什么要用 Service,而不用 Thread 呢,因为用 Thread 是很方便的,比起 Service 也方便多了,下面我详细的来解释一下. 1). Thread:Thre ...
- 创建android 模拟器并在cmd中打开
因为在运行monkeyrunner之前必须先运行相应的模拟器或连接真机,否则monkeyrunner无法连接到设备,运行模拟器有两种方法:1.通过eclipse中执行模拟器 2.在CMD中通过命令调用 ...
- [转]swift 学习资源 大集合
今天看到了一个swift的学习网站,里面收集了很多学习资源 [转自http://blog.csdn.net/sqc3375177/article/details/29206779] Swift 介绍 ...
- HTML5 appcache
参考http://www.zation.me/2013/05/28/build_offline_mobile_web_app.html 他的事件总结的比较好 checking:客户端正在检查manif ...
- 帝国cms语句调用
帝国cms系统,灵动标签,有着非常强大的数据调用功能.这里为广大菜鸟站长普及一下. 我们来看这段代码. [e:loop={,,}] <li>·<a target="_bla ...
- QT 自动获取可用串口
本来想直接用Settings来获取的,但是串口信息类似 "\Device\Serial0",死活获取不了,用了转义.反斜杠还是获取不到,所以就放弃了,网上好像也没有获取成功的. 所 ...
- lightoj Again Array Queries
1100 - Again Array Queries PDF (English) Statistics Forum Time Limit: 3 second(s) Memory Limit: 32 ...
- JavaBean讲解 规范
一:JavaBean定义 JavaBean是一种可重复使用.跨平台的软件组件.JavaBean可分为两种:一种是有用户界面(UI,User Interface)的JavaBean,例如NetBean中 ...
- Java对象序列化与反序列化一 JSON
Java对象序列化与反序列化一 JSON 1. 依赖库 jackson-all-1.6.1.jar 2. 代码 public class Student { private String nam ...