[数据结构与算法]第一届ACC联赛 决赛 C(什么换根dp)

题目
题意: 给定一棵n个点的有向树，选中一个点为根，使得到达个点的距离和最小。对于无法到达的点，可以将边进行翻转，求最少的翻转次数。
思路: 先建树嘛，给定的边边权为0，反边边权为1。随便选中一点root进行dfs，求出以root为根的最小代价。之后进行换根dp，这个换根还算简单的，连我都能想出来。。。
对于u的儿子j，如果根从u换成了j。
假设u到j的边权为1，换根以后j不需要代价即可到达u，所有u能到的点j也能到。f[j] = f[u]
假设u到j的边权为0，j到u边权则为1，需要1的代价从j到u。f[j] = f[u] + 1.
时间复杂度: O(n)
代码:

// Problem: 翻转树边
// Contest: AcWing
// URL: https://www.acwing.com/problem/content/4384/
// Memory Limit: 256 MB
// Time Limit: 1000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 2e5+10;
const int M = 1e6+10;
const int mod = 1e9+7;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
int n,m,k,T;
int h[N],e[N<<1],ne[N<<1],w[N<<1],idx = 0;
void add(int a,int b,int c)
{
	e[idx] = b,w[idx] = c,ne[idx] = h[a],h[a] = idx++;
}
int f[N]; //以u为根的最小代价
void dfs(int cur,int fa)
{
	for(int i=h[cur];~i;i=ne[i])
	{
		int j = e[i];
		if(j==fa) continue;
		dfs(j,cur);
		f[cur] += f[j] + w[i];
	}
}
void dfs2(int cur,int fa)
{
	for(int i=h[cur];~i;i=ne[i])
	{
		int j = e[i];
		if(j==fa) continue;
		if(w[i]==1) f[j] = f[cur] - 1;
		else f[j] = f[cur] + 1;
		dfs2(j,cur);
	}
}
int ans = 1e9;
void solve()
{
   mem(h,-1); idx = 0;
   read(n);
   for(int i=0;i<n-1;++i)
   {
   	  int x,y; read(x),read(y);
   	  add(x,y,0),add(y,x,1);
   }
   dfs(1,0);
   dfs2(1,0);
   for(int i=1;i<=n;++i) ans = min(ans,f[i]);
   // for(int i=1;i<=n;++i) cout<<i<<":"<<f[i]<<"\n";
   write(ans); puts("");
   for(int i=1;i<=n;++i)
   {
   	  if(f[i] == ans) write(i),putchar(' ');
   }
}
signed main(void)
{  
   T = 1;
   // OldTomato; cin>>T;
   // read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}