【三分法】题集

---

Line belt

思路： 三分套三分。当AB上的离开点E固定时，当进入点F从C到D时，总时间总是递减，或先递减后递增。对于每个三分的E，三分F求极小值点。

这里有一个规律：

以D为起点，DC为方向作极坐标系。过E作EL垂直于DC于点L。设|DF|=x，则时间为t(x)

• 当F从正无穷远接近D至L时，我们易知函数t是递减的。
• 当F从L接近D至D时：
• 1.当外部速度比线上速度更快或相等时，显然函数t是递减的。
• 2.反之，函数t先递减后递增（不会证）。

我们可以保证t是单调或单谷的，E确定后就可以三分点F求极（最）小值了。

---

更新（2021.10.21）：发现了一篇文章，和本题模型是一样的。

• 初中数学几何最值问题之“胡不归”问题
• 初中几何最值——胡不归问题详解
• 函数下凹性质由图易得。

AC代码：

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;

const int N=1e5+10;
typedef long long LL;
const double eps = 1e-5;

struct point{
    double x, y;
    void input(){
        scanf("%lf %lf", &x, &y);
    }
};
point operator + (point a, point b) { return {a.x + b.x, a.y + b.y}; }
point operator - (point a, point b) { return {a.x - b.x, a.y - b.y}; }
point operator * (point a, double b) { return {a.x * b, a.y * b}; }
point operator / (point a, double b) { return {a.x / b, a.y / b}; }
double get_dis(point a, point b) { return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y)); }
double length(point a) { return sqrt(a.x * a.x + a.y * a.y); }


point a, b, c, d;
double p, q, r;



double in_for(point e) // 返回从AB上的点e为起点离开，到终点D的最短距离
{
    double le = 0, ri = 1;
    point v = d - c;
    while(ri - le > eps)
    {
        double mid1 = le + (ri - le) / 3, mid2 = ri - (ri - le) / 3;
        point v1 = c + v * mid1, v2 = c + v * mid2;
        if(get_dis(e, v1) / r + get_dis(v1, d) / q < get_dis(e, v2) / r + get_dis(v2, d) / q) ri = mid2;
        else le = mid1;
    }
    return get_dis(e, c + v * ri) / r + get_dis(c + v * ri, d) / q;
}

double out_for()
{
    double le = 0, ri = 1;
    point v = b - a;
    while(ri - le > eps)
    {
        double mid1 = le + (ri - le) / 3, mid2 = ri - (ri - le) / 3;
        point v1 = a + v * mid1, v2 = a + v * mid2;
        if(in_for(v1) + get_dis(a, v1) / p < in_for(v2) + get_dis(a, v2) / p) ri = mid2;
        else le = mid1;
    }
    return in_for(a + v * ri) + get_dis(a, a + v * ri) / p;
}



int main()
{
    int T; scanf("%d", &T);
    while(T--)
    {
        a.input(); b.input(); c.input(); d.input();
        scanf("%lf %lf %lf", &p, &q, &r);

        printf("%.2f\n", out_for());
    }

    system("pause");
    return 0;
}
// 136.602540

---

Party All the Time

这道题是浮点数三分板子，感觉有一些可拓展结论，就写一下。

思路： 在坐标轴上找一个浮点坐标x，使得 $f(x)=\sum |s[i]-x{|}^3\ast w[i]$ 最小。易知 $f(x)$ 是单谷凹函数，三分即可。

拓展：

• 对于任意的 $k\in [1,\infty )$ ,当$s[i]\in R,w[i]>=0$ ，函数$f(x)=\sum |s[i]-x{|}^k\ast w[i]$ 都是单谷凹函数。

拓展不会证，正确性不能保证。拍了几千组数据，应该没问题。

AC代码：

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;

const int N=5e4+10;
typedef long long LL;
const double eps = 1e-8;

double s[N];
double w[N];
int n;

double f(double k)
{
    double sum = 0;
    for(int i=1; i<=n; i++)
        sum += fabs(s[i] - k) * (s[i] - k) * (s[i] - k) * w[i];
    return sum;
}

int main()
{
    int T; scanf("%d", &T);
    for(int _=1; _<=T; _++)
    {
        scanf("%d", &n);
        double l = 2e6, r = -2e6;
        for(int i=1; i<=n; i++)
        {
            scanf("%lf %lf", &s[i], &w[i]);
            l = min(l, s[i]); r = max(r, s[i]);
        }

        while(r - l > eps)
        {
            double mid1 = l + (r - l) / 3, mid2 = r - (r - l) / 3;
            if(f(mid1) > f(mid2)) l = mid1;
            else r = mid2;
        }
        printf("Case #%d: %.0f\n", _, f(r));
    }

    return 0;
}

---

感觉好多三分的题，结果函数的极值唯一性都很难证，都是感觉出来的，离谱。

---

【整数三分】CF1355E - Restorer Distance

题解：三分专题 几道练习题_forward77-CSDN博客

思路：以最终高度为自变量，总花费就是下凹函数。操作一二等价于操作三，比较一下A+R和M得到砌墙策略。三分。

AC代码：

#include<iostream>
#include<cstdio>
using namespace std;
 
const int N=1e5+10;
typedef long long LL;
 
LL h[N];
LL a, b, c, n;
 
LL f(LL k)
{
    LL res = 0, up = 0, down = 0;
    for(int i=1; i<=n; i++)
        if(h[i] < k) up += k - h[i];
        else down += h[i] - k;
 
    if(a + b <= c) res = a * up + b * down;
    else{
        res = min(up, down) * c;
        if(up < down) res += (down - up) * b;
        else res += (up - down) * a;
    }
    return res;
}
 
int main()
{
    scanf("%lld %lld %lld %lld", &n, &a, &b, &c);
    LL l = 2e9, r = -2e9;
    for(int i=1; i<=n; i++)
    {
        scanf("%lld", &h[i]);
        l = min(l, h[i]); r = max(r, h[i]);
    }
    while(l < r)
    {
        LL mid1 = l + (r - l) / 3, mid2 = r - (r - l) / 3;
        if(f(mid1) < f(mid2)) r = mid2 - 1;
        else l = mid1 + 1;
        // cout<<"# "<<l<<" "<<r<<endl;
    }
    printf("%lld", f(r));
 
    return 0;
}

---

2021 ICPC济南 D Arithmetic Sequence

题意：给你一个序列，每次操作对某个元素++或- -，求最少操作次数使其变成等差数列。

思路：猜测与公差 $d$ 有关的代价函数是凹函数，三分公差即可。$f(d)$ 的计算方法是：$d$ 已知，把相邻元素的差分抵消掉，即使$a_1-=d\ast 1,a_2-=d\ast 2,a_3-=d\ast 3...$，以此类推。问题转化为：求一个数 $x$ ，使得 $a_{0\cdots n-1}$ 到 $x$ 的距离之和最短，也就是经典的货仓选址问题。

注意：1. $f(d)$ 可能会爆 long long，记得开 __int128。

2.三分不优化的话就不能使用排序了，找中位数要使用 nth_element() 函数。链接：nth_element()函数用法

AC代码：

#include <bits/stdc++.h>

using namespace std;

const int N = 2e5 + 10;
typedef long long LL;

LL a[N];
__int128 back[N];
int n;

__int128 ans(__int128 mid)
{
	__int128 sum = 0;
	for(int i=0; i<n; i++) back[i] = a[i];
	for(int i=1; i<n; i++) back[i] -= mid * i;
	nth_element(back, back + (n - 1) / 2, back + n);
	__int128 m = back[(n - 1) / 2];
	for(int i=0; i<n; i++){
		__int128 del = back[i] - m;
		if(del < 0) del = -del;
		sum += del;
	}
	return sum;
}


int main()
{
	scanf("%d", &n);
	LL l = 1e18, r = -1e18;
	LL d;
	
	for(int i=0; i<n;i++) scanf("%lld", &a[i]);
	for(int i=1; i<n; i++) d = a[i] - a[i - 1], l = min(l, d), r = max(r, d);
	
	while(l < r)
	{
		LL mid1 = l + (r - l) / 3;
		LL mid2 = r - (r - l) / 3;
		
		if(ans(mid1) < ans(mid2)) r = mid2 - 1;
		else l = mid1 + 1;
	}
	
	LL anss = (LL)ans(r);
	
	printf("%lld", anss);
	
	return 0;
}