标签：ch 剖分 int CSPS 树链 2019 联考 dp define

因为版权原因，不放题面。

幸福

Solution

只会写 \(70\) 分的找规律代码，可以考虑每个斐波那契数对卷积的贡献，发现贡献是类似于倒斐波那契数列，就可以 \(O(n)\) 做了，有 \(70\) 分。

std 是化简一下式子，化简结果为

\[F_n = F_{n-1} + F_{n-2} + f_n\]

发现这是一个递推式，用矩阵乘法加速即可。但是我不会矩阵乘法，先咕了。不会矩阵乘法也可以用退出来的递推式拿 \(70\) 分。

Code

#include <bits/stdc++.h>
using namespace std;
#define re register
#define F first
#define S second
typedef long long ll;
typedef pair<int, int> P;
const int N = 1e7 + 5, mod = 998244353;
const int INF = 0x3f3f3f3f;
int read() {
    int x = 0, f = 0; char ch = 0;
    while (!isdigit(ch)) f |= ch == '-', ch = getchar();
    while (isdigit(ch)) x = (x << 3) + (x << 1) + (ch ^ 48), ch = getchar();
    return f ? -x : x;
}
ll fb[N], cnt[N]; 
ll ans = 0;
int main(){
    int n = read();
    fb[0] = fb[1] = 1;
    for (int i = 2; i <= n; i++) fb[i] = (fb[i - 1] + fb[i - 2]) % mod;
    cnt[n] = 1;
    for (int i = n - 1; i >= 0; i--) cnt[i] = (cnt[i + 1] + cnt[i + 2] + 1) % mod;
    ll ans = 0;
    for (int i = 0; i <= n; i++) ans = (ans + fb[i] * cnt[i]) % mod;
    printf("%lld\n", ans);
    return 0;
}

树链剖分

Solution

可以用 lca + 树上差分将所有的路径转换为边权，那么我们要找到边权之和最大的树链剖分。因为这是一棵无根树，所以要进行换根 dp，分为两次树形 dp。不过可以暴力枚举根，所以有 \(65\) 分是送的。

用 \(f_i\) 表示以任意点为根的树下，以 \(i\) 为根的子树的树链剖分。可以选 \(1\) 为根。这一步可以用树形 dp 搞定。

考虑如何换根，设 \(F_i\) 为以 \(i\) 为根的树的树链剖分。如果把点 \(y\) 作为根，它的父亲为 \(x\)，他在 \(x\) 的树链剖分中的儿子是 \(z\)，从 \(1\) 开始，从上往下进行树形 dp，其中 \(F_1 = f_1\)。换根需要注意换根之后还要满足每个节点的儿子只选一个， 如果在 \(x\) 的树链剖分中，\(y\) 是 \(x\) 的儿子，那么换根后就出现了 \(x\) 和 \(z\) 都是 \(y\) 的儿子的情况，违背了树链剖分的要求，所以我们只能将 \(x\) 连向 \(y\) 改到 \(x\) 的所有出边中仅次于连向 \(y\) 的边的次优边，或者不选择 \(z\)，而把 \(x\) 的最优边和次优边加入 \(F_y\) 中。如果没有这个问题，那么就选择 \(x\) 和 \(y\) 的最优边即可。对于刚刚的最优边次优边等问题，可以在第一次树形 dp 中预处理出每个点连向它所有儿子中，最大边权的和次大边权。

总而言之，一道好的树论题。

Code

#include <bits/stdc++.h>
using namespace std;
#define re register
#define F first
#define S second
typedef long long ll;
typedef pair<int, int> P;
const int N = 1e6 + 5, M = 1e6 + 5;
const int INF = 0x3f3f3f3f;
int read() {
    int x = 0, f = 0; char ch = 0;
    while (!isdigit(ch)) f |= ch == '-', ch = getchar();
    while (isdigit(ch)) x = (x << 3) + (x << 1) + (ch ^ 48), ch = getchar();
    return f ? -x : x;
}
struct edge{
    int to, nxt, cost;
}e[M];
int head[N], tot;
void addedge(int x, int y){
    e[++tot].to = y; e[tot].nxt = head[x]; head[x] = tot;
}
int dep[N], fa[N][30], lg[N];
ll f[N], F[N]; int mx1[N], mx2[N], val[N];
int n, m; ll sum, ans;
void dfs1(int x, int pre){
    dep[x] = dep[pre] + 1; fa[x][0] = pre;
    for (int i = 1; (1 << i) <= dep[x]; i++) fa[x][i] = fa[fa[x][i - 1]][i - 1];
    for (int i = head[x]; i; i = e[i].nxt) if (e[i].to != pre) dfs1(e[i].to, x);
}
int lca(int x, int y){
    if (dep[x] < dep[y]) swap(x, y);
    while (dep[x] > dep[y]) x = fa[x][lg[dep[x] - dep[y]] - 1];
    if(x == y) return x;
    for (int k = lg[dep[x]] - 1; k >= 0; k--) if (fa[x][k] != fa[y][k]) x = fa[x][k], y = fa[y][k];
    return fa[x][0];
}
void dfs2(int x, int pre){
    for (int i = head[x]; i; i = e[i].nxt){
        int y = e[i].to;
        if (y != pre){
            dfs2(y, x); val[x] += val[y]; f[x] += f[y];
            if (mx1[x] <= val[y]) mx2[x] = mx1[x], mx1[x] = val[y]; 
            else if (mx2[x] <= val[y]) mx2[x] = val[y];
        }
    }
    f[x] += mx1[x];
}
void dfs3(int x, int pre){
    for (int i = head[x]; i; i = e[i].nxt){
        int y = e[i].to;
        if (mx1[x] <= val[x]) mx2[x] = mx1[x], mx1[x] = val[x];
        else if (mx2[x] <= val[x]) mx2[x] = val[x];
        if (y != pre){
            if (mx1[x] == val[y]) F[y] = F[x] - mx1[x] + mx2[x] - mx1[y] + max(mx1[x], mx1[y]);
            else F[y] = F[x] - mx1[y] + max(mx1[y], val[y]); 
            dfs3(y, x);
        }
    }
}
int main(){
    int n = read(), m = read();
    for (int i = 1; i <= n; i++) lg[i] = lg[i - 1] + (1 << lg[i - 1] == i);
    for (int i = 1; i < n; i++){
        int x = read(), y = read();
        addedge(x, y); addedge(y, x);
    }
    dfs1(1, 0);
    for (int i = 1; i <= m; i++){
        int x = read(), y = read();
        int k = lca(x, y);
        val[x]++, val[y]++, val[k] -= 2;
        sum += dep[x] + dep[y] - dep[k] * 2;
    }
    dfs2(1, 0); F[1] = f[1]; dfs3(1, 0);
    for (int i = 1; i <= n; i++) ans = max(ans, F[i]);  
    printf("%lld\n", sum - ans);
    return 0;
}

小 E 和小 F strikes again

好像是一个乱搞 dp，咕咕。

标签：ch,剖分,int,CSPS,树链,2019,联考,dp,define
来源： https://www.cnblogs.com/lyfoi/p/11616104.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。