1. 程式人生 > >hdu 6070 Dirt Ratio(分數規劃)

hdu 6070 Dirt Ratio(分數規劃)

分享 int check def space pen void hid ati

題目鏈接:hdu 6070 Dirt Ratio

題意:

給你n個數,讓你找一段區間[l,r],使得[l,r]中不同的數的個數size/(r-l+1)最小。

題解:

Claris官方題解:

技術分享

技術分享
 1 #include<bits/stdc++.h>
 2 #define F(i,a,b) for(int i=(a);i<=(b);++i)
 3 using namespace std;
 4 
 5 const int N=6e4+7;
 6 int t,n,a[N],lazy[N*4],la[N];
 7 double T[N*4],eps=1e-4;
 8 
 9 void build(double
v,int l=1,int r=n,int rt=1) 10 { 11 lazy[rt]=0; 12 if(l==r){T[rt]=l*v;return;} 13 int mid=l+r>>1; 14 build(v,l,mid,rt<<1),build(v,mid+1,r,rt<<1|1); 15 T[rt]=min(T[rt<<1],T[rt<<1|1]); 16 } 17 18 void PD(int rt) 19 { 20 lazy[rt<<1]+=lazy[rt],lazy[rt<<1
|1]+=lazy[rt]; 21 T[rt<<1]+=lazy[rt],T[rt<<1|1]+=lazy[rt],lazy[rt]=0; 22 } 23 24 void update(int L,int R,int v,int l=1,int r=n,int rt=1) 25 { 26 if(L<=l&&r<=R){lazy[rt]+=v,T[rt]+=v;return;} 27 int mid=l+r>>1; 28 if(lazy[rt])PD(rt); 29 if(L<=mid)update(L,R,v,l,mid,rt<<1
); 30 if(R>mid)update(L,R,v,mid+1,r,rt<<1|1); 31 T[rt]=min(T[rt<<1],T[rt<<1|1]); 32 } 33 34 double ask(int L,int R,int l=1,int r=n,int rt=1) 35 { 36 if(L<=l&&r<=R)return T[rt]; 37 int mid=l+r>>1; 38 if(lazy[rt])PD(rt);double ans=1e12; 39 if(L<=mid)ans=min(ans,ask(L,R,l,mid,rt<<1)); 40 if(R>mid)ans=min(ans,ask(L,R,mid+1,r,rt<<1|1)); 41 return ans; 42 } 43 44 int check(double mid) 45 { 46 build(mid),memset(la,0,sizeof(la)); 47 F(r,1,n) 48 { 49 update(la[a[r]]+1,r,1); 50 la[a[r]]=r; 51 if(ask(1,r)<=mid*(r+1))return 1; 52 } 53 return 0; 54 } 55 56 int main(){ 57 scanf("%d",&t); 58 while(t--) 59 { 60 scanf("%d",&n); 61 F(i,1,n)scanf("%d",a+i); 62 double l=0,r=1,mid; 63 F(ic,1,15) 64 { 65 if(r-l<eps)break; 66 mid=(l+r)/2; 67 if(check(mid))r=mid;else l=mid; 68 } 69 printf("%f\n",mid); 70 } 71 return 0; 72 }
View Code

hdu 6070 Dirt Ratio(分數規劃)