QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#259078 | #7790. 最短路求和 | zhouhuanyi | 10 | 39ms | 40892kb | C++14 | 5.0kb | 2023-11-20 16:21:06 | 2023-11-20 16:21:08 |
Judging History
answer
#include<iostream>
#include<cstdio>
#include<vector>
#include<map>
#include<queue>
#include<algorithm>
#define N 100000
#define M 8000
using namespace std;
const int inf=(int)(1e9);
int read()
{
char c=0;
int sum=0;
while (c<'0'||c>'9') c=getchar();
while ('0'<=c&&c<='9') sum=sum*10+c-'0',c=getchar();
return sum;
}
struct reads
{
int num,data,snum;
};
struct node
{
int x,y;
bool operator < (const node &t)const
{
return x!=t.x?x<t.x:y<t.y;
}
};
node st[N+1];
struct rds
{
int num,data;
bool operator < (const rds &t)const
{
return data>t.data;
}
};
int n,m,lg[N+1],tfa[N+1],num[N+1],length,leng,tong[N+1],fa[N+1],nt[N+1],nt2[N+1],rst[N+1],sz[N+1],szt[N+1],dque[N+1],top;
long long depth[N+1],delta[N+1],delta2[N+1],ans,dis[M+1][M+1];
bool used[N+1],vis[N+1],vst[N+1],dst[N+1],vist[N+1];
vector<reads>E[N+1];
vector<reads>ES[N+1];
vector<int>v[N+1];
vector<int>cnt[N+1];
vector<int>cnt2[N+1];
map<node,int>p;
void add(int x,int y,int z,int w)
{
E[x].push_back((reads){y,z,w}),E[y].push_back((reads){x,z,w});
return;
}
void add2(int x,int y,int z)
{
ES[x].push_back((reads){y,z}),ES[y].push_back((reads){x,z});
return;
}
void dfs(int x)
{
used[x]=1;
for (int i=0;i<E[x].size();++i)
if (!used[E[x][i].num])
depth[E[x][i].num]=depth[x]+E[x][i].data,fa[E[x][i].num]=x,vis[E[x][i].snum]=1,dfs(E[x][i].num);
return;
}
void dfs2(int x)
{
int cnt=0,y=0;
dst[x]=vst[x];
for (int i=0;i<E[x].size();++i)
if (fa[E[x][i].num]==x)
{
dfs2(E[x][i].num),dst[x]|=dst[E[x][i].num],cnt+=dst[E[x][i].num];
if (dst[E[x][i].num]) y=E[x][i].num;
}
if (cnt>=2) vst[x]=1;
else nt[x]=nt[y];
if (vst[x]) nt[x]=nt2[x]=x;
return;
}
void dfs3(int x)
{
if (dst[x]) rst[x]=x;
for (int i=0;i<E[x].size();++i)
if (fa[E[x][i].num]==x)
{
if (dst[E[x][i].num]) nt2[E[x][i].num]=nt2[x];
rst[E[x][i].num]=rst[x],dfs3(E[x][i].num);
}
return;
}
void dijkstra(int x)
{
int top;
priority_queue<rds>q;
for (int i=1;i<=length;++i) dis[x][i]=inf;
q.push((rds){x,0}),dis[x][x]=0;
while (!q.empty())
{
top=q.top().num,q.pop();
if (used[top]) continue;
used[top]=1;
for (int i=0;i<E[top].size();++i)
if (dis[x][E[top][i].num]>dis[x][top]+E[top][i].data)
dis[x][E[top][i].num]=dis[x][top]+E[top][i].data,q.push((rds){E[top][i].num,dis[x][E[top][i].num]});
}
return;
}
bool cmp(int x,int y)
{
return delta2[x]-delta[x]>delta2[y]-delta[y];
}
void dfs4(int x)
{
dque[++top]=x,vist[x]=szt[x]=1;
for (int i=0;i<E[x].size();++i)
if (!vist[E[x][i].num]&&nt[x]==nt[E[x][i].num]&&nt2[x]==nt2[E[x][i].num])
tfa[E[x][i].num]=E[x][i].data,dfs4(E[x][i].num),szt[x]+=szt[E[x][i].num];
return;
}
int main()
{
int x,y,z,ps;
long long d1,d2,d3,d4,ds1,ds2;
for (int i=2;i<=N;++i) lg[i]=lg[i>>1]+1;
n=read(),m=read();
for (int i=1;i<=m;++i) x=read(),y=read(),z=read(),add(x,y,z,i);
dfs(1);
for (int i=1;i<=n;++i)
for (int j=0;j<E[i].size();++j)
if (!vis[E[i][j].snum])
vst[i]=1;
dfs2(1),dfs3(1);
for (int i=1;i<=n;++i)
{
nt[i]=nt[rst[i]],nt2[i]=nt2[rst[i]];
if (nt[i]) delta[i]=depth[i]-depth[nt[i]]+depth[i]-depth[rst[i]];
if (nt2[i]) delta2[i]=depth[nt2[i]]-depth[i]+depth[i]-depth[rst[i]];
}
for (int i=1;i<=n;++i)
if (vst[i])
tong[++length]=i,num[i]=length;
for (int i=1;i<=n;++i)
for (int j=0;j<E[i].size();++j)
if (!vis[E[i][j].snum])
add2(num[i],num[E[i][j].num],E[i][j].data);
for (int i=1;i<=length;++i)
if (fa[tong[i]])
add2(i,num[rst[fa[tong[i]]]],depth[i]-depth[rst[fa[tong[i]]]]);
for (int i=0;i<=length;++i) dis[0][i]=dis[i][0]=inf;
for (int i=1;i<=length;++i) dijkstra(i);
for (int i=1;i<=n;++i)
{
if (p.find((node){nt[i],nt2[i]})==p.end()) st[++leng]=(node){nt[i],nt2[i]},p[(node){nt[i],nt2[i]}]=leng;
v[p[(node){nt[i],nt2[i]}]].push_back(i);
}
for (int i=1;i<=length;++i)
{
sort(v[i].begin(),v[i].end(),cmp),cnt[i].resize(v[i].size()),cnt2[i].resize(v[i].size());
for (int j=0;j<v[i].size();++j) cnt[i][j]=(j?cnt[i][j-1]:0)+delta[v[i][j]];
for (int j=(int)(v[i].size())-1;j>=0;--j) cnt2[i][j]=(j+1!=v[i].size()?cnt2[i][j+1]:0)+delta2[v[i][j]];
}
for (int i=1;i<=leng;++i)
for (int j=i+1;j<=leng;++j)
if (!v[i].empty()&&!v[j].empty())
{
d1=dis[st[i].x][st[j].x],d2=dis[st[i].x][st[j].y],d3=dis[st[i].y][st[j].x],d4=dis[st[i].y][st[j].y];
for (int k=0;k<v[i].size();++k)
{
ds1=min(delta[v[i][k]]+d1,delta2[v[i][k]]+d3),ds2=min(delta[v[i][k]]+d2,delta2[v[i][k]]+d4),ps=-1;
for (int t=lg[v[j].size()];t>=0;--t)
if (ps+(1<<t)<v[j].size()&&ds1+delta[v[j][ps+(1<<t)]]<ds2+delta2[v[j][ps+(1<<t)]])
ps+=(1<<t);
if (ps!=-1) ans+=1ll*ds1*(ps+1)+cnt[j][ps];
if (ps!=v[j].size()) ans+=1ll*ds2*(v[j].size()-ps)+cnt2[j][ps+1];
}
}
for (int i=1;i<=n;++i)
if (!vis[i])
{
top=0,dfs4(i);
for (int j=2;j<=top;++j) ans+=1ll*(szt[dque[1]]-szt[dque[j]])*szt[dque[j]]*tfa[dque[j]];
}
printf("%lld\n",ans);
return 0;
}
詳細信息
Subtask #1:
score: 0
Wrong Answer
Test #1:
score: 0
Wrong Answer
time: 0ms
memory: 40892kb
input:
300 1300 90 125 9397 157 77 3704 197 112 8218 152 235 1702 271 107 5600 117 92 1401 104 61 2242 127 230 1471 91 116 2740 29 127 4326 151 78 2569 273 241 7487 170 115 3100 152 171 2504 193 95 5921 30 281 1309 285 262 6462 100 265 8151 200 90 277 237 151 1123 231 219 974 238 176 2239 89 147 2256 233 2...
output:
-98914298877
result:
wrong answer 1st lines differ - expected: '324731073', found: '-98914298877'
Subtask #2:
score: 10
Accepted
Test #5:
score: 10
Accepted
time: 39ms
memory: 33864kb
input:
100000 99999 54625 54626 7146 20763 20764 300 41530 41531 9968 37448 37449 7434 81056 81055 700 27731 27730 8783 12408 12409 514 90652 90653 99 84104 84105 2524 83093 83094 195 17757 17756 2560 81925 81926 8935 14220 14219 9619 25516 25515 5883 89413 89412 275 46936 46937 3997 82755 82754 2775 53080...
output:
834269687204155387
result:
ok single line: '834269687204155387'
Test #6:
score: 0
Accepted
time: 36ms
memory: 25052kb
input:
100000 99999 13706 21290 3420 3037 78334 3887 94743 35121 9291 67873 91038 345 48348 12825 56 25237 56325 19 44215 92806 6788 40110 98929 5038 43250 87034 907 19698 18774 44 79406 51075 9523 79992 15613 4062 91111 66707 1595 1223 12300 3924 65613 22546 9008 24856 20394 393 14915 86273 5876 39594 160...
output:
4573680940298584
result:
ok single line: '4573680940298584'
Subtask #3:
score: 0
Runtime Error
Test #7:
score: 0
Runtime Error
input:
100000 100000 35241 48789 5098 4546 39869 6127 31415 22834 6026 25703 1952 6807 86143 78951 3421 34193 9615 4329 31012 98959 1664 81244 37874 3542 600 74315 9939 91066 57088 2111 5064 33313 9799 78834 28718 1133 41687 82171 4214 44801 87500 4238 40150 73606 5172 17787 30281 3718 52715 82529 4419 924...
output:
result:
Subtask #4:
score: 0
Skipped
Dependency #2:
100%
Accepted
Dependency #3:
0%
Subtask #5:
score: 0
Skipped
Dependency #4:
0%
Subtask #6:
score: 0
Skipped
Dependency #1:
0%