이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
using namespace std;
#define lg long long
#define fastio ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0);
const int N = 5e4+5, M = 305, K = (5e4+5)/(305)+1000;
vector<int> adj[N];
int loc[N], stop[N], stops = 0;
long long dist[N][K], dist2[N][K], sto[K], doub[N], sing[N], ans2[N], dist3[N][2], gh[N], inc[N];
bitset<K> vis[N], visd;
int main()
{
fastio;
int n, m, k;
cin >> n >> m >> k;
for(int i = 0; i < m; i++)
{
int u, v;
cin >> u >> v;
adj[u].push_back(v);
adj[v].push_back(u);
}
string s;
cin >> s;
vector<int> st;
for(int i = 1; i <= n; i++)
{
stop[i] = (s[i-1] == '1');
stops += stop[i];
if(stop[i]) st.push_back(i);
}
for(int i = 0; i < k; i++)
{
cin >> loc[i];
}
if(!stops)
{
memset(ans2, 0x3f, sizeof(ans2));
queue<int>q;
visd = 0;
visd[loc[0]] = 1;
ans2[loc[0]] = 0;
q.push(loc[0]);
while(q.size())
{
int u = q.front();
q.pop();
if(stop[u] && u != loc[0]) continue;
for(auto it : adj[u])
{
if(ans2[it] > ans2[u]+1)
{
ans2[it] = ans2[u]+1;
q.push(it);
}
}
}
for(int i = 1; i <= n; i++) cout << ans2[i] << '\n';
return 0;
}
if(k > M)
{
memset(dist, 0x3f, sizeof(dist));
memset(dist2, 0x3f, sizeof(dist2));
stops = min(stops, K-5);
queue<int> q[stops+3];
q[0].push(loc[0]);
dist[loc[0]][0] = 0;
vis[loc[0]][0] = true;
for(int i = 0; i <= stops+2; i++)
{
int its = 0;
while(q[i].size())
{
int u = q[i].front();
q[i].pop();
vis[u][i] = false;
int newOg = i;
if(stop[u] && (its || i))
{
newOg++;
}
if(newOg > stops+2) continue;
for(auto it : adj[u])
{
if(dist[it][newOg] > dist[u][i]+1)
{
if(!vis[it][newOg]) q[newOg].push(it);
vis[it][newOg] = true;
dist[it][newOg] = dist[u][i]+1;
}
}
its++;
}
}
for(int i = 1; i <= n; i++) vis[i] = 0, dist2[i][0] = 0;
for(auto it : st) q[0].push(it), vis[it][0] = true;
for(int i = 0; i <= stops+2; i++)
{
while(q[i].size())
{
int u = q[i].front();
q[i].pop();
vis[u][i] = false;
int newOg = i;
if(stop[u])
{
newOg++;
}
if(newOg > stops+2) continue;
for(auto it : adj[u])
{
if(dist2[it][newOg] > dist2[u][i]+1)
{
if(!vis[it][newOg]) q[newOg].push(it);
vis[it][newOg] = true;
dist2[it][newOg] = dist2[u][i]+1;
}
}
}
}
for(int i = 1; i < k; i++)
{
for(int j = 0; j <= stops+2; j++)
{
if(dist2[loc[i]][j] > 1e18)
{
sto[j] = -1;
break;
}
sto[j] += dist2[loc[i]][j];
}
}
for(int i = 1; i <= n; i++)
{
long long ans = dist[i][0];
for(int j = 1; j <= stops+2; j++)
{
if(sto[j] == -1) continue;
ans = min(ans, dist[i][j]+sto[j]);
}
cout << ans << '\n';
}
return 0;
}
memset(sing, 0x3f, sizeof(sing));
memset(doub, 0x3f, sizeof(doub));
queue<int> q;
for(auto it : st)
{
q.push(it);
sing[it] = 0;
}
while(q.size())
{
int u = q.front();
q.pop();
for(auto it : adj[u])
{
if(sing[it] > sing[u]+1)
{
sing[it] = sing[u]+1;
q.push(it);
}
}
}
for(int i = 1; i <= n; i++)
{
sing[i] = 2*(!sing[i])+sing[i];
}
for(int i = 1; i < k; i++) gh[0] += sing[loc[i]];
deque<int> dq;
for(auto it : st)
{
for(auto it2 : adj[it])
{
if(stop[it2])
{
dq.push_back(it);
doub[it] = 0;
break;
}
}
}
while(dq.size())
{
int u = dq.front();
dq.pop_front();
for(auto it : adj[u])
{
if(doub[it] > doub[u]+(!stop[u]))
{
doub[it] = doub[u]+(!stop[u]);
if(stop[u]) dq.push_front(it);
else dq.push_back(it);
}
}
}
bool bad = false;
for(int i = 1; i < k; i++)
{
if(doub[loc[i]] >= 1e18)
{
bad = true;
break;
}
inc[i] = doub[loc[i]]-sing[loc[i]];
}
if(!bad)
{
sort(inc+1, inc+k);
for(int i = 1; i < k; i++)
{
gh[i] = gh[i-1]+inc[i];
}
}
memset(ans2, 0x3f, sizeof(ans2));
visd = 0;
visd[loc[0]] = 1;
ans2[loc[0]] = 0;
q.push(loc[0]);
while(q.size())
{
int u = q.front();
q.pop();
if(stop[u] && u != loc[0]) continue;
for(auto it : adj[u])
{
if(ans2[it] > ans2[u]+1)
{
ans2[it] = ans2[u]+1;
q.push(it);
}
}
}
for(int its = 0; its < k; its++)
{
if(bad && its) break;
memset(dist3, 0x3f, sizeof(dist3));
dist3[loc[0]][0] = 0;
vis[loc[0]][0] = true;
priority_queue<array<lg, 3>> pq;
pq.push({0, loc[0], 0});
while(pq.size())
{
int u = pq.top()[1];
int ty = pq.top()[2];
vis[u][ty] = 0;
int stopp = (u != loc[0] && stop[u]);
int newty = (ty||stopp);
pq.pop();
for(auto it : adj[u])
{
if(dist3[it][newty] > dist3[u][ty]+(stopp ? 2*(k-1-its)+its+1 : 1))
{
dist3[it][newty] = dist3[u][ty]+(stopp ? 2*(k-1-its)+its+1 : 1);
if(!vis[it][newty])
{
pq.push({-dist3[it][newty], it, newty});
}
vis[it][newty] = true;
}
}
}
for(int i = 1; i <= n; i++)
{
ans2[i] = min(ans2[i], dist3[i][1]+gh[its]-2*(k-1-its));
}
}
for(int i = 1; i <= n; i++) cout << ans2[i] << '\n';
return 0;
}
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |