洛谷 P1351 联合权值

洛谷 P1351 联合权值

---

题目

题目描述

无向连通图G 有n 个点，n - 1 条边。点从1 到n 依次编号，编号为 i 的点的权值为W i ，每条边的长度均为1 。图上两点( u , v ) 的距离定义为u 点到v 点的最短距离。对于图G 上的点对( u, v) ，若它们的距离为2 ，则它们之间会产生Wu×Wv 的联合权值。

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

输入输出格式

输入格式：
输入文件名为link .in。

第一行包含1 个整数n 。

接下来n - 1 行，每行包含 2 个用空格隔开的正整数u 、v ，表示编号为 u 和编号为v 的点之间有边相连。

最后1 行，包含 n 个正整数，每两个正整数之间用一个空格隔开，其中第 i 个整数表示图G 上编号为i 的点的权值为W i 。

输出格式：
输出文件名为link .out 。

输出共1 行，包含2 个整数，之间用一个空格隔开，依次为图G 上联合权值的最大值

和所有联合权值之和。由于所有联合权值之和可能很大，输出它时要对10007 取余。

输入输出样例

输入样例#1：
5
1 2
2 3
3 4
4 5
1 5 2 3 10
输出样例#1：
20 74
说明

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

其联合权值分别为2 、15、2 、20、15、20。其中最大的是20，总和为74。

【数据说明】

对于30% 的数据，1 < n≤ 100 ；

对于60% 的数据，1 < n≤ 2000；

对于100%的数据，1 < n≤ 200 , 000 ，0 < wi≤ 10, 000 。

---

题解

因为题目说只有距离为2的两点才有联合权值，那么我们可以知道一个点所连接的所有点之间都有联合权值，那么我们就可以把每一条边建成两条单向边，然后进行排序在逐个做每个点所连接的所有点互相之间的联合权值即可

优化：求所有的乘积和可以用前缀和来优化

注意，取模是若是负数要将它加上tt后再取模

---

代码（Pascal）

const tt=10007;
var n:longint;
    maxp,ans:int64;
    x,y:array[0..400005]of longint;
    p,q:array[0..200005]of longint;
 procedure init;
  var i:longint;
   begin
     readln(n);
     for i:=1 to n-1 do
      begin
        readln(x[i*2-1],y[i*2-1]);
        x[i*2]:=y[i*2-1];
        y[i*2]:=x[i*2-1];
       end;
     for i:=1 to n do
      read(p[i]);
    end;
 procedure qsort(l,r:longint);
  var i,j,mid,t:longint;
   begin
     i:=l;j:=r;
     mid:=x[random(r-l+1)+l];
     repeat
       while x[i]<mid do inc(i);
       while x[j]>mid do dec(j);
       if i<=j then begin
                      t:=x[i];x[i]:=x[j];x[j]:=t;
                      t:=y[i];y[i]:=y[j];y[j]:=t;
                      inc(i);dec(j);
                     end;
      until i>j;
     if i<r then qsort(i,r);
     if l<j then qsort(l,j);
    end;
 function max(x,y:int64):int64;
  begin
    if x>y then exit(x) else exit(y);
   end;
 procedure main;
  var i,j,sum,tot,max1,max2:longint;
   begin
     qsort(1,2*(n-1));
     sum:=0;tot:=0;ans:=0;maxp:=0;max1:=0;max2:=0;
     for i:=1 to 2*(n-1) do
      if x[i]=x[i-1] then
       begin
         inc(tot);
         q[tot]:=p[y[i]];
         sum:=(sum+q[tot])mod tt;
         if q[tot]>max1 then begin max2:=max1;max1:=q[tot];end else
         if q[tot]>max2 then max2:=q[tot];
        end
      else begin
             for j:=1 to tot do
              begin
                ans:=(ans+(((sum-q[j]) mod tt)*(q[j] mod tt)) mod tt) 

mod tt;
               end;
             maxp:=max(max1*max2,maxp);
             tot:=1;
             q[tot]:=p[y[i]];
             sum:=q[tot];
             max1:=p[y[i]];max2:=0;
            end;
    end;
 procedure print;
  begin
    write(maxp,' ',(ans+tt)mod tt);
   end;
 begin
   init;
   main;
   print;
  end.