#include <bits/stdc++.h>
#define ll long long
#define pb push_back
#define s second
#define f first
#define pii pair <int,int>
#define left h*2,l,(l + r)/2
#define right h*2+1,(l + r)/2 + 1,r
// #define int ll
using namespace std;
const int N = 2e5 + 5,inf = 1e18;
int arr[N];
ll dist[8][N],dd[37][N],ddd[N];
ll DD[37][N],val[37][N];
vector <pii> v[N];
signed main (){
ios_base::sync_with_stdio(0),cin.tie(NULL),cout.tie(NULL);
int n,k,m;
cin>>n>>k>>m;
for (int i = 1; i <= k; i++)
cin >> arr[i];
for (int i =1; i <= m; i++){
int a,b,c;
cin>>a>>b>>c;
v[a].pb({b,c});
v[b].pb({a,c});
}
for (int i = 1; i <= k; i++){
for (int j = 1; j <= n; j++)
dist[i][j] = inf;
priority_queue<pair <ll,int>> pq;
dist[i][arr[i]] = 0;
pq.push({0,arr[i]});
while (pq.size()){
ll D = -pq.top().f; int x = pq.top().s;
pq.pop();
if (dist[i][x] != D) continue;
for (auto [to,w]: v[x]){
if (dist[i][to] > D + w){
dist[i][to] = D + w;
pq.push({-dist[i][to],to});
}
}
}
}
ll ans=inf;
for (int i = 1; i <= n; i++){
ll cur=0;
for (int j = 1; j <= k; j++)
cur += dist[j][i];
ans=min(ans,cur);
}
if (k < 4){
cout<<ans<<endl;
exit(0);
}
for (int mask1=0;mask1<(1<<k); mask1++)
for(int i=1;i<=n;i++)
DD[mask1][i] = inf;
for (int mask2=0;mask2<(1<<k); mask2++){
for (int i = 1; i <= n; i++){
for (int j=0;j<k;j++){
if (((mask2>>j) & 1)) val[mask2][i] += (dist[j + 1][i]);
}
}
}
for (int mask1 = 0; mask1 < (1<<k); mask1++){
for (int i = 1; i <= n; i++){
dd[mask1][i] = 0;
for (int j = 0; j < k; j++){
if (((mask1>>j) & 1)) dd[mask1][i] += dist[j + 1][i];
}
}
priority_queue <pair <ll,int>> pq;
for (int i = 1; i <= n; i++)
pq.push({-dd[mask1][i],i});
while (pq.size()){
ll D = -pq.top().f; int x = pq.top().s;
pq.pop();
if (dd[mask1][x] != D) continue;
for (auto [to,w]: v[x]){
if (dd[mask1][to] > D + w){
dd[mask1][to] = D + w;
pq.push({-dd[mask1][to],to});
}
}
}
for (int mask2 = 0; mask2 < (1<<k); mask2++){
int y = (mask1 & mask2);
if (y > 0) continue;
for (int i = 1; i <= n; i++){
ll kk = dd[mask1][i] + val[mask2][i];
DD[(mask1 | mask2)][i] = min(DD[(mask1 | mask2)][i],kk);
}
}
}
for (int mask1=0;mask1<(1<<k); mask1++){
for (int i =1 ; i <= n; i++){
ddd[i] = DD[mask1][i];
}
priority_queue <pair <ll,int>> pq;
for (int i = 1; i <= n; i++)
pq.push({-ddd[i],i});
while (pq.size()){
ll D = -pq.top().f; int x = pq.top().s;
pq.pop();
if (ddd[x] != D) continue;
for (auto [to,w]: v[x]){
if (ddd[to] > D + w){
ddd[to] = D + w;
pq.push({-ddd[to],to});
}
}
}
for (int i = 1; i <= n; i++){
for (int j =0;j < k; j++){
if (((mask1>>j) & 1)) continue;
ddd[i] += (dist[j + 1][i]);
}
ans=min(ans,ddd[i]);
}
}
cout<<ans;
}
Compilation message
cities.cpp:12:29: warning: overflow in conversion from 'double' to 'int' changes value from '1.0e+18' to '2147483647' [-Woverflow]
12 | const int N = 2e5 + 5,inf = 1e18;
| ^~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
12892 KB |
Output is correct |
2 |
Correct |
1 ms |
12892 KB |
Output is correct |
3 |
Correct |
2 ms |
12892 KB |
Output is correct |
4 |
Correct |
6 ms |
68268 KB |
Output is correct |
5 |
Incorrect |
7 ms |
74840 KB |
Output isn't correct |
6 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
60 ms |
24660 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
12 ms |
74588 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
476 ms |
98296 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
1054 ms |
135240 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |