D2. Mocha and Diana (Hard Version)

至于D1,由于范围是1000,我们直接枚举所有的边,看看能不能加上去就行,复杂度是\(O(n^2logn)\).至于\(n\)到了\(1e5\),就要重新考虑解法。

考虑到树的边数是\(n-1\)。也就是说我们枚举的大多数边都是无效的,这个时候就要考虑我们能不能只连有效的边,对于两个连通块而言,显然在这两个连通块之间我们只需要连一条边就可以了。还有一般这种两个图的题,一般都是在一个图中操作,在另一个图中进行合不合法的判定。这里我们以第一个图作为操作的对象,以第二个图作为判定的依据。首先想到的就是我们先将两个图都划分连通块,然后枚举第一个图任意两个连通块x,y,考虑两个连通块的点,若存在u属于x,v属于y并且u和v在图二中不属于同一个连通块,那么u和v就可以相连,这两个连通块就可以联通。并且两个连通块不能连的充要条件是这两个连通块内的所有点在图二中都属于同一个连通块内。枚举任意两个连通块显然复杂度是过不去的,我们能不能固定一个连通块(比如说连通块1),然后尝试将其他连通块与它联通,但这样会不会少连一些边?比如连通块x与连通块y都不能与1联通,但这两个能不能联通呢?考虑下刚才提到的充要条件,这两个连通块不能和1联通,说明这两个连通块内的所有点和1的所有点在图二中都属于同一个连通块内。那这两个连通块的点在图二也都属于同一个连通块内啊!肯定不能连啊。好了,那我们固定一号连通块,尝试将其他的所有连通块与它联通。显然我们找点的时候我们肯定不能枚举所有的点来判断,既然要求找到一个点对即可。我们考虑保存下当前一号连通块内的所有点在图二中的连通块编号,以及该编号中的其中一个点(方便输出方案)。当下一个连通块来的时候我们直接枚举这个连通块内的所有点来判断是否有没有符合要求的点对。并且当可以联通时,我们还需要在图二中连边,合并两个连通块,这个可以用并查集实现,但有时这两个连通块(符合要求的点对在图二中的两个连通块)都在我们维护的一号连通块内(图一中),这就需要我们维护的东西支持删除操作,这个我用了set,不想再动脑了。你以为这就结束了吗?不不不....我们发现当一号连通块内对应的在图二中的连通块个数越多,越容易与其它连通块联通,为了避免一些恶心情况(因为当前1号连通块对应的图二中的连通块个数为1,而失败联通,但之后有其他的连通块对应图二有更多的连通块,调整下顺序就能联通原本不能联通的),我们将图一中的连通块按照在图二中对应地连通块的个数排序,个数多着在前,这样我们就先尽可能的增加我们维护的一号图在二号图中的连通块的个数,之后就能尽可能多的接纳其他的连通块。总复杂度\(O(nlogn)\),多少次想放弃,但还是想坚持下去....

查看代码

  1. //不等,不问,不犹豫,不回头.
  2. #include
  3. #define _ 0
  4. #define ls p<<1
  5. #define db double
  6. #define rs p<<1|1
  7. #define P 1000000007
  8. #define ll long long
  9. #define INF 1000000000
  10. #define get(x) x=read()
  11. #define PLI pair
  12. #define PII pair
  13. #define ull unsigned long long
  14. #define put(x) printf("%d\n",x)
  15. #define putl(x) printf("%lld\n",x)
  16. #define rep(x,y,z) for(int x=y;x<=z;++x)
  17. #define fep(x,y,z) for(int x=y;x>=z;--x)
  18. #define go(x) for(int i=link[x],y=a[i].y;i;y=a[i=a[i].next].y)
  19. using namespace std;
  20. const int N=1e5+10;
  21. int n,m1,m2,f[3][N],vis[N],num,bx[N],by[N],id[N];
  22. struct wy{int len,v;}a[N];
  23. vectorv[N];
  24. mapmp;
  25. sets;
  26. inline int read()

  27. {

  28. int x=0,ff=1;

  29. char ch=getchar();

  30. while(!isdigit(ch)) {if(ch=='-') ff=-1;ch=getchar();}

  31. while(isdigit(ch)) {x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}

  32. return x*ff;

  33. }
  34. inline int getf(int op,int x){return f[op][x]==x?x:f[op][x]=getf(op,f[op][x]);}
  35. inline void init()

  36. {

  37. get(n);get(m1);get(m2);

  38. rep(i,1,n) f[1][i]=f[2][i]=i;

  39. rep(i,1,m1)

  40. {

  41. int get(x),get(y);

  42. int t1=getf(1,x),t2=getf(1,y);

  43. if(t1!=t2) f[1][t1]=t2;

  44. }

  45. rep(i,1,m2)

  46. {

  47. int get(x),get(y);

  48. int t1=getf(2,x),t2=getf(2,y);

  49. if(t1!=t2) f[2][t1]=t2;

  50. }

  51. }
  52. inline bool cmp(wy a,wy b) {return a.len>b.len;}
  53. inline void prework()

  54. {

  55. rep(i,1,n) f[1][i]=getf(1,i),f[2][i]=getf(2,i);

  56. rep(i,1,n)

  57. {

  58. if(!vis[f[1][i]])

  59. {

  60. vis[f[1][i]]=1;

  61. id[f[1][i]]=++num;

  62. }

  63. a[id[f[1][i]]].len++;

  64. a[id[f[1][i]]].v=f[1][i];

  65. }

  66. sort(a+1,a+num+1,cmp);

  67. rep(i,1,num) id[a[i].v]=i;

  68. rep(i,1,n) v[id[f[1][i]]].push_back(i);

  69. memset(vis,0,sizeof(vis));

  70. for(auto x:v[1])

  71. {

  72. if(vis[f[2][x]]) continue;

  73. vis[f[2][x]]=1;

  74. s.insert(f[2][x]);

  75. mp[f[2][x]]=x;

  76. }

  77. }
  78. inline void solve()

  79. {

  80. int ans=0;

  81. rep(i,2,num)//枚举从第二个连通块开始的所有连通块,尝试与第一个连通块合并。

  82. {

  83. bool flag=false;

  84. for(auto x:s)//枚举每个连通块

  85. {

  86. if(flag) break;

  87. for(auto y:v[i])

  88. {

  89. int t=getf(2,y);

  90. if(t!=x&&!flag)

  91. {

  92. bx[++ans]=mp[x];

  93. by[ans]=y;

  94. f[2][t]=x;

  95. if(s.find(t)!=s.end()) s.erase(t);

  96. flag=true;

  97. }

  98. else if(flag&&!vis[t])

  99. {

  100. vis[t]=1;

  101. s.insert(t);

  102. mp[t]=y;

  103. }

  104. }

  105. }

  106. }

  107. put(ans);

  108. rep(i,1,ans) printf("%d %d\n",bx[i],by[i]);

  109. }
  110. int main()

  111. {

  112. // freopen("1.in","r",stdin);

  113. //freopen("sol.out","w",stdout);

  114. init();

  115. prework();

  116. solve();

  117. return (0_0);

  118. }

  119. //以吾之血,铸吾最后的亡魂.
  120.  
  121.  

Codeforces Round #738 (Div. 2) D2题解的更多相关文章

  1. # Codeforces Round #529(Div.3)个人题解

    Codeforces Round #529(Div.3)个人题解 前言: 闲来无事补了前天的cf,想着最近刷题有点点怠惰,就直接一场cf一场cf的刷算了,以后的题解也都会以每场的形式写出来 A. Re ...

  2. Codeforces Round #557 (Div. 1) 简要题解

    Codeforces Round #557 (Div. 1) 简要题解 codeforces A. Hide and Seek 枚举起始位置\(a\),如果\(a\)未在序列中出现,则对答案有\(2\ ...

  3. Codeforces Round #738 (Div. 2)

    Codeforces Round #738 (Div. 2) 跳转链接 A. Mocha and Math 题目大意 有一个长度为\(n\)的数组 可以进行无数次下面的操作,问操作后数组中的最大值的最 ...

  4. Codeforces Round #540 (Div. 3) 部分题解

    Codeforces Round #540 (Div. 3) 题目链接:https://codeforces.com/contest/1118 题目太多啦,解释题意都花很多时间...还有事情要做,就选 ...

  5. Codeforces Round #527 (Div. 3) ABCDEF题解

    Codeforces Round #527 (Div. 3) 题解 题目总链接:https://codeforces.com/contest/1092 A. Uniform String 题意: 输入 ...

  6. Codeforces Round #538 (Div. 2) (A-E题解)

    Codeforces Round #538 (Div. 2) 题目链接:https://codeforces.com/contest/1114 A. Got Any Grapes? 题意: 有三个人, ...

  7. Codeforces Round #531 (Div. 3) ABCDEF题解

    Codeforces Round #531 (Div. 3) 题目总链接:https://codeforces.com/contest/1102 A. Integer Sequence Dividin ...

  8. Codeforces Round #499 (Div. 1)部分题解(B,C,D)

    Codeforces Round #499 (Div. 1) 这场本来想和同学一起打\(\rm virtual\ contest\)的,结果有事耽搁了,之后又陆陆续续写了些,就综合起来发一篇题解. B ...

  9. Codeforces Round #350 (Div. 2) D2. Magic Powder - 2

    题目链接: http://codeforces.com/contest/670/problem/D2 题解: 二分答案. #include<iostream> #include<cs ...

随机推荐

  1. 一文看懂String类中的常用方法

    1.int length(): 返回字符串的长度: return value.length 2.char charAt(int index): 返回某索引处的字符return value[index] ...

  2. Jmeter扩展组件开发(8) - 函数助手扩展开发demo

    前提条件 1.pom文件引用ApacheJMeter_functions包 <dependency> <groupId>org.apache.jmeter</groupI ...

  3. git 报错 gitThere is no tracking information for the current branch. Please specify which branch you w

    新建本地分支后将本地分支推送到远程库, 使用git pull 或者 git push 的时候报错gitThere is no tracking information for the current ...

  4. promise对象总结

    一.Promise是异步编程的一种解决方案,它是一个容器,里面保存着某个未来才会结束的事件(通常是一个异步操作)的结果.从语法上说,Promise 是一个对象,从它可以获取异步操作的消息.Promis ...

  5. First Linux Centos 7.2 rpm 安装MySQL 5.7

    服务器需要换python环境,手贱重装了,今天凑巧需要测试数据库,花了一个小时搞了一下MySQL安装. 1.删除原有Mariadb 说明:目前centos默认的MySQL是Mariadb,由于习惯了M ...

  6. 大型项目源码集合「GitHub 热点速览 v.21.39」

    作者:HelloGitHub-小鱼干 代码,尤其是优雅规范的代码,一直都是学习编程技巧的捷径.虽然有实用的代码小片段,能拯救当前业务的燃眉之急,但是真要去提升自己的技能还是得从大型的项目,尤其是有一定 ...

  7. CSS 小技巧 | 一行代码实现头像与国旗的融合

    到国庆了,大家都急着给祖国母亲庆生. 每年每到此时,微信朋友圈就会流行起给头像装饰上国旗,而今年又流行这款: emm,很不错. 那么,将一张国旗图片与我们的头像,快速得到想要的头像,使用 CSS 如何 ...

  8. JAVA 150道笔试题知识点整理

    JAVA 笔试题 整理了几天才整理的题目,都是在笔试或者面试碰到的,好好理解消化下,对你会有帮助,祝你找工作顺利,收到满意的 offer . 1.Java 基础知识 1.1 Java SE 语法 &a ...

  9. Spring源码阅读一

    引导: 众所周知,阅读spring源码最开始的就是去了解spring bean的生命周期:bean的生命周期是怎么样的呢,见图知意: 大致流程: 首先后通过BeanDefinitionReader读取 ...

  10. Go语言核心36讲(Go语言基础知识一)--学习笔记

    01 | 工作区和GOPATH 从 Go 1.5 版本的自举(即用 Go 语言编写程序来实现 Go 语言自身),到 Go 1.7 版本的极速 GC(也称垃圾回收器),再到 2018 年 2 月发布的 ...