This documentation is automatically generated by online-judge-tools/verification-helper
View the Project on GitHub mugen1337/procon
#define PROBLEM "https://judge.yosupo.jp/problem/tree_diameter" #include "../template.hpp" #include "../Graph2/TreeDiameter.hpp" signed main(){ int n;cin>>n; Graph<ll> g(n); g.read(n-1,0,true,false); auto [D,P]=TreeDiameter(g); cout<<D<<" "<<P.size()<<endl; cout<<P<<endl; return 0; }
#line 1 "test/yosupo_TreeDiameter.test.cpp" #define PROBLEM "https://judge.yosupo.jp/problem/tree_diameter" #line 1 "template.hpp" #include<bits/stdc++.h> using namespace std; #define ALL(x) begin(x),end(x) #define rep(i,n) for(int i=0;i<(n);i++) #define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl; #define mod 1000000007 using ll=long long; const int INF=1000000000; const ll LINF=1001002003004005006ll; int dx[]={1,0,-1,0},dy[]={0,1,0,-1}; // ll gcd(ll a,ll b){return b?gcd(b,a%b):a;} template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;} template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;} struct IOSetup{ IOSetup(){ cin.tie(0); ios::sync_with_stdio(0); cout<<fixed<<setprecision(12); } } iosetup; template<typename T> ostream &operator<<(ostream &os,const vector<T>&v){ for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" "); return os; } template<typename T> istream &operator>>(istream &is,vector<T>&v){ for(T &x:v)is>>x; return is; } #line 4 "test/yosupo_TreeDiameter.test.cpp" #line 1 "Graph2/GraphTemplate.hpp" // graph template // ref : https://ei1333.github.io/library/graph/graph-template.hpp template<typename T=int> struct Edge{ int from,to; T w; int idx; Edge()=default; Edge(int from,int to,T w=1,int idx=-1):from(from),to(to),w(w),idx(idx){} operator int() const{return to;} }; template<typename T=int> struct Graph{ vector<vector<Edge<T>>> g; int V,E; Graph()=default; Graph(int n):g(n),V(n),E(0){} int size(){ return (int)g.size(); } void resize(int k){ g.resize(k); V=k; } inline const vector<Edge<T>> &operator[](int k)const{ return (g.at(k)); } inline vector<Edge<T>> &operator[](int k){ return (g.at(k)); } void add_directed_edge(int from,int to,T cost=1){ g[from].emplace_back(from,to,cost,E++); } void add_edge(int from,int to,T cost=1){ g[from].emplace_back(from,to,cost,E); g[to].emplace_back(to,from,cost,E++); } void read(int m,int pad=-1,bool weighted=false,bool directed=false){ for(int i=0;i<m;i++){ int u,v;cin>>u>>v; u+=pad,v+=pad; T w=T(1); if(weighted) cin>>w; if(directed) add_directed_edge(u,v,w); else add_edge(u,v,w); } } }; #line 2 "Graph2/TreeDiameter.hpp" template<typename T> pair<T,vector<int>> TreeDiameter(const Graph<T> &g){ assert(g.V==g.E+1); vector<T> dis(g.V); vector<int> pre(g.V,-1); dis[0]=0; function<void(int,int)> dfs=[&g,&dfs,&dis,&pre](int par,int cur){ for(auto &e:g[cur])if(par!=e.to){ dis[e.to]=dis[cur]+e.w; pre[e.to]=cur; dfs(cur,e.to); } }; dfs(-1,0); int l=max_element(ALL(dis))-begin(dis); dis[l]=0; pre.assign(g.V,-1); dfs(-1,l); int r=max_element(ALL(dis))-begin(dis); vector<int> path; for(int i=r;i!=-1;i=pre[i]) path.push_back(i); return make_pair(dis[r],path); } #line 6 "test/yosupo_TreeDiameter.test.cpp" signed main(){ int n;cin>>n; Graph<ll> g(n); g.read(n-1,0,true,false); auto [D,P]=TreeDiameter(g); cout<<D<<" "<<P.size()<<endl; cout<<P<<endl; return 0; }