hdu6070 Dirt Ratio 二分+线段树

/**
题目：hdu6070 Dirt Ratio
链接：http://acm.hdu.edu.cn/showproblem.php?pid=6070
题意：给定n个数，求1.0*x/y最小是多少。x表示一段区间内不同数字的个数，y表示区间长度。
思路：二分+线段树
二分答案x/y。 找一段区间满足 size(l,r)/(r-l+1) <= mid  , size(l,r)表示[l,r]内不同数的个数。

size(l,r)<=mid(r-l+1) => size(l,r)+mid*l<=mid*(r+1);

用线段树维护size(l,r)+mid*l的最小值。每一个节点存储size(l,r)+mid*l,枚举r，更新当前a[r]到不包含该数字的节点值+1

*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
#define lson L,m,rt<<1
#define rson m+1,R,rt<<1|1
const double eps = 1e-;
const int N = 6e4+;
double add[N<<], sum[N<<];
int last[N], a[N];
void pushup(int rt)
{
    sum[rt] = min(sum[rt<<],sum[rt<<|]);
}
void pushdown(int rt)
{
    sum[rt<<] += add[rt];
    sum[rt<<|] += add[rt];
    add[rt<<] += add[rt];
    add[rt<<|] += add[rt];
    add[rt] = ;
}
void update(int l,int r,int L,int R,int rt,double d)
{
    if(l<=L&&R<=r){
        sum[rt]+=d;
        add[rt]+=d;
        return ;
    }
    if(add[rt]>eps) pushdown(rt);
    int m = (L+R)/;
    if(r<=m) update(l,r,lson,d);
    else if(l>m) update(l,r,rson,d);
    else{
        update(l,r,lson,d);
        update(l,r,rson,d);
    }
    pushup(rt);
}
double query(int l,int r,int L,int R,int rt)
{
    if(l<=L&&R<=r){
        return sum[rt];
    }
    if(add[rt]>eps) pushdown(rt);
    int m = (L+R)/;
    double mis;
    if(r<=m) mis = query(l,r,lson);
    else if(l>m) mis = query(l,r,rson);
    else mis = min(query(l,r,lson),query(l,r,rson));
    pushup(rt);
    return mis;
}
int main()
{
    //freopen("C:\\Users\\accqx\\Desktop\\in.txt","r",stdin);
    int T, n;
    cin>>T;
    while(T--)
    {
        scanf("%d",&n);
        for(int i = ; i <= n; i++) scanf("%d",&a[i]);
        double hi = , lo = , mid;
        double temp;
        while(hi-lo>eps){
            mid = (hi+lo)/;
            memset(add, , sizeof add);
            memset(sum, , sizeof sum);
            memset(last, , sizeof last);
            int flag = ;
            for(int i = ; i <= n; i++){
                update(last[a[i]]+,i,,n,,);
                update(i,i,,n,,mid*i);
                temp = query(,i,,n,);
                if(temp<=mid*(i+)){
                    flag = ; break;
                }
                last[a[i]] = i;
            }
            if(flag) hi = mid;
            else lo = mid;
        }
        //cout<<mid<<endl;
        printf("%.10f\n",mid);
    }
    return ;
}