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 73 74 72 75 63 74 20 55 6e 69 P;....struct Uni
01a0: 6f 6e 46 69 6e 64 0d 0a 7b 0d 0a 09 76 65 63 74 onFind..{...vect
01b0: 6f 72 3c 69 6e 74 3e 20 75 66 2c 20 73 7a 3b 0d or<int> uf, sz;.
01c0: 0a 09 69 6e 74 20 6e 63 3b 0d 0a 0d 0a 09 55 6e ..int nc;.....Un
01d0: 69 6f 6e 46 69 6e 64 28 69 6e 74 20 4e 29 3a 20 ionFind(int N):
01e0: 75 66 28 4e 29 2c 20 73 7a 28 4e 2c 31 29 2c 20 uf(N), sz(N,1),
01f0: 6e 63 28 4e 29 0d 0a 09 09 7b 20 66 6f 72 28 69 nc(N)....{ for(i
0200: 6e 74 20 69 3d 30 3b 20 69 3c 4e 3b 20 2b 2b 69 nt i=0; i<N; ++i
0210: 29 20 75 66 5b 69 5d 20 3d 20 69 3b 20 7d 0d 0a ) uf[i] = i; }..
0220: 09 69 6e 74 20 73 69 7a 65 28 29 0d 0a 09 09 7b .int size()....{
0230: 20 72 65 74 75 72 6e 20 6e 63 3b 20 7d 0d 0a 09 return nc; }...
0240: 69 6e 74 20 73 69 7a 65 28 69 6e 74 20 61 29 0d int size(int a).
0250: 0a 09 09 7b 20 72 65 74 75 72 6e 20 73 7a 5b 46 ...{ return sz[F
0260: 69 6e 64 28 61 29 5d 3b 20 7d 0d 0a 09 69 6e 74 ind(a)]; }...int
0270: 20 46 69 6e 64 28 69 6e 74 20 61 29 0d 0a 09 09 Find(int a)....
0280: 7b 20 72 65 74 75 72 6e 20 75 66 5b 61 5d 3d 3d { return uf[a]==
0290: 61 20 3f 20 61 20 3a 20 75 66 5b 61 5d 3d 46 69 a ? a : uf[a]=Fi
02a0: 6e 64 28 75 66 5b 61 5d 29 3b 20 7d 0d 0a 09 62 nd(uf[a]); }...b
02b0: 6f 6f 6c 20 55 6e 69 6f 6e 28 69 6e 74 20 61 2c ool Union(int a,
02c0: 20 69 6e 74 20 62 29 0d 0a 09 09 7b 0d 0a 09 09 int b)....{....
02d0: 09 61 20 3d 20 46 69 6e 64 28 61 29 3b 0d 0a 09 .a = Find(a);...
02e0: 09 09 62 20 3d 20 46 69 6e 64 28 62 29 3b 0d 0a ..b = Find(b);..
02f0: 09 09 09 69 66 28 20 61 20 21 3d 20 62 20 29 0d ...if( a != b ).
0300: 0a 09 09 09 7b 0d 0a 09 09 09 09 69 66 28 20 73 ....{......if( s
0310: 7a 5b 61 5d 20 3e 3d 20 73 7a 5b 62 5d 20 29 20 z[a] >= sz[b] )
0320: 73 77 61 70 28 61 2c 20 62 29 3b 0d 0a 09 09 09 swap(a, b);.....
0330: 09 75 66 5b 61 5d 20 3d 20 62 3b 0d 0a 09 09 09 .uf[a] = b;.....
0340: 09 73 7a 5b 62 5d 20 2b 3d 20 73 7a 5b 61 5d 3b .sz[b] += sz[a];
0350: 0d 0a 09 09 09 09 2d 2d 6e 63 3b 0d 0a 09 09 09 ......--nc;.....
0360: 7d 0d 0a 09 09 09 72 65 74 75 72 6e 20 28 61 21 }.....return (a!
0370: 3d 62 29 3b 0d 0a 09 09 7d 0d 0a 7d 3b 0d 0a 0d =b);....}..};...
0380: 0a 63 6c 61 73 73 20 45 67 61 6c 69 74 61 72 69 .class Egalitari
0390: 61 6e 69 73 6d 32 20 7b 20 70 75 62 6c 69 63 3a anism2 { public:
03a0: 0d 0a 09 69 6e 74 20 4e 3b 0d 0a 09 76 65 63 74 ...int N;...vect
03b0: 6f 72 3c 76 65 63 74 6f 72 3c 4c 4c 3e 3e 20 64 or<vector<LL>> d
03c0: 32 74 61 62 6c 65 3b 0d 0a 0d 0a 09 74 79 70 65 2table;.....type
03d0: 64 65 66 20 74 75 70 6c 65 3c 4c 4c 2c 20 69 6e def tuple<LL, in
03e0: 74 2c 20 69 6e 74 3e 20 65 64 67 65 3b 0d 0a 0d t, int> edge;...
03f0: 0a 09 64 6f 75 62 6c 65 20 6d 69 6e 53 74 64 65 ..double minStde
0400: 76 28 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 v(vector <int> x
0410: 2c 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 , vector <int> y
0420: 29 0d 0a 09 7b 0d 0a 09 09 4e 20 3d 20 78 2e 73 )...{....N = x.s
0430: 69 7a 65 28 29 3b 0d 0a 09 09 64 32 74 61 62 6c ize();....d2tabl
0440: 65 2e 72 65 73 69 7a 65 28 4e 2c 20 76 65 63 74 e.resize(N, vect
0450: 6f 72 3c 4c 4c 3e 28 4e 29 29 3b 0d 0a 09 09 66 or<LL>(N));....f
0460: 6f 72 28 69 6e 74 20 76 3d 30 3b 20 76 3c 4e 3b or(int v=0; v<N;
0470: 20 2b 2b 76 29 0d 0a 09 09 66 6f 72 28 69 6e 74 ++v)....for(int
0480: 20 75 3d 30 3b 20 75 3c 4e 3b 20 2b 2b 75 29 0d u=0; u<N; ++u).
0490: 0a 09 09 09 64 32 74 61 62 6c 65 5b 76 5d 5b 75 ....d2table[v][u
04a0: 5d 20 3d 20 28 78 5b 76 5d 2d 78 5b 75 5d 29 2a ] = (x[v]-x[u])*
04b0: 28 78 5b 76 5d 2d 78 5b 75 5d 29 20 2b 20 28 79 (x[v]-x[u]) + (y
04c0: 5b 76 5d 2d 79 5b 75 5d 29 2a 28 79 5b 76 5d 2d [v]-y[u])*(y[v]-
04d0: 79 5b 75 5d 29 3b 0d 0a 0d 0a 09 09 76 65 63 74 y[u]);......vect
04e0: 6f 72 3c 65 64 67 65 3e 20 65 64 67 65 73 3b 0d or<edge> edges;.
04f0: 0a 09 09 66 6f 72 28 69 6e 74 20 76 3d 30 3b 20 ...for(int v=0;
0500: 76 3c 4e 3b 20 2b 2b 76 29 0d 0a 09 09 66 6f 72 v<N; ++v)....for
0510: 28 69 6e 74 20 75 3d 76 2b 31 3b 20 75 3c 4e 3b (int u=v+1; u<N;
0520: 20 2b 2b 75 29 0d 0a 09 09 09 65 64 67 65 73 2e ++u).....edges.
0530: 65 6d 70 6c 61 63 65 5f 62 61 63 6b 28 64 32 74 emplace_back(d2t
0540: 61 62 6c 65 5b 76 5d 5b 75 5d 2c 20 76 2c 20 75 able[v][u], v, u
0550: 29 3b 0d 0a 09 09 73 6f 72 74 28 65 64 67 65 73 );....sort(edges
0560: 2e 62 65 67 69 6e 28 29 2c 20 65 64 67 65 73 2e .begin(), edges.
0570: 65 6e 64 28 29 29 3b 0d 0a 0d 0a 09 09 55 6e 69 end());......Uni
0580: 6f 6e 46 69 6e 64 20 75 66 28 4e 29 3b 0d 0a 09 onFind uf(N);...
0590: 09 72 65 74 75 72 6e 20 73 71 72 74 28 72 65 63 .return sqrt(rec
05a0: 28 30 2c 20 75 66 2c 20 65 64 67 65 73 2c 20 30 (0, uf, edges, 0
05b0: 2c 20 30 29 29 3b 0d 0a 09 7d 0d 0a 0d 0a 09 64 , 0));...}.....d
05c0: 6f 75 62 6c 65 20 72 65 63 28 69 6e 74 20 69 2c ouble rec(int i,
05d0: 20 55 6e 69 6f 6e 46 69 6e 64 26 20 75 66 2c 20 UnionFind& uf,
05e0: 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c 65 64 67 const vector<edg
05f0: 65 3e 26 20 65 64 67 65 73 2c 20 64 6f 75 62 6c e>& edges, doubl
0600: 65 20 73 78 2c 20 64 6f 75 62 6c 65 20 73 78 32 e sx, double sx2
0610: 29 0d 0a 09 7b 0d 0a 09 09 69 66 28 75 66 2e 73 )...{....if(uf.s
0620: 69 7a 65 28 29 20 3d 3d 20 31 29 0d 0a 09 09 09 ize() == 1).....
0630: 72 65 74 75 72 6e 20 73 78 32 2f 28 4e 2d 31 29 return sx2/(N-1)
0640: 20 2d 20 73 78 2a 73 78 2f 28 4e 2d 31 29 2f 28 - sx*sx/(N-1)/(
0650: 4e 2d 31 29 3b 0d 0a 09 09 69 66 28 69 3d 3d 65 N-1);....if(i==e
0660: 64 67 65 73 2e 73 69 7a 65 28 29 29 0d 0a 09 09 dges.size())....
0670: 09 72 65 74 75 72 6e 20 31 65 2b 33 30 30 3b 0d .return 1e+300;.
0680: 0a 0d 0a 09 09 76 65 63 74 6f 72 3c 69 6e 74 3e .....vector<int>
0690: 20 63 61 6e 64 3b 0d 0a 09 09 69 6e 74 20 65 3d cand;....int e=
06a0: 69 3b 0d 0a 09 09 66 6f 72 28 3b 20 65 3c 65 64 i;....for(; e<ed
06b0: 67 65 73 2e 73 69 7a 65 28 29 20 26 26 20 67 65 ges.size() && ge
06c0: 74 3c 30 3e 28 65 64 67 65 73 5b 69 5d 29 3d 3d t<0>(edges[i])==
06d0: 67 65 74 3c 30 3e 28 65 64 67 65 73 5b 65 5d 29 get<0>(edges[e])
06e0: 3b 20 2b 2b 65 29 0d 0a 09 09 09 69 66 28 75 66 ; ++e).....if(uf
06f0: 2e 46 69 6e 64 28 67 65 74 3c 31 3e 28 65 64 67 .Find(get<1>(edg
0700: 65 73 5b 65 5d 29 29 20 21 3d 20 75 66 2e 46 69 es[e])) != uf.Fi
0710: 6e 64 28 67 65 74 3c 32 3e 28 65 64 67 65 73 5b nd(get<2>(edges[
0720: 65 5d 29 29 29 0d 0a 09 09 09 09 63 61 6e 64 2e e])))......cand.
0730: 70 75 73 68 5f 62 61 63 6b 28 65 29 3b 0d 0a 0d push_back(e);...
0740: 0a 09 09 69 66 28 63 61 6e 64 2e 65 6d 70 74 79 ...if(cand.empty
0750: 28 29 29 0d 0a 09 09 09 72 65 74 75 72 6e 20 72 ()).....return r
0760: 65 63 28 65 2c 20 75 66 2c 20 65 64 67 65 73 2c ec(e, uf, edges,
0770: 20 73 78 2c 20 73 78 32 29 3b 0d 0a 0d 0a 09 09 sx, sx2);......
0780: 64 6f 75 62 6c 65 20 62 65 73 74 20 3d 20 31 65 double best = 1e
0790: 2b 33 30 30 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 +300;....for(int
07a0: 20 65 3a 20 63 61 6e 64 29 20 7b 0d 0a 09 09 09 e: cand) {.....
07b0: 55 6e 69 6f 6e 46 69 6e 64 20 75 66 32 20 3d 20 UnionFind uf2 =
07c0: 75 66 3b 0d 0a 09 09 09 4c 4c 20 20 64 20 3d 20 uf;.....LL d =
07d0: 67 65 74 3c 30 3e 28 65 64 67 65 73 5b 65 5d 29 get<0>(edges[e])
07e0: 3b 0d 0a 09 09 09 69 6e 74 20 75 20 3d 20 67 65 ;.....int u = ge
07f0: 74 3c 31 3e 28 65 64 67 65 73 5b 65 5d 29 3b 0d t<1>(edges[e]);.
0800: 0a 09 09 09 69 6e 74 20 76 20 3d 20 67 65 74 3c ....int v = get<
0810: 32 3e 28 65 64 67 65 73 5b 65 5d 29 3b 0d 0a 09 2>(edges[e]);...
0820: 09 09 75 66 32 2e 55 6e 69 6f 6e 28 75 2c 20 76 ..uf2.Union(u, v
0830: 29 3b 0d 0a 09 09 09 62 65 73 74 20 3d 20 6d 69 );.....best = mi
0840: 6e 28 62 65 73 74 2c 20 72 65 63 28 65 2b 31 2c n(best, rec(e+1,
0850: 20 75 66 32 2c 20 65 64 67 65 73 2c 20 73 78 2b uf2, edges, sx+
0860: 73 71 72 74 28 64 29 2c 20 73 78 32 2b 64 29 29 sqrt(d), sx2+d))
0870: 3b 0d 0a 09 09 7d 0d 0a 09 09 72 65 74 75 72 6e ;....}....return
0880: 20 62 65 73 74 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a best;...}..};..
0890: 0d 0a 2f 2f 20 42 45 47 49 4e 20 43 55 54 20 48 ..// BEGIN CUT H
08a0: 45 52 45 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 ERE..#include <c
08b0: 74 69 6d 65 3e 0d 0a 64 6f 75 62 6c 65 20 73 74 time>..double st
08c0: 61 72 74 5f 74 69 6d 65 3b 20 73 74 72 69 6e 67 art_time; string
08d0: 20 74 69 6d 65 72 28 29 0d 0a 20 7b 20 6f 73 74 timer().. { ost
08e0: 72 69 6e 67 73 74 72 65 61 6d 20 6f 73 3b 20 6f ringstream os; o
08f0: 73 20 3c 3c 20 22 20 28 22 20 3c 3c 20 69 6e 74 s << " (" << int
0900: 28 28 63 6c 6f 63 6b 28 29 2d 73 74 61 72 74 5f ((clock()-start_
0910: 74 69 6d 65 29 2f 43 4c 4f 43 4b 53 5f 50 45 52 time)/CLOCKS_PER
0920: 5f 53 45 43 2a 31 30 30 30 29 20 3c 3c 20 22 20 _SEC*1000) << "
0930: 6d 73 65 63 29 22 3b 20 72 65 74 75 72 6e 20 6f msec)"; return o
0940: 73 2e 73 74 72 28 29 3b 20 7d 0d 0a 74 65 6d 70 s.str(); }..temp
0950: 6c 61 74 65 3c 74 79 70 65 6e 61 6d 65 20 54 3e late<typename T>
0960: 20 6f 73 74 72 65 61 6d 26 20 6f 70 65 72 61 74 ostream& operat
0970: 6f 72 3c 3c 28 6f 73 74 72 65 61 6d 26 20 6f 73 or<<(ostream& os
0980: 2c 20 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c 54 , const vector<T
0990: 3e 26 20 76 29 0d 0a 20 7b 20 6f 73 20 3c 3c 20 >& v).. { os <<
09a0: 22 7b 20 22 3b 0d 0a 20 20 20 66 6f 72 28 74 79 "{ ";.. for(ty
09b0: 70 65 6e 61 6d 65 20 76 65 63 74 6f 72 3c 54 3e pename vector<T>
09c0: 3a 3a 63 6f 6e 73 74 5f 69 74 65 72 61 74 6f 72 ::const_iterator
09d0: 20 69 74 3d 76 2e 62 65 67 69 6e 28 29 3b 20 69 it=v.begin(); i
09e0: 74 21 3d 76 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 t!=v.end(); ++it
09f0: 29 0d 0a 20 20 20 6f 73 20 3c 3c 20 27 5c 22 27 ).. os << '\"'
0a00: 20 3c 3c 20 2a 69 74 20 3c 3c 20 27 5c 22 27 20 << *it << '\"'
0a10: 3c 3c 20 28 69 74 2b 31 3d 3d 76 2e 65 6e 64 28 << (it+1==v.end(
0a20: 29 20 3f 20 22 22 20 3a 20 22 2c 20 22 29 3b 20 ) ? "" : ", ");
0a30: 6f 73 20 3c 3c 20 22 20 7d 22 3b 20 72 65 74 75 os << " }"; retu
0a40: 72 6e 20 6f 73 3b 20 7d 0d 0a 76 6f 69 64 20 76 rn os; }..void v
0a50: 65 72 69 66 79 5f 63 61 73 65 28 63 6f 6e 73 74 erify_case(const
0a60: 20 64 6f 75 62 6c 65 26 20 45 78 70 65 63 74 65 double& Expecte
0a70: 64 2c 20 63 6f 6e 73 74 20 64 6f 75 62 6c 65 26 d, const double&
0a80: 20 52 65 63 65 69 76 65 64 29 20 7b 0d 0a 20 62 Received) {.. b
0a90: 6f 6f 6c 20 6f 6b 20 3d 20 28 61 62 73 28 45 78 ool ok = (abs(Ex
0aa0: 70 65 63 74 65 64 20 2d 20 52 65 63 65 69 76 65 pected - Receive
0ab0: 64 29 20 3c 20 31 65 2d 39 29 3b 0d 0a 20 69 66 d) < 1e-9);.. if
0ac0: 28 6f 6b 29 20 63 65 72 72 20 3c 3c 20 22 50 41 (ok) cerr << "PA
0ad0: 53 53 45 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 SSED" << timer()
0ae0: 20 3c 3c 20 65 6e 64 6c 3b 20 20 65 6c 73 65 20 << endl; else
0af0: 7b 20 63 65 72 72 20 3c 3c 20 22 46 41 49 4c 45 { cerr << "FAILE
0b00: 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c D" << timer() <<
0b10: 20 65 6e 64 6c 3b 0d 0a 20 63 65 72 72 20 3c 3c endl;.. cerr <<
0b20: 20 22 5c 74 6f 3a 20 5c 22 22 20 3c 3c 20 45 78 "\to: \"" << Ex
0b30: 70 65 63 74 65 64 20 3c 3c 20 27 5c 22 27 20 3c pected << '\"' <
0b40: 3c 20 65 6e 64 6c 20 3c 3c 20 22 5c 74 78 3a 20 < endl << "\tx:
0b50: 5c 22 22 20 3c 3c 20 52 65 63 65 69 76 65 64 20 \"" << Received
0b60: 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 3b << '\"' << endl;
0b70: 20 7d 20 7d 0d 0a 23 64 65 66 69 6e 65 20 43 41 } }..#define CA
0b80: 53 45 28 4e 29 20 7b 63 65 72 72 20 3c 3c 20 22 SE(N) {cerr << "
0b90: 54 65 73 74 20 43 61 73 65 20 23 22 20 3c 3c 20 Test Case #" <<
0ba0: 4e 20 3c 3c 20 22 2e 2e 2e 22 20 3c 3c 20 66 6c N << "..." << fl
0bb0: 75 73 68 3b 20 73 74 61 72 74 5f 74 69 6d 65 3d ush; start_time=
0bc0: 63 6c 6f 63 6b 28 29 3b 0d 0a 23 64 65 66 69 6e clock();..#defin
0bd0: 65 20 45 4e 44 09 20 76 65 72 69 66 79 5f 63 61 e END. verify_ca
0be0: 73 65 28 5f 2c 20 45 67 61 6c 69 74 61 72 69 61 se(_, Egalitaria
0bf0: 6e 69 73 6d 32 28 29 2e 6d 69 6e 53 74 64 65 76 nism2().minStdev
0c00: 28 78 2c 20 79 29 29 3b 7d 0d 0a 69 6e 74 20 6d (x, y));}..int m
0c10: 61 69 6e 28 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 ain(){....CASE(0
0c20: 29 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b )...int x_[] = {
0c30: 30 2c 30 2c 31 2c 31 7d 3b 0d 0a 09 20 20 76 65 0,0,1,1};... ve
0c40: 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c ctor <int> x(x_,
0c50: 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 x_+sizeof(x_)/s
0c60: 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 izeof(*x_)); ...
0c70: 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 30 2c 31 2c int y_[] = {0,1,
0c80: 30 2c 31 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 0,1};... vector
0c90: 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79 5f 2b <int> y(y_, y_+
0ca0: 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a 65 6f sizeof(y_)/sizeo
0cb0: 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 f(*y_)); ...doub
0cc0: 6c 65 20 5f 20 3d 20 30 2e 30 3b 20 0d 0a 45 4e le _ = 0.0; ..EN
0cd0: 44 0d 0a 43 41 53 45 28 31 29 0d 0a 09 69 6e 74 D..CASE(1)...int
0ce0: 20 78 5f 5b 5d 20 3d 20 7b 30 2c 30 2c 30 7d 3b x_[] = {0,0,0};
0cf0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 ... vector <int
0d00: 3e 20 78 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f > x(x_, x_+sizeo
0d10: 66 28 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 5f f(x_)/sizeof(*x_
0d20: 29 29 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 )); ...int y_[]
0d30: 3d 20 7b 30 2c 39 2c 31 30 7d 3b 0d 0a 09 20 20 = {0,9,10};...
0d40: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79 vector <int> y(y
0d50: 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 _, y_+sizeof(y_)
0d60: 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d /sizeof(*y_)); .
0d70: 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 30 2e 35 ..double _ = 0.5
0d80: 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 32 29 ; ..END..CASE(2)
0d90: 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b 31 ...int x_[] = {1
0da0: 32 2c 34 36 2c 38 31 2c 35 36 7d 3b 0d 0a 09 20 2,46,81,56};...
0db0: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 vector <int> x(
0dc0: 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f x_, x_+sizeof(x_
0dd0: 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 )/sizeof(*x_));
0de0: 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 30 ...int y_[] = {0
0df0: 2c 34 35 2c 32 2c 36 37 7d 3b 0d 0a 09 20 20 76 ,45,2,67};... v
0e00: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f ector <int> y(y_
0e10: 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f , y_+sizeof(y_)/
0e20: 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a sizeof(*y_)); ..
0e30: 09 64 6f 75 62 6c 65 20 5f 20 3d 20 36 2e 31 30 .double _ = 6.10
0e40: 32 37 39 39 39 37 31 33 32 30 38 37 32 3b 20 0d 2799971320872; .
0e50: 0a 45 4e 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09 .END..CASE(3)...
0e60: 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b 30 2c 30 2c int x_[] = {0,0,
0e70: 30 2c 30 2c 30 2c 30 2c 30 7d 3b 0d 0a 09 20 20 0,0,0,0,0};...
0e80: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 vector <int> x(x
0e90: 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 _, x_+sizeof(x_)
0ea0: 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d /sizeof(*x_)); .
0eb0: 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 30 2c ..int y_[] = {0,
0ec0: 32 2c 33 2c 39 2c 31 30 2c 31 35 2c 31 36 7d 3b 2,3,9,10,15,16};
0ed0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 ... vector <int
0ee0: 3e 20 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f > y(y_, y_+sizeo
0ef0: 66 28 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f f(y_)/sizeof(*y_
0f00: 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 )); ...double _
0f10: 3d 20 30 2e 39 34 32 38 30 39 30 34 31 35 38 32 = 0.942809041582
0f20: 30 36 31 37 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 0617; ..END..CAS
0f30: 45 28 34 29 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 E(4)...int x_[]
0f40: 3d 20 7b 31 36 37 30 35 33 2c 20 35 33 36 37 37 = {167053, 53677
0f50: 30 2c 20 2d 35 39 30 34 30 31 2c 20 35 30 37 30 0, -590401, 5070
0f60: 34 37 2c 20 33 35 30 31 37 38 2c 20 2d 32 37 34 47, 350178, -274
0f70: 35 32 33 2c 20 2d 35 38 34 36 37 39 2c 20 2d 37 523, -584679, -7
0f80: 36 36 37 39 35 2c 20 2d 36 36 34 31 37 37 2c 20 66795, -664177,
0f90: 32 36 37 37 35 37 2c 20 2d 32 39 31 38 35 36 2c 267757, -291856,
0fa0: 20 2d 37 36 35 35 34 37 2c 20 36 30 34 38 30 31 -765547, 604801
0fb0: 2c 20 2d 36 38 32 39 32 32 2c 20 2d 34 30 34 35 , -682922, -4045
0fc0: 39 30 2c 20 34 36 38 30 30 31 2c 20 36 30 37 39 90, 468001, 6079
0fd0: 32 35 2c 20 35 30 33 38 34 39 2c 20 2d 34 39 39 25, 503849, -499
0fe0: 36 39 39 2c 20 2d 37 39 38 36 33 37 7d 3b 0d 0a 699, -798637};..
0ff0: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
1000: 78 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 x(x_, x_+sizeof(
1010: 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 x_)/sizeof(*x_))
1020: 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 ; ...int y_[] =
1030: 7b 2d 31 32 33 39 36 2c 20 2d 36 36 30 39 38 2c {-12396, -66098,
1040: 20 2d 35 36 38 34 33 2c 20 32 30 32 37 30 2c 20 -56843, 20270,
1050: 38 31 35 31 30 2c 20 2d 32 33 32 39 34 2c 20 31 81510, -23294, 1
1060: 30 34 32 33 2c 20 32 34 30 30 37 2c 20 2d 32 34 0423, 24007, -24
1070: 33 34 33 2c 20 2d 32 31 35 38 37 2c 20 2d 36 33 343, -21587, -63
1080: 31 38 2c 20 2d 37 33 39 36 2c 20 2d 36 38 36 32 18, -7396, -6862
1090: 32 2c 20 35 36 33 30 34 2c 20 2d 38 35 36 38 30 2, 56304, -85680
10a0: 2c 20 2d 31 34 38 39 30 2c 20 2d 33 38 33 37 33 , -14890, -38373
10b0: 2c 20 2d 32 35 34 37 37 2c 20 2d 33 38 32 34 30 , -25477, -38240
10c0: 2c 20 31 31 37 33 36 7d 3b 0d 0a 09 20 20 76 65 , 11736};... ve
10d0: 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f 2c ctor <int> y(y_,
10e0: 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f 73 y_+sizeof(y_)/s
10f0: 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 izeof(*y_)); ...
1100: 64 6f 75 62 6c 65 20 5f 20 3d 20 34 30 30 35 36 double _ = 40056
1110: 2e 39 35 39 34 36 34 35 31 36 37 38 3b 20 0d 0a .95946451678; ..
1120: 45 4e 44 0d 0a 43 41 53 45 28 35 29 0d 0a 09 69 END..CASE(5)...i
1130: 6e 74 20 78 5f 5b 5d 20 3d 20 7b 2d 33 30 36 38 nt x_[] = {-3068
1140: 38 30 2c 20 31 36 39 34 38 30 2c 20 2d 35 35 38 80, 169480, -558
1150: 34 30 34 2c 20 2d 31 39 33 39 32 35 2c 20 36 35 404, -193925, 65
1160: 34 34 34 34 2c 20 2d 33 30 30 32 34 37 2c 20 2d 4444, -300247, -
1170: 34 35 36 34 32 30 2c 20 2d 31 31 39 34 33 36 2c 456420, -119436,
1180: 20 2d 36 32 30 39 32 30 2c 20 2d 34 37 30 30 31 -620920, -47001
1190: 38 2c 20 2d 39 31 34 32 37 32 2c 20 2d 36 39 31 8, -914272, -691
11a0: 32 35 36 2c 20 2d 34 39 34 31 38 2c 20 2d 32 31 256, -49418, -21
11b0: 30 35 34 2c 20 36 30 33 33 37 33 2c 20 2d 32 33 054, 603373, -23
11c0: 36 35 36 2c 20 38 39 31 36 39 31 2c 20 32 35 38 656, 891691, 258
11d0: 39 38 36 2c 20 2d 34 35 33 37 39 33 2c 20 2d 37 986, -453793, -7
11e0: 38 32 39 34 30 7d 3b 0d 0a 09 20 20 76 65 63 74 82940};... vect
11f0: 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 78 or <int> x(x_, x
1200: 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69 7a _+sizeof(x_)/siz
1210: 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e eof(*x_)); ...in
1220: 74 20 79 5f 5b 5d 20 3d 20 7b 2d 37 37 33 31 38 t y_[] = {-77318
1230: 2c 20 2d 36 33 32 36 32 39 2c 20 2d 33 34 34 39 , -632629, -3449
1240: 34 32 2c 20 2d 33 36 31 37 30 36 2c 20 31 39 31 42, -361706, 191
1250: 39 38 32 2c 20 33 34 39 34 32 34 2c 20 36 37 36 982, 349424, 676
1260: 35 33 36 2c 20 31 36 36 31 32 34 2c 20 32 39 31 536, 166124, 291
1270: 33 34 32 2c 20 2d 32 36 38 39 36 38 2c 20 31 38 342, -268968, 18
1280: 38 32 36 32 2c 20 2d 35 33 37 39 35 33 2c 20 2d 8262, -537953, -
1290: 37 30 34 33 32 2c 20 31 35 36 38 30 33 2c 20 31 70432, 156803, 1
12a0: 36 36 31 37 34 2c 20 33 34 35 31 32 38 2c 20 35 66174, 345128, 5
12b0: 38 36 31 34 2c 20 2d 36 37 31 37 34 37 2c 20 35 8614, -671747, 5
12c0: 30 38 32 36 35 2c 20 39 32 33 32 34 7d 3b 0d 0a 08265, 92324};..
12d0: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
12e0: 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 y(y_, y_+sizeof(
12f0: 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 y_)/sizeof(*y_))
1300: 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 ; ...double _ =
1310: 33 36 38 37 39 2e 31 35 31 32 37 36 33 34 32 39 36879.1512763429
1320: 3b 20 0d 0a 45 4e 44 0d 0a 2f 2a 0d 0a 43 41 53 ; ..END../*..CAS
1330: 45 28 36 29 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 E(6)...int x_[]
1340: 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c = ;... vector <
1350: 69 6e 74 3e 20 78 28 78 5f 2c 20 78 5f 2b 73 69 int> x(x_, x_+si
1360: 7a 65 6f 66 28 78 5f 29 2f 73 69 7a 65 6f 66 28 zeof(x_)/sizeof(
1370: 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 79 5f *x_)); ...int y_
1380: 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f [] = ;... vecto
1390: 72 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79 5f r <int> y(y_, y_
13a0: 2b 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a 65 +sizeof(y_)/size
13b0: 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 64 6f 75 of(*y_)); ...dou
13c0: 62 6c 65 20 5f 20 3d 20 3b 20 0d 0a 45 4e 44 0d ble _ = ; ..END.
13d0: 0a 43 41 53 45 28 37 29 0d 0a 09 69 6e 74 20 78 .CASE(7)...int x
13e0: 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 _[] = ;... vect
13f0: 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 78 or <int> x(x_, x
1400: 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69 7a _+sizeof(x_)/siz
1410: 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e eof(*x_)); ...in
1420: 74 20 79 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 t y_[] = ;... v
1430: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f ector <int> y(y_
1440: 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f , y_+sizeof(y_)/
1450: 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a sizeof(*y_)); ..
1460: 09 64 6f 75 62 6c 65 20 5f 20 3d 20 3b 20 0d 0a .double _ = ; ..
1470: 45 4e 44 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f 20 45 END..*/..}..// E
1480: 4e 44 20 43 55 54 20 48 45 52 45 0d 0a ND CUT HERE..