HDU 6070 Dirt Ratio 二分+线段树

    题意是说给出一组数，在这一组数中找一个区间，使这个区间内AC题数/提交数的值最小。（假设每个区间内最后一次提交都是AC的）

    因为保证每个区间都是有题AC的，所以这里可以直接简化为区间内题目种数/区间长度，但在如何更快的处理题目种数的时候一直没能找到合适的处理办法，看了好多博客才能明白这个整体的思路= =（感觉让我自己想怎么都想不出要二分答案....）

    首先假设答案为k，区间的题目种数为size,这样就可以得到一个式子size/(r-l+1)=k，但是要直接处理左边的值没有什么好的办法，所以将式子转化一下，可以变成size=(r-l+1)*k，化简成这样还不够，因为这样还是必须要枚举每个区间，所以再进一步转化，改为size+l*k=(r+1)*k，这样的话在处理的时候可以只枚举r，在查询size的时候存入线段树种来进行更新与查询，因为放入了线段树，所以就可以将l*k看作是一个定值，在初始化的时候将线段树的sum直接初始化为l*k再进行计算就可以不用再管l*k的值了。

    在线段树更新的时候，每次读入数据的时候记录下他的位置，然后只需去更新他的上个位置+1到现在的位置权值+1就可以确定这个区间的种数了。

    但上面的还只是假设，k的值我们是不知道的，但我们知道答案k的区间一定是在[0,1]中的，所以这里可以采用二分的方法，将初始的式子写成size/(r-l+1)<=k，题目中说误差为10^-4，所以只要跑一个循环控制一下精度，然后对于k的值进行二分处理，不断缩小范围，直到条件成立就可以得出结果了。

    （PS：这个题好像用的是分数规划的思想，但从来没见过，貌似知道这个思想的都很快想到二分来做了= =改天再去看一下这个东西）

    （再PS：代码用的是mid写的k....结合上面的看可能有点乱，见谅见谅= =）

    下面AC代码：

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<climits>using namespace std;const double acc=0.00001;int a[100005];int pre[100005];double mid;int n;struct node{    int l;    int r;    double sum;    double mark;}s[300005];void pushup(int k){    s[k].sum=min(s[k<<1].sum,s[(k<<1)+1].sum);}void pushdown(int k){    if(s[k].mark)    {        s[k<<1].mark+=s[k].mark;        s[(k<<1)+1].mark+=s[k].mark;        s[k<<1].sum+=s[k].mark;        s[(k<<1)+1].sum+=s[k].mark;        s[k].mark=0;    }}void build(int l,int r,int k){    s[k].l=l;    s[k].r=r;    s[k].mark=0;    if(l==r)    {        s[k].mark=0;        s[k].sum=mid*l;        return ;    }    int mid1=(l+r)/2;    build(l,mid1,k<<1);    build(mid1+1,r,(k<<1)+1);    pushup(k);}double query(int l,int r,int k){    if(l<=s[k].l&&r>=s[k].r)    {        return s[k].sum;    }    pushdown(k);    int mid3=(s[k].l+s[k].r)/2;    double res=INT_MAX;    if(l<=mid3)        res=min(res,query(l,r,(k<<1)));    if(r>mid3)        res=min(res,query(l,r,(k<<1)+1));    pushup(k);    return res;}void update(int l,int r,int c,int k){    if(l<=s[k].l&&r>=s[k].r)    {        s[k].mark+=c;        s[k].sum+=c;        return ;    }    pushdown(k);    int mid2=(s[k].l+s[k].r)/2;    if(l<=mid2)        update(l,r,c,k<<1);    if(r>mid2)        update(l,r,c,(k<<1)+1);    pushup(k);}int check(){    int i;    double k;    memset(pre,0,sizeof(pre));    build(1,n,1);    for(i=1;i<=n;i++)    {        k=mid*(i+1.0);        update(pre[a[i]]+1,i,1,1);        if(query(1,i,1)<=k)            return 1;        pre[a[i]]=i;    }    return 0;}int main(){    int T;    int i,j;    double ans;    double l,r;    scanf("%d",&T);    while(T--)    {        scanf("%d",&n);        for(i=1;i<=n;i++)        {            scanf("%d",&a[i]);        }        l=0;        r=1;        ans=1;        while(r-l>=acc)        {            mid=(l+r)/2.0;            if(check())            {                ans=mid;                r=mid-acc;            }            else                l=mid+acc;        }        printf("%.10f\n",ans);    }    return 0;}