经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » C » 查看文章
POJ1659 Frogs' Neighborhood(Havel–Hakimi定理)
来源:cnblogs  作者:自为风月马前卒  时间:2018/10/15 9:27:06  对本文有异议

题意

题目链接

\(T\)组数据,给出\(n\)个点的度数,问是否可以构造出一个简单图

Sol

Havel–Hakimi定理:

  • 给定一串有限多个非负整数组成的序列,是否存在一个简单图使得其度数列恰为这个序列。

\(S=(d_1,d_2,\dots,d_n)\)为有限多个非负整数组成的非递增序列。 S可简单图化当且仅当有穷序列\(S’=(d_2-1,d_3-1,...,d(d_1+1)-1,d(d_1+2),...,d_n)\)只含有非负整数且是可简单图化的。

最后判断一下是否都是零就好了

感觉这个算法。。就是个贪心吧。。

当然判断这类问题的可行性还有另外一种方法:Erd?s–Gallai定理

\(S=(d_1,d_2,...,d_n)\)为有限多个非负整数组成的非递增序列。\(S\)可简单图化当且仅当这些数字的和为偶数,并且

\(\sum_{i = 1}^k d_i \leqslant k(k - 1) + \sum_{i = k + 1}^n min(d_i, k)\)

对所有\(1 \leqslant k \leqslant n\)都成立

不过这个好像没办法输出方案??。。。

  1. #include<cstdio>
  2. #include<iostream>
  3. #include<cstring>
  4. #include<algorithm>
  5. #define Pair pair<int, int>
  6. #define MP(x, y) make_pair(x, y)
  7. #define fi first
  8. #define se second
  9. using namespace std;
  10. const int MAXN = 1e5 + 10, INF = 1e9 + 7;
  11. inline int read() {
  12. char c = getchar(); int x = 0, f = 1;
  13. while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
  14. while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
  15. return x * f;
  16. }
  17. int T, N, reach[101][101], sum = 0;
  18. Pair a[MAXN];
  19. void init() {
  20. memset(reach, 0, sizeof(reach));
  21. sum = 0;
  22. }
  23. int main() {
  24. // freopen("a.in", "r", stdin);
  25. T = read();
  26. while(T--) {
  27. init();
  28. N = read();
  29. for(int i = 1; i <= N; i++) a[i] = MP(read(), i), sum += a[i].fi;
  30. if(sum % 2 != 0) {puts("NO\n"); continue;}
  31. bool f = 0;
  32. for(int i = 1; i <= N; i++) {
  33. sort(a + i, a + N + 1, greater<Pair>());
  34. if(a[i].fi <= 0) continue;
  35. for(int j = i + 1; j <= i + a[i].fi; j++) a[j].fi -= 1, reach[a[i].se][a[j].se] = 1, reach[a[j].se][a[i].se] = 1;
  36. a[i].fi = 0;
  37. }
  38. for(int i = 1; i <= N; i++) if(a[i].fi != 0) {puts("NO\n"); f = 1; break;}
  39. if(f) continue;
  40. puts("YES");
  41. for(int i = 1; i <= N; i++, puts(""))
  42. for(int j = 1; j <= N; j++)
  43. printf("%d ", reach[i][j]);
  44. puts("");
  45. }
  46. }
  47. /*
  48. 1
  49. 6
  50. 4 3 1 4 2 0
  51. */
 友情链接:直通硅谷  点职佳  北美留学生论坛

本站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号