当前位置:   article > 正文

WQS二分—简介及基本应用

wqs二分

WQS二分是啥?搞了这么久我连它的全名是什么都一头雾水,姑且称其为忘情水二分吧。
先推荐一个博客 https://www.cnblogs.com/CreeperLKF/p/9045491.html
再来膜一个dalao https://www.cnblogs.com/AKCqhzdy/

理论

忘情水二分的模型给出一些选物品的条件,不同的组合法会有不同的价值,要求刚好选m个,最大化(最小化)总价值。
设选x个物品时的最优价值为g(x),以x为横坐标,价值为纵坐标,那么g(x)会呈现出一个上凸包的形状。
我们要求得显然是g(m),但是并不好直接求,因为这个形状我们只是大概知道而已。
那么考虑构造一些直线来切这个凸包,从而求出这个凸包的形状,进而找到g(m)的值。
我们枚举一条斜率为mid的直线通过一次O(N)的求解,可以找到这条直线的最高位置g(x),即找到恰好相切的截距。在求解过程中我们还可以顺便得到x,也就是说我们得到了凸壳上的一个点(x,g(x))。
关键不在于这个g(x),而在于x。这个x如果是小于m的,那么我们要考虑下次枚举时加大x,也就是说要加大斜率;同理,如果x大于了m,那么就要使下一次枚举的x变小,就是要减小斜率。相当于问题就是找到一个直线的斜率,使得这条直线与凸包的切点横坐标为x。横坐标很显然满足二分性,同样的,因为凸包的斜率单调递减,即斜率也是可以二分的。这两个都满足二分性,所以最终结果可以二分。
二分直到x=m时,结束二分,此时计算的ans才是真正的ans。设截距为f(x),因为g(x)=kx+f(x),所以f(x)=g(x)-kx,那么ans=g(x)-k*x。

例题

bzoj2654 tree

这是一道最最最基础的忘情水二分了。
二分一个值mid,白边所有白边加上mid之后做一个最小生成树,注意统计最小生成树中白边的条数tot和最小生成树的边权和sum。这个tot相当于x,sum相当于g(x)。所以当x大于need是减小斜率,注意考虑斜率与mid的关系(mid越大,选的白边越少,也就是斜率越小),实际上是要增大mid;反之同理。
在考虑一个细节,如果当x=x0时,白边数大于need,但x=x0+1时白边数小于need,这说明有在x=x0时,有很多白边与黑边的权值相同,此时应当减少一下白边数量。

代码

  1. #include<cstdio>
  2. #include<cstring>
  3. #include<algorithm>
  4. using namespace std;
  5. const int INF=1<<30;
  6. const int MAXN=50010,MAXM=100010;
  7.  
  8. int v,e,need;
  9.  
  10. struct E{int x,y,v,c;}a[MAXM];
  11. bool cmp(E a1,E a2)
  12. {
  13.     if(a1.v!=a2.v) return a1.v<a2.v;
  14.     return a1.c<a2.c;//debug
  15. }
  16.  
  17. int fa[MAXN];
  18. int findfa(int x){return x==fa[x]?fa[x]:findfa(fa[x]);}
  19.  
  20. int mintree;
  21. int kruskal(int mid)
  22. {
  23.     for(int i=1;i<=e;i++) if(a[i].c==0) a[i].v+=mid;
  24.     for(int i=1;i<=v;i++) fa[i]=i;
  25.     sort(a+1,a+e+1,cmp);//debug
  26.     mintree=0;int re=0;
  27.     for(int i=1,cnt=1;i<=e;i++)
  28.     {
  29.         int fx=findfa(a[i].x),fy=findfa(a[i].y);
  30.         if(fx!=fy)
  31.         {
  32.             fa[fx]=fy;
  33.             mintree+=a[i].v;
  34.             if(a[i].c==0) re++;
  35.             cnt++;if(cnt==v) break;
  36.         }
  37.     }
  38.     for(int i=1;i<=e;i++) if(a[i].c==0) a[i].v-=mid;
  39.     return re;
  40. }
  41.  
  42. int main()
  43. {
  44.     scanf("%d%d%d",&v,&e,&need);
  45.     for(int i=1;i<=e;i++) scanf("%d%d%d%d",&a[i].x,&a[i].y,&a[i].v,&a[i].c),a[i].x++,a[i].y++;
  46.      
  47.     int l=-105,r=105,ans;//???
  48.     while(l<=r)
  49.     {
  50.         int mid=l+r>>1;
  51.         if(kruskal(mid)>=need) ans=mintree-mid*need,l=mid+1;
  52.         else r=mid-1;
  53.     }
  54.     printf("%d\n",ans);
  55.     return 0;
  56. }

下面两篇没有那么模板了,我们另外来讲。

bzoj5311 贞鱼 题解链接

洛谷T51924 忘情 题解链接

总结

WQS二分一般是用来优化那种斜率满足递减性的DP,如果同时在不要求满足说要选多少个时能够比较快的的求出一个结果,还可以降维的话,那么WQS二分就是一个优选。

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/Monodyee/article/detail/436713
推荐阅读
相关标签
  

闽ICP备14008679号