Submission #927089

#TimeUsernameProblemLanguageResultExecution timeMemory
927089socpite한자 끝말잇기 (JOI14_kanji)C++14
100 / 100
129 ms14612 KiB
#include "Annalib.h" #include<bits/stdc++.h> using namespace std; const long long INF = LLONG_MAX/2; const int maxn = 305; namespace anna{ vector<pair<int, int>> g[maxn]; vector<long long> shortest_path(int N, int s, long long C[]){ vector<long long> d(N, INF); vector<int> vis(N, 0); d[s] = 0; for(int i = 1; i <= N; i++){ long long best = INF; int x = -1; for(int j = 0; j < N; j++){ if(vis[j])continue; if(d[j] < best){ best = d[j]; x = j; } } if(x == -1)break; vis[x] = 1; for(auto v: g[x]) if(d[v.first] > C[v.second] + best)d[v.first] = C[v.second] + best; } return d; } void solve(int K, int Q, vector<long long>& forgor_C, vector<vector<long long>>& cost_to_T){ // eval(i, x) = forgor(x) + cost_to_T vector<vector<int>> queries(K+1); for(int i = 0; i < Q; i++)queries[0].push_back(i); unsigned long long total_case = 1; vector<pair<int, int>> cc; for(int b = 1; b <= K; b++){ for(int a = 0; a < b; a++){ vector<int> a_wins; for(auto v: queries[a]){ if(forgor_C[a] + cost_to_T[a][v] < forgor_C[b] + cost_to_T[b][v])a_wins.push_back(v); else queries[b].push_back(v); } int cases = queries[a].size() + 1, num = a_wins.size(); total_case *= cases; queries[a] = a_wins; cc.push_back({num, cases}); } } // cout << cc.size() << endl; unsigned long long crr = 0; while(!cc.empty()){ crr = crr*cc.back().second; crr = crr+cc.back().first; cc.pop_back(); } while(total_case > 1){ Tap(crr&1); total_case = (total_case+1)>>1; crr>>=1; } } } void Anna(int N, int M, int A[], int B[], long long C[], int Q, int S[], int T[], int K, int U[]) { vector<long long> forgor_C(K+1, 0); vector<vector<long long>> cost_to_T(K+1, vector<long long>(Q)); for(int i = 0; i < K; i++){ forgor_C[i] = C[U[i]]; } for(int i = 0; i < K; i++)C[U[i]] = INF; for(int i = 0; i < M; i++){ anna::g[A[i]].push_back({B[i], i}); } for(int i = 0; i < K; i++){ auto d = anna::shortest_path(N, B[U[i]], C); for(int j = 0; j < Q; j++)cost_to_T[i][j] = d[T[j]]; } for(int i = 0; i < Q; i++){ auto d = anna::shortest_path(N, S[i], C); for(int j = 0; j < K; j++)cost_to_T[j][i] += d[A[U[j]]]; cost_to_T[K][i] = d[T[i]]; } anna::solve(K, Q, forgor_C, cost_to_T); // eval(i, x) = forgor(x) + cost_to_A + cost_to_T }
#include "Brunolib.h" #include<bits/stdc++.h> using namespace std; const int maxn = 305; const long long INF = 1e16*maxn; namespace bruno { vector<pair<int, int>> g[maxn]; vector<long long> shortest_path(int N, int s, long long C[]){ vector<long long> d(N, INF); vector<int> vis(N, 0); d[s] = 0; for(int i = 1; i <= N; i++){ long long best = INF; int x = -1; for(int j = 0; j < N; j++){ if(vis[j])continue; if(d[j] < best){ best = d[j]; x = j; } } if(x == -1)break; vis[x] = 1; for(auto v: g[x]) if(d[v.first] > C[v.second] + best)d[v.first] = C[v.second] + best; } return d; } vector<int> path[60]; vector<int> solve_query(int N, int s, int t, long long C[]){ vector<long long> d(N, INF); vector<int> vis(N, 0); vector<pair<int, int>> trace(N); d[s] = 0; for(int i = 1; i <= N; i++){ long long best = INF; int x = -1; for(int j = 0; j < N; j++){ if(vis[j])continue; if(d[j] < best){ best = d[j]; x = j; } } if(x == -1)break; vis[x] = 1; for(auto v: g[x]) if(d[v.first] > C[v.second] + best){ d[v.first] = C[v.second] + best; trace[v.first] = {x, v.second}; } } vector<int> ans; while(t != s){ ans.push_back(trace[t].second); t = trace[t].first; } reverse(ans.begin(), ans.end()); return ans; } vector<vector<int>> solve(int K, int Q, vector<vector<long long>>& cost_to_T, int L, int X[]){ unsigned long long crr = 0; for(int i = 0; i < L; i++)crr ^= (1ULL*X[i])<<i; vector<vector<int>> queries(K+1); for(int i = 0; i < Q; i++)queries[0].push_back(i); for(int b = 1; b <= K; b++){ for(int a = 0; a < b; a++){ auto cmp = [&a, &b, &cost_to_T](int lhs, int rhs){ return cost_to_T[a][lhs] - cost_to_T[b][lhs] < cost_to_T[a][rhs] - cost_to_T[b][rhs]; }; int cases = queries[a].size() + 1; sort(queries[a].begin(), queries[a].end(), cmp); while(queries[a].size() > crr%cases){ queries[b].push_back(queries[a].back()); queries[a].pop_back(); } crr/=cases; } } return queries; } } void Bruno(int N, int M, int A[], int B[], long long C[], int Q, int S[], int T[], int K, int U[], int L, int X[]) { vector<vector<long long>> cost_to_T(K+1, vector<long long>(Q)); for(int i = 0; i < K; i++)C[U[i]] = INF; for(int i = 0; i < M; i++){ bruno::g[A[i]].push_back({B[i], i}); } for(int i = 0; i < K; i++){ auto d = bruno::shortest_path(N, B[U[i]], C); for(int j = 0; j < Q; j++)cost_to_T[i][j] = d[T[j]]; } for(int i = 0; i < Q; i++){ auto d = bruno::shortest_path(N, S[i], C); for(int j = 0; j < K; j++)cost_to_T[j][i] += d[A[U[j]]]; cost_to_T[K][i] = d[T[i]]; }; auto qtypes = bruno::solve(K, Q, cost_to_T, L, X); for(int i = 0; i <= K; i++){ if(i < K)C[U[i]] = 0; for(auto v: qtypes[i]){ bruno::path[v] = bruno::solve_query(N, S[v], T[v], C); if(i == K)for(auto ele: bruno::path[v]){ for(int j = 0; j < K; j++)assert(ele != U[j]); } } if(i < K)C[U[i]] = INF; } for(int i = 0; i < Q; i++){ for(auto v: bruno::path[i])Answer(v); Answer(-1); } }
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...