传送门:CF
题目描述:
题目较长,此处暂略
输入:
5
0 2 3 4 5
2 0 4 5 6
3 4 0 6 7
4 5 6 0 8
5 6 7 8 0
输出:
26
4 5
3 4
2 3
2 1
赛场上被较复杂的题面给哄住了(主要是英文题很难快速读懂),导致没有下定决心好好去做
对于此题,我们最终得出的结论是求一颗最大生成树,接下来分析一下得出这个结论的过程
对于 1 1 1~ n n n的n个节点,我们先假设最开始合并的三个点为 i , k , j i,k,j i,k,j,并且递增.所以对于这三个点,我们如果想合并成一个点,可以先 i 与 k 再 { i , k } 与 j i与k再\{i,k\}与j i与k再{i,k}与j 或者 i 与 j 再 { i , j } 与 k i与j再\{i,j\}与k i与j再{i,j}与k 或者 k 与 j 再 { j , k } 与 i k与j再\{j,k\}与i k与j再{j,k}与i 一共有三种合并方式,然后我们就会发现不论是那种合并方式,在这几个点的合并过程中,我们只用到了 P [ i ] [ j ] , P [ i ] [ k ] , P [ k ] [ j ] 以及他们对应的相等数 P[i][j],P[i][k],P[k][j]以及他们对应的相等数 P[i][j],P[i][k],P[k][j]以及他们对应的相等数,并且不论是那种方式,我们这三个点最终会合并成一个大点,这个大点的所有与其他点的 p o p u l a r i t y popularity popularity值等与这三个点取一个 m a x max max(也就是说只要我们三个最终是合并的,那么这三个合并的过程与我们的大点没有关系,最终的大点是一个确定值).那么对于上述的合并来说,我们想取一个最大值,我们会发现这个最大值其实就是上述的三个值中最大的两个.为什么呢,下面仔细证明一下
为了方便起见,我们直接设三个值为
a
,
b
,
c
a,b,c
a,b,c分别代表<A,B>,<B,C>,<A,C>
那么我们有值1=a+max{b,c} 因为假设我们合并了A,B,那么此时对于我们的C来说,我们此时合并的贡献应该是max{AC,BC}也就是{C,B}
类似的有值2=b+max{a,c},值3=c+max{a,b},并且我们最终的答案应该为:
max{值1,值2,值3},不难发现,这个式子等价于
a
,
b
,
c
a,b,c
a,b,c中取最大的两个值
类似的,当我们有四个点时,显然是取与这四个点有关的几个值中最大的三个(可以自己手画一下,结合最大生成树,很好理解)
当我们有n个点时,取与n个点有关的几个值中最大的n-1个
其实,在我们进行合并的过程中, 我们不难发现这就是将所有点变成一个联通分量的过程,并且合并的花费相当于两者的P值,并且对于外部的一个点想和这个分量联通时,我们取得是这个连通分量中的所有点中的最优点(这个最优点恰好等价于我们的P值取max合并过程).所以这就是一个求最大生成树的一道题目
下面是具体的代码部分:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define root 1,n,1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
inline ll read() {
ll x=0,w=1;char ch=getchar();
for(;ch>'9'||ch<'0';ch=getchar()) if(ch=='-') w=-1;
for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
return x*w;
}
#define maxn 1000000
const double eps=1e-8;
#define int_INF 0x3f3f3f3f
#define ll_INF 0x3f3f3f3f3f3f3f3f
int n;
struct Edge{
int u,v,w;
}edge[maxn];
int mp[1010][1010];
struct P{
int a,b;
};
vector<P>p;
bool cmp(Edge a,Edge b) {
return a.w>b.w;
}
int fa[maxn];
int find(int x) {
if(x==fa[x]) return x;
else return fa[x]=find(fa[x]);
}
int main() {
n=read();int cnt=0;
for(int i=1;i<=n;i++) fa[i]=i;
for(int i=1;i<=n;i++) {
for(int j=1;j<=n;j++) {
int w=read();
edge[++cnt]={i,j,w};
}
}
int ans=0;
sort(edge+1,edge+cnt+1,cmp);
for(int i=1;i<=cnt;i++) {
int u=edge[i].u,v=edge[i].v,w=edge[i].w;
if(find(u)!=find(v)) {
ans+=w;
p.push_back({find(u),find(v)});
fa[max(find(u),find(v))]=min(find(u),find(v));
}
}
cout<<ans<<endl;
for(int i=0;i<p.size();i++) {
printf("%d %d\n",p[i].a,p[i].b);
}
return 0;
}