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 63 73 74 72 69 6e 67 3e 0d 0a 75 73 69 6e 67 <cstring>..using
0150: 20 6e 61 6d 65 73 70 61 63 65 20 73 74 64 3b 0d namespace std;.
0160: 0a 74 79 70 65 64 65 66 20 6c 6f 6e 67 20 6c 6f .typedef long lo
0170: 6e 67 20 4c 4c 3b 0d 0a 74 79 70 65 64 65 66 20 ng LL;..typedef
0180: 63 6f 6d 70 6c 65 78 3c 64 6f 75 62 6c 65 3e 20 complex<double>
0190: 43 4d 50 3b 0d 0a 0d 0a 74 65 6d 70 6c 61 74 65 CMP;....template
01a0: 3c 74 79 70 65 6e 61 6d 65 20 54 3e 0d 0a 63 6c <typename T>..cl
01b0: 61 73 73 20 49 64 47 65 6e 0d 0a 7b 0d 0a 09 6d ass IdGen..{...m
01c0: 61 70 3c 54 2c 20 69 6e 74 3e 20 76 32 69 64 5f ap<T, int> v2id_
01d0: 3b 0d 0a 09 76 65 63 74 6f 72 3c 54 3e 20 20 20 ;...vector<T>
01e0: 69 64 32 76 5f 3b 0d 0a 70 75 62 6c 69 63 3a 0d id2v_;..public:.
01f0: 0a 09 69 6e 74 20 76 32 69 64 28 63 6f 6e 73 74 ..int v2id(const
0200: 20 54 26 20 76 29 20 7b 0d 0a 09 09 69 66 28 20 T& v) {....if(
0210: 21 76 32 69 64 5f 2e 63 6f 75 6e 74 28 76 29 20 !v2id_.count(v)
0220: 29 20 7b 20 76 32 69 64 5f 5b 76 5d 20 3d 20 73 ) { v2id_[v] = s
0230: 69 7a 65 28 29 3b 20 69 64 32 76 5f 2e 70 75 73 ize(); id2v_.pus
0240: 68 5f 62 61 63 6b 28 76 29 3b 20 7d 0d 0a 09 09 h_back(v); }....
0250: 72 65 74 75 72 6e 20 76 32 69 64 5f 5b 76 5d 3b return v2id_[v];
0260: 0d 0a 09 7d 0d 0a 09 63 6f 6e 73 74 20 54 26 20 ...}...const T&
0270: 69 64 32 76 28 69 6e 74 20 69 29 20 63 6f 6e 73 id2v(int i) cons
0280: 74 20 7b 20 72 65 74 75 72 6e 20 69 64 32 76 5f t { return id2v_
0290: 5b 69 5d 3b 20 7d 0d 0a 09 69 6e 74 20 73 69 7a [i]; }...int siz
02a0: 65 28 29 20 63 6f 6e 73 74 20 7b 20 72 65 74 75 e() const { retu
02b0: 72 6e 20 69 64 32 76 5f 2e 73 69 7a 65 28 29 3b rn id2v_.size();
02c0: 20 7d 0d 0a 7d 3b 0d 0a 0d 0a 74 65 6d 70 6c 61 }..};....templa
02d0: 74 65 3c 74 79 70 65 6e 61 6d 65 20 56 65 72 74 te<typename Vert
02e0: 2c 20 74 79 70 65 6e 61 6d 65 20 43 6f 73 74 2c , typename Cost,
02f0: 20 74 79 70 65 6e 61 6d 65 20 46 6c 6f 77 2c 20 typename Flow,
0300: 69 6e 74 20 4e 56 3d 32 35 36 3e 0d 0a 63 6c 61 int NV=256>..cla
0310: 73 73 20 4d 69 6e 43 6f 73 74 46 6c 6f 77 0d 0a ss MinCostFlow..
0320: 7b 0d 0a 09 49 64 47 65 6e 3c 56 65 72 74 3e 20 {...IdGen<Vert>
0330: 69 64 67 65 6e 3b 0d 0a 0d 0a 09 76 65 63 74 6f idgen;.....vecto
0340: 72 3c 69 6e 74 3e 20 47 5b 4e 56 5d 3b 0d 0a 09 r<int> G[NV];...
0350: 46 6c 6f 77 20 46 5b 4e 56 5d 5b 4e 56 5d 3b 0d Flow F[NV][NV];.
0360: 0a 09 43 6f 73 74 20 43 5b 4e 56 5d 5b 4e 56 5d ..Cost C[NV][NV]
0370: 3b 0d 0a 0d 0a 70 75 62 6c 69 63 3a 0d 0a 09 76 ;....public:...v
0380: 6f 69 64 20 61 64 64 45 64 67 65 28 20 56 65 72 oid addEdge( Ver
0390: 74 20 73 5f 2c 20 56 65 72 74 20 74 5f 2c 20 43 t s_, Vert t_, C
03a0: 6f 73 74 20 63 2c 20 46 6c 6f 77 20 66 20 29 0d ost c, Flow f ).
03b0: 0a 09 7b 0d 0a 09 09 69 6e 74 20 73 20 3d 20 69 ..{....int s = i
03c0: 64 67 65 6e 2e 76 32 69 64 28 73 5f 29 2c 20 74 dgen.v2id(s_), t
03d0: 20 3d 20 69 64 67 65 6e 2e 76 32 69 64 28 74 5f = idgen.v2id(t_
03e0: 29 3b 0d 0a 09 09 47 5b 73 5d 2e 70 75 73 68 5f );....G[s].push_
03f0: 62 61 63 6b 28 74 29 3b 0d 0a 09 09 47 5b 74 5d back(t);....G[t]
0400: 2e 70 75 73 68 5f 62 61 63 6b 28 73 29 3b 0d 0a .push_back(s);..
0410: 09 09 43 5b 73 5d 5b 74 5d 20 3d 20 63 3b 0d 0a ..C[s][t] = c;..
0420: 09 09 43 5b 74 5d 5b 73 5d 20 3d 20 2d 63 3b 0d ..C[t][s] = -c;.
0430: 0a 09 09 46 5b 73 5d 5b 74 5d 20 3d 20 66 3b 0d ...F[s][t] = f;.
0440: 0a 09 09 46 5b 74 5d 5b 73 5d 20 3d 20 30 3b 0d ...F[t][s] = 0;.
0450: 0a 09 7d 0d 0a 0d 0a 09 70 61 69 72 3c 43 6f 73 ..}.....pair<Cos
0460: 74 2c 20 46 6c 6f 77 3e 20 63 61 6c 63 28 20 56 t, Flow> calc( V
0470: 65 72 74 20 73 5f 2c 20 56 65 72 74 20 74 5f 20 ert s_, Vert t_
0480: 29 0d 0a 09 7b 0d 0a 09 09 63 6f 6e 73 74 20 69 )...{....const i
0490: 6e 74 20 4e 3d 69 64 67 65 6e 2e 73 69 7a 65 28 nt N=idgen.size(
04a0: 29 2c 20 53 3d 69 64 67 65 6e 2e 76 32 69 64 28 ), S=idgen.v2id(
04b0: 73 5f 29 2c 20 54 3d 69 64 67 65 6e 2e 76 32 69 s_), T=idgen.v2i
04c0: 64 28 74 5f 29 3b 0d 0a 09 09 73 74 61 74 69 63 d(t_);....static
04d0: 20 63 6f 6e 73 74 20 43 6f 73 74 20 43 4f 53 54 const Cost COST
04e0: 5f 49 4e 46 20 3d 20 31 65 2b 33 30 30 3b 20 2f _INF = 1e+300; /
04f0: 2f 20 21 21 45 44 49 54 20 48 45 52 45 21 21 0d / !!EDIT HERE!!.
0500: 0a 09 09 73 74 61 74 69 63 20 63 6f 6e 73 74 20 ...static const
0510: 46 6c 6f 77 20 46 4c 4f 57 5f 49 4e 46 20 3d 20 Flow FLOW_INF =
0520: 30 78 37 66 66 66 66 66 66 66 3b 0d 0a 0d 0a 09 0x7fffffff;.....
0530: 09 43 6f 73 74 20 74 6f 74 61 6c 5f 63 6f 73 74 .Cost total_cost
0540: 20 3d 20 30 3b 0d 0a 09 09 46 6c 6f 77 20 74 6f = 0;....Flow to
0550: 74 61 6c 5f 66 6c 6f 77 20 3d 20 30 3b 0d 0a 09 tal_flow = 0;...
0560: 09 76 65 63 74 6f 72 3c 43 6f 73 74 3e 20 64 69 .vector<Cost> di
0570: 73 74 28 4e 2c 20 30 29 3b 20 2f 2f 20 44 69 73 st(N, 0); // Dis
0580: 74 61 6e 63 65 20 66 72 6f 6d 20 53 20 3a 20 69 tance from S : i
0590: 6e 69 74 69 61 6c 6c 79 20 75 6e 6b 6e 6f 77 6e nitially unknown
05a0: 2e 0d 0a 09 09 66 6f 72 28 3b 3b 29 0d 0a 09 09 .....for(;;)....
05b0: 7b 0d 0a 09 09 09 2f 2f 20 44 69 6a 6b 73 74 72 {.....// Dijkstr
05c0: 61 20 3a 20 66 69 6e 64 20 74 68 65 20 22 73 68 a : find the "sh
05d0: 6f 72 74 65 73 74 20 70 61 74 68 22 20 66 72 6f ortest path" fro
05e0: 6d 20 53 20 74 6f 20 54 20 77 72 74 20 43 5b 5d m S to T wrt C[]
05f0: 5b 5d 2e 0d 0a 09 09 09 2f 2f 20 20 20 43 5b 5d []......// C[]
0600: 5b 5d 20 63 61 6e 20 62 65 20 3c 30 20 73 6f 20 [] can be <0 so
0610: 77 65 20 6d 75 73 74 20 62 65 20 63 61 72 65 66 we must be caref
0620: 75 6c 2e 20 49 6e 73 74 65 61 64 20 6f 66 20 63 ul. Instead of c
0630: 6f 6d 70 75 74 69 6e 67 20 74 68 65 20 73 68 6f omputing the sho
0640: 72 74 65 73 74 20 70 61 74 68 20 64 69 72 65 63 rtest path direc
0650: 74 6c 79 2c 0d 0a 09 09 09 2f 2f 20 20 20 77 65 tly,.....// we
0660: 20 63 6f 6d 70 75 74 65 20 74 68 65 20 69 6e 63 compute the inc
0670: 72 65 61 73 65 20 28 22 64 65 6c 74 61 22 29 20 rease ("delta")
0680: 66 72 6f 6d 20 74 68 65 20 73 68 6f 72 74 65 73 from the shortes
0690: 74 20 70 61 74 68 20 69 6e 20 74 68 65 20 70 72 t path in the pr
06a0: 65 76 69 6f 75 73 20 69 74 65 72 61 74 69 6f 6e evious iteration
06b0: 2e 0d 0a 09 09 09 2f 2f 20 20 20 53 69 6e 63 65 ......// Since
06c0: 20 73 68 6f 72 74 65 73 74 20 70 61 74 68 20 63 shortest path c
06d0: 61 6e 6e 6f 74 20 64 65 63 72 65 61 73 65 2c 20 annot decrease,
06e0: 64 65 6c 74 61 20 69 73 20 61 6c 77 61 79 73 20 delta is always
06f0: 3e 3d 30 20 77 68 65 6e 20 74 72 61 76 65 72 73 >=0 when travers
0700: 69 6e 67 20 65 64 67 65 73 2e 0d 0a 09 09 09 2f ing edges....../
0710: 2f 20 20 20 53 6d 61 6c 6c 65 73 74 20 64 65 6c / Smallest del
0720: 74 61 20 69 6d 70 6c 69 65 73 20 73 6d 61 6c 6c ta implies small
0730: 65 73 74 20 64 69 73 74 5b 54 5d 2b 64 65 6c 74 est dist[T]+delt
0740: 61 5b 54 5d 2e 0d 0a 09 09 09 76 65 63 74 6f 72 a[T]......vector
0750: 3c 43 6f 73 74 3e 20 64 65 6c 74 61 28 4e 2c 20 <Cost> delta(N,
0760: 43 4f 53 54 5f 49 4e 46 29 3b 20 64 65 6c 74 61 COST_INF); delta
0770: 5b 53 5d 20 3d 20 30 3b 0d 0a 09 09 09 76 65 63 [S] = 0;.....vec
0780: 74 6f 72 3c 69 6e 74 3e 20 20 70 72 65 76 28 4e tor<int> prev(N
0790: 2c 20 2d 31 29 3b 0d 0a 0d 0a 09 09 09 74 79 70 , -1);.......typ
07a0: 65 64 65 66 20 70 61 69 72 3c 20 43 6f 73 74 2c edef pair< Cost,
07b0: 20 70 61 69 72 3c 69 6e 74 2c 20 69 6e 74 3e 20 pair<int, int>
07c0: 3e 20 63 65 64 67 65 3b 0d 0a 09 09 09 70 72 69 > cedge;.....pri
07d0: 6f 72 69 74 79 5f 71 75 65 75 65 3c 20 63 65 64 ority_queue< ced
07e0: 67 65 2c 20 76 65 63 74 6f 72 3c 63 65 64 67 65 ge, vector<cedge
07f0: 3e 2c 20 67 72 65 61 74 65 72 3c 63 65 64 67 65 >, greater<cedge
0800: 3e 20 3e 20 51 3b 0d 0a 09 09 09 51 2e 70 75 73 > > Q;.....Q.pus
0810: 68 28 20 63 65 64 67 65 28 30 2c 20 6d 61 6b 65 h( cedge(0, make
0820: 5f 70 61 69 72 28 53 2c 20 53 29 29 20 29 3b 0d _pair(S, S)) );.
0830: 0a 09 09 09 77 68 69 6c 65 28 20 21 51 2e 65 6d ....while( !Q.em
0840: 70 74 79 28 29 20 29 20 7b 0d 0a 09 09 09 09 63 pty() ) {......c
0850: 6f 6e 73 74 20 63 65 64 67 65 20 65 20 3d 20 51 onst cedge e = Q
0860: 2e 74 6f 70 28 29 3b 20 51 2e 70 6f 70 28 29 3b .top(); Q.pop();
0870: 0d 0a 09 09 09 09 63 6f 6e 73 74 20 69 6e 74 20 ......const int
0880: 75 5f 70 72 65 76 20 3d 20 65 2e 73 65 63 6f 6e u_prev = e.secon
0890: 64 2e 66 69 72 73 74 3b 0d 0a 09 09 09 09 63 6f d.first;......co
08a0: 6e 73 74 20 69 6e 74 20 75 20 3d 20 65 2e 73 65 nst int u = e.se
08b0: 63 6f 6e 64 2e 73 65 63 6f 6e 64 3b 0d 0a 09 09 cond.second;....
08c0: 09 09 69 66 28 20 70 72 65 76 5b 75 5d 20 3e 3d ..if( prev[u] >=
08d0: 20 30 20 29 20 2f 2f 20 76 69 73 69 74 65 64 0d 0 ) // visited.
08e0: 0a 09 09 09 09 09 63 6f 6e 74 69 6e 75 65 3b 0d ......continue;.
08f0: 0a 09 09 09 09 70 72 65 76 5b 75 5d 20 3d 20 75 .....prev[u] = u
0900: 5f 70 72 65 76 3b 0d 0a 0d 0a 09 09 09 09 66 6f _prev;........fo
0910: 72 28 69 6e 74 20 69 3d 30 3b 20 69 3c 47 5b 75 r(int i=0; i<G[u
0920: 5d 2e 73 69 7a 65 28 29 3b 20 2b 2b 69 29 20 7b ].size(); ++i) {
0930: 0d 0a 09 09 09 09 09 63 6f 6e 73 74 20 69 6e 74 .......const int
0940: 20 20 76 20 3d 20 47 5b 75 5d 5b 69 5d 3b 0d 0a v = G[u][i];..
0950: 09 09 09 09 09 63 6f 6e 73 74 20 43 6f 73 74 20 .....const Cost
0960: 76 5f 64 65 6c 74 61 20 3d 20 64 69 73 74 5b 75 v_delta = dist[u
0970: 5d 2b 64 65 6c 74 61 5b 75 5d 2b 43 5b 75 5d 5b ]+delta[u]+C[u][
0980: 76 5d 20 2d 20 64 69 73 74 5b 76 5d 3b 0d 0a 09 v] - dist[v];...
0990: 09 09 09 09 69 66 28 20 46 5b 75 5d 5b 76 5d 3e ....if( F[u][v]>
09a0: 30 20 26 26 20 64 65 6c 74 61 5b 76 5d 3e 76 5f 0 && delta[v]>v_
09b0: 64 65 6c 74 61 20 29 0d 0a 09 09 09 09 09 09 51 delta )........Q
09c0: 2e 70 75 73 68 28 20 63 65 64 67 65 28 64 65 6c .push( cedge(del
09d0: 74 61 5b 76 5d 3d 76 5f 64 65 6c 74 61 2c 20 6d ta[v]=v_delta, m
09e0: 61 6b 65 5f 70 61 69 72 28 75 2c 76 29 29 20 29 ake_pair(u,v)) )
09f0: 3b 0d 0a 09 09 09 09 7d 0d 0a 09 09 09 7d 0d 0a ;......}.....}..
0a00: 0d 0a 09 09 09 2f 2f 20 49 66 20 54 20 69 73 20 .....// If T is
0a10: 75 6e 72 65 61 63 68 61 62 6c 65 2c 20 66 69 6e unreachable, fin
0a20: 69 73 68 65 64 2e 0d 0a 09 09 09 69 66 28 20 70 ished......if( p
0a30: 72 65 76 5b 54 5d 20 3c 20 30 20 29 0d 0a 09 09 rev[T] < 0 )....
0a40: 09 09 62 72 65 61 6b 3b 0d 0a 0d 0a 09 09 09 2f ..break;......./
0a50: 2f 20 55 70 64 61 74 65 20 74 68 65 20 64 69 73 / Update the dis
0a60: 74 61 6e 63 65 20 74 61 62 6c 65 2e 0d 0a 09 09 tance table.....
0a70: 09 66 6f 72 28 69 6e 74 20 75 3d 30 3b 20 75 3c .for(int u=0; u<
0a80: 4e 3b 20 2b 2b 75 29 0d 0a 09 09 09 09 69 66 28 N; ++u)......if(
0a90: 20 64 65 6c 74 61 5b 75 5d 20 21 3d 20 43 4f 53 delta[u] != COS
0aa0: 54 5f 49 4e 46 20 29 0d 0a 09 09 09 09 09 64 69 T_INF ).......di
0ab0: 73 74 5b 75 5d 20 2b 3d 20 64 65 6c 74 61 5b 75 st[u] += delta[u
0ac0: 5d 3b 0d 0a 0d 0a 09 09 09 2f 2f 20 48 6f 77 20 ];.......// How
0ad0: 6d 75 63 68 20 77 61 74 65 72 20 63 61 6e 20 66 much water can f
0ae0: 6c 6f 77 20 6f 6e 20 74 68 65 20 6d 69 6e 2d 63 low on the min-c
0af0: 6f 73 74 20 70 61 74 68 3f 0d 0a 09 09 09 46 6c ost path?.....Fl
0b00: 6f 77 20 66 20 3d 20 46 4c 4f 57 5f 49 4e 46 3b ow f = FLOW_INF;
0b10: 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 75 3d 54 .....for(int u=T
0b20: 3b 20 75 21 3d 53 3b 20 75 3d 70 72 65 76 5b 75 ; u!=S; u=prev[u
0b30: 5d 29 0d 0a 09 09 09 09 66 20 3d 20 6d 69 6e 28 ])......f = min(
0b40: 66 2c 20 46 5b 70 72 65 76 5b 75 5d 5d 5b 75 5d f, F[prev[u]][u]
0b50: 29 3b 0d 0a 0d 0a 09 09 09 2f 2f 20 52 75 6e 20 );.......// Run
0b60: 74 68 65 20 66 6c 6f 77 20 61 73 20 6d 75 63 68 the flow as much
0b70: 20 61 73 20 70 6f 73 73 69 62 6c 65 0d 0a 09 09 as possible....
0b80: 09 74 6f 74 61 6c 5f 66 6c 6f 77 20 2b 3d 20 66 .total_flow += f
0b90: 3b 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 75 3d ;.....for(int u=
0ba0: 54 3b 20 75 21 3d 53 3b 20 75 3d 70 72 65 76 5b T; u!=S; u=prev[
0bb0: 75 5d 29 20 7b 0d 0a 09 09 09 09 74 6f 74 61 6c u]) {......total
0bc0: 5f 63 6f 73 74 20 20 20 20 2b 3d 20 66 20 2a 20 _cost += f *
0bd0: 43 5b 70 72 65 76 5b 75 5d 5d 5b 75 5d 3b 0d 0a C[prev[u]][u];..
0be0: 09 09 09 09 46 5b 70 72 65 76 5b 75 5d 5d 5b 75 ....F[prev[u]][u
0bf0: 5d 20 2d 3d 20 66 3b 0d 0a 09 09 09 09 46 5b 75 ] -= f;......F[u
0c00: 5d 5b 70 72 65 76 5b 75 5d 5d 20 2b 3d 20 66 3b ][prev[u]] += f;
0c10: 0d 0a 09 09 09 7d 0d 0a 09 09 7d 0d 0a 09 09 72 .....}....}....r
0c20: 65 74 75 72 6e 20 6d 61 6b 65 5f 70 61 69 72 28 eturn make_pair(
0c30: 74 6f 74 61 6c 5f 63 6f 73 74 2c 20 74 6f 74 61 total_cost, tota
0c40: 6c 5f 66 6c 6f 77 29 3b 0d 0a 09 7d 0d 0a 7d 3b l_flow);...}..};
0c50: 0d 0a 0d 0a 63 6c 61 73 73 20 46 6f 78 43 61 72 ....class FoxCar
0c60: 64 47 61 6d 65 20 7b 20 70 75 62 6c 69 63 3a 0d dGame { public:.
0c70: 0a 09 64 6f 75 62 6c 65 20 74 68 65 4d 61 78 50 ..double theMaxP
0c80: 72 6f 70 6f 72 74 69 6f 6e 28 76 65 63 74 6f 72 roportion(vector
0c90: 20 3c 64 6f 75 62 6c 65 3e 20 70 69 6c 65 41 2c <double> pileA,
0ca0: 20 76 65 63 74 6f 72 20 3c 64 6f 75 62 6c 65 3e vector <double>
0cb0: 20 70 69 6c 65 42 2c 20 69 6e 74 20 6b 29 20 0d pileB, int k) .
0cc0: 0a 09 7b 0d 0a 09 09 64 6f 75 62 6c 65 20 4c 3d ..{....double L=
0cd0: 31 2c 20 52 3d 35 30 3b 0d 0a 09 09 66 6f 72 28 1, R=50;....for(
0ce0: 69 6e 74 20 69 3d 30 3b 20 69 3c 35 30 3b 20 2b int i=0; i<50; +
0cf0: 2b 69 29 0d 0a 09 09 09 28 70 6f 73 73 69 62 6c +i).....(possibl
0d00: 65 28 70 69 6c 65 41 2c 20 70 69 6c 65 42 2c 20 e(pileA, pileB,
0d10: 6b 2c 20 28 4c 2b 52 29 2f 32 29 20 3f 20 4c 20 k, (L+R)/2) ? L
0d20: 3a 20 52 29 20 3d 20 28 4c 2b 52 29 2f 32 3b 0d : R) = (L+R)/2;.
0d30: 0a 09 09 72 65 74 75 72 6e 20 4c 3b 0d 0a 09 7d ...return L;...}
0d40: 0d 0a 0d 0a 09 65 6e 75 6d 20 54 61 67 20 7b 53 .....enum Tag {S
0d50: 72 63 2c 20 4c 65 66 74 2c 20 52 69 67 68 74 2c rc, Left, Right,
0d60: 20 54 61 72 67 65 74 2c 20 47 6f 61 6c 7d 3b 0d Target, Goal};.
0d70: 0a 09 62 6f 6f 6c 20 70 6f 73 73 69 62 6c 65 28 ..bool possible(
0d80: 76 65 63 74 6f 72 20 3c 64 6f 75 62 6c 65 3e 20 vector <double>
0d90: 70 69 6c 65 41 2c 20 76 65 63 74 6f 72 20 3c 64 pileA, vector <d
0da0: 6f 75 62 6c 65 3e 20 70 69 6c 65 42 2c 20 69 6e ouble> pileB, in
0db0: 74 20 6b 2c 20 64 6f 75 62 6c 65 20 52 29 0d 0a t k, double R)..
0dc0: 09 7b 0d 0a 09 09 4d 69 6e 43 6f 73 74 46 6c 6f .{....MinCostFlo
0dd0: 77 3c 70 61 69 72 3c 54 61 67 2c 69 6e 74 3e 2c w<pair<Tag,int>,
0de0: 20 64 6f 75 62 6c 65 2c 20 69 6e 74 3e 20 6d 63 double, int> mc
0df0: 66 3b 0d 0a 0d 0a 09 09 66 6f 72 28 69 6e 74 20 f;......for(int
0e00: 61 3d 30 3b 20 61 3c 70 69 6c 65 41 2e 73 69 7a a=0; a<pileA.siz
0e10: 65 28 29 3b 20 2b 2b 61 29 0d 0a 09 09 09 6d 63 e(); ++a).....mc
0e20: 66 2e 61 64 64 45 64 67 65 28 20 6d 61 6b 65 5f f.addEdge( make_
0e30: 70 61 69 72 28 53 72 63 2c 30 29 2c 20 6d 61 6b pair(Src,0), mak
0e40: 65 5f 70 61 69 72 28 4c 65 66 74 2c 61 29 2c 20 e_pair(Left,a),
0e50: 30 2e 30 2c 20 31 20 29 3b 0d 0a 09 09 66 6f 72 0.0, 1 );....for
0e60: 28 69 6e 74 20 61 3d 30 3b 20 61 3c 70 69 6c 65 (int a=0; a<pile
0e70: 41 2e 73 69 7a 65 28 29 3b 20 2b 2b 61 29 0d 0a A.size(); ++a)..
0e80: 09 09 66 6f 72 28 69 6e 74 20 62 3d 30 3b 20 62 ..for(int b=0; b
0e90: 3c 70 69 6c 65 42 2e 73 69 7a 65 28 29 3b 20 2b <pileB.size(); +
0ea0: 2b 62 29 0d 0a 09 09 7b 0d 0a 09 09 09 64 6f 75 +b)....{.....dou
0eb0: 62 6c 65 20 78 5f 20 3d 20 70 69 6c 65 41 5b 61 ble x_ = pileA[a
0ec0: 5d 2b 70 69 6c 65 42 5b 62 5d 2c 20 79 5f 20 3d ]+pileB[b], y_ =
0ed0: 20 70 69 6c 65 41 5b 61 5d 2a 70 69 6c 65 42 5b pileA[a]*pileB[
0ee0: 62 5d 3b 0d 0a 09 09 09 64 6f 75 62 6c 65 20 78 b];.....double x
0ef0: 20 3d 20 6d 61 78 28 78 5f 2c 79 5f 29 2c 20 79 = max(x_,y_), y
0f00: 20 3d 20 6d 69 6e 28 78 5f 2c 79 5f 29 3b 0d 0a = min(x_,y_);..
0f10: 09 09 09 6d 63 66 2e 61 64 64 45 64 67 65 28 20 ...mcf.addEdge(
0f20: 6d 61 6b 65 5f 70 61 69 72 28 4c 65 66 74 2c 61 make_pair(Left,a
0f30: 29 2c 20 6d 61 6b 65 5f 70 61 69 72 28 52 69 67 ), make_pair(Rig
0f40: 68 74 2c 62 29 2c 20 52 2a 79 2d 78 2b 31 30 30 ht,b), R*y-x+100
0f50: 30 30 2e 30 2c 20 31 20 29 3b 0d 0a 09 09 7d 0d 00.0, 1 );....}.
0f60: 0a 09 09 66 6f 72 28 69 6e 74 20 62 3d 30 3b 20 ...for(int b=0;
0f70: 62 3c 70 69 6c 65 42 2e 73 69 7a 65 28 29 3b 20 b<pileB.size();
0f80: 2b 2b 62 29 0d 0a 09 09 09 6d 63 66 2e 61 64 64 ++b).....mcf.add
0f90: 45 64 67 65 28 20 6d 61 6b 65 5f 70 61 69 72 28 Edge( make_pair(
0fa0: 52 69 67 68 74 2c 62 29 2c 20 6d 61 6b 65 5f 70 Right,b), make_p
0fb0: 61 69 72 28 54 61 72 67 65 74 2c 30 29 2c 20 30 air(Target,0), 0
0fc0: 2e 30 2c 20 31 20 29 3b 0d 0a 09 09 6d 63 66 2e .0, 1 );....mcf.
0fd0: 61 64 64 45 64 67 65 28 20 6d 61 6b 65 5f 70 61 addEdge( make_pa
0fe0: 69 72 28 54 61 72 67 65 74 2c 30 29 2c 20 6d 61 ir(Target,0), ma
0ff0: 6b 65 5f 70 61 69 72 28 47 6f 61 6c 2c 30 29 2c ke_pair(Goal,0),
1000: 20 30 2e 30 2c 20 6b 20 29 3b 0d 0a 0d 0a 09 09 0.0, k );......
1010: 72 65 74 75 72 6e 20 6d 63 66 2e 63 61 6c 63 28 return mcf.calc(
1020: 20 6d 61 6b 65 5f 70 61 69 72 28 53 72 63 2c 30 make_pair(Src,0
1030: 29 2c 20 6d 61 6b 65 5f 70 61 69 72 28 47 6f 61 ), make_pair(Goa
1040: 6c 2c 30 29 20 29 2e 66 69 72 73 74 20 3c 3d 20 l,0) ).first <=
1050: 31 30 30 30 30 2e 30 2a 6b 3b 0d 0a 09 7d 0d 0a 10000.0*k;...}..
1060: 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45 47 49 4e 20 43 };....// BEGIN C
1070: 55 54 20 48 45 52 45 0d 0a 23 69 6e 63 6c 75 64 UT HERE..#includ
1080: 65 20 3c 63 74 69 6d 65 3e 0d 0a 64 6f 75 62 6c e <ctime>..doubl
1090: 65 20 73 74 61 72 74 5f 74 69 6d 65 3b 20 73 74 e start_time; st
10a0: 72 69 6e 67 20 74 69 6d 65 72 28 29 0d 0a 20 7b ring timer().. {
10b0: 20 6f 73 74 72 69 6e 67 73 74 72 65 61 6d 20 6f ostringstream o
10c0: 73 3b 20 6f 73 20 3c 3c 20 22 20 28 22 20 3c 3c s; os << " (" <<
10d0: 20 69 6e 74 28 28 63 6c 6f 63 6b 28 29 2d 73 74 int((clock()-st
10e0: 61 72 74 5f 74 69 6d 65 29 2f 43 4c 4f 43 4b 53 art_time)/CLOCKS
10f0: 5f 50 45 52 5f 53 45 43 2a 31 30 30 30 29 20 3c _PER_SEC*1000) <
1100: 3c 20 22 20 6d 73 65 63 29 22 3b 20 72 65 74 75 < " msec)"; retu
1110: 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d 0d 0a rn os.str(); }..
1120: 74 65 6d 70 6c 61 74 65 3c 74 79 70 65 6e 61 6d template<typenam
1130: 65 20 54 3e 20 6f 73 74 72 65 61 6d 26 20 6f 70 e T> ostream& op
1140: 65 72 61 74 6f 72 3c 3c 28 6f 73 74 72 65 61 6d erator<<(ostream
1150: 26 20 6f 73 2c 20 63 6f 6e 73 74 20 76 65 63 74 & os, const vect
1160: 6f 72 3c 54 3e 26 20 76 29 0d 0a 20 7b 20 6f 73 or<T>& v).. { os
1170: 20 3c 3c 20 22 7b 20 22 3b 0d 0a 20 20 20 66 6f << "{ ";.. fo
1180: 72 28 74 79 70 65 6e 61 6d 65 20 76 65 63 74 6f r(typename vecto
1190: 72 3c 54 3e 3a 3a 63 6f 6e 73 74 5f 69 74 65 72 r<T>::const_iter
11a0: 61 74 6f 72 20 69 74 3d 76 2e 62 65 67 69 6e 28 ator it=v.begin(
11b0: 29 3b 20 69 74 21 3d 76 2e 65 6e 64 28 29 3b 20 ); it!=v.end();
11c0: 2b 2b 69 74 29 0d 0a 20 20 20 6f 73 20 3c 3c 20 ++it).. os <<
11d0: 27 5c 22 27 20 3c 3c 20 2a 69 74 20 3c 3c 20 27 '\"' << *it << '
11e0: 5c 22 27 20 3c 3c 20 28 69 74 2b 31 3d 3d 76 2e \"' << (it+1==v.
11f0: 65 6e 64 28 29 20 3f 20 22 22 20 3a 20 22 2c 20 end() ? "" : ",
1200: 22 29 3b 20 6f 73 20 3c 3c 20 22 20 7d 22 3b 20 "); os << " }";
1210: 72 65 74 75 72 6e 20 6f 73 3b 20 7d 0d 0a 76 6f return os; }..vo
1220: 69 64 20 76 65 72 69 66 79 5f 63 61 73 65 28 63 id verify_case(c
1230: 6f 6e 73 74 20 64 6f 75 62 6c 65 26 20 45 78 70 onst double& Exp
1240: 65 63 74 65 64 2c 20 63 6f 6e 73 74 20 64 6f 75 ected, const dou
1250: 62 6c 65 26 20 52 65 63 65 69 76 65 64 29 20 7b ble& Received) {
1260: 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 61 62 .. bool ok = (ab
1270: 73 28 45 78 70 65 63 74 65 64 20 2d 20 52 65 63 s(Expected - Rec
1280: 65 69 76 65 64 29 20 3c 20 31 65 2d 39 29 3b 0d eived) < 1e-9);.
1290: 0a 20 69 66 28 6f 6b 29 20 63 65 72 72 20 3c 3c . if(ok) cerr <<
12a0: 20 22 50 41 53 53 45 44 22 20 3c 3c 20 74 69 6d "PASSED" << tim
12b0: 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 20 20 65 er() << endl; e
12c0: 6c 73 65 20 7b 20 63 65 72 72 20 3c 3c 20 22 46 lse { cerr << "F
12d0: 41 49 4c 45 44 22 20 3c 3c 20 74 69 6d 65 72 28 AILED" << timer(
12e0: 29 20 3c 3c 20 65 6e 64 6c 3b 0d 0a 20 63 65 72 ) << endl;.. cer
12f0: 72 20 3c 3c 20 22 5c 74 6f 3a 20 5c 22 22 20 3c r << "\to: \"" <
1300: 3c 20 45 78 70 65 63 74 65 64 20 3c 3c 20 27 5c < Expected << '\
1310: 22 27 20 3c 3c 20 65 6e 64 6c 20 3c 3c 20 22 5c "' << endl << "\
1320: 74 78 3a 20 5c 22 22 20 3c 3c 20 52 65 63 65 69 tx: \"" << Recei
1330: 76 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 ved << '\"' << e
1340: 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 64 65 66 69 6e ndl; } }..#defin
1350: 65 20 43 41 53 45 28 4e 29 20 7b 63 65 72 72 20 e CASE(N) {cerr
1360: 3c 3c 20 22 54 65 73 74 20 43 61 73 65 20 23 22 << "Test Case #"
1370: 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e 22 20 3c << N << "..." <
1380: 3c 20 66 6c 75 73 68 3b 20 73 74 61 72 74 5f 74 < flush; start_t
1390: 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 0d 0a 23 64 ime=clock();..#d
13a0: 65 66 69 6e 65 20 45 4e 44 09 20 76 65 72 69 66 efine END. verif
13b0: 79 5f 63 61 73 65 28 5f 2c 20 46 6f 78 43 61 72 y_case(_, FoxCar
13c0: 64 47 61 6d 65 28 29 2e 74 68 65 4d 61 78 50 72 dGame().theMaxPr
13d0: 6f 70 6f 72 74 69 6f 6e 28 70 69 6c 65 41 2c 20 oportion(pileA,
13e0: 70 69 6c 65 42 2c 20 6b 29 29 3b 7d 0d 0a 69 6e pileB, k));}..in
13f0: 74 20 6d 61 69 6e 28 29 7b 0d 0a 0d 0a 43 41 53 t main(){....CAS
1400: 45 28 30 29 0d 0a 09 64 6f 75 62 6c 65 20 70 69 E(0)...double pi
1410: 6c 65 41 5f 5b 5d 20 3d 20 7b 31 2c 20 32 2c 20 leA_[] = {1, 2,
1420: 33 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 3};... vector <
1430: 64 6f 75 62 6c 65 3e 20 70 69 6c 65 41 28 70 69 double> pileA(pi
1440: 6c 65 41 5f 2c 20 70 69 6c 65 41 5f 2b 73 69 7a leA_, pileA_+siz
1450: 65 6f 66 28 70 69 6c 65 41 5f 29 2f 73 69 7a 65 eof(pileA_)/size
1460: 6f 66 28 2a 70 69 6c 65 41 5f 29 29 3b 20 0d 0a of(*pileA_)); ..
1470: 09 64 6f 75 62 6c 65 20 70 69 6c 65 42 5f 5b 5d .double pileB_[]
1480: 20 3d 20 7b 34 2c 20 35 2c 20 36 7d 3b 0d 0a 09 = {4, 5, 6};...
1490: 20 20 76 65 63 74 6f 72 20 3c 64 6f 75 62 6c 65 vector <double
14a0: 3e 20 70 69 6c 65 42 28 70 69 6c 65 42 5f 2c 20 > pileB(pileB_,
14b0: 70 69 6c 65 42 5f 2b 73 69 7a 65 6f 66 28 70 69 pileB_+sizeof(pi
14c0: 6c 65 42 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 69 leB_)/sizeof(*pi
14d0: 6c 65 42 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 6b leB_)); ...int k
14e0: 20 3d 20 32 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 = 2; ...double
14f0: 5f 20 3d 20 31 2e 37 36 39 32 33 30 37 36 39 32 _ = 1.7692307692
1500: 33 30 39 39 34 38 3b 20 0d 0a 45 4e 44 0d 0a 43 309948; ..END..C
1510: 41 53 45 28 31 29 0d 0a 09 64 6f 75 62 6c 65 20 ASE(1)...double
1520: 70 69 6c 65 41 5f 5b 5d 20 3d 20 7b 31 2e 32 33 pileA_[] = {1.23
1530: 34 2c 20 35 2e 36 37 38 2c 20 39 2e 30 31 32 2c 4, 5.678, 9.012,
1540: 20 33 2e 34 35 36 2c 20 37 2e 38 39 7d 3b 0d 0a 3.456, 7.89};..
1550: 09 20 20 76 65 63 74 6f 72 20 3c 64 6f 75 62 6c . vector <doubl
1560: 65 3e 20 70 69 6c 65 41 28 70 69 6c 65 41 5f 2c e> pileA(pileA_,
1570: 20 70 69 6c 65 41 5f 2b 73 69 7a 65 6f 66 28 70 pileA_+sizeof(p
1580: 69 6c 65 41 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 ileA_)/sizeof(*p
1590: 69 6c 65 41 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 ileA_)); ...doub
15a0: 6c 65 20 70 69 6c 65 42 5f 5b 5d 20 3d 20 7b 32 le pileB_[] = {2
15b0: 2e 33 34 35 2c 20 36 2e 37 38 39 2c 20 39 2e 38 .345, 6.789, 9.8
15c0: 37 36 2c 20 35 2e 34 33 32 2c 20 31 2e 30 31 32 76, 5.432, 1.012
15d0: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 64 };... vector <d
15e0: 6f 75 62 6c 65 3e 20 70 69 6c 65 42 28 70 69 6c ouble> pileB(pil
15f0: 65 42 5f 2c 20 70 69 6c 65 42 5f 2b 73 69 7a 65 eB_, pileB_+size
1600: 6f 66 28 70 69 6c 65 42 5f 29 2f 73 69 7a 65 6f of(pileB_)/sizeo
1610: 66 28 2a 70 69 6c 65 42 5f 29 29 3b 20 0d 0a 09 f(*pileB_)); ...
1620: 69 6e 74 20 6b 20 3d 20 33 3b 20 0d 0a 09 64 6f int k = 3; ...do
1630: 75 62 6c 65 20 5f 20 3d 20 34 2e 31 35 39 34 32 uble _ = 4.15942
1640: 34 34 32 30 30 37 39 35 38 36 3b 20 0d 0a 45 4e 4420079586; ..EN
1650: 44 0d 0a 43 41 53 45 28 32 29 0d 0a 09 64 6f 75 D..CASE(2)...dou
1660: 62 6c 65 20 70 69 6c 65 41 5f 5b 5d 20 3d 20 7b ble pileA_[] = {
1670: 31 2c 20 31 2e 31 2c 20 31 2e 32 2c 20 31 2e 33 1, 1.1, 1.2, 1.3
1680: 2c 20 31 2e 34 2c 20 31 2e 35 7d 3b 0d 0a 09 20 , 1.4, 1.5};...
1690: 20 76 65 63 74 6f 72 20 3c 64 6f 75 62 6c 65 3e vector <double>
16a0: 20 70 69 6c 65 41 28 70 69 6c 65 41 5f 2c 20 70 pileA(pileA_, p
16b0: 69 6c 65 41 5f 2b 73 69 7a 65 6f 66 28 70 69 6c ileA_+sizeof(pil
16c0: 65 41 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 69 6c eA_)/sizeof(*pil
16d0: 65 41 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 eA_)); ...double
16e0: 20 70 69 6c 65 42 5f 5b 5d 20 3d 20 7b 35 2c 20 pileB_[] = {5,
16f0: 31 30 2c 20 31 35 2c 20 32 30 2c 20 32 35 2c 20 10, 15, 20, 25,
1700: 33 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 30};... vector
1710: 3c 64 6f 75 62 6c 65 3e 20 70 69 6c 65 42 28 70 <double> pileB(p
1720: 69 6c 65 42 5f 2c 20 70 69 6c 65 42 5f 2b 73 69 ileB_, pileB_+si
1730: 7a 65 6f 66 28 70 69 6c 65 42 5f 29 2f 73 69 7a zeof(pileB_)/siz
1740: 65 6f 66 28 2a 70 69 6c 65 42 5f 29 29 3b 20 0d eof(*pileB_)); .
1750: 0a 09 69 6e 74 20 6b 20 3d 20 32 3b 20 0d 0a 09 ..int k = 2; ...
1760: 64 6f 75 62 6c 65 20 5f 20 3d 20 31 2e 33 39 37 double _ = 1.397
1770: 32 36 30 32 37 33 39 37 32 36 38 32 37 3b 20 0d 2602739726827; .
1780: 0a 45 4e 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09 .END..CASE(3)...
1790: 64 6f 75 62 6c 65 20 70 69 6c 65 41 5f 5b 5d 20 double pileA_[]
17a0: 3d 20 7b 38 35 2e 33 30 32 2c 20 39 32 2e 37 39 = {85.302, 92.79
17b0: 38 2c 20 37 36 2e 38 31 33 2c 20 33 37 2e 39 39 8, 76.813, 37.99
17c0: 34 2c 20 33 36 2e 37 33 37 2c 20 39 38 2e 36 35 4, 36.737, 98.65
17d0: 39 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 9};... vector <
17e0: 64 6f 75 62 6c 65 3e 20 70 69 6c 65 41 28 70 69 double> pileA(pi
17f0: 6c 65 41 5f 2c 20 70 69 6c 65 41 5f 2b 73 69 7a leA_, pileA_+siz
1800: 65 6f 66 28 70 69 6c 65 41 5f 29 2f 73 69 7a 65 eof(pileA_)/size
1810: 6f 66 28 2a 70 69 6c 65 41 5f 29 29 3b 20 0d 0a of(*pileA_)); ..
1820: 09 64 6f 75 62 6c 65 20 70 69 6c 65 42 5f 5b 5d .double pileB_[]
1830: 20 3d 20 7b 31 33 2e 33 35 32 2c 20 37 2e 33 30 = {13.352, 7.30
1840: 39 34 2c 20 35 34 2e 37 36 31 2c 20 38 2e 32 37 94, 54.761, 8.27
1850: 30 36 2c 20 36 33 2e 32 32 33 2c 20 33 37 2e 34 06, 63.223, 37.4
1860: 38 36 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 86};... vector
1870: 3c 64 6f 75 62 6c 65 3e 20 70 69 6c 65 42 28 70 <double> pileB(p
1880: 69 6c 65 42 5f 2c 20 70 69 6c 65 42 5f 2b 73 69 ileB_, pileB_+si
1890: 7a 65 6f 66 28 70 69 6c 65 42 5f 29 2f 73 69 7a zeof(pileB_)/siz
18a0: 65 6f 66 28 2a 70 69 6c 65 42 5f 29 29 3b 20 0d eof(*pileB_)); .
18b0: 0a 09 69 6e 74 20 6b 20 3d 20 33 3b 20 0d 0a 09 ..int k = 3; ...
18c0: 64 6f 75 62 6c 65 20 5f 20 3d 20 33 33 2e 35 38 double _ = 33.58
18d0: 36 30 33 38 38 39 38 33 36 31 37 35 3b 20 0d 0a 603889836175; ..
18e0: 45 4e 44 0d 0a 43 41 53 45 28 34 29 0d 0a 09 64 END..CASE(4)...d
18f0: 6f 75 62 6c 65 20 70 69 6c 65 41 5f 5b 5d 20 3d ouble pileA_[] =
1900: 20 7b 31 2e 30 2c 20 32 2e 30 2c 20 33 2e 30 2c {1.0, 2.0, 3.0,
1910: 20 34 2e 30 2c 20 35 2e 30 2c 20 36 2e 30 2c 20 4.0, 5.0, 6.0,
1920: 37 2e 30 2c 20 38 2e 30 2c 20 39 2e 30 2c 20 31 7.0, 8.0, 9.0, 1
1930: 30 2e 30 2c 20 31 31 2e 30 2c 20 31 32 2e 30 2c 0.0, 11.0, 12.0,
1940: 20 31 33 2e 30 2c 20 31 34 2e 30 2c 20 31 35 2e 13.0, 14.0, 15.
1950: 30 2c 20 31 36 2e 30 2c 20 31 37 2e 30 2c 20 31 0, 16.0, 17.0, 1
1960: 38 2e 30 2c 20 31 39 2e 30 2c 20 32 30 2e 30 2c 8.0, 19.0, 20.0,
1970: 20 32 31 2e 30 2c 20 32 32 2e 30 2c 20 32 33 2e 21.0, 22.0, 23.
1980: 30 2c 20 32 34 2e 30 2c 20 32 35 2e 30 2c 20 32 0, 24.0, 25.0, 2
1990: 36 2e 30 2c 20 32 37 2e 30 2c 20 32 38 2e 30 2c 6.0, 27.0, 28.0,
19a0: 20 32 39 2e 30 2c 20 33 30 2e 30 2c 20 33 31 2e 29.0, 30.0, 31.
19b0: 30 2c 20 33 32 2e 30 2c 20 33 33 2e 30 2c 20 33 0, 32.0, 33.0, 3
19c0: 34 2e 30 2c 20 33 35 2e 30 2c 20 33 36 2e 30 2c 4.0, 35.0, 36.0,
19d0: 20 33 37 2e 30 2c 20 33 38 2e 30 2c 20 33 39 2e 37.0, 38.0, 39.
19e0: 30 2c 20 34 30 2e 30 2c 20 34 31 2e 30 2c 20 34 0, 40.0, 41.0, 4
19f0: 32 2e 30 2c 20 34 33 2e 30 2c 20 34 34 2e 30 2c 2.0, 43.0, 44.0,
1a00: 20 34 35 2e 30 2c 20 34 36 2e 30 2c 20 34 37 2e 45.0, 46.0, 47.
1a10: 30 2c 20 34 38 2e 30 2c 20 34 39 2e 30 2c 20 35 0, 48.0, 49.0, 5
1a20: 30 2e 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 0.0};... vector
1a30: 20 3c 64 6f 75 62 6c 65 3e 20 70 69 6c 65 41 28 <double> pileA(
1a40: 70 69 6c 65 41 5f 2c 20 70 69 6c 65 41 5f 2b 73 pileA_, pileA_+s
1a50: 69 7a 65 6f 66 28 70 69 6c 65 41 5f 29 2f 73 69 izeof(pileA_)/si
1a60: 7a 65 6f 66 28 2a 70 69 6c 65 41 5f 29 29 3b 20 zeof(*pileA_));
1a70: 0d 0a 09 64 6f 75 62 6c 65 20 70 69 6c 65 42 5f ...double pileB_
1a80: 5b 5d 20 3d 20 7b 31 2e 30 2c 20 32 2e 30 2c 20 [] = {1.0, 2.0,
1a90: 33 2e 30 2c 20 34 2e 30 2c 20 35 2e 30 2c 20 36 3.0, 4.0, 5.0, 6
1aa0: 2e 30 2c 20 37 2e 30 2c 20 38 2e 30 2c 20 39 2e .0, 7.0, 8.0, 9.
1ab0: 30 2c 20 31 30 2e 30 2c 20 31 31 2e 30 2c 20 31 0, 10.0, 11.0, 1
1ac0: 32 2e 30 2c 20 31 33 2e 30 2c 20 31 34 2e 30 2c 2.0, 13.0, 14.0,
1ad0: 20 31 35 2e 30 2c 20 31 36 2e 30 2c 20 31 37 2e 15.0, 16.0, 17.
1ae0: 30 2c 20 31 38 2e 30 2c 20 31 39 2e 30 2c 20 32 0, 18.0, 19.0, 2
1af0: 30 2e 30 2c 20 32 31 2e 30 2c 20 32 32 2e 30 2c 0.0, 21.0, 22.0,
1b00: 20 32 33 2e 30 2c 20 32 34 2e 30 2c 20 32 35 2e 23.0, 24.0, 25.
1b10: 30 2c 20 32 36 2e 30 2c 20 32 37 2e 30 2c 20 32 0, 26.0, 27.0, 2
1b20: 38 2e 30 2c 20 32 39 2e 30 2c 20 33 30 2e 30 2c 8.0, 29.0, 30.0,
1b30: 20 33 31 2e 30 2c 20 33 32 2e 30 2c 20 33 33 2e 31.0, 32.0, 33.
1b40: 30 2c 20 33 34 2e 30 2c 20 33 35 2e 30 2c 20 33 0, 34.0, 35.0, 3
1b50: 36 2e 30 2c 20 33 37 2e 30 2c 20 33 38 2e 30 2c 6.0, 37.0, 38.0,
1b60: 20 33 39 2e 30 2c 20 34 30 2e 30 2c 20 34 31 2e 39.0, 40.0, 41.
1b70: 30 2c 20 34 32 2e 30 2c 20 34 33 2e 30 2c 20 34 0, 42.0, 43.0, 4
1b80: 34 2e 30 2c 20 34 35 2e 30 2c 20 34 36 2e 30 2c 4.0, 45.0, 46.0,
1b90: 20 34 37 2e 30 2c 20 34 38 2e 30 2c 20 34 39 2e 47.0, 48.0, 49.
1ba0: 30 2c 20 35 30 2e 30 7d 3b 0d 0a 09 20 20 76 65 0, 50.0};... ve
1bb0: 63 74 6f 72 20 3c 64 6f 75 62 6c 65 3e 20 70 69 ctor <double> pi
1bc0: 6c 65 42 28 70 69 6c 65 42 5f 2c 20 70 69 6c 65 leB(pileB_, pile
1bd0: 42 5f 2b 73 69 7a 65 6f 66 28 70 69 6c 65 42 5f B_+sizeof(pileB_
1be0: 29 2f 73 69 7a 65 6f 66 28 2a 70 69 6c 65 42 5f )/sizeof(*pileB_
1bf0: 29 29 3b 20 0d 0a 09 69 6e 74 20 6b 20 3d 20 35 )); ...int k = 5
1c00: 30 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 0; ...double _ =
1c10: 20 31 36 2e 38 34 36 39 33 38 37 37 35 35 31 30 16.846938775510
1c20: 32 30 33 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 203; ..END..CASE
1c30: 28 35 29 0d 0a 09 64 6f 75 62 6c 65 20 70 69 6c (5)...double pil
1c40: 65 41 5f 5b 5d 20 3d 20 7b 31 2e 30 2c 20 32 2e eA_[] = {1.0, 2.
1c50: 30 2c 20 33 2e 30 2c 20 34 2e 30 2c 20 35 2e 30 0, 3.0, 4.0, 5.0
1c60: 2c 20 36 2e 30 2c 20 37 2e 30 2c 20 38 2e 30 2c , 6.0, 7.0, 8.0,
1c70: 20 39 2e 30 2c 20 31 30 2e 30 2c 20 31 31 2e 30 9.0, 10.0, 11.0
1c80: 2c 20 31 32 2e 30 2c 20 31 33 2e 30 2c 20 31 34 , 12.0, 13.0, 14
1c90: 2e 30 2c 20 31 35 2e 30 2c 20 31 36 2e 30 2c 20 .0, 15.0, 16.0,
1ca0: 31 37 2e 30 2c 20 31 38 2e 30 2c 20 31 39 2e 30 17.0, 18.0, 19.0
1cb0: 2c 20 32 30 2e 30 2c 20 32 31 2e 30 2c 20 32 32 , 20.0, 21.0, 22
1cc0: 2e 30 2c 20 32 33 2e 30 2c 20 32 34 2e 30 2c 20 .0, 23.0, 24.0,
1cd0: 32 35 2e 30 2c 20 32 36 2e 30 2c 20 32 37 2e 30 25.0, 26.0, 27.0
1ce0: 2c 20 32 38 2e 30 2c 20 32 39 2e 30 2c 20 33 30 , 28.0, 29.0, 30
1cf0: 2e 30 2c 20 33 31 2e 30 2c 20 33 32 2e 30 2c 20 .0, 31.0, 32.0,
1d00: 33 33 2e 30 2c 20 33 34 2e 30 2c 20 33 35 2e 30 33.0, 34.0, 35.0
1d10: 2c 20 33 36 2e 30 2c 20 33 37 2e 30 2c 20 33 38 , 36.0, 37.0, 38
1d20: 2e 30 2c 20 33 39 2e 30 2c 20 34 30 2e 30 2c 20 .0, 39.0, 40.0,
1d30: 34 31 2e 30 2c 20 34 32 2e 30 2c 20 34 33 2e 30 41.0, 42.0, 43.0
1d40: 2c 20 34 34 2e 30 2c 20 34 35 2e 30 2c 20 34 36 , 44.0, 45.0, 46
1d50: 2e 30 2c 20 34 37 2e 30 2c 20 34 38 2e 30 2c 20 .0, 47.0, 48.0,
1d60: 34 39 2e 30 2c 20 35 30 2e 30 7d 3b 0d 0a 09 20 49.0, 50.0};...
1d70: 20 76 65 63 74 6f 72 20 3c 64 6f 75 62 6c 65 3e vector <double>
1d80: 20 70 69 6c 65 41 28 70 69 6c 65 41 5f 2c 20 70 pileA(pileA_, p
1d90: 69 6c 65 41 5f 2b 73 69 7a 65 6f 66 28 70 69 6c ileA_+sizeof(pil
1da0: 65 41 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 69 6c eA_)/sizeof(*pil
1db0: 65 41 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 eA_)); ...double
1dc0: 20 70 69 6c 65 42 5f 5b 5d 20 3d 20 7b 35 31 2e pileB_[] = {51.
1dd0: 30 2c 20 35 32 2e 30 2c 20 35 33 2e 30 2c 20 35 0, 52.0, 53.0, 5
1de0: 34 2e 30 2c 20 35 35 2e 30 2c 20 35 36 2e 30 2c 4.0, 55.0, 56.0,
1df0: 20 35 37 2e 30 2c 20 35 38 2e 30 2c 20 35 39 2e 57.0, 58.0, 59.
1e00: 30 2c 20 36 30 2e 30 2c 20 36 31 2e 30 2c 20 36 0, 60.0, 61.0, 6
1e10: 32 2e 30 2c 20 36 33 2e 30 2c 20 36 34 2e 30 2c 2.0, 63.0, 64.0,
1e20: 20 36 35 2e 30 2c 20 36 36 2e 30 2c 20 36 37 2e 65.0, 66.0, 67.
1e30: 30 2c 20 36 38 2e 30 2c 20 36 39 2e 30 2c 20 37 0, 68.0, 69.0, 7
1e40: 30 2e 30 2c 20 37 31 2e 30 2c 20 37 32 2e 30 2c 0.0, 71.0, 72.0,
1e50: 20 37 33 2e 30 2c 20 37 34 2e 30 2c 20 37 35 2e 73.0, 74.0, 75.
1e60: 30 2c 20 37 36 2e 30 2c 20 37 37 2e 30 2c 20 37 0, 76.0, 77.0, 7
1e70: 38 2e 30 2c 20 37 39 2e 30 2c 20 38 30 2e 30 2c 8.0, 79.0, 80.0,
1e80: 20 38 31 2e 30 2c 20 38 32 2e 30 2c 20 38 33 2e 81.0, 82.0, 83.
1e90: 30 2c 20 38 34 2e 30 2c 20 38 35 2e 30 2c 20 38 0, 84.0, 85.0, 8
1ea0: 36 2e 30 2c 20 38 37 2e 30 2c 20 38 38 2e 30 2c 6.0, 87.0, 88.0,
1eb0: 20 38 39 2e 30 2c 20 39 30 2e 30 2c 20 39 31 2e 89.0, 90.0, 91.
1ec0: 30 2c 20 39 32 2e 30 2c 20 39 33 2e 30 2c 20 39 0, 92.0, 93.0, 9
1ed0: 34 2e 30 2c 20 39 35 2e 30 2c 20 39 36 2e 30 2c 4.0, 95.0, 96.0,
1ee0: 20 39 37 2e 30 2c 20 39 38 2e 30 2c 20 39 39 2e 97.0, 98.0, 99.
1ef0: 30 2c 20 31 30 30 2e 30 7d 3b 0d 0a 09 20 20 76 0, 100.0};... v
1f00: 65 63 74 6f 72 20 3c 64 6f 75 62 6c 65 3e 20 70 ector <double> p
1f10: 69 6c 65 42 28 70 69 6c 65 42 5f 2c 20 70 69 6c ileB(pileB_, pil
1f20: 65 42 5f 2b 73 69 7a 65 6f 66 28 70 69 6c 65 42 eB_+sizeof(pileB
1f30: 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 69 6c 65 42 _)/sizeof(*pileB
1f40: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 6b 20 3d 20 _)); ...int k =
1f50: 35 30 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 50; ...double _
1f60: 3d 20 32 31 2e 31 32 38 31 34 34 31 38 36 39 36 = 21.12814418696
1f70: 37 37 31 37 3b 20 0d 0a 45 4e 44 0d 0a 7d 0d 0a 7717; ..END..}..
1f80: 2f 2f 20 45 4e 44 20 43 55 54 20 48 45 52 45 0d // END CUT HERE.
1f90: 0a 0d 0a ...