洛谷 P1351 [NOIP 2014 提高组] 联合权值-普及/提高-

P1351 [NOIP 2014 提高组] 联合权值

题目背景

NOIP2014 提高组 D1T2

题目描述

无向连通图 GGGnnn 个点,n−1n-1n1 条边。点从 111nnn 依次编号,编号为 iii 的点的权值为 WiW_iWi,每条边的长度均为 111。图上两点 (u,v)(u, v)(u,v) 的距离定义为 uuu 点到 vvv 点的最短距离。对于图 GGG 上的点对 (u,v)(u, v)(u,v),若它们的距离为 222,则它们之间会产生 Wv×WuW_v \times W_uWv×Wu 的联合权值。

请问图 GGG 上所有可产生联合权值的有序点对中,联合权值最大的是多少?所有联合权值之和是多少?

输入格式

第一行包含 111 个整数 nnn

接下来 n−1n-1n1 行,每行包含 222 个用空格隔开的正整数 u,vu,vu,v,表示编号为 uuu 和编号为 vvv 的点之间有边相连。

最后 111 行,包含 nnn 个正整数,每两个正整数之间用一个空格隔开,其中第 iii 个整数表示图 GGG 上编号为 iii 的点的权值为 WiW_iWi

输出格式

输出共 111 行,包含 222 个整数,之间用一个空格隔开,依次为图 GGG 上联合权值的最大值和所有联合权值之和。由于所有联合权值之和可能很大,输出它时要对 100071000710007 取余。

输入输出样例 #1

输入 #1

5  
1 2  
2 3
3 4  
4 5  
1 5 2 3 10

输出 #1

20 74

说明/提示

样例解释

本例输入的图如上所示,距离为 222 的有序点对有(1,3)(1,3)(1,3)(2,4)(2,4)(2,4)(3,1)(3,1)(3,1) 、$(3,5) 、、(4,2)$ 、$(5,3) $。

其联合权值分别为 2,15,2,20,15,202,15,2,20,15,202,15,2,20,15,20。其中最大的是 202020,总和为 747474

数据说明

  • 对于 30%30\%30% 的数据,1<n≤1001 < n \leq 1001<n100
  • 对于 60%60\%60% 的数据,1<n≤20001 < n \leq 20001<n2000
  • 对于 100%100\%100% 的数据,1<n≤2×1051 < n \leq 2\times 10^51<n2×1050<Wi≤100000 < W_i \leq 100000<Wi10000

保证一定存在可产生联合权值的有序点对。

solution

对于一个节点 u ,它的所有相邻顶点之间的距离和为 2,找到其中最大的两个 w[i] 和 w[j],w[i] * w[j] 为局部最大值。对于 要算 ∑i,j,i≠j(w[i]∗w[j])\sum_{i,j,i \neq j }(w[i] * w[j])i,j,i=j(w[i]w[j]) ,由于点太多,不能两两计算,而是求它们的和 S 与平方和 T。则
S2−T=∑i,j,i≠j(w[i]∗w[j])S^2 - T = \sum_{i,j,i \neq j }(w[i] * w[j]) S2T=i,j,i=j(w[i]w[j])

代码

#include <iostream>
#include "bit"
#include "vector"
#include "unordered_set"
#include "set"
#include "queue"
#include "algorithm"
#include "bitset"
#include "cstring"

using namespace std;

/*
 * 贪心:在连接前面端点的情况下,下一站尽可能远
 */


const int N = 2e5 + 1, M = 10007;

int n, a[N];
vector<int> e[N];

int main() {
    cin >> n;
    for (int i = 1; i < n; i++) {
        int x, y;
        scanf("%d %d", &x, &y);
        e[x].push_back(y);
        e[y].push_back(x);
    }

    for (int i = 1; i <= n; i++) cin >> a[i];

    int sum = 0, Max = 0;
    for (int i = 1; i <= n; i++) {
        int s = 0, t = 0, _1st = 0, _2nd = 0;
        for (int u: e[i]) {
            s = (s + a[u]) % M;
            if(a[u] > _1st) _2nd = _1st, _1st = a[u];
            else if(a[u] > _2nd) _2nd = a[u];
            t = (t + a[u] * a[u]) % M;
        }
        Max = max(Max, _1st * _2nd);
        s = (s * s - t + M) % M;
        sum = (sum + s) % M;
    }

    cout << Max << " " << sum;
}

结果

在这里插入图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

智趣代码实验室

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值