QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#235884 | #5313. Please Save Pigeland | ushg8877# | WA | 460ms | 51428kb | C++20 | 3.3kb | 2023-11-03 11:58:49 | 2023-11-03 11:58:50 |
Judging History
answer
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define MP make_pair
mt19937 rnd(time(0));
const int MAXN=5e5+5;
int siz[MAXN];ll sum[MAXN],rep[MAXN],f[MAXN];bool mark[MAXN];
// f 表示内部差的 gcd
int n,k;
int eu[MAXN],ev[MAXN],ew[MAXN];ll ans=5e18;
vector<int> edg[MAXN];
int gcd(int x,int y){
if(!x||!y) return x|y;
else if(x%y==0) return y;
else return gcd(y,x%y);
}
void dfs(int u,int fa){
rep[u]=-1;siz[u]=0;sum[u]=0;
if(mark[u]){
siz[u]=1;rep[u]=0;
}
for(int i:edg[u]){
int v=(eu[i]==u?ev[i]:eu[i]);
if(v==fa) continue;
dfs(v,u);
siz[u]+=siz[v];sum[u]+=sum[v]+siz[v]*ew[i];
f[u]=gcd(f[u],f[v]);
if(rep[v]!=-1){
if(rep[u]!=-1) f[u]=gcd(f[u],abs(rep[u]-(rep[v]+ew[i])));
else rep[u]=rep[v]+ew[i];
}
}
// cout<<"get "<<u<<" "<<f[u]<<' '<<rep[u]<<endl;
return;
}
void root(int u,int fa){
rep[u]=-1;siz[u]=0;sum[u]=0;f[u]=0;
// cout<<"Init "<<u<<" "<<sum[u]<<' '<<siz[u]<<' '<<gcd(rep[u],f[u])<<endl;
// for(int i=1;i<=n;i++) cout<<rep[i]<<' ';cout<<endl;
// for(int i=1;i<=n;i++) cout<<f[i]<<' ';cout<<endl;
// for(int i=1;i<=n;i++) cout<<siz[i]<<' ';cout<<endl;
// for(int i=1;i<=n;i++) cout<<sum[i]<<' ';cout<<endl;
if(mark[u]){
siz[u]=1;rep[u]=0;
}
for(int i:edg[u]){
int v=(eu[i]==u?ev[i]:eu[i]);
siz[u]+=siz[v];sum[u]+=sum[v]+siz[v]*ew[i];
f[u]=gcd(f[u],f[v]);
if(rep[v]!=-1){
if(rep[u]!=-1) f[u]=gcd(f[u],abs(rep[u]-(rep[v]+ew[i])));
else rep[u]=rep[v]+ew[i];
}
}
ans=min(ans,sum[u]/gcd(rep[u],f[u]));
// cout<<"Get "<<u<<" "<<sum[u]<<' '<<siz[u]<<' '<<gcd(rep[u],f[u])<<endl;
// for(int i=1;i<=n;i++) cout<<rep[i]<<' ';cout<<endl;
// for(int i=1;i<=n;i++) cout<<f[i]<<' ';cout<<endl;
// for(int i=1;i<=n;i++) cout<<siz[i]<<' ';cout<<endl;
// for(int i=1;i<=n;i++) cout<<sum[i]<<' ';cout<<endl;
int e=edg[u].size();
vector<ll> pf(e),sf(e),pr(e),sr(e);
rep[u]=-1;f[u]=0;
if(mark[u]) rep[u]=0;
for(int id=0;id<e;id++){
int i=edg[u][id];
int v=(eu[i]==u?ev[i]:eu[i]);
f[u]=gcd(f[u],f[v]);
if(rep[v]!=-1){
if(rep[u]!=-1) f[u]=gcd(f[u],abs(rep[u]-(rep[v]+ew[i])));
else rep[u]=rep[v]+ew[i];
}
pf[id]=f[u],pr[id]=rep[u];
}
rep[u]=-1;f[u]=0;
if(mark[u]) rep[u]=0;
for(int id=e-1;id>=0;id--){
int i=edg[u][id];
int v=(eu[i]==u?ev[i]:eu[i]);
f[u]=gcd(f[u],f[v]);
if(rep[v]!=-1){
if(rep[u]!=-1) f[u]=gcd(f[u],abs(rep[u]-(rep[v]+ew[i])));
else rep[u]=rep[v]+ew[i];
}
sf[id]=f[u],sr[id]=rep[u];
}
for(int id=0;id<e;id++){
int i=edg[u][id],v=(eu[i]==u?ev[i]:eu[i]);
if(v==fa) continue;
ll _sizv=siz[v],_sizu=siz[u],_sumu=sum[u],_sumv=sum[v];
siz[u]=_sizu-_sizv;
sum[u]=_sumu-_sumv-ew[i]*_sizv;
// cout<<"Change "<<v<<' '<<siz[u]<<' '<<_sizu<<' '<<_sizv<<endl;
f[u]=gcd(id?pf[id-1]:0,id<e-1?sf[id+1]:0);
if(id>0&&pr[id-1]!=-1&&id<e-1&&sr[id+1]!=-1) f[u]=gcd(f[u],abs(pr[id-1]-sr[id+1]));
rep[u]=(id>0&&pr[id-1]!=-1?pr[id-1]:(id<e-1?sr[id+1]:-1));
root(v,u);
siz[v]=_sizv;siz[u]=_sizu;sum[u]=_sumu;sum[v]=_sumv;
}
return;
}
int main(){
ios::sync_with_stdio(false);
cin>>n>>k;
if(k==1){cout<<"0"<<endl;return 0;}
while(k--){int u;cin>>u;mark[u]=true;}
for(int i=1;i<n;i++){
cin>>eu[i]>>ev[i]>>ew[i];
edg[eu[i]].push_back(i);
edg[ev[i]].push_back(i);
}
dfs(1,0);root(1,0);
cout<<ans*2<<endl;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 29096kb
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: 0ms
memory: 28872kb
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: 4ms
memory: 18148kb
input:
1 1 1
output:
0
result:
ok 1 number(s): "0"
Test #4:
score: 0
Accepted
time: 0ms
memory: 18860kb
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: 27380kb
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: 29288kb
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: 3ms
memory: 28456kb
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: 27424kb
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: 0
Accepted
time: 460ms
memory: 51392kb
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:
193870600
result:
ok 1 number(s): "193870600"
Test #10:
score: 0
Accepted
time: 447ms
memory: 51428kb
input:
500000 150000 1 4 6 12 19 20 21 24 30 34 35 37 39 40 41 45 47 48 49 51 60 65 71 72 75 80 83 84 90 93 94 96 97 99 105 113 115 125 128 132 135 136 138 139 142 145 146 150 152 159 164 167 179 180 186 188 189 190 195 200 204 205 208 210 211 212 214 217 223 226 230 233 238 243 247 260 263 264 268 269 271...
output:
8719582
result:
ok 1 number(s): "8719582"
Test #11:
score: 0
Accepted
time: 4ms
memory: 28808kb
input:
1000 10 334 271 592 707 890 566 590 450 169 129 784 868 7940645 969 691 6692155 864 76 6691895 853 343 2376735 12 897 7467840 783 560 3342995 908 372 952130 5 843 3055195 531 64 7912240 313 701 2196620 246 52 9831250 808 788 6735060 414 340 168805 278 276 6271820 97 247 2379595 490 459 701440 422 45...
output:
15519516
result:
ok 1 number(s): "15519516"
Test #12:
score: 0
Accepted
time: 0ms
memory: 29084kb
input:
1000 10 935 536 379 686 787 66 829 374 369 957 913 887 7565610 118 846 8734750 109 698 7065147 156 813 2228752 769 869 9612951 887 99 4384470 253 494 6908549 655 124 4406482 242 567 5548735 385 694 8703600 705 894 798435 930 179 9333610 346 83 9417442 621 768 953905 490 586 7780155 149 897 7212005 1...
output:
23133220
result:
ok 1 number(s): "23133220"
Test #13:
score: -100
Wrong Answer
time: 3ms
memory: 27436kb
input:
1000 500 474 953 613 122 585 861 627 461 231 466 197 437 807 645 724 259 221 65 312 263 441 895 643 507 159 279 257 801 577 546 160 706 674 840 433 761 921 818 40 862 535 443 396 75 311 836 959 57 986 603 179 133 110 775 990 269 539 858 153 124 529 979 423 241 561 911 151 223 570 275 586 969 49 390 ...
output:
541794
result:
wrong answer 1st numbers differ - expected: '1414024', found: '541794'