結果
問題 | No.1333 Squared Sum |
ユーザー | penguinman |
提出日時 | 2020-11-03 07:03:29 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 351 ms / 2,000 ms |
コード長 | 2,681 bytes |
コンパイル時間 | 3,312 ms |
コンパイル使用メモリ | 182,332 KB |
実行使用メモリ | 55,368 KB |
最終ジャッジ日時 | 2023-09-29 14:23:51 |
合計ジャッジ時間 | 16,736 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
4,376 KB |
testcase_01 | AC | 1 ms
4,380 KB |
testcase_02 | AC | 2 ms
4,380 KB |
testcase_03 | AC | 313 ms
30,544 KB |
testcase_04 | AC | 330 ms
30,572 KB |
testcase_05 | AC | 324 ms
30,536 KB |
testcase_06 | AC | 317 ms
30,524 KB |
testcase_07 | AC | 320 ms
30,532 KB |
testcase_08 | AC | 323 ms
30,508 KB |
testcase_09 | AC | 332 ms
30,632 KB |
testcase_10 | AC | 328 ms
30,488 KB |
testcase_11 | AC | 329 ms
30,500 KB |
testcase_12 | AC | 324 ms
30,612 KB |
testcase_13 | AC | 278 ms
55,368 KB |
testcase_14 | AC | 349 ms
42,940 KB |
testcase_15 | AC | 351 ms
52,480 KB |
testcase_16 | AC | 1 ms
4,380 KB |
testcase_17 | AC | 2 ms
4,376 KB |
testcase_18 | AC | 1 ms
4,376 KB |
testcase_19 | AC | 1 ms
4,376 KB |
testcase_20 | AC | 2 ms
4,376 KB |
testcase_21 | AC | 1 ms
4,380 KB |
testcase_22 | AC | 2 ms
4,380 KB |
testcase_23 | AC | 2 ms
4,380 KB |
testcase_24 | AC | 2 ms
4,380 KB |
testcase_25 | AC | 1 ms
4,380 KB |
testcase_26 | AC | 335 ms
54,856 KB |
testcase_27 | AC | 319 ms
43,728 KB |
testcase_28 | AC | 332 ms
54,300 KB |
testcase_29 | AC | 255 ms
55,352 KB |
testcase_30 | AC | 86 ms
14,568 KB |
testcase_31 | AC | 42 ms
9,560 KB |
testcase_32 | AC | 143 ms
19,584 KB |
testcase_33 | AC | 104 ms
16,144 KB |
testcase_34 | AC | 221 ms
25,220 KB |
testcase_35 | AC | 133 ms
19,804 KB |
testcase_36 | AC | 70 ms
13,012 KB |
testcase_37 | AC | 71 ms
13,592 KB |
testcase_38 | AC | 86 ms
15,100 KB |
testcase_39 | AC | 160 ms
21,936 KB |
testcase_40 | AC | 198 ms
31,964 KB |
testcase_41 | AC | 194 ms
31,960 KB |
testcase_42 | AC | 213 ms
31,968 KB |
testcase_43 | AC | 199 ms
31,920 KB |
ソースコード
#include<bits/stdc++.h> using std::cin; using std::cout; #define endl "\n" using std::vector; using ll=long long; //宣言、modpow const int mod=1e9+7; vector<vector<int>> edge,weight; vector<ll> dp,subtree,sum; vector<bool> flag; ll ans=0; ll modpow(ll x,ll y){ ll ret=1; while(y){ if(y&1){ ret*=x; ret%=mod; } x*=x; x%=mod; y/=2; } return ret; } //dfs void dfs(int now){ flag[now]=1; for(int i=0;i<edge[now].size();i++){ int next=edge[now][i]; ll w=weight[now][i]; if(flag[next]) continue; dfs(next); subtree[now]+=subtree[next]; dp[now]+=dp[next]+subtree[next]*w%mod*w%mod+sum[next]*w%mod*2%mod; dp[now]%=mod; sum[now]+=sum[next]+w*subtree[next]%mod; sum[now]%=mod; } } //rerooting void reroot(int now){ flag[now]=0; ans+=dp[now]; ans%=mod; for(int i=0;i<edge[now].size();i++){ int next=edge[now][i]; ll w=weight[now][i]; if(!flag[next]) continue; ll dp2=dp[now]-dp[next]-subtree[next]*w%mod*w%mod-sum[next]*w%mod*2%mod; dp2%=mod; if(dp2<0) dp2+=mod; ll sum2=sum[now]-sum[next]-w*subtree[next]%mod; sum2%=mod; if(sum2<0) sum2+=mod; ll subtree2=subtree[now]-subtree[next]; dp[next]+=dp2+subtree2*w%mod*w%mod+sum2*w%mod*2%mod; dp[next]%=mod; sum[next]+=sum2+w*subtree2%mod; sum[next]%=mod; subtree[next]+=subtree2; reroot(next); } } struct Union_Find{ int N; vector<int> par; Union_Find(int n):N(n){ par.resize(N); for(int i=0;i<N;i++) par[i]=i; } int root(int x){ if(par[x]==x) return x; return par[x]=root(par[x]); } void unite(int x,int y){ x=root(x); y=root(y); if(x==y) return; par[x]=y; } }; const int max=2e5,inf=1e9; //main int main(){ std::ios::sync_with_stdio(false); std::cin.tie(nullptr); //resize、input int N; cin>>N; assert(1<=N&&N<=max); edge.resize(N); weight.resize(N); dp.resize(N); subtree.resize(N,1); sum.resize(N); flag.resize(N); Union_Find tree(N); for(int i=1;i<N;i++){ int x,y,z; cin>>x>>y>>z; assert(1<=x&&x<=N&&1<=y&&y<=N&&1<=z&&z<=inf); assert(tree.root(x-1)!=tree.root(y-1)); tree.unite(x-1,y-1); edge[x-1].push_back(y-1); edge[y-1].push_back(x-1); weight[x-1].push_back(z); weight[y-1].push_back(z); } //求値、output dfs(0); reroot(0); ans*=modpow(2,mod-2); ans%=mod; cout<<ans<<endl; }