经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » C++ » 查看文章
BZOJ4602: [Sdoi2016]齿轮(并查集 启发式合并)
来源:cnblogs  作者:自为风月马前卒  时间:2018/12/7 9:33:31  对本文有异议

题意

题目链接

Sol

和cc的一道题很像啊

对于初始的\(N\)个点,每加一条限制实际上就是合并了两个联通块。

那么我们预处理出\(val[i]\)表示的是\(i\)节点所在的联通块根节点转了\(1\)圈,该节点会转多少圈

并查集维护联通性以及联通块大小,直接启发式合并就可以了

跑的好像还挺快

  1. #include<bits/stdc++.h>
  2. #define siz(v) ((int)v.size())
  3. using namespace std;
  4. const int MAXN = 2001;
  5. const double eps = 1e-9;
  6. inline int read() {
  7. char c = getchar(); int x = 0, f = 1;
  8. while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
  9. while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
  10. return x * f;
  11. }
  12. int N, M, fa[MAXN], siz[MAXN];
  13. double val[MAXN];
  14. vector<int> v[MAXN];
  15. void init() {
  16. for(int i = 1; i <= N; i++) fa[i] = i, v[i].clear(), siz[i] = 1, val[i] = 1;
  17. }
  18. int find(int x) {
  19. return fa[x] == x ? fa[x] : fa[x] = find(fa[x]);
  20. }
  21. bool dcmp(double x) {
  22. return fabs(x) < eps;
  23. }
  24. void mem(int x, int fa, double va) {
  25. for(int i = 0; i < siz(v[x]); i++) {
  26. int to = v[x][i];
  27. if(to == fa) continue;
  28. val[to] *= va;
  29. mem(to, x, va);
  30. }
  31. }
  32. bool solve() {
  33. N = read(); M = read(); bool flag = 0;
  34. init();
  35. for(int i = 1; i <= M; i++) {
  36. int x = read(), y = read(), fx = find(x), fy = find(y); double vx = read(), vy = read();
  37. if(flag) continue;
  38. if(fx == fy) {if(!(dcmp(val[x] / val[y] - vx / vy))) flag = 1; continue;}
  39. if(siz[fx] > siz[fy]) swap(x, y), swap(vx, vy), swap(fx, fy);
  40. val[fx] = vx / vy * val[y] / val[x];
  41. mem(fx, 0, val[fx]); fa[fx] = fy; v[y].push_back(fx);
  42. siz[fy] += siz[fx]; siz[fx] = 0;
  43. }
  44. return flag ^ 1;
  45. }
  46. int main() {
  47. #ifndef ONLINE_JUDGE
  48. freopen("a.in", "r", stdin);
  49. #endif
  50. int T = read();
  51. for(int i = 1; i <= T; i++) {
  52. printf("Case #%d: ", i);
  53. puts(solve() ? "Yes" : "No");
  54. }
  55. return 0;
  56. }
  57. /*
  58. 2
  59. 3 3
  60. 1 2 3 5
  61. 2 3 5 -7
  62. 1 3 3 7
  63. 3 3
  64. 1 2 3 5
  65. 2 3 5 -7
  66. 1 3 3 -7
  67. */
 友情链接:直通硅谷  点职佳  北美留学生论坛

本站QQ群:前端 618073944 | Java 606181507 | Python 626812652 | C/C++ 612253063 | 微信 634508462 | 苹果 692586424 | C#/.net 182808419 | PHP 305140648 | 运维 608723728

W3xue 的所有内容仅供测试,对任何法律问题及风险不承担任何责任。通过使用本站内容随之而来的风险与本站无关。
关于我们  |  意见建议  |  捐助我们  |  报错有奖  |  广告合作、友情链接(目前9元/月)请联系QQ:27243702 沸活量
皖ICP备17017327号-2 皖公网安备34020702000426号