HDU 6070 (线段树)(统计颜色)
HDU 6070 Partition
Problem :
给一段长度为n的序列,要求找出一段区间,使得这段区间的数字种类除以区间长度最小。输出最后的答案即可。(n <= 60000)(9s时限)
Solution :
显然,答案是0~1中的一个数字,可以很自然的想到二分答案的做法。假设目前二分到的答案为mid,那么需要判断
\]
其中cnt(l,r)为l到r这个区间内的数字种类。变化一下式子可以得到:
\]
通过枚举有端点r,使用线段树维护左边的式子,每当右端点r向右移动1时,所影响的区间为r到对应颜色上一次出现的位置,区间整体加1就行了。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <string>
#include <cmath>
#include <algorithm>
#include <vector>
#include <map>
#include <queue>
using namespace std;
#define eps 1e-10
const int N = 1e5 + 8;
int a[N], pre[N];
int n;
struct Segment_Tree
{
double tag[N << 2];
double lazy[N << 2];
void pushup(int rt)
{
int l = rt << 1, r = rt << 1 | 1;
tag[rt] = min(tag[l], tag[r]);
}
void pushdown(int rt)
{
int l = rt << 1, r = rt << 1 | 1;
if (lazy[rt])
{
tag[l] += lazy[rt];
tag[r] += lazy[rt];
lazy[l] += lazy[rt];
lazy[r] += lazy[rt];
lazy[rt] = 0;
}
}
void build(int l, int r, int rt, double x)
{
tag[rt] = lazy[rt] = 0;
if (l == r)
{
tag[rt] = l * x;
return;
}
int m = l + r >> 1;
build(l, m, rt << 1, x);
build(m + 1, r, rt << 1 | 1, x);
pushup(rt);
}
void update(int L, int R, int val, int l, int r, int rt)
{
if (L <= l && r <= R)
{
tag[rt] = tag[rt] + val;
lazy[rt] += val;
return;
}
pushdown(rt);
int m = l + r >> 1;
if (L <= m) update(L, R, val, l, m, rt << 1);
if (m < R) update(L, R, val, m + 1, r, rt << 1 | 1);
pushup(rt);
}
double query(int L, int R, int l, int r, int rt)
{
if (L <= l && r <= R)
{
return tag[rt];
}
pushdown(rt);
int m = l + r >> 1;
double ans = 1e12;
if (L <= m) ans = min(ans, query(L, R, l, m, rt << 1));
if (m < R) ans = min(ans, query(L, R, m + 1, r, rt << 1 | 1));
return ans;
}
}T;
void init()
{
cin >> n;
for (int i = 1; i <= n; ++i) cin >> a[i];
}
int sgn(double x)
{
if (fabs(x) < eps) return 0;
if (x > 0) return 1;
return -1;
}
bool check(double mid)
{
for (int i = 1; i <= n; ++i) pre[i] = 0;
T.build(1, n, 1, mid);
for (int i = 1; i <= n; ++i)
{
T.update(pre[a[i]] + 1, i, 1, 1, n, 1);
pre[a[i]] = i;
if (sgn(mid * (i + 1) - T.query(1, i, 1, n, 1) >= 0)) return 1;
}
return 0;
}
void solve()
{
double l = 0, r = 1;
while (l + eps < r)
{
double mid = (l + r) / 2;
if (check(mid)) r = mid; else l = mid;
}
printf("%.6f\n", l);
}
int main()
{
cin.sync_with_stdio(0);
int T; cin >> T;
for (int cas = 1; cas <= T; ++cas)
{
init();
solve();
}
}
HDU 6070 (线段树)(统计颜色)的更多相关文章
- HDU 6070 线段树
题意:求AC率,x/y 的最小值,x是区间数字的种类数,y是区间的长度. 分析: 二分答案比率.ans, 动态插入结点,一些区间的size会发生变化,是那些前面暂时没有新的结点的区间 size + 1 ...
- HDU 6035 (虚树)(统计颜色)
HDU 6035 Colorful Tree Problem : 给一棵树,每个结点有一种颜色,定义每条路径的权值为这条路径上颜色的种数,询问所有路径(C(n,2)条)的权值之和. Solution ...
- hdu 5877 线段树(2016 ACM/ICPC Asia Regional Dalian Online)
Weak Pair Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Total ...
- hdu 3974 线段树 将树弄到区间上
Assign the task Time Limit: 15000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) ...
- hdu 3436 线段树 一顿操作
Queue-jumpers Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) To ...
- hdu 3397 线段树双标记
Sequence operation Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Othe ...
- hdu 4578 线段树(标记处理)
Transformation Time Limit: 15000/8000 MS (Java/Others) Memory Limit: 65535/65536 K (Java/Others) ...
- hdu 4533 线段树(问题转化+)
威威猫系列故事——晒被子 Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others) Tot ...
- hdu 2871 线段树(各种操作)
Memory Control Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) T ...
随机推荐
- CF385C Bear and Prime Numbers
思路: 需要对埃氏筛法的时间复杂度有正确的认识(O(nlog(log(n)))),我都以为肯定超时了,结果能过. 实现: #include <bits/stdc++.h> using na ...
- OCP 11g 第一章练习
练习 1-1 研究所在环境的DBMS 这是一个书面练习,没有具体的解决方案. 确定自己所在环境使用的应用程序, 应用服务器 , 和数据库. 然后集中精力研究数据库, 体验一下数据库的规模和忙碌程度. ...
- Oracle Storage in Action : 删除物理数据文件
ALTER TABLESPACE XH_DM DROP DATAFILE 'F:\DEV\ORACLE-DATA\ORCL\XH_DM_1.DBF'; SQL> ALTER TABLESPACE ...
- Android-Emulator使用
1.查看当前android支持的avd版本 2.创建Emulator avd: android create avd -n magicyu -t 2 -n后面接需要创建avd的名字,-t后面接需要创建 ...
- ie 导出不行,不兼容问题,或只出现后缀文件无法识别
// 下载模板 @RequestMapping("/download") @ResponseBody public ResponseEntity<byte[]> dow ...
- 洛谷 P1006 传纸条
题目描述 小渊和小轩是好朋友也是同班同学,他们在一起总有谈不完的话题.一次素质拓展活动中,班上同学安排做成一个m行n列的矩阵,而小渊和小轩被安排在矩阵对角线的两端,因此,他们就无法直接交谈了.幸运的是 ...
- example - 在这里插入一句话的简介
总览 (SYNOPSIS) example [options] arguments 描述 (DESCRIPTION) 在这里插入描述 man9 应当是 “内核文档” 但是由于内核文档一般不以 man ...
- dockerfile 的最佳实践
Dockerfile 编写nginx容器 [root@mast nginx]# cat Dockerfile FROM centos MAINTAINER zhaoruidong RUN yum -y ...
- pooling需要注意的一个地方
max pooling 在不同的 depth 上是分开执行的,且不需要参数控制.也就是说,pooling之后,feature map的维度不会改变
- QT +菜单栏和工具栏
#include "mainwindow.h" #include <QMenuBar>//菜单栏需要的头文件 #include <QMenu>//菜单 #i ...