23.7.18 杭电暑期多校1部分题解

Belief_yfly

已于 2023-07-26 09:36:53 修改

阅读量67

点赞数

分类专栏：题解文章标签：换根 LCA 扩展欧几里得

于 2023-07-23 12:20:09 首次发布

本文链接：https://blog.csdn.net/weixin_44346868/article/details/131878223

版权

题解专栏收录该内容

40 篇文章 1 订阅

订阅专栏

1012 - Play on Tree

题目大意

给定一棵树，Alice和Bob轮流操作，每次操作可以删去一个节点以及它的所有子树，删掉根节点的人输，根节点随机的情况下，先手有多少概率必胜

解题思路

典型的博弈，当前点的 $s g$ 代表删掉它所有子树的游戏，那么叶子节点 $s g$ 为 $0$ ，当前点的 $s g$ 为子节点的 $s g + 1$ 的异或和，先求一遍以 $1$ 为根节点的树的是 $s g$ ，在进行一遍换根即可，遍历当前边时将父节点的 $s g$ 除去子节点的 $s g$ ，然后再将子节点的 $s g$ 叠加上父节点的 $s g$ ，返回时还原即可

code

#include <bits/stdc++.h>
using namespace std;
const int N = 2e5 + 9;
const int MOD = 1e9 + 7;
struct lol {int x, y;} e[N << 1];
int t, n, top[N], ans, num, sg[N];
int mul(int a, int b) {return 1ll * a * b % MOD;}
void ein(int x, int y) {
    e[++ ans].x = top[x];
    e[ans].y = y;
    top[x] = ans;
}
int pw(int a, int b) {
    int res = 1;
    while (b) {
        if (b & 1) res = mul(res, a);
        a = mul(a, a);
        b >>= 1;
    }
    return res;
}
void dfs(int x, int fa) {
    for (int i = top[x]; i; i = e[i].x) {
        int y = e[i].y;
        if (y == fa) continue;
        dfs(y, x);
        sg[x] ^= (sg[y] + 1);
    }
}
void dfs1(int x, int fa) {
    if (sg[x]) ++ num;
    for (int i = top[x]; i; i = e[i].x) {
        int y = e[i].y;
        if (y == fa) continue;
        sg[x] ^= (sg[y] + 1);
        sg[y] ^= (sg[x] + 1);
        dfs1(y, x);
        sg[y] ^= (sg[x] + 1);
        sg[x] ^= (sg[y] + 1);
    }
}
int main() {
    scanf("%d", &t);
    while (t --) {
        scanf("%d", &n);
        for (int i = 1, u, v; i < n; ++ i)
            scanf("%d%d", &u, &v), ein(u, v), ein(v, u);
        dfs(1, 0);
        dfs1(1, 0);
        printf("%d\n", mul(num, pw(n, MOD - 2)));
        ans = num = 0;
        for (int i = 1; i <= n; ++ i) top[i] = sg[i] = 0;
    }
    return 0;
}

1001 - Hide-And-Seek Game

题目大意

给定一棵大小为 $n$ 的树，有 $m$ 组查询，每组查询包含 $s_a,\space t_a,\space s_b,\space t_b$ 分别表示 $a$ 的起点和终点， $b$ 的起点和终点， $a$ 和 $b$ 将从起点出发，每秒跑一个点地折返跑，问最早相遇的点（或者不能相遇）

解题思路

先将两个人的路径标记寻找重合的点，对于每个可能的点， $a$ 到达的时间和 $b$ 到达的时间都可以用两个函数来表示，即 $2*len_a*t1\pm b1$ 和 $2*len_b*t2\pm b2$ （ $len_a,\space len_b$ 为 $a,\space b$ 路径的长度， $t1,\space t2$ 为未知数， $b1,\space b2$ 为当前点到 $sa,\space sb$ 的距离，函数的结果就是到达当前点的时间），用二者相等写出四个等式，化简后用扩展欧几里德算法分别求解，并且记录最小值和产生最小值的点，再处理亿点点细节就可以啦

code

#include <bits/stdc++.h>
using namespace std;
const int N = 3e3 + 9;
struct lol {int x, y;} e[N << 1];
int t, n, m, top[N], ans, sa, ta, sb, tb, f[N], dep[N], p[N], sum, k1, k2, id;
int exgcd(int a, int b, int &x, int &y) {//扩展欧几里德
    if (!b) {x = 1, y = 0; return a;}
    int d = exgcd(b, a % b, x, y);
    int z = x; x = y, y = z - y * (a / b);
    return d;
}
void ein(int x, int y) {//建图
    e[++ ans].x = top[x];
    e[ans].y = y;
    top[x] = ans;
}
void dfs(int x, int fa) {//预处理求深度和父节点，方便求lca和路径长度
    f[x] = fa; dep[x] = dep[fa] + 1;
    for (int i = top[x]; i; i = e[i].x) {
        int y = e[i].y;
        if (y == fa) continue;
        dfs(y, x);
    }
}
int getlca(int x, int y) {//求lca
    if (dep[x] < dep[y]) swap(x, y);
    while (dep[x] > dep[y]) x = f[x];
    while (x != y) x = f[x], y = f[y];
    return x;
}
void bj(int x, int rt) {while (x != rt) ++ p[x], x = f[x];}//标记路径上的点
void tj(int x, int rt) {//统计答案并维护最小值和产生最小值的点
    do {
        if (p[x] <= 1) {x = f[x]; continue;}
        int nl1 = getlca(sa, x), nl2 = getlca(sb, x);
        int b1 = dep[sa] + dep[x] - 2 * dep[nl1];
        int b2 = dep[sb] + dep[x] - 2 * dep[nl2];
        int x0, y0, val;
        int d = exgcd(k1, k2, x0, y0);
        //化简后应该k1和k2中有一个是负数，我不知道怎么搞，就把负号放到xx和yy上处理了，有另外方法可以简单一点
        //注意一定要把特解单独存下来，不然后面一直用的时候会出现用非当前函数的解计算导致出错
        val = b2 - b1;
        if (val % d == 0) {
            int xx = x0 * (val / d), yy = y0 * (val / d);
            xx = (xx % (k2 / d) + (k2 / d)) % (k2 / d);//求最小非负整数解
            yy = (val - xx * k1) / k2;
            if (yy <= 0) if (k1 * xx + b1 < sum)//xx为非负整数，yy就必须是非正整数
                sum = k1 * xx + b1, id = x;
        }
        val = -val;//等是两边移项，k1和k2符号交换，即xx和yy符号交换
        if (val % d == 0) {
            int xx = x0 * (val / d), yy = y0 * (val / d);
            yy = (yy % (k1 / d) + (k1 / d)) % (k1 / d);
            xx = (val - yy * k2) / k1;
            if (xx <= 0) if (k2 * yy + b2 < sum)
                sum = k2 * yy + b2, id = x;
        }
        val = b2 + b1 - k1;//可以发现xx和yy有时不能为零，我的处理是让它变成xx+1和yy+1
        if (val % d == 0) {
            xx = x0 * (val / d), yy = y0 * (val / d);
            xx = (xx % (k2 / d) + (k2 / d)) % (k2 / d);
            yy = (val - xx * k1) / k2;
            if (yy <= 0) if (k1 * (xx + 1) - b1 < sum)
                sum = k1 * (xx + 1) - b1, id = x;
        }
        val = -val;
        if (val % d == 0) {
            int xx = x0 * (val / d), yy = y0 * (val / d);
            yy = (yy % (k1 / d) + (k1 / d)) % (k1 / d);
            xx = (val - yy * k2) / k1;
            if (xx <= 0) if (k2 * yy + b2 < sum)
                sum = k2 * yy + b2, id = x;
        }
        val = -b2 - b1 + k2;
        if (val % d == 0) {
            int xx = x0 * (val / d), yy = y0 * (val / d);
            xx = (xx % (k2 / d) + (k2 / d)) % (k2 / d);
            yy = (val - xx * k1) / k2;
            if (yy <= 0) if (k1 * xx + b1 < sum)
                sum = k1 * xx + b1, id = x;
        }
        val = -val;
        if (val % d == 0) {
            int xx = x0 * (val / d), yy = y0 * (val / d);
            yy = (yy % (k1 / d) + (k1 / d)) % (k1 / d);
            xx = (val - yy * k2) / k1;
            if (xx <= 0) if (k2 * (yy + 1) - b2 < sum)
                sum = k2 * (yy + 1) - b2, id = x;
        }
        val = -b2 + b1 + k2 - k1;
        if (val % d == 0) {
            int xx = x0 * (val / d), yy = y0 * (val / d);
            xx = (xx % (k2 / d) + (k2 / d)) % (k2 / d);
            yy = (val - xx * k1) / k2;
            if (yy <= 0) if (k1 * (xx + 1) - b1 < sum)
                sum = k1 * (xx + 1) - b1, id = x;
        }
        val = -val;
        if (val % d == 0) {
            int xx = x0 * (val / d), yy = y0 * (val / d);
            yy = (yy % (k1 / d) + (k1 / d)) % (k1 / d);
            xx = (val - yy * k2) / k1;
            if (xx <= 0) if (k2 * (yy + 1) - b2 < sum)
                sum = k2 * (yy + 1) - b2, id = x;
        }
        x = f[x];
    } while (x != rt && x != f[rt]);
}
int main() {
    scanf("%d", &t);
    while (t --) {
        scanf("%d%d", &n, &m);
        for (int i = 1, u, v; i < n; ++ i)
            scanf("%d%d", &u, &v), ein(u, v), ein(v, u);
        dfs(1, 0);
        for (int i = 1; i <= m; ++ i) {
            sum = 1e9;
            scanf("%d%d%d%d", &sa, &ta, &sb, &tb);
            int la = getlca(sa, ta), lb = getlca(sb, tb);
            bj(sa, la); bj(ta, la); ++ p[la];//lca单独处理防止多次计算，下方同理
            bj(sb, lb); bj(tb, lb); ++ p[lb];
            k1 = 2 * (dep[sa] + dep[ta] - 2 * dep[la]);//即2*len_a
            k2 = 2 * (dep[sb] + dep[tb] - 2 * dep[lb]);
            tj(sa, la); tj(ta, la); tj(la, la);//只用在某一方的路径上统计
            if (sum == 1e9) printf("-1\n");
            else printf("%d\n", id);
            for (int j = 1; j <= n; ++ j) p[j] = 0;//记得每个查询清空
        }
        ans = 0;//记得每组数据清空
        for (int i = 1; i <= n; ++ i) top[i] = 0;
    }
    return 0;
}

Belief_yfly

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
23.7.18 杭电暑期多校1部分题解

的距离，函数的结果就是到达当前点的时间），用二者相等写出四个等式，化简后用扩展欧几里德算法分别求解，并且记录最小值和产生最小值的点，再处理亿点点细节就可以啦。将从起点出发，每秒跑一个点地折返跑，问最早相遇的点（或者不能相遇）先将两个人的路径标记寻找重合的点，对于每个可能的点，到达的时间都可以用两个函数来表示，即。组查询，每组查询包含。
复制链接

扫一扫