CF1394C 计算几何

题目大意：

$n$个字符串，每个字符串仅包含$B，N$两种字符，可以对$BN-string$进行如下操作：

• 移走一个字符
• 移走$“BN”$,$“NB”$的子串
• 添加一个字符$B$，或$N$在$s$的末尾
• 添加$“BN”$或者$“NB”$在$s$的末尾

定义两个字符串相似的条件：$B$与$N$的个数分别都相等

且定义$dist(s,t)$代表$s$与$t$相似的最小操作次数

找到一个字符串$t$，使得$ma{x}_{i=1}^{n}dist(s_i,t)$最小

解题思路：（思路来源：WYXkk）

首先队对以上的操作进行分析：

• 假设字符串$t$的$B$，$N$都大于$s$，那么此时最好先对$t$进行第二种操作，再执行第一种操作，而第二种操作必须是子串，所以我们可以将$t$变成类似形式：$B...BN...N$
• 假设字符串$t$的$B$，$N$都小于$s$，那么此时先对$t$执行第四种操作，再执行第三种操作
• 最后这种情况，就是分别执行第一种，第三种操作即可

通过以上分析，操作可以变成如下：

• 移走一个字符或添加一个字符
• 添加/删除一个$B$以及一个$N$

对于最大值最小这类问题，一般都是二分答案，那么现在就是如何检验

假设二分的答案是$l$

• 可以将每个字符串的$B$的个数以及$N$的个数，分别当作横纵坐标画在坐标轴上

那么每个点移动的方向就如下：

最后范围就是：

所以只要对所有点这样的图形判断是否有交集即可

• 对于这种非常规则的图形，很多可以不用半平面交来判断

先写出方程：
$$\{\begin{array}{c}x\le x[i]+l\\ x\ge x[i]-l\\ y\le y[i]+l\\ y\ge y[i]-l\\ x-y\le x[i]-y[i]+l\\ x-y\ge x[i]-y[i]-l\end{array}$$

• 对于以上方程分别对每个范围的最小值取$max$，以及最大值取$min$

判断是否图形有交集：（$xmin$为$x$的最小范围，其余同理，设$zmin$为斜的最小距离）

• $xmin>xmax$ 或者 $ymin>ymax$ 或者$zmin>zmax$即交集为空

• 因为$zmin\le x-y\le zmax$，而$x-y$又有范围：$xmin-ymax\le x-y\le xmax-ymin$，所以还要判断两个区间是否有交集

最后得到一个答案$l$，还要找出有多少个$B$以及$N$

• 因为$zmin\le x-y\le zmax$，所以$zmin+ymin\le x\le zmax+ymax$，而$x$又有$xmin\le x\le xmax$，取交集中的一个数$x_0$就好
• 已知$x_0$以及$zmin\le x-y\le zmax$，所以$x_0-zmax\le y\le x0-zmin$，而$y$又有$ymin\le y\le ymax$，取交集即可

AC代码：

#include <bits/stdc++.h>
#define ft first
#define sd second
#define pb push_back
#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0) //不能跟puts混用
#define seteps(N) fixed << setprecision(N)
#define endl "\n"
const int maxn = 5e5 + 10;
using namespace std;
typedef long long ll;
typedef double db;
typedef pair<int, int> pii;
const ll mod = 1e9 + 7;
const int inf = 0x3f3f3f3f;
int n, x[maxn], y[maxn];
bool check(int l) {
    int xmin = -inf, xmax = inf, ymin = -inf, ymax = inf, zmin = -inf, zmax = inf;
    for (int i = 1; i <= n; i++) {
        xmin = max(xmin, x[i] - l);
        xmax = min(xmax, x[i] + l);
        ymin = max(ymin, y[i] - l);
        ymax = min(ymax, y[i] + l);
        zmin = max(zmin, x[i] - y[i] - l);
        zmax = min(zmax, x[i] - y[i] + l);
    }
    xmin = max(xmin, 0); ymin = max(ymin, 0);
    if (xmin > xmax || ymin > ymax || zmin > zmax) return 0;
    int zzmin = xmin - ymax, zzmax = xmax - ymin;
    if (zzmax < zmin || zmax < zzmin) return 0;
    return 1;
}
int main() {
    cin >> n;
    string s;
    for (int i = 1; i <= n; i++) {
        cin >> s;
        for (int j = 0; j < s.size(); j++)
            (s[j] == 'B') ? x[i]++ : y[i]++;
    }
    int lc = 0, rc = 1e6, mc, res;
    while (lc <= rc) {
        mc = (lc + rc) >> 1;
        if (check(mc)) rc = mc - 1, res = mc;
        else lc = mc + 1;
    }
    cout << res << endl;
    int xmin = -inf, xmax = inf, ymin = -inf, ymax = inf, zmin = -inf, zmax = inf;
    for (int i = 1; i <= n; i++) {
        xmin = max(xmin, x[i] - res);
        xmax = min(xmax, x[i] + res);
        ymin = max(ymin, y[i] - res);
        ymax = min(ymax, y[i] + res);
        zmin = max(zmin, x[i] - y[i] - res);
        zmax = min(zmax, x[i] - y[i] + res);
    }
    int x0 = min(xmax, zmax + ymax), y0 = min(x0 - zmin, ymax);
    for (int i = 1; i <= x0 + y0; i++)
        if (i <= x0) cout << 'B';
        else cout << 'N';
    cout << endl;
    return 0;
}