Hex Artifact Content
Not logged in

Artifact 541b463b7581621e32a29657bc67241305661ae2:


0000: 23 69 6e 63 6c 75 64 65 20 3c 69 6f 73 74 72 65  #include <iostre
0010: 61 6d 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 73  am>..#include <s
0020: 73 74 72 65 61 6d 3e 0d 0a 23 69 6e 63 6c 75 64  stream>..#includ
0030: 65 20 3c 69 6f 6d 61 6e 69 70 3e 0d 0a 23 69 6e  e <iomanip>..#in
0040: 63 6c 75 64 65 20 3c 76 65 63 74 6f 72 3e 0d 0a  clude <vector>..
0050: 23 69 6e 63 6c 75 64 65 20 3c 73 74 72 69 6e 67  #include <string
0060: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 6d 61 70  >..#include <map
0070: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 73 65 74  >..#include <set
0080: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 61 6c 67  >..#include <alg
0090: 6f 72 69 74 68 6d 3e 0d 0a 23 69 6e 63 6c 75 64  orithm>..#includ
00a0: 65 20 3c 6e 75 6d 65 72 69 63 3e 0d 0a 23 69 6e  e <numeric>..#in
00b0: 63 6c 75 64 65 20 3c 69 74 65 72 61 74 6f 72 3e  clude <iterator>
00c0: 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 66 75 6e 63  ..#include <func
00d0: 74 69 6f 6e 61 6c 3e 0d 0a 23 69 6e 63 6c 75 64  tional>..#includ
00e0: 65 20 3c 63 6f 6d 70 6c 65 78 3e 0d 0a 23 69 6e  e <complex>..#in
00f0: 63 6c 75 64 65 20 3c 71 75 65 75 65 3e 0d 0a 23  clude <queue>..#
0100: 69 6e 63 6c 75 64 65 20 3c 73 74 61 63 6b 3e 0d  include <stack>.
0110: 0a 23 69 6e 63 6c 75 64 65 20 3c 63 6d 61 74 68  .#include <cmath
0120: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 61 73  >..#include <cas
0130: 73 65 72 74 3e 0d 0a 23 69 6e 63 6c 75 64 65 20  sert>..#include 
0140: 3c 74 75 70 6c 65 3e 0d 0a 75 73 69 6e 67 20 6e  <tuple>..using n
0150: 61 6d 65 73 70 61 63 65 20 73 74 64 3b 0d 0a 74  amespace std;..t
0160: 79 70 65 64 65 66 20 6c 6f 6e 67 20 6c 6f 6e 67  ypedef long long
0170: 20 4c 4c 3b 0d 0a 74 79 70 65 64 65 66 20 63 6f   LL;..typedef co
0180: 6d 70 6c 65 78 3c 64 6f 75 62 6c 65 3e 20 43 4d  mplex<double> CM
0190: 50 3b 0d 0a 0d 0a 63 6c 61 73 73 20 4d 61 6b 69  P;....class Maki
01a0: 6e 67 52 65 67 75 6c 61 72 47 72 61 70 68 20 7b  ngRegularGraph {
01b0: 20 70 75 62 6c 69 63 3a 0d 0a 09 76 65 63 74 6f   public:...vecto
01c0: 72 20 3c 69 6e 74 3e 20 61 64 64 45 64 67 65 73  r <int> addEdges
01d0: 28 69 6e 74 20 6e 2c 20 76 65 63 74 6f 72 20 3c  (int n, vector <
01e0: 69 6e 74 3e 20 78 2c 20 76 65 63 74 6f 72 20 3c  int> x, vector <
01f0: 69 6e 74 3e 20 79 29 0d 0a 09 7b 0d 0a 09 09 76  int> y)...{....v
0200: 65 63 74 6f 72 3c 69 6e 74 3e 20 64 65 67 28 6e  ector<int> deg(n
0210: 29 3b 0d 0a 09 09 66 6f 72 20 28 69 6e 74 20 69  );....for (int i
0220: 20 3d 20 30 3b 20 69 20 3c 20 78 2e 73 69 7a 65   = 0; i < x.size
0230: 28 29 3b 20 2b 2b 69 29 0d 0a 09 09 09 64 65 67  (); ++i).....deg
0240: 5b 78 5b 69 5d 5d 2b 2b 2c 20 64 65 67 5b 79 5b  [x[i]]++, deg[y[
0250: 69 5d 5d 2b 2b 3b 0d 0a 0d 0a 09 09 73 65 74 3c  i]]++;......set<
0260: 69 6e 74 3e 20 64 30 2c 20 64 31 3b 0d 0a 09 09  int> d0, d1;....
0270: 66 6f 72 20 28 69 6e 74 20 76 20 3d 20 30 3b 20  for (int v = 0; 
0280: 76 20 3c 20 6e 3b 20 2b 2b 76 29 0d 0a 09 09 09  v < n; ++v).....
0290: 69 66 20 28 64 65 67 5b 76 5d 20 3d 3d 20 30 29  if (deg[v] == 0)
02a0: 20 64 30 2e 69 6e 73 65 72 74 28 76 29 3b 0d 0a   d0.insert(v);..
02b0: 09 09 09 65 6c 73 65 20 69 66 20 28 64 65 67 5b  ...else if (deg[
02c0: 76 5d 20 3d 3d 20 31 29 20 64 31 2e 69 6e 73 65  v] == 1) d1.inse
02d0: 72 74 28 76 29 3b 0d 0a 09 09 6d 61 70 3c 69 6e  rt(v);....map<in
02e0: 74 2c 20 69 6e 74 3e 20 64 31 65 3b 0d 0a 09 09  t, int> d1e;....
02f0: 66 6f 72 20 28 69 6e 74 20 69 20 3d 20 30 3b 20  for (int i = 0; 
0300: 69 20 3c 20 78 2e 73 69 7a 65 28 29 3b 20 2b 2b  i < x.size(); ++
0310: 69 29 0d 0a 09 09 09 69 66 20 28 64 31 2e 63 6f  i).....if (d1.co
0320: 75 6e 74 28 78 5b 69 5d 29 20 26 26 20 64 31 2e  unt(x[i]) && d1.
0330: 63 6f 75 6e 74 28 79 5b 69 5d 29 29 20 7b 0d 0a  count(y[i])) {..
0340: 09 09 09 09 64 31 65 5b 78 5b 69 5d 5d 20 3d 20  ....d1e[x[i]] = 
0350: 79 5b 69 5d 3b 0d 0a 09 09 09 09 64 31 65 5b 79  y[i];......d1e[y
0360: 5b 69 5d 5d 20 3d 20 78 5b 69 5d 3b 0d 0a 09 09  [i]] = x[i];....
0370: 09 09 64 31 2e 65 72 61 73 65 28 78 5b 69 5d 29  ..d1.erase(x[i])
0380: 3b 0d 0a 09 09 09 09 64 31 2e 65 72 61 73 65 28  ;......d1.erase(
0390: 79 5b 69 5d 29 3b 0d 0a 09 09 09 7d 0d 0a 0d 0a  y[i]);.....}....
03a0: 09 09 23 64 65 66 69 6e 65 20 42 41 44 20 28 28  ..#define BAD ((
03b0: 28 64 30 2e 73 69 7a 65 28 29 3d 3d 31 7c 7c 64  (d0.size()==1||d
03c0: 30 2e 73 69 7a 65 28 29 3d 3d 32 29 20 26 26 20  0.size()==2) && 
03d0: 64 31 2e 65 6d 70 74 79 28 29 20 26 26 20 64 31  d1.empty() && d1
03e0: 65 2e 65 6d 70 74 79 28 29 29 20 7c 7c 20 28 64  e.empty()) || (d
03f0: 30 2e 65 6d 70 74 79 28 29 20 26 26 20 64 31 2e  0.empty() && d1.
0400: 65 6d 70 74 79 28 29 20 26 26 20 64 31 65 2e 73  empty() && d1e.s
0410: 69 7a 65 28 29 3d 3d 32 20 26 26 20 64 31 65 2e  ize()==2 && d1e.
0420: 62 65 67 69 6e 28 29 2d 3e 66 69 72 73 74 3d 3d  begin()->first==
0430: 64 31 65 2e 72 62 65 67 69 6e 28 29 2d 3e 73 65  d1e.rbegin()->se
0440: 63 6f 6e 64 29 29 0d 0a 09 09 69 66 28 42 41 44  cond))....if(BAD
0450: 29 0d 0a 09 09 09 72 65 74 75 72 6e 20 76 65 63  ).....return vec
0460: 74 6f 72 3c 69 6e 74 3e 28 31 2c 20 2d 31 29 3b  tor<int>(1, -1);
0470: 0d 0a 0d 0a 09 09 76 65 63 74 6f 72 3c 62 6f 6f  ......vector<boo
0480: 6c 3e 20 75 73 65 64 28 6e 2a 6e 29 3b 0d 0a 09  l> used(n*n);...
0490: 09 66 6f 72 20 28 69 6e 74 20 69 20 3d 20 30 3b  .for (int i = 0;
04a0: 20 69 20 3c 20 78 2e 73 69 7a 65 28 29 3b 20 2b   i < x.size(); +
04b0: 2b 69 29 0d 0a 09 09 09 75 73 65 64 5b 78 5b 69  +i).....used[x[i
04c0: 5d 20 2a 20 6e 20 2b 20 79 5b 69 5d 5d 20 3d 20  ] * n + y[i]] = 
04d0: 74 72 75 65 3b 0d 0a 0d 0a 09 09 76 65 63 74 6f  true;......vecto
04e0: 72 3c 69 6e 74 3e 20 61 6e 73 3b 0d 0a 09 09 66  r<int> ans;....f
04f0: 6f 72 28 69 6e 74 20 76 3d 30 3b 20 76 3c 6e 3b  or(int v=0; v<n;
0500: 20 2b 2b 76 29 0d 0a 09 09 66 6f 72 28 69 6e 74   ++v)....for(int
0510: 20 75 3d 76 2b 31 3b 20 75 3c 6e 3b 20 2b 2b 75   u=v+1; u<n; ++u
0520: 29 20 69 66 28 21 75 73 65 64 5b 76 2a 6e 2b 75  ) if(!used[v*n+u
0530: 5d 29 0d 0a 09 09 7b 0d 0a 09 09 09 69 66 20 28  ])....{.....if (
0540: 64 30 2e 65 6d 70 74 79 28 29 20 26 26 20 64 31  d0.empty() && d1
0550: 2e 65 6d 70 74 79 28 29 20 26 26 20 64 31 65 2e  .empty() && d1e.
0560: 65 6d 70 74 79 28 29 29 0d 0a 09 09 09 09 62 72  empty())......br
0570: 65 61 6b 3b 0d 0a 0d 0a 09 09 09 62 6f 6f 6c 20  eak;.......bool 
0580: 62 61 64 20 3d 20 74 72 75 65 3b 0d 0a 09 09 09  bad = true;.....
0590: 69 66 20 28 64 30 2e 63 6f 75 6e 74 28 76 29 29  if (d0.count(v))
05a0: 20 7b 0d 0a 09 09 09 09 64 30 2e 65 72 61 73 65   {......d0.erase
05b0: 28 76 29 3b 0d 0a 09 09 09 09 69 66 20 28 64 30  (v);......if (d0
05c0: 2e 63 6f 75 6e 74 28 75 29 29 20 7b 0d 0a 09 09  .count(u)) {....
05d0: 09 09 09 64 30 2e 65 72 61 73 65 28 75 29 3b 0d  ...d0.erase(u);.
05e0: 0a 09 09 09 09 09 64 31 65 5b 76 5d 20 3d 20 75  ......d1e[v] = u
05f0: 3b 0d 0a 09 09 09 09 09 64 31 65 5b 75 5d 20 3d  ;.......d1e[u] =
0600: 20 76 3b 0d 0a 09 09 09 09 09 69 66 20 28 62 61   v;.......if (ba
0610: 64 3d 42 41 44 29 20 7b 0d 0a 09 09 09 09 09 09  d=BAD) {........
0620: 64 30 2e 69 6e 73 65 72 74 28 75 29 3b 0d 0a 09  d0.insert(u);...
0630: 09 09 09 09 09 64 31 65 2e 65 72 61 73 65 28 76  .....d1e.erase(v
0640: 29 3b 0d 0a 09 09 09 09 09 09 64 31 65 2e 65 72  );........d1e.er
0650: 61 73 65 28 75 29 3b 0d 0a 09 09 09 09 09 7d 0d  ase(u);.......}.
0660: 0a 09 09 09 09 7d 0d 0a 09 09 09 09 65 6c 73 65  .....}......else
0670: 20 69 66 20 28 64 31 2e 63 6f 75 6e 74 28 75 29   if (d1.count(u)
0680: 29 20 7b 0d 0a 09 09 09 09 09 64 31 2e 65 72 61  ) {.......d1.era
0690: 73 65 28 75 29 3b 0d 0a 09 09 09 09 09 64 31 2e  se(u);.......d1.
06a0: 69 6e 73 65 72 74 28 76 29 3b 0d 0a 09 09 09 09  insert(v);......
06b0: 09 69 66 20 28 62 61 64 20 3d 20 42 41 44 29 20  .if (bad = BAD) 
06c0: 7b 0d 0a 09 09 09 09 09 09 64 31 2e 69 6e 73 65  {........d1.inse
06d0: 72 74 28 75 29 3b 0d 0a 09 09 09 09 09 09 64 31  rt(u);........d1
06e0: 2e 65 72 61 73 65 28 76 29 3b 0d 0a 09 09 09 09  .erase(v);......
06f0: 09 7d 0d 0a 09 09 09 09 7d 0d 0a 09 09 09 09 65  .}......}......e
0700: 6c 73 65 20 69 66 20 28 64 31 65 2e 63 6f 75 6e  lse if (d1e.coun
0710: 74 28 75 29 29 20 7b 0d 0a 09 09 09 09 09 69 6e  t(u)) {.......in
0720: 74 20 75 74 20 3d 20 64 31 65 5b 75 5d 3b 20 64  t ut = d1e[u]; d
0730: 31 65 2e 65 72 61 73 65 28 75 29 3b 0d 0a 09 09  1e.erase(u);....
0740: 09 09 09 64 31 2e 69 6e 73 65 72 74 28 76 29 3b  ...d1.insert(v);
0750: 0d 0a 09 09 09 09 09 69 66 20 28 62 61 64 20 3d  .......if (bad =
0760: 20 42 41 44 29 20 7b 0d 0a 09 09 09 09 09 09 64   BAD) {........d
0770: 31 65 5b 75 5d 20 3d 20 75 74 3b 0d 0a 09 09 09  1e[u] = ut;.....
0780: 09 09 09 64 31 2e 65 72 61 73 65 28 76 29 3b 0d  ...d1.erase(v);.
0790: 0a 09 09 09 09 09 7d 0d 0a 09 09 09 09 7d 0d 0a  ......}......}..
07a0: 09 09 09 09 69 66 20 28 62 61 64 29 0d 0a 09 09  ....if (bad)....
07b0: 09 09 09 64 30 2e 69 6e 73 65 72 74 28 76 29 3b  ...d0.insert(v);
07c0: 0d 0a 09 09 09 7d 0d 0a 09 09 09 65 6c 73 65 20  .....}.....else 
07d0: 69 66 20 28 64 31 2e 63 6f 75 6e 74 28 76 29 29  if (d1.count(v))
07e0: 20 7b 0d 0a 09 09 09 09 64 31 2e 65 72 61 73 65   {......d1.erase
07f0: 28 76 29 3b 0d 0a 09 09 09 09 69 66 20 28 64 30  (v);......if (d0
0800: 2e 63 6f 75 6e 74 28 75 29 29 20 7b 0d 0a 09 09  .count(u)) {....
0810: 09 09 09 64 30 2e 65 72 61 73 65 28 75 29 3b 0d  ...d0.erase(u);.
0820: 0a 09 09 09 09 09 64 31 2e 69 6e 73 65 72 74 28  ......d1.insert(
0830: 75 29 3b 0d 0a 09 09 09 09 09 69 66 20 28 62 61  u);.......if (ba
0840: 64 20 3d 20 42 41 44 29 20 7b 0d 0a 09 09 09 09  d = BAD) {......
0850: 09 09 64 30 2e 69 6e 73 65 72 74 28 75 29 3b 0d  ..d0.insert(u);.
0860: 0a 09 09 09 09 09 09 64 31 2e 65 72 61 73 65 28  .......d1.erase(
0870: 75 29 3b 0d 0a 09 09 09 09 09 7d 0d 0a 09 09 09  u);.......}.....
0880: 09 7d 0d 0a 09 09 09 09 65 6c 73 65 20 69 66 20  .}......else if 
0890: 28 64 31 2e 63 6f 75 6e 74 28 75 29 29 20 7b 0d  (d1.count(u)) {.
08a0: 0a 09 09 09 09 09 64 31 2e 65 72 61 73 65 28 75  ......d1.erase(u
08b0: 29 3b 0d 0a 09 09 09 09 09 69 66 20 28 62 61 64  );.......if (bad
08c0: 20 3d 20 42 41 44 29 20 7b 0d 0a 09 09 09 09 09   = BAD) {.......
08d0: 09 64 31 2e 69 6e 73 65 72 74 28 75 29 3b 0d 0a  .d1.insert(u);..
08e0: 09 09 09 09 09 7d 0d 0a 09 09 09 09 7d 0d 0a 09  .....}......}...
08f0: 09 09 09 65 6c 73 65 20 69 66 20 28 64 31 65 2e  ...else if (d1e.
0900: 63 6f 75 6e 74 28 75 29 29 20 7b 0d 0a 09 09 09  count(u)) {.....
0910: 09 09 69 6e 74 20 75 74 20 3d 20 64 31 65 5b 75  ..int ut = d1e[u
0920: 5d 3b 20 64 31 65 2e 65 72 61 73 65 28 75 29 3b  ]; d1e.erase(u);
0930: 0d 0a 09 09 09 09 09 69 66 20 28 62 61 64 20 3d  .......if (bad =
0940: 20 42 41 44 29 20 7b 0d 0a 09 09 09 09 09 09 64   BAD) {........d
0950: 31 65 5b 75 5d 20 3d 20 75 74 3b 0d 0a 09 09 09  1e[u] = ut;.....
0960: 09 09 7d 0d 0a 09 09 09 09 7d 0d 0a 09 09 09 09  ..}......}......
0970: 69 66 20 28 62 61 64 29 0d 0a 09 09 09 09 09 64  if (bad).......d
0980: 31 2e 69 6e 73 65 72 74 28 76 29 3b 0d 0a 09 09  1.insert(v);....
0990: 09 7d 0d 0a 09 09 09 65 6c 73 65 20 69 66 20 28  .}.....else if (
09a0: 64 31 65 2e 63 6f 75 6e 74 28 76 29 29 20 7b 0d  d1e.count(v)) {.
09b0: 0a 09 09 09 09 69 6e 74 20 76 74 20 3d 20 64 31  .....int vt = d1
09c0: 65 5b 76 5d 3b 20 20 64 31 65 2e 65 72 61 73 65  e[v];  d1e.erase
09d0: 28 76 29 3b 0d 0a 09 09 09 09 69 66 20 28 64 30  (v);......if (d0
09e0: 2e 63 6f 75 6e 74 28 75 29 29 20 7b 0d 0a 09 09  .count(u)) {....
09f0: 09 09 09 64 30 2e 65 72 61 73 65 28 75 29 3b 0d  ...d0.erase(u);.
0a00: 0a 09 09 09 09 09 64 31 2e 69 6e 73 65 72 74 28  ......d1.insert(
0a10: 75 29 3b 0d 0a 09 09 09 09 09 69 66 20 28 62 61  u);.......if (ba
0a20: 64 20 3d 20 42 41 44 29 20 7b 0d 0a 09 09 09 09  d = BAD) {......
0a30: 09 09 64 30 2e 69 6e 73 65 72 74 28 75 29 3b 0d  ..d0.insert(u);.
0a40: 0a 09 09 09 09 09 09 64 31 2e 65 72 61 73 65 28  .......d1.erase(
0a50: 75 29 3b 0d 0a 09 09 09 09 09 7d 0d 0a 09 09 09  u);.......}.....
0a60: 09 7d 0d 0a 09 09 09 09 65 6c 73 65 20 69 66 20  .}......else if 
0a70: 28 64 31 2e 63 6f 75 6e 74 28 75 29 29 20 7b 0d  (d1.count(u)) {.
0a80: 0a 09 09 09 09 09 64 31 2e 65 72 61 73 65 28 75  ......d1.erase(u
0a90: 29 3b 0d 0a 09 09 09 09 09 69 66 20 28 62 61 64  );.......if (bad
0aa0: 20 3d 20 42 41 44 29 20 7b 0d 0a 09 09 09 09 09   = BAD) {.......
0ab0: 09 64 31 2e 69 6e 73 65 72 74 28 75 29 3b 0d 0a  .d1.insert(u);..
0ac0: 09 09 09 09 09 7d 0d 0a 09 09 09 09 7d 0d 0a 09  .....}......}...
0ad0: 09 09 09 65 6c 73 65 20 69 66 20 28 64 31 65 2e  ...else if (d1e.
0ae0: 63 6f 75 6e 74 28 75 29 29 20 7b 0d 0a 09 09 09  count(u)) {.....
0af0: 09 09 69 6e 74 20 75 74 20 3d 20 64 31 65 5b 75  ..int ut = d1e[u
0b00: 5d 3b 20 64 31 65 2e 65 72 61 73 65 28 75 29 3b  ]; d1e.erase(u);
0b10: 0d 0a 09 09 09 09 09 69 66 20 28 62 61 64 20 3d  .......if (bad =
0b20: 20 42 41 44 29 20 7b 0d 0a 09 09 09 09 09 09 64   BAD) {........d
0b30: 31 65 5b 75 5d 20 3d 20 75 74 3b 0d 0a 09 09 09  1e[u] = ut;.....
0b40: 09 09 7d 0d 0a 09 09 09 09 7d 0d 0a 09 09 09 09  ..}......}......
0b50: 69 66 20 28 62 61 64 29 0d 0a 09 09 09 09 09 64  if (bad).......d
0b60: 31 65 5b 76 5d 20 3d 20 76 74 3b 0d 0a 09 09 09  1e[v] = vt;.....
0b70: 7d 0d 0a 0d 0a 09 09 09 69 66 28 21 62 61 64 29  }.......if(!bad)
0b80: 0d 0a 09 09 09 09 61 6e 73 2e 70 75 73 68 5f 62  ......ans.push_b
0b90: 61 63 6b 28 76 29 2c 20 61 6e 73 2e 70 75 73 68  ack(v), ans.push
0ba0: 5f 62 61 63 6b 28 75 29 3b 0d 0a 09 09 7d 0d 0a  _back(u);....}..
0bb0: 0d 0a 09 09 72 65 74 75 72 6e 20 61 6e 73 3b 0d  ....return ans;.
0bc0: 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45  ..}..};....// BE
0bd0: 47 49 4e 20 43 55 54 20 48 45 52 45 0d 0a 23 69  GIN CUT HERE..#i
0be0: 6e 63 6c 75 64 65 20 3c 63 74 69 6d 65 3e 0d 0a  nclude <ctime>..
0bf0: 64 6f 75 62 6c 65 20 73 74 61 72 74 5f 74 69 6d  double start_tim
0c00: 65 3b 20 73 74 72 69 6e 67 20 74 69 6d 65 72 28  e; string timer(
0c10: 29 0d 0a 20 7b 20 6f 73 74 72 69 6e 67 73 74 72  ).. { ostringstr
0c20: 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c 20 22 20  eam os; os << " 
0c30: 28 22 20 3c 3c 20 69 6e 74 28 28 63 6c 6f 63 6b  (" << int((clock
0c40: 28 29 2d 73 74 61 72 74 5f 74 69 6d 65 29 2f 43  ()-start_time)/C
0c50: 4c 4f 43 4b 53 5f 50 45 52 5f 53 45 43 2a 31 30  LOCKS_PER_SEC*10
0c60: 30 30 29 20 3c 3c 20 22 20 6d 73 65 63 29 22 3b  00) << " msec)";
0c70: 20 72 65 74 75 72 6e 20 6f 73 2e 73 74 72 28 29   return os.str()
0c80: 3b 20 7d 0d 0a 74 65 6d 70 6c 61 74 65 3c 74 79  ; }..template<ty
0c90: 70 65 6e 61 6d 65 20 54 3e 20 6f 73 74 72 65 61  pename T> ostrea
0ca0: 6d 26 20 6f 70 65 72 61 74 6f 72 3c 3c 28 6f 73  m& operator<<(os
0cb0: 74 72 65 61 6d 26 20 6f 73 2c 20 63 6f 6e 73 74  tream& os, const
0cc0: 20 76 65 63 74 6f 72 3c 54 3e 26 20 76 29 0d 0a   vector<T>& v)..
0cd0: 20 7b 20 6f 73 20 3c 3c 20 22 7b 20 22 3b 0d 0a   { os << "{ ";..
0ce0: 20 20 20 66 6f 72 28 74 79 70 65 6e 61 6d 65 20     for(typename 
0cf0: 76 65 63 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73 74  vector<T>::const
0d00: 5f 69 74 65 72 61 74 6f 72 20 69 74 3d 76 2e 62  _iterator it=v.b
0d10: 65 67 69 6e 28 29 3b 20 69 74 21 3d 76 2e 65 6e  egin(); it!=v.en
0d20: 64 28 29 3b 20 2b 2b 69 74 29 0d 0a 20 20 20 6f  d(); ++it)..   o
0d30: 73 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 2a 69 74  s << '\"' << *it
0d40: 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 28 69 74 2b   << '\"' << (it+
0d50: 31 3d 3d 76 2e 65 6e 64 28 29 20 3f 20 22 22 20  1==v.end() ? "" 
0d60: 3a 20 22 2c 20 22 29 3b 20 6f 73 20 3c 3c 20 22  : ", "); os << "
0d70: 20 7d 22 3b 20 72 65 74 75 72 6e 20 6f 73 3b 20   }"; return os; 
0d80: 7d 0d 0a 76 6f 69 64 20 76 65 72 69 66 79 5f 63  }..void verify_c
0d90: 61 73 65 28 63 6f 6e 73 74 20 76 65 63 74 6f 72  ase(const vector
0da0: 20 3c 69 6e 74 3e 26 20 45 78 70 65 63 74 65 64   <int>& Expected
0db0: 2c 20 63 6f 6e 73 74 20 76 65 63 74 6f 72 20 3c  , const vector <
0dc0: 69 6e 74 3e 26 20 52 65 63 65 69 76 65 64 29 20  int>& Received) 
0dd0: 7b 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 45  {.. bool ok = (E
0de0: 78 70 65 63 74 65 64 20 3d 3d 20 52 65 63 65 69  xpected == Recei
0df0: 76 65 64 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63  ved);.. if(ok) c
0e00: 65 72 72 20 3c 3c 20 22 50 41 53 53 45 44 22 20  err << "PASSED" 
0e10: 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e  << timer() << en
0e20: 64 6c 3b 20 20 65 6c 73 65 20 7b 20 63 65 72 72  dl;  else { cerr
0e30: 20 3c 3c 20 22 46 41 49 4c 45 44 22 20 3c 3c 20   << "FAILED" << 
0e40: 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b  timer() << endl;
0e50: 0d 0a 20 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a  .. cerr << "\to:
0e60: 20 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20 3c   " << Expected <
0e70: 3c 20 65 6e 64 6c 20 3c 3c 20 22 5c 74 78 3a 20  < endl << "\tx: 
0e80: 22 20 3c 3c 20 52 65 63 65 69 76 65 64 20 3c 3c  " << Received <<
0e90: 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 64 65 66   endl; } }..#def
0ea0: 69 6e 65 20 43 41 53 45 28 4e 29 20 7b 63 65 72  ine CASE(N) {cer
0eb0: 72 20 3c 3c 20 22 54 65 73 74 20 43 61 73 65 20  r << "Test Case 
0ec0: 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e 22  #" << N << "..."
0ed0: 20 3c 3c 20 66 6c 75 73 68 3b 20 73 74 61 72 74   << flush; start
0ee0: 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 0d 0a  _time=clock();..
0ef0: 23 64 65 66 69 6e 65 20 45 4e 44 09 20 76 65 72  #define END. ver
0f00: 69 66 79 5f 63 61 73 65 28 5f 2c 20 4d 61 6b 69  ify_case(_, Maki
0f10: 6e 67 52 65 67 75 6c 61 72 47 72 61 70 68 28 29  ngRegularGraph()
0f20: 2e 61 64 64 45 64 67 65 73 28 6e 2c 20 78 2c 20  .addEdges(n, x, 
0f30: 79 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 6e 28  y));}..int main(
0f40: 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 29 0d 0a 09  ){....CASE(0)...
0f50: 69 6e 74 20 6e 20 3d 20 36 3b 20 0d 0a 09 69 6e  int n = 6; ...in
0f60: 74 20 78 5f 5b 5d 20 3d 20 7b 30 2c 32 7d 3b 0d  t x_[] = {0,2};.
0f70: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e  ..  vector <int>
0f80: 20 78 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66   x(x_, x_+sizeof
0f90: 28 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29  (x_)/sizeof(*x_)
0fa0: 29 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d  ); ...int y_[] =
0fb0: 20 7b 32 2c 33 7d 3b 0d 0a 09 20 20 76 65 63 74   {2,3};...  vect
0fc0: 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79  or <int> y(y_, y
0fd0: 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a  _+sizeof(y_)/siz
0fe0: 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 69 6e  eof(*y_)); ...in
0ff0: 74 20 5f 5f 5b 5d 20 3d 20 7b 30 2c 20 31 2c 20  t __[] = {0, 1, 
1000: 31 2c 20 34 2c 20 33 2c 20 35 2c 20 34 2c 20 35  1, 4, 3, 5, 4, 5
1010: 20 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c   };...  vector <
1020: 69 6e 74 3e 20 5f 28 5f 5f 2c 20 5f 5f 2b 73 69  int> _(__, __+si
1030: 7a 65 6f 66 28 5f 5f 29 2f 73 69 7a 65 6f 66 28  zeof(__)/sizeof(
1040: 2a 5f 5f 29 29 3b 20 0d 0a 45 4e 44 0d 0a 43 41  *__)); ..END..CA
1050: 53 45 28 31 29 0d 0a 09 69 6e 74 20 6e 20 3d 20  SE(1)...int n = 
1060: 34 3b 20 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d  4; ...int x_[] =
1070: 20 7b 31 2c 32 2c 31 7d 3b 0d 0a 09 20 20 76 65   {1,2,1};...  ve
1080: 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c  ctor <int> x(x_,
1090: 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73   x_+sizeof(x_)/s
10a0: 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09  izeof(*x_)); ...
10b0: 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 32 2c 33 2c  int y_[] = {2,3,
10c0: 33 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  3};...  vector <
10d0: 69 6e 74 3e 20 79 28 79 5f 2c 20 79 5f 2b 73 69  int> y(y_, y_+si
10e0: 7a 65 6f 66 28 79 5f 29 2f 73 69 7a 65 6f 66 28  zeof(y_)/sizeof(
10f0: 2a 79 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 5f 5f  *y_)); ...int __
1100: 5b 5d 20 3d 20 7b 2d 31 20 7d 3b 0d 0a 09 20 20  [] = {-1 };...  
1110: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 5f 28 5f  vector <int> _(_
1120: 5f 2c 20 5f 5f 2b 73 69 7a 65 6f 66 28 5f 5f 29  _, __+sizeof(__)
1130: 2f 73 69 7a 65 6f 66 28 2a 5f 5f 29 29 3b 20 0d  /sizeof(*__)); .
1140: 0a 45 4e 44 0d 0a 43 41 53 45 28 32 29 0d 0a 09  .END..CASE(2)...
1150: 69 6e 74 20 6e 20 3d 20 33 3b 20 0d 0a 09 76 65  int n = 3; ...ve
1160: 63 74 6f 72 20 3c 69 6e 74 3e 20 78 3b 20 0d 0a  ctor <int> x; ..
1170: 09 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 3b  .vector <int> y;
1180: 20 0d 0a 09 69 6e 74 20 5f 5f 5b 5d 20 3d 20 7b   ...int __[] = {
1190: 30 2c 20 31 2c 20 30 2c 20 32 2c 20 31 2c 20 32  0, 1, 0, 2, 1, 2
11a0: 20 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c   };...  vector <
11b0: 69 6e 74 3e 20 5f 28 5f 5f 2c 20 5f 5f 2b 73 69  int> _(__, __+si
11c0: 7a 65 6f 66 28 5f 5f 29 2f 73 69 7a 65 6f 66 28  zeof(__)/sizeof(
11d0: 2a 5f 5f 29 29 3b 20 0d 0a 45 4e 44 0d 0a 43 41  *__)); ..END..CA
11e0: 53 45 28 33 29 0d 0a 09 69 6e 74 20 6e 20 3d 20  SE(3)...int n = 
11f0: 35 3b 20 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d  5; ...int x_[] =
1200: 20 7b 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72   {0};...  vector
1210: 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 78 5f 2b   <int> x(x_, x_+
1220: 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69 7a 65 6f  sizeof(x_)/sizeo
1230: 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20  f(*x_)); ...int 
1240: 79 5f 5b 5d 20 3d 20 7b 34 7d 3b 0d 0a 09 20 20  y_[] = {4};...  
1250: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79  vector <int> y(y
1260: 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29  _, y_+sizeof(y_)
1270: 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d  /sizeof(*y_)); .
1280: 0a 09 69 6e 74 20 5f 5f 5b 5d 20 3d 20 7b 30 2c  ..int __[] = {0,
1290: 20 31 2c 20 31 2c 20 32 2c 20 32 2c 20 33 2c 20   1, 1, 2, 2, 3, 
12a0: 33 2c 20 34 20 7d 3b 0d 0a 09 20 20 76 65 63 74  3, 4 };...  vect
12b0: 6f 72 20 3c 69 6e 74 3e 20 5f 28 5f 5f 2c 20 5f  or <int> _(__, _
12c0: 5f 2b 73 69 7a 65 6f 66 28 5f 5f 29 2f 73 69 7a  _+sizeof(__)/siz
12d0: 65 6f 66 28 2a 5f 5f 29 29 3b 20 0d 0a 45 4e 44  eof(*__)); ..END
12e0: 0d 0a 43 41 53 45 28 34 29 0d 0a 09 69 6e 74 20  ..CASE(4)...int 
12f0: 6e 20 3d 20 35 3b 20 0d 0a 09 69 6e 74 20 78 5f  n = 5; ...int x_
1300: 5b 5d 20 3d 20 7b 32 7d 3b 0d 0a 09 20 20 76 65  [] = {2};...  ve
1310: 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c  ctor <int> x(x_,
1320: 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73   x_+sizeof(x_)/s
1330: 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09  izeof(*x_)); ...
1340: 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 33 7d 3b 0d  int y_[] = {3};.
1350: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e  ..  vector <int>
1360: 20 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66   y(y_, y_+sizeof
1370: 28 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29  (y_)/sizeof(*y_)
1380: 29 3b 20 0d 0a 09 69 6e 74 20 5f 5f 5b 5d 20 3d  ); ...int __[] =
1390: 20 7b 30 2c 20 31 2c 20 30 2c 20 32 2c 20 31 2c   {0, 1, 0, 2, 1,
13a0: 20 34 2c 20 33 2c 20 34 20 7d 3b 0d 0a 09 20 20   4, 3, 4 };...  
13b0: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 5f 28 5f  vector <int> _(_
13c0: 5f 2c 20 5f 5f 2b 73 69 7a 65 6f 66 28 5f 5f 29  _, __+sizeof(__)
13d0: 2f 73 69 7a 65 6f 66 28 2a 5f 5f 29 29 3b 20 0d  /sizeof(*__)); .
13e0: 0a 45 4e 44 0d 0a 2f 2a 0d 0a 43 41 53 45 28 35  .END../*..CASE(5
13f0: 29 0d 0a 09 69 6e 74 20 6e 20 3d 20 3b 20 0d 0a  )...int n = ; ..
1400: 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 3b 0d 0a 09  .int x_[] = ;...
1410: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78    vector <int> x
1420: 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78  (x_, x_+sizeof(x
1430: 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b  _)/sizeof(*x_));
1440: 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 3b   ...int y_[] = ;
1450: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
1460: 3e 20 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f  > y(y_, y_+sizeo
1470: 66 28 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f  f(y_)/sizeof(*y_
1480: 29 29 3b 20 0d 0a 09 69 6e 74 20 5f 5f 5b 5d 20  )); ...int __[] 
1490: 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  = ;...  vector <
14a0: 69 6e 74 3e 20 5f 28 5f 5f 2c 20 5f 5f 2b 73 69  int> _(__, __+si
14b0: 7a 65 6f 66 28 5f 5f 29 2f 73 69 7a 65 6f 66 28  zeof(__)/sizeof(
14c0: 2a 5f 5f 29 29 3b 20 0d 0a 45 4e 44 0d 0a 43 41  *__)); ..END..CA
14d0: 53 45 28 36 29 0d 0a 09 69 6e 74 20 6e 20 3d 20  SE(6)...int n = 
14e0: 3b 20 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20  ; ...int x_[] = 
14f0: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e  ;...  vector <in
1500: 74 3e 20 78 28 78 5f 2c 20 78 5f 2b 73 69 7a 65  t> x(x_, x_+size
1510: 6f 66 28 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 78  of(x_)/sizeof(*x
1520: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d  _)); ...int y_[]
1530: 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20   = ;...  vector 
1540: 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79 5f 2b 73  <int> y(y_, y_+s
1550: 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a 65 6f 66  izeof(y_)/sizeof
1560: 28 2a 79 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 5f  (*y_)); ...int _
1570: 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74  _[] = ;...  vect
1580: 6f 72 20 3c 69 6e 74 3e 20 5f 28 5f 5f 2c 20 5f  or <int> _(__, _
1590: 5f 2b 73 69 7a 65 6f 66 28 5f 5f 29 2f 73 69 7a  _+sizeof(__)/siz
15a0: 65 6f 66 28 2a 5f 5f 29 29 3b 20 0d 0a 45 4e 44  eof(*__)); ..END
15b0: 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f 20 45 4e 44 20  ..*/..}..// END 
15c0: 43 55 54 20 48 45 52 45 0d 0a                    CUT HERE..