hdu6070Dirt Ratio 多校题 套路二分

比赛中我看了一眼题目就觉得是二分的套路，跟miaom说，结果发现miaom开始碎觉

miaom醒来以后表示这是道凸包合并（%%%）

我&wzf2000：那您快写啊

miaom：我不会写啊

莫名其妙的忘记了最初二分的思路

赛后看题解，果然是套路移动右端点

假设要mid能达到，那么一定有一组lr满足d(l,r)/(r-l+1)<=mid （d为区间不同数的种数）

d(l,r)+l*mid<=(r+1)*mid

在右端点右移1的时候只有一个后缀的d会改变，于是能轻松维护所有左端点对应的值

 #include <bits/stdc++.h>
 #define mi (l+r>>1)
 #define eps 0.00001
 #define INF 2000000000
 using namespace std;
 int n,T;
 int a[],b[];
 double tr[],flag[];
 void add(int now,int l,int r,int x,int y,double z)
 {
     if(l==x && r==y)
     {
         tr[now]+=z;
         flag[now]+=z;
         return;
     }
     if(x<=mi) add(now<<,l,mi,x,min(mi,y),z);
     if(y>mi) add(now<<|,mi+,r,max(mi+,x),y,z);
     tr[now]=min(tr[now<<],tr[now<<|])+flag[now];
 }
 double que(int now,int l,int r,int x,int y)
 {
     if(l==x && r==y) return tr[now];
     double ret=INF;
     if(x<=mi) ret=min(ret,que(now<<,l,mi,x,min(mi,y)));
     if(y>mi) ret=min(ret,que(now<<|,mi+,r,max(mi+,x),y));
     return ret+flag[now];
 }
 bool check(double mid)
 {
     for(int i=;i<=n;i++)
         b[i]=;
     for(int i=;i<=*n;i++)
         tr[i]=,flag[i]=;
     for(int i=;i<=n;i++)
     {
         add(,,n,i,i,mid*i);
         add(,,n,b[a[i]]+,i,);
         b[a[i]]=i;
         if(que(,,n,,i)<=mid*(i+))
             return ;
     }
     return ;
 }
 int main()
 {
     for(scanf("%d",&T);T;T--)
     {
         scanf("%d",&n);
         for(int i=;i<=n;i++)
             scanf("%d",&a[i]);
         double l=,r=;
         for(double mid=(l+r)/;r-l>eps;mid=(l+r)/)
         if(check(mid)) r=mid;else l=mid;
         printf("%.5f",l);
     }
     return ;
 }