[关闭]
@ysner 2018-05-15T21:31:30.000000Z 字数 1824 阅读 2060

[CTSC2018]混合果汁

主席树 二分


题面

戳我

解析

热烈庆祝中国最难OI赛事CTSC au线达到NOI水平
显然,我们可以二分出美味值的下界。
然后发现各个价格的果汁好像并没有什么区别。
那么贪心尽可能取其小者就好。
求前小?
主席树轻松解决。

  1. #include<iostream>
  2. #include<cstdio>
  3. #include<cstdlib>
  4. #include<cstring>
  5. #include<cmath>
  6. #include<algorithm>
  7. #define ll long long
  8. #define re register
  9. #define il inline
  10. #define ls x<<1
  11. #define rs x<<1|1
  12. #define max(a,b) (a>b?a:b)
  13. #define fp(i,a,b) for(re int i=a;i<=b;i++)
  14. #define fq(i,a,b) for(re int i=a;i>=b;i--)
  15. using namespace std;
  16. const int N=1e5+100;
  17. ll n,m,tot,mxd,mxp,p,v,rt[N*100];
  18. struct fru
  19. {
  20. ll d,p,v;
  21. bool operator < (const fru &o) const {return d<o.d;}
  22. }a[N];
  23. struct chi
  24. {
  25. ll p,v;
  26. }b[N];
  27. struct tree
  28. {
  29. ll l,r,sum,num;
  30. }t[N*100];
  31. il ll gi()
  32. {
  33. re ll x=0,t=1;
  34. re char ch=getchar();
  35. while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
  36. if(ch=='-') t=-1,ch=getchar();
  37. while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
  38. return x*t;
  39. }
  40. il void Modify(re ll &x,re ll y,re ll l,re ll r,re ll p,re ll v)
  41. {
  42. t[x=++tot]=t[y];
  43. t[x].num+=v;
  44. t[x].sum+=v*p;
  45. re ll mid=l+r>>1;
  46. if(l==r) return;
  47. if(p<=mid) Modify(t[x].l,t[y].l,l,mid,p,v);
  48. else Modify(t[x].r,t[y].r,mid+1,r,p,v);
  49. }
  50. il int Query(re ll x,re ll y,re ll l,re ll r,re ll p,re ll v)
  51. {
  52. re ll num=t[t[x].l].num-t[t[y].l].num,sum=t[t[x].l].sum-t[t[y].l].sum;
  53. if(l==r)
  54. {
  55. re ll l1=t[x].num-t[y].num,l2=t[x].sum-t[y].sum;
  56. if(l1<v) return 0;
  57. return p>=(l2/l1)*v;
  58. }
  59. re ll mid=l+r>>1;
  60. if(num>=v)
  61. {
  62. if(p>=sum) return 1;
  63. return Query(t[x].l,t[y].l,l,mid,p,v);
  64. }
  65. else
  66. {
  67. if(p<sum) return 0;
  68. return Query(t[x].r,t[y].r,mid+1,r,p-sum,v-num);
  69. }
  70. }
  71. il bool check(re ll x)
  72. {
  73. re ll l=1,r=n;
  74. while(l<r)
  75. {
  76. re ll mid=l+r>>1;
  77. if(a[mid].d>=x) r=mid;
  78. else l=mid+1;
  79. }
  80. return Query(rt[n],rt[l-1],1,mxp,p,v);
  81. }
  82. il ll solve()
  83. {
  84. re ll l=1,r=mxd;
  85. if(!check(1)) return -1;
  86. while(l<r)
  87. {
  88. re ll mid=l+r+1>>1;
  89. if(check(mid)) l=mid;
  90. else r=mid-1;
  91. }
  92. return l;//我们取的是最小值,所以要用l
  93. }
  94. int main()
  95. {
  96. n=gi();m=gi();
  97. fp(i,1,n) a[i].d=gi(),a[i].p=gi(),a[i].v=gi(),mxd=max(mxd,a[i].d),mxp=max(mxp,a[i].p);
  98. sort(a+1,a+1+n);
  99. fp(i,1,n) Modify(rt[i],rt[i-1],1,mxp,a[i].p,a[i].v);
  100. fp(i,1,m)
  101. {
  102. p=gi(),v=gi();
  103. printf("%lld\n",solve());
  104. }
  105. return 0;
  106. }
添加新批注
在作者公开此批注前,只有你和作者可见。
回复批注