【知识小结】凸优化

算法简述

凸优化是一个神奇的东西，通过加权，二分权值，把数量限制去掉，从而优化DP，或者其他很难满足的限制

这里有讲解和证明
wqs的论文

这个东西bytedance冬令营时就学了，但是没有完全掌握。今天再翻以前的题，竟然没有立即反应过来。
赶快复习，然后写几份代码！
学习算法不能偷懒，一定要实现！还要定期复习！否则会忘的！

例题

CF739E Gosha is hunting
首先这道题可以费用流。只需要把两个同时选多出的贡献：ui * pi在第二个流量的费用中减去。
然后费用流可以用贪心和数据结构维护，CF官方题解，然而似乎很麻烦，不想看
这道题的神奇之处是wqs二分
详细题解

一定要注意精度！
一点点的差距都会导致多选一个物品，必须使用eps，并且保证cost稍微多一点（使用二分后的r，用l会挂！），否则会爆炸
最后因为此题胶囊肯定用完更好，就直接贡献+a * ca+b * cb,不要用最后一次DP的选的个数。二分的时候输出一下选的个数就会发现一直在不规则跳动（a，a + 1），好像不会收敛。因为精度有很大问题！
整数二分不靠谱，至少！1e9+7这个质数还不够大，好像直接小数二分50次更好一些
WA了无数发的血的教训！

#include<bits/stdc++.h>
using namespace std;

#define rep(i,l,r) for(register int i = l ; i <= r ; i++)
#define repd(i,r,l) for(register int i = r ; i >= l ; i--)
#define rvc(i,S) for(register int i = 0 ; i < (int)S.size() ; i++)
#define rvcd(i,S) for(register int i = ((int)S.size()) - 1 ; i >= 0 ; i--)
#define fore(i,x)for (register int i = head[x] ; i ; i = e[i].next)
#define forup(i,l,r) for (register int i = l ; i <= r ; i += lowbit(i))
#define fordown(i,id) for (register int i = id ; i ; i -= lowbit(i))
#define pb push_back
#define prev prev_
#define stack stack_
#define mp make_pair
#define fi first
#define se second
#define lowbit(x) ((x)&(-(x)))
#define eps 1e-10

typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef pair<bool,double> pr;

const double inf = 2e8;
const int N = 3e6 + 10;
const int maxn = 100020;
const ll mod = 1e9 + 7;

double p[maxn],u[maxn],f[maxn];
int n,a,b,numa[maxn],numb[maxn];

bool check2(double ca,double cb){
   
	rep(i,1,n) f[i] = -inf , numa[i] = numb[i] = 0;
	f[0] = 0;
	rep(i,1,n){
   
		f[i] = f[i - 1] , numa[i] = numa[i - 1] , numb[i] = numb[i - 1];
		if ( f[i - 1] + p[i] - ca - eps > f[i] ){
   
		   	f[i] = f[i - 1] + p[i] - ca , numa[i] = numa[i - 1] + 1