QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#368048 | #7790. 最短路求和 | NATURAL6 | 0 | 23ms | 22088kb | C++14 | 4.2kb | 2024-03-26 19:09:10 | 2024-03-26 19:09:10 |
Judging History
answer
#include<bits/stdc++.h>
using namespace std;
inline int qread()
{
int a=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)){(a*=10)+=(ch^48);ch=getchar();}
return a*f;
}
int n,m,cnt,tot,U[200010],V[200010],W[200010],deg[100010],siz[100010],del[100010],key[100010],id[100010];
int fr[3010],to[3010],C[3010],vis[100010],le[100010][2],dis[3010][3010];
long long ans,sum[100010][2],s[100010];
queue<int>que;
vector< pair<int,int> >e[100010],E[100010];
vector<int>S[3010];
inline int dfs(int rt,int ST,int da,int len)
{
vis[rt]=1;
if(da&&key[rt])
{
E[key[rt]].emplace_back(make_pair(key[ST],len));
E[key[ST]].emplace_back(make_pair(key[rt],len));
++tot;
fr[tot]=rt;to[tot]=ST;C[tot]=len;
return tot;
}
int p=0;
for(pair<int,int>i:e[rt])
{
if(i.first==da)continue;
if(!key[i.first]&&vis[i.first])continue;
p=dfs(i.first,ST,rt,len+i.second);
if(p&&rt^ST)
{
S[p].emplace_back(rt);
le[rt][0]=len;
le[rt][1]=C[p]-len;
}
}
return p;
}
struct node
{
int id,val;
}cx;
inline bool operator <(node x,node y){return x.val>y.val;}
priority_queue<node>qu;
inline void dij(int st,int *dis)
{
memset(dis,63,(cnt+1)<<2);
dis[st]=0;qu.push((node){st,0});
while(!qu.empty())
{
cx=qu.top();qu.pop();
if(cx.val^dis[cx.id])continue;
for(pair<int,int>i:e[cx.id])if(dis[i.first]>cx.val+i.second)qu.push((node){i.first,dis[i.first]=cx.val+i.second});
}
return ;
}
int main()
{
n=qread(),m=qread();
for(int i=1,u,v,w;i<=m;++i)
{
++deg[U[i]=u=qread()],++deg[V[i]=v=qread()],W[i]=w=qread();
e[u].emplace_back(make_pair(v,w));
e[v].emplace_back(make_pair(u,w));
}
for(int i=1;i<=n;++i)
{
siz[i]=1;
if(deg[i]<=1)que.push(i);
}
while(!que.empty())
{
++cnt;
int p=que.front();que.pop();
del[p]=1;
for(pair<int,int>i:e[p])
{
if(--deg[i.first]==1)que.push(i.first);
siz[i.first]+=siz[p];
ans+=1ll*i.second*siz[p]*(n-siz[p]);
}
}
if(cnt==n)
{
printf("%lld\n",ans);
return 0;
}
for(int i=1;i<=n;++i)e[i].clear();
for(int i=1;i<=m;++i)if(!del[U[i]]&&!del[V[i]])
{
e[U[i]].emplace_back(make_pair(V[i],W[i]));
e[V[i]].emplace_back(make_pair(U[i],W[i]));
}
cnt=0;
for(int i=1;i<=n;++i)if(deg[i]>=3)
{
key[i]=++cnt;
id[cnt]=i;
}
if(!cnt)
{
for(int i=1;i<=n;++i)if(!del[i])
{
key[i]=++cnt;
id[cnt]=i;
break;
}
}
for(int i=1;i<=cnt;++i)dfs(id[i],id[i],0,0);
for(int i=1;i<=cnt;++i)dij(i,dis[i]);
for(int i=1;i<=cnt;++i)for(int j=i+1;j<=cnt;++j)ans+=1ll*dis[i][j]*siz[id[i]]*siz[id[j]];
for(int A=1,l,r,sz;A<=tot;++A)if(!S[A].empty())
{
l=fr[A],r=to[A];sz=S[A].size();
for(int i=0;i<sz;++i)
{
sum[S[A][i]][0]=1ll*siz[S[A][i]]*le[S[A][i]][0];
sum[S[A][i]][1]=1ll*siz[S[A][i]]*le[S[A][i]][1];
s[S[A][i]]=siz[S[A][i]];
}
for(int i=1;i<sz;++i)
{
sum[S[A][i]][0]+=sum[S[A][i-1]][0];
sum[S[A][i]][1]+=sum[S[A][i-1]][1];
s[S[A][i]]+=s[S[A][i-1]];
}
for(int i=1,p=-1;i<sz;++i)
{
while(p+1<i&&le[S[A][i]][0]-le[S[A][p+1]][0]>=le[S[A][i]][1]+dis[key[l]][key[r]]+le[S[A][p+1]][0])++p;
while(p>=0&&le[S[A][i]][0]-le[S[A][p]][0]<le[S[A][i]][1]+dis[key[l]][key[r]]+le[S[A][p]][0])--p;
ans+=1ll*siz[S[A][i]]*(1ll*le[S[A][i]][0]*s[S[A][i-1]]-sum[S[A][i-1]][0]);
if(p>=0)
{
ans+=1ll*siz[S[A][i]]*(1ll*(le[S[A][i]][1]+dis[key[l]][key[r]]-le[S[A][i]][0])*s[S[A][p]]);
ans+=1ll*siz[S[A][i]]*(1ll*(sum[S[A][p]][0]+sum[S[A][p]][0]));
}
}
for(int i=1;i<=cnt;++i)for(int j:S[A])ans+=1ll*siz[j]*siz[id[i]]*min(dis[i][key[l]]+le[j][0],dis[i][key[r]]+le[j][1]);
for(int B=1,L,R,SZ,p,disL,disR;B<A;++B)if(!S[B].empty())
{
L=fr[B],R=to[B];
SZ=S[B].size(),p=0;
for(int i:S[A])
{
disL=min(dis[key[L]][key[l]]+le[i][0],dis[key[L]][key[r]]+le[i][1]);
disR=min(dis[key[R]][key[l]]+le[i][0],dis[key[R]][key[r]]+le[i][1]);
while(p>=1&&le[S[B][p-1]][0]+disL>=le[S[B][p-1]][1]+disR)--p;
while(p+1<=SZ&&le[S[B][p]][0]+disL<le[S[B][p]][1]+disR)++p;
ans+=1ll*siz[i]*(sum[S[B][SZ-1]][1]+disR*s[S[B][SZ-1]]);
if(p)ans+=1ll*siz[i]*(sum[S[B][p-1]][0]-sum[S[B][p-1]][1]+(disL-disR)*s[S[B][p-1]]);
}
}
}
printf("%lld\n",ans);
return 0;
}
詳細信息
Subtask #1:
score: 0
Wrong Answer
Test #1:
score: 0
Wrong Answer
time: 14ms
memory: 19392kb
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:
326072691
result:
wrong answer 1st lines differ - expected: '324731073', found: '326072691'
Subtask #2:
score: 0
Wrong Answer
Test #5:
score: 0
Wrong Answer
time: 18ms
memory: 17976kb
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:
1668554693120732537
result:
wrong answer 1st lines differ - expected: '834269687204155387', found: '1668554693120732537'
Subtask #3:
score: 0
Wrong Answer
Test #7:
score: 0
Wrong Answer
time: 23ms
memory: 22088kb
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:
9617450304999034
result:
wrong answer 1st lines differ - expected: '3317259529562659', found: '9617450304999034'
Subtask #4:
score: 0
Skipped
Dependency #2:
0%
Subtask #5:
score: 0
Skipped
Dependency #4:
0%
Subtask #6:
score: 0
Skipped
Dependency #1:
0%