QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#368053#7790. 最短路求和NATURAL60 26ms22584kbC++144.2kb2024-03-26 19:17:062024-03-26 19:17:07

Judging History

This is the latest submission verdict.

  • [2024-03-26 19:17:07]
  • Judged
  • Verdict: 0
  • Time: 26ms
  • Memory: 22584kb
  • [2024-03-26 19:17:06]
  • Submitted

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: 21132kb

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:

323986820

result:

wrong answer 1st lines differ - expected: '324731073', found: '323986820'

Subtask #2:

score: 0
Wrong Answer

Test #5:

score: 0
Wrong Answer
time: 15ms
memory: 16812kb

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: 26ms
memory: 22584kb

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%