QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#279722#5313. Please Save Pigelandzzuqy#WA 420ms54228kbC++112.9kb2023-12-09 00:25:412023-12-09 00:25:41

Judging History

你现在查看的是最新测评结果

  • [2023-12-09 00:25:41]
  • 评测
  • 测评结果:WA
  • 用时:420ms
  • 内存:54228kb
  • [2023-12-09 00:25:41]
  • 提交

answer

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll read()
{
    ll x;cin>>x;
    return x;
}
ll lowbit(ll x)
{
    return x&(-x);
}
const ll N=5e5+10;
ll n,k,a[N],c[N],gcd[N*4];
ll head[N],tot,L[N],R[N],vis[N],D[N],cnt[N],f[N];
ll SUM[N];
struct edge
{
	ll y,v,next;
}e[N*2];
void add0(ll l,ll r,ll v)
{
	while(l<=n)
    {
        c[l]+=v;
        l+=lowbit(l);
    }
    r++;
    while(r<=n)
    {
        c[r]-=v;
        r+=lowbit(r);
    }
}
ll ask0(ll x)
{
    ll t=0;
    while(x)
    {
        t=t+c[x];
        x=x-lowbit(x);
    }
    return t;
}
void add(ll x,ll l,ll r,ll d,ll v)
{
    if(l==r)
    {
        gcd[x]=gcd[x]+v;
        return ;
    }
    ll mid=(l+r)/2;
    if(d<=mid)
        add(x*2,l,mid,d,v);
    else
        add(x*2+1,mid+1,r,d,v);
    gcd[x]=__gcd(gcd[x*2],gcd[x*2+1]);
}
void e_add(ll x,ll y,ll v)
{
	tot++;
	e[tot]={y,v,head[x]};
	head[x]=tot;
}
void dfs(ll x,ll fa,ll d)
{
	if(vis[x])
	{
		tot++;
		L[x]=R[x]=tot;
		a[tot]=d;
		cnt[x]=1;
	}
	else
		L[x]=1e9,R[x]=-1;
	for(ll i=head[x];i;i=e[i].next)
	{
		if(e[i].y==fa)
			continue;
		dfs(e[i].y,x,d+e[i].v);
		L[x]=min(L[x],L[e[i].y]);
		R[x]=max(R[x],R[e[i].y]);
		cnt[x]+=cnt[e[i].y];
		f[x]=f[x]+f[e[i].y]+e[i].v*cnt[e[i].y];
	}
}
void ADD(ll l,ll r,ll v)
{
	add0(l,r,v);
	if(l!=1)
		add(1,2,k,l,v);
	if(r!=k)
		add(1,2,k,r+1,-v);
}
void dfs1(ll x,ll fa)
{
	D[x]=abs(__gcd(gcd[1],ask0(1)));
	for(ll i=head[x];i;i=e[i].next)
	{
		if(e[i].y==fa)
			continue;
		if(L[e[i].y]!=1000000000)
		{
			if(R[e[i].y]!=k)
				ADD(R[e[i].y]+1,k,e[i].v);
			if(L[e[i].y]!=1)
				ADD(1,L[e[i].y]-1,e[i].v);
			ADD(L[e[i].y],R[e[i].y],-e[i].v);
		}
		dfs1(e[i].y,x);
		if(L[e[i].y]!=1000000000)
		{
			if(R[e[i].y]!=k)
				ADD(R[e[i].y]+1,k,-e[i].v);
			if(L[e[i].y]!=1)
				ADD(1,L[e[i].y]-1,-e[i].v);
			ADD(L[e[i].y],R[e[i].y],e[i].v);
		}
	}
}
void dfs2(ll x,ll fa,ll sum)
{
	SUM[x]=sum+f[x];
	for(ll i=head[x];i;i=e[i].next)
	{
		if(e[i].y==fa)
			continue;
		f[x]=f[x]-f[e[i].y]-e[i].v*cnt[e[i].y];
		dfs2(e[i].y,x,sum+f[x]+(k-cnt[e[i].y])*e[i].v);
		f[x]=f[x]+f[e[i].y]+e[i].v*cnt[e[i].y];
	}
}
int main()
{
    // freopen("1.in","r",stdin);
    // freopen("2.out","w",stdout);
    n=read();k=read();
    if(k==1)
    {
    	cout<<"0\n";
    	return 0;
    }
    if(k==2)
    {
    	cout<<"1\n";
    	return 0;
    }
    for(ll i=1;i<=k;i++){
    	int x=read();
    	assert(vis[x]==0);
    	vis[x]=1;

    }
    for(ll i=1;i<n;i++)
    {
    	ll x=read(),y=read(),v=read();
    	e_add(x,y,v);
    	e_add(y,x,v);
    }
    tot=0;
    dfs(1,0,0);
    for(ll i=1;i<=k;i++)
    {
    	add0(i,i,a[i]);
    	if(i!=1)
    		add(1,2,k,i,a[i]);
    	if(i!=k)
    		add(1,2,k,i+1,-a[i]);
    }
    dfs1(1,0);
    dfs2(1,0,0);
    ll ans=SUM[1]*2/D[1];
    for(ll i=1;i<=n;i++)
    	ans=min(ans,SUM[i]*2/D[i]);
    cout<<ans;
}

详细

Test #1:

score: 100
Accepted
time: 2ms
memory: 9572kb

input:

5 3
3 4 5
1 2 2
2 3 4
2 5 4
3 4 6

output:

8

result:

ok 1 number(s): "8"

Test #2:

score: 0
Accepted
time: 1ms
memory: 7724kb

input:

10 3
1 7 10
7 6 3
1 8 3
3 6 3
8 6 2
4 1 1
10 6 4
2 8 3
9 10 3
5 10 3

output:

24

result:

ok 1 number(s): "24"

Test #3:

score: 0
Accepted
time: 0ms
memory: 3392kb

input:

1 1
1

output:

0

result:

ok 1 number(s): "0"

Test #4:

score: 0
Accepted
time: 0ms
memory: 3572kb

input:

100000 1
79187
72704 72659 15
32741 43496 10
21580 97447 17
55758 36700 21
32116 3643 14
60460 58764 12
75894 50624 7
58295 49393 22
43733 17210 1
58093 68769 15
1086 58916 17
25632 37710 11
49555 92976 8
32547 27060 18
84896 12811 1
3196 1242 16
18870 78236 14
2414 7945 12
48745 15399 1
17648 83791...

output:

0

result:

ok 1 number(s): "0"

Test #5:

score: 0
Accepted
time: 0ms
memory: 7572kb

input:

100 10
3 27 33 45 48 72 76 91 92 100
66 98 4
70 17 2
28 59 4
26 25 3
77 92 1
40 61 2
11 27 2
85 35 1
57 26 1
68 99 4
50 84 1
20 82 3
31 39 1
71 7 4
54 55 4
60 26 4
56 61 2
15 66 3
95 53 2
8 60 4
21 82 1
18 81 2
29 73 3
94 4 1
10 4 4
86 43 1
62 41 1
45 57 1
25 66 3
69 89 2
14 53 3
27 92 1
42 98 4
13 ...

output:

200

result:

ok 1 number(s): "200"

Test #6:

score: 0
Accepted
time: 0ms
memory: 7516kb

input:

100 90
75 17 78 84 52 69 54 24 74 35 58 51 72 7 21 70 93 15 60 87 13 40 23 92 99 53 27 22 91 46 56 86 61 19 44 98 50 28 14 12 55 64 30 80 95 38 18 43 31 89 20 16 8 65 63 79 59 34 97 25 2 11 67 71 29 9 37 76 77 26 39 68 32 62 90 10 85 49 42 45 96 83 94 3 6 100 81 57 88 47
67 65 1
43 78 4
98 71 3
71 2...

output:

1798

result:

ok 1 number(s): "1798"

Test #7:

score: 0
Accepted
time: 1ms
memory: 7656kb

input:

1000 10
111 240 479 530 572 583 644 652 753 869
121 923 2
886 685 4
446 284 4
352 250 1
540 485 2
72 154 4
522 693 3
664 917 4
792 941 3
132 832 4
709 186 3
509 114 2
824 978 2
216 265 2
138 570 1
498 959 4
434 222 1
803 693 1
253 677 4
172 463 3
383 978 2
718 959 3
369 421 4
568 454 4
256 938 1
6 1...

output:

226

result:

ok 1 number(s): "226"

Test #8:

score: 0
Accepted
time: 0ms
memory: 7632kb

input:

1000 957
233 514 228 739 827 85 840 175 766 807 19 276 549 611 145 511 895 121 116 525 280 431 810 629 990 509 542 324 241 801 849 506 178 176 49 528 221 742 444 513 111 505 442 794 107 392 291 674 298 803 198 927 738 590 706 804 860 512 421 618 697 516 335 420 418 288 544 694 330 776 104 510 621 47...

output:

32602

result:

ok 1 number(s): "32602"

Test #9:

score: -100
Wrong Answer
time: 420ms
memory: 54228kb

input:

500000 4
182462 188845 259396 281751
456733 79213 9204078
395954 45205 3919968
454058 310013 734433
433648 435834 3887333
448797 138275 9946222
385528 63721 3037094
44276 184047 1799127
169565 81666 3752583
459111 229807 5534913
374868 374333 8627923
476055 408523 2692999
445258 424229 3038119
92885...

output:

486956092

result:

wrong answer 1st numbers differ - expected: '193870600', found: '486956092'