2021牛客多校4

E-Tree Xor（线段树，区间交）

题意：

给定一棵树，已知：每个点$u$ 的 权值范围 $[l_i,r_i]$ ，每条边$<u,v>$对应两点权值的异或值 $W_u\oplus W_v$

求合法的 $w_{1\dots n}$ 的解数。

思路：

因为是一棵树，所以，如果确定一个点的权值，就可以确定所有点的权值。

令 $w[1]=0$ 可以得到一组解（不一定满足权值范围要求）$w[1\dots n]$

如果 $w[1]=a$ ，那么其他所有点的权值都要 $\oplus a$ 。

那么对于 $\forall i\in [1,n],l_i\le w[i]\oplus a\le r_i$ ，即 $a\in ([l_i,r_i]\oplus w[i])$ 。对 $n$ 个这样的区间，求交集中点的个数，即 合法 $a$ 的数量 。

但是因为一个区间异或一个数后，它就不是一个完整的区间了，会变成多个区间。

有结论：

$[...0000，...1111]\oplus xor$ ：该区间（$...$ 表示相同的前缀）异或一个数后，得到的结果依然是在一个区间中，

并且区间的后 $k$ 位依旧是 $0000->1111$。

原区间为 $[L,R]$ ，满足前 $n-k$ 位为相同前缀，后 $k$ 位是 $00..00->11..11$。那么得到的新区间过程：

$\text{异或的数位 a}$
$prel=L\&(((1<<n-1)-1)\oplus ((1<<k)-1))$
$prer=R\&(((1<<n-1)-1)\oplus ((1<<k)-1))$
$pree=a\&(((1<<n-1)-1)\oplus ((1<<k)-1))$

新区间：$[prel\oplus pree,(prel\oplus pree)+(1<<k)-1]$

一个区间，满足前缀相同，后缀从 $00\dots 0$ 到 $11\dots 1$ 。不就是线段树上一个节点所代表的区间吗

线段树一个节点代表的区间异或一个数得到的结果还是一个区间，那么就可以将 $[l_i,r_i]$ 分成 $logW$个区间，每个小区间异或一个数得到一个区间。

将所有的区间分组统计出来，然后就是 $n$ 组区间 求交集：将区间排序，用扫描线的方式从前往后扫一遍，就能得到交集长度。

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e5 + 10;
int n;
ll l[N], r[N];
ll w[N];
ll a[N]; // 一组特解
int vis[N];

struct edge{
    int v;
    ll cost;
    int nxt;
}es[N << 1];
int head[N], tot;
void init(int _n) {
    tot = 0;
    for (int i = 0; i <= _n; i++)
        head[i] = -1, vis[i] = 0;
}
void add_edge(int u, int v, ll cost) {
    es[tot].v = v;
    es[tot].cost = cost;
    es[tot].nxt = head[u];
    head[u] = tot++;
}

void dfs(int u, int fa, ll ww) { // 找到一组解，并不一定要满足区间要求
    a[u] = ww;
    vis[u] = 1;
    for (int i = head[u]; i != -1; i = es[i].nxt) {
        int v = es[i].v;
        if (vis[v] || v == fa) continue;
        dfs(v, u, ww ^ es[i].cost);
    }
}

vector<pair<int, int> > seg;

void modify(int s, int t, int l, int r, int dep, int xx) {
    if (dep < 0) return;
    if (l > t || r < s) return;
    if (s <= l && r <= t) {
        int pl = l & (((1 << 30) - 1) ^ ((1 << dep) - 1)); // 取出 L 的前 30 - dep 位的值
//        int pr = r & (((1 << 30) - 1) ^ ((1 << dep) - 1)); // 取出 R 的前 30 - dep 位的值
        int pxor = xx & (((1 << 30) - 1) ^ ((1 << dep) - 1)); // 取出异或值 xx 的前 30 - dep 位的值
        seg.push_back({pl ^ pxor, pl ^ pxor + (1 << dep) - 1}); // 右区间直接加上 (1 << dep) - 1 即可
//        pair<int, int> now = *seg.rbegin();
//        cout << now.first << " " << now.second << "\n";
        return;
    }
    int mid = (l + r) >> 1;
    modify(s, t, l, mid, dep - 1, xx);
    modify(s, t, mid + 1, r, dep - 1, xx);
}

int get_ans() {
    vector<pair<int, int> > p;
    for (auto it : seg) {
        p.push_back({it.first, -1});
        p.push_back({it.second, 1});
    }
    int ans = 0;
    sort(p.begin(), p.end());
    int len = p.size();
    int cnt = 0; // 当前包含了多少个区间
    for (int i = 0; i < len; i++) {
//        cout << p[i].first << " " << p[i].second << "\n";
        cnt -= p[i].second;
        if (cnt == n) { // 已经包含了 n 个区间
            if (i + 1 < len) ans += p[i + 1].first - p[i].first + 1;
        }
    }
    return ans;
}

int main() {
    scanf("%d", &n);
    for (int i = 0; i < n; i++) {
        scanf("%d%d", &l[i], &r[i]);
    }
    init(n);
    for (int i = 1; i < n; i++) {
        int u, v;
        ll cost;
        scanf("%d%d%lld", &u, &v, &cost);
        add_edge(u - 1, v - 1, cost);
        add_edge(v - 1, u - 1, cost);
    }
    dfs(0, -1, 0);
//    for (int i = 0; i < n; i++)
//        printf("%lld ", a[i]);
//    printf("\n");
    seg.clear();
    for (int i = 0; i < n; i++) {
        modify(l[i], r[i], 0, (1 << 30) - 1, 30, a[i]);
//        cout << "\n";
    }
    printf("%d\n", get_ans());
    return 0;
}

/*
4
0 7
1 6
2 5
3 4
1 2 0
1 3 7
2 4 6
*/

J-Average（二分）

题意：

给一个 $n\times m$ 的矩阵 $W$，$W_{i,j}=a_i+b_j$

要求一个子矩阵 宽 $\ge x$ , 高 $\ge y$ 的平均值最大。

思路：

将题目要求的平均值用公式表示出来：
$$\begin{gathered} \frac{y{\sum }_{i=lx}^{rx}b[i]+x{\sum }_{i=ly}^{ry}a[i]}{(rx-lx)(ry-ly)} \\ =\frac{{\sum }_{i=lx}^{rx}b[i]}{rx-lx}+\frac{{\sum }_{i=ly}^{ry}a[i]}{ry-ly} \end{gathered}$$
然后就是求一维上，长度 $>=x$ 的 最大平均值。

二分答案，$check$ 函数实现：将当前数组全部减去 $mid$ ，再判断有无长度 $>=x$ 的字段和 $>=0$ 。

代码：

代码摸了，放一份队友的。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const double eps = 1e-12;
const int maxn = 2e5 + 5;
const ll INF = 2e6;
double a[maxn], b[maxn], sum[maxn];
int n, m, x, y;
int dcmp(double x) {
    if (fabs(x) < eps) return 0;
    return x > 0 ? 1 : -1;
}
int check1(double mid) {
    for(int i = 1; i <= n; i++) {
        sum[i] = a[i] - mid + sum[i - 1];
    }
    double mi = 2e6;
    for(int i = x; i <= n; i++) {
        mi = min(mi, sum[i - x]);
        if(sum[i] - mi >= eps) return true;
    }
    return false;
}
int check2(double mid) {
    for(int i = 1; i <= m; i++) {
        sum[i] = b[i] - mid + sum[i - 1];
    }
    double mi = 2e6;
    for(int i = y; i <= m; i++) {
        mi = min(mi, sum[i - y]);
        if(sum[i] - mi >= eps) return true;
    }
    return false;
}
int main() {
    scanf("%d%d%d%d", &n, &m, &x, &y);
    double ans1 = 0, ans2 = 0;
    double l = 0, r = 0;
    for(int i = 1; i <= n; i++){
        scanf("%lf", &a[i]);
        r = max(r, a[i]);
    }
    for(int i = 0; i < 2000; i++){
        double mid = (l + r) / 2;
        if(check1(mid)) {
            l = mid;
            ans1 = mid;
        } else {
            r = mid;
        }
    }
    l = 0, r = 0;
    for(int i = 1; i <= m; i++) {
        scanf("%lf", &b[i]);
        r = max(r, b[i]);
    }
    for(int i = 0; i < 2000; i++){
        double mid = (l + r) / 2;
        if(check2(mid)) {
            l = mid;
            ans2 = mid;
        } else {
            r = mid;
        }
    }
    printf("%.10f", ans1 + ans2);
}