D2. Mocha and Diana (Hard Version)

参考博客了博客 同时也参考了cf评论中humbertoyusta的题解

然后详细的对它的最后一步进行了一个证明

先说下题意:
给你两个结点数相同的2个森林它有n个点,你每次同时给两个森林的同一对结点连边,使得这两个森林最后还是两个森林(森林:1~n棵树组成),问你最多加哪几条边。

解析下:
你能加的边数 = min(n - 1 - m1, n - 1 - m2)
因为加到最后两个森林的中至少会有一个变成一棵树,这棵树结点数是n,那么它的边数就是 n - 1。你每次选择哪一个结点去加边实际上不影响你最多可以加几条边,因为极限是其中1个森林变成一棵树也就是只有一个联通块,所以不管你怎么加边,只要两个森林都是有两个或两个以上联通块你总能找出一个连边的方法(可以自己去证明下) 从而达到只有一个联通块的状态(只有一棵树),而这个过程中森林的边数的变化还是从max(m1,m2)到 n - 1,也就是min(n - 1 - m1, n - 1 - m2);所以不论过程如何,最后可以加的最多的边数是不变的。

既然如此不妨先尽可能的给一个点加边,从而减小联通块数量

在保证森林1和森林2中的所有树还是树的前提下。你可以试着对每一个结点x建立一个(1,x)的边。

假设这是原图
在这里插入图片描述

进行上述操作后就是这样了
在这里插入图片描述



之后你会发现对结点1已经没有可以操作的余地了,你只能考虑除了结点1的其它联通块了。然后你就会发现这两个森林除1号结点联通块外的其它联通块之间都是互不联通的,比如上图森林1的5号结点和森林2的,2,3结点在各自森林里不连通,6号结点也和森林2的2,3结点在各自森林里不连通。

所以直接从两个森林除了1号结点的联通块的其它联通块中各自选出一个两两建边连接起来就行了

那么以上发现具有一般性吗?证明如下:

定义2个集合分别为s1,s2分别表示森林1和森林2中除结点编号1所在的树以外其它树的根结点编号,也就是各自森林除了结点1外的联通块集合。

如上图,对应的s1,s2分别为下图所示
在这里插入图片描述

那么它有如下两个性质
下面证明中树就是联通块

①这两个集合内不存在相同结点。反证: 如果有,因为s1,s2集合中的点和结点1不联通,所以这两个相同结点x的树一定可以各自建立一条边(1,x)与1号结点所对应的树联通,那么根据对这两个集合的定义,既然结点x和结点1联通了,就自然不会出现在集合s1,s2中。

如下图
在这里插入图片描述

②对于s2的任意一个根结点b对应在森林1的b结点一定在1号树上。反证:如果不在,那么就说明s1的树a上有一个b号结点它与森林1的1号树不联通,森林2中的b也和森林2的1号树不联通,那么他们一定可以通过之前的操作和建立边(b,1)从而与1号结点所对应的树联通,既然结点b和结点1联通了,就自然不会出现在集合s1,s2中。

如下图
在这里插入图片描述



再陈述下上面得两个性质
①s1中结点编号和s2中得结点编号一定都不相同。②集合s2中的任何一个结点一定不在集合s1中的任何一个结点对应的树上,而是在对应森林里的1号结点的树上,反之亦然。也就是集合s1,和s2上的任意两点在各自森林都不联通。所以可以大胆的对这两个点加边,并且使得结点与各自森林的1号树联通。



所以你每次可以在集合s1,s2中各自取出一个结点,放心的对森林1和森林2各自将这两点连接。

如下图
在这里插入图片描述

时间复杂度分析:
①每个森林对结点1的连边操作(1,x) ,用并查集维护联通块 O(nlogn)
②求出除了1号结点所在联通块的其它联通块O(n)

代码:

#include<bits/stdc++.h>
using namespace std;
const int N = 1e5 + 10;
int ans[N][2];
int fa1[N];
int fa2[N];
int s1[N],s2[N];

int find(int fa[],int x){
    return x == fa[x] ? x : fa[x] = find(fa,fa[x]);
}

void merge(int fa[],int x,int y){
    int fx = find(fa,x);
    int fy = find(fa,y);
    fa[fx] = fy;
}

void init(int fa[],int n){
    for(int i = 1 ;i <= n; ++i)
        fa[i] = i;
}

int main(){
   /*  freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);   */
    int n,m1,m2;    scanf("%d %d %d",&n,&m1,&m2);
    init(fa1,n);
    init(fa2,n);
    for(int i = 1; i <= m1; ++i){
        int u,v;    scanf("%d %d",&u,&v);
        merge(fa1,u,v);
    }
    for(int i = 1; i <= m2; ++i){
        int u,v;    scanf("%d %d",&u,&v);
        merge(fa2,u,v);
    }
    int ct = 0;
    for(int i = 2; i <= n; ++i){
        if(find(fa1,1) != find(fa1,i) && find(fa2,1) != find(fa2,i))
        {
            ans[++ct][0] = 1;
            ans[ct][1] = i;
            merge(fa1,i,1);
            merge(fa2,i,1);
        }
    }
    int sz1 = 0,sz2 = 0;
    for(int i = 2 ; i <= n; ++i)
    {
        if(fa1[i] == i && fa1[i] != find(fa1,1))
            s1[sz1++] = i;
        if(fa2[i] == i && fa2[i] != find(fa2,1))
            s2[sz2++] = i;
    }
    
    for( int i = 0,j = 0; i < sz1 && j < sz2; ++i, ++j ){
         ans[++ct][0] = s1[i];
         ans[ct][1] = s2[j];
     }

    printf("%d\n",ct);
    for(int i = 1; i <= ct; ++i)
        printf("%d %d\n",ans[i][0],ans[i][1]);
    return 0;
}
  • 16
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值