//#pragma GCC optimize("Ofast") //#pragma GCC target("avx") //#undef LOCAL #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using uint = unsigned int; using ll = long long; using ull = unsigned long long; constexpr ll TEN(int n) { return (n == 0) ? 1 : 10 * TEN(n - 1); } template using V = vector; template using VV = V>; #include struct Scanner { int fd = -1; char line[(1 << 15) + 1]; size_t st = 0, ed = 0; void reread() { memmove(line, line + st, ed - st); ed -= st; st = 0; ed += ::read(fd, line + ed, (1 << 15) - ed); line[ed] = '\0'; } bool succ() { while (true) { if (st == ed) { reread(); if (st == ed) return false; } while (st != ed && isspace(line[st])) st++; if (st != ed) break; } if (ed - st <= 50) { bool sep = false; for (size_t i = st; i < ed; i++) { if (isspace(line[i])) { sep = true; break; } } if (!sep) reread(); } return true; } template ::value, int> = 0> bool read_single(T& ref) { if (!succ()) return false; while (true) { size_t sz = 0; while (st + sz < ed && !isspace(line[st + sz])) sz++; ref.append(line + st, sz); st += sz; if (!sz || st != ed) break; reread(); } return true; } template ::value, int> = 0> bool read_single(T& ref) { if (!succ()) return false; bool neg = false; if (line[st] == '-') { neg = true; st++; } ref = T(0); while (isdigit(line[st])) { ref = 10 * ref + (line[st++] & 0xf); } if (neg) ref = -ref; return true; } template bool read_single(V& ref) { for (auto& d : ref) { if (!read_single(d)) return false; } return true; } void read() {} template void read(H& h, T&... t) { bool f = read_single(h); assert(f); read(t...); } Scanner(FILE* fp) : fd(fileno(fp)) {} }; struct Printer { public: template void write() {} template void write(const H& h, const T&... t) { if (F) write_single(' '); write_single(h); write(t...); } template void writeln(const T&... t) { write(t...); write_single('\n'); } Printer(FILE* _fp) : fp(_fp) {} ~Printer() { flush(); } private: static constexpr size_t SIZE = 1 << 15; FILE* fp; char line[SIZE], small[50]; size_t pos = 0; void flush() { fwrite(line, 1, pos, fp); pos = 0; } void write_single(const char& val) { if (pos == SIZE) flush(); line[pos++] = val; } template ::value, int> = 0> void write_single(T val) { if (pos > (1 << 15) - 50) flush(); if (val == 0) { write_single('0'); return; } if (val < 0) { write_single('-'); val = -val; // todo min } size_t len = 0; while (val) { small[len++] = char(0x30 | (val % 10)); val /= 10; } for (size_t i = 0; i < len; i++) { line[pos + i] = small[len - 1 - i]; } pos += len; } void write_single(const string& s) { for (char c : s) write_single(c); } void write_single(const char* s) { size_t len = strlen(s); for (size_t i = 0; i < len; i++) write_single(s[i]); } template void write_single(const V& val) { auto n = val.size(); for (size_t i = 0; i < n; i++) { if (i) write_single(' '); write_single(val[i]); } } }; template struct MinDist { V dist; V from; }; template MinDist mindist(const VV& g, int s, D inf = numeric_limits::max()) { int n = (int)g.size(); V dist = V(n, inf); V from = V(n); struct P { D key; int to; bool operator<(P r) const { return key > r.key; } }; priority_queue

q; q.push(P{0, s}); dist[s] = D(0); while (!q.empty()) { P p = q.top(); q.pop(); if (dist[p.to] < p.key) continue; for (E e : g[p.to]) { if (p.key + e.dist < dist[e.to]) { dist[e.to] = p.key + e.dist; from[e.to] = p.to; q.push(P{dist[e.to], e.to}); } } } return MinDist{dist, from}; } Scanner sc = Scanner(stdin); Printer pr = Printer(stdout); using D = double; struct E { int to; D dist; }; ostream& operator<<(ostream& os, const E& e) { return os << "E(" << e.to << ", " << e.dist << ")"; }; using B = bitset<2020>; struct ST { B used; int pos; D now, min; bool operator<(const ST& r) const { return min < r.min; } bool operator>(const ST& r) const { return min > r.min; } }; #include #include #include struct Random { private: // Use xoshiro256** // Refereces: http://xoshiro.di.unimi.it/xoshiro256starstar.c static uint64_t rotl(const uint64_t x, int k) { return (x << k) | (x >> (64 - k)); } std::array s; uint64_t next() { const uint64_t result_starstar = rotl(s[1] * 5, 7) * 9; const uint64_t t = s[1] << 17; s[2] ^= s[0]; s[3] ^= s[1]; s[1] ^= s[2]; s[0] ^= s[3]; s[2] ^= t; s[3] = rotl(s[3], 45); return result_starstar; } // random choice from [0, upper] uint64_t next(uint64_t upper) { if (!(upper & (upper + 1))) { // b = 00..0011..11 return next() & upper; } int lg = 63 - __builtin_clzll(upper); uint64_t mask = (lg == 63) ? ~0ULL : (1ULL << (lg + 1)) - 1; while (true) { uint64_t r = next() & mask; if (r <= upper) return r; } } public: Random(uint64_t seed = 0) { // Use splitmix64 // Reference: http://xoshiro.di.unimi.it/splitmix64.c for (int i = 0; i < 4; i++) { uint64_t z = (seed += 0x9e3779b97f4a7c15); z = (z ^ (z >> 30)) * 0xbf58476d1ce4e5b9; z = (z ^ (z >> 27)) * 0x94d049bb133111eb; s[i] = z ^ (z >> 31); } } // random choice from [lower, upper] template T uniform(T lower, T upper) { assert(lower <= upper); return T(lower + next(uint64_t(upper - lower))); } bool uniform_bool() { return uniform(0, 1) == 1; } double uniform01() { uint64_t v = next(1ULL << 63); return double(v) / (1ULL << 63); } template std::pair uniform_pair(T lower, T upper) { assert(upper - lower >= 1); T a, b; do { a = uniform(lower, upper); b = uniform(lower, upper); } while (a == b); if (a > b) std::swap(a, b); return {a, b}; } // generate random lower string that length = n std::string lower_string(size_t n) { std::string res = ""; for (size_t i = 0; i < n; i++) { res += uniform('a', 'z'); } return res; } // random shuffle template void shuffle(Iter first, Iter last) { if (first == last) return; // Reference and edit: // cpprefjp - C++日本語リファレンス // (https://cpprefjp.github.io/reference/algorithm/shuffle.html) int len = 1; for (auto it = first + 1; it != last; it++) { len++; int j = uniform(0, len - 1); if (j != len - 1) iter_swap(it, first + j); } } // generate random permutation that length = n template std::vector perm(size_t n) { std::vector idx(n); std::iota(idx.begin(), idx.end(), T(0)); shuffle(idx.begin(), idx.end()); return idx; } template std::vector choice(size_t n, T lower, T upper) { assert(n <= upper - lower + 1); std::set res; while (res.size() < n) res.insert(uniform(lower, upper)); return {res.begin(), res.end()}; } }; Random& global_gen() { static Random gen; return gen; } Random get_random_gen() { return Random(chrono::steady_clock::now().time_since_epoch().count()); } Random& global_runtime_gen() { static Random gen = get_random_gen(); return gen; } int main() { int n, m, k; scanf("%d %d %d", &n, &m, &k); assert(2 <= n && n <= 2000); assert(1 <= m && m <= 2000); assert(1 <= k && k <= 10); //n = 1000; //m = 2000; //k = 10; int s, t; scanf("%d %d", &s, &t); s--; t--; //s = 0; t = n - 1; V x(n), y(n); for (int i = 0; i < n; i++) { scanf("%lf %lf", &(x[i]), &(y[i])); //x[i] = global_gen().uniform(-1000, 1000); //y[i] = global_gen().uniform(-1000, 1000); x[i] += global_gen().uniform01() / TEN(9); y[i] += global_gen().uniform01() / TEN(9); } VV g(n); for (int i = 0; i < m; i++) { int a, b; scanf("%d %d", &a, &b); a--; b--; /*do { a = global_gen().uniform(0, n - 1); b = global_gen().uniform(0, n - 1); } while (a == b);*/ D di = hypot(abs(x[a] - x[b]), abs(y[a] - y[b])); g[a].push_back({b, di}); g[b].push_back({a, di}); } auto dist = mindist(g, s).dist; priority_queue, greater> que; { B fi = B(); fi[s] = true; que.push(ST{fi, s, 0, dist[t]}); } V ans; while (que.size()) { ST st = que.top(); que.pop(); if (st.pos == t) { ans.push_back(st.now); if (int(ans.size()) == k) break; } for (auto e: g[st.pos]) { if (st.used[e.to]) continue; ST nst = st; nst.used[e.to] = true; nst.pos = e.to; nst.now += e.dist; VV ng(n); for (int i = 0; i < n; i++) { for (auto ne: g[i]) { if (nst.used[ne.to]) continue; ng[i].push_back(ne); } } nst.min = nst.now + mindist(ng, nst.pos).dist[t]; if (nst.min < 1e30) que.push(nst); } } int l = int(ans.size()); for (int i = 0; i < l; i++) { printf("%.20lf\n", ans[i]); } for (int i = l; i < k; i++) { printf("-1\n"); } return 0; }