結果
問題 | No.399 動的な領主 |
ユーザー |
![]() |
提出日時 | 2016-07-15 22:48:11 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 169 ms / 2,000 ms |
コード長 | 1,613 bytes |
コンパイル時間 | 1,562 ms |
コンパイル使用メモリ | 162,280 KB |
実行使用メモリ | 24,192 KB |
最終ジャッジ日時 | 2024-11-07 19:07:39 |
合計ジャッジ時間 | 4,628 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 19 |
ソースコード
#include <bits/stdc++.h>using namespace std;typedef signed long long ll;#undef _P#define _P(...) (void)printf(__VA_ARGS__)#define FOR(x,to) for(x=0;x<(to);x++)#define FORR(x,arr) for(auto& x:arr)#define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++)#define ALL(a) (a.begin()),(a.end())#define ZERO(a) memset(a,0,sizeof(a))#define MINUS(a) memset(a,0xff,sizeof(a))//-------------------------------------------------------int N;vector<int> E[101010];int P[21][200005],D[200005];ll tot[200005];ll ret;void dfs(int cur) {ITR(it,E[cur]) if(*it!=P[0][cur]) D[*it]=D[cur]+1, P[0][*it]=cur, dfs(*it);}int lca(int a,int b) {int ret=0,i,aa=a,bb=b;if(D[aa]>D[bb]) swap(aa,bb);for(i=19;i>=0;i--) if(D[bb]-D[aa]>=1<<i) bb=P[i][bb];for(i=19;i>=0;i--) if(P[i][aa]!=P[i][bb]) aa=P[i][aa], bb=P[i][bb];return (aa==bb)?aa:P[0][aa]; // vertex}void dfs2(int cur,int pre) {FORR(e,E[cur]) if(e!=pre) dfs2(e,cur);if(pre>=0) tot[pre] += tot[cur];ret += tot[cur]*(tot[cur]+1)/2;}void solve() {int i,j,k,l,r,x,y; string s;cin>>N;FOR(i,N-1) {cin>>x>>y;E[x-1].push_back(y-1);E[y-1].push_back(x-1);}dfs(0);FOR(i,19) FOR(x,N) P[i+1][x]=P[i][P[i][x]];cin>>i;while(i--) {cin>>x>>y;x--,y--;int lc=lca(x,y);tot[x]++;tot[y]++;tot[lc]--;if(lc!=0) tot[P[0][lc]]--;}dfs2(0,-1);cout<<ret<<endl;}int main(int argc,char** argv){string s;int i;if(argc==1) ios::sync_with_stdio(false), cin.tie(0);FOR(i,argc-1) s+=argv[i+1],s+='\n';FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin);solve(); return 0;}