首页 > 代码库 > BZOJ1131: [POI2008]Sta

BZOJ1131: [POI2008]Sta

1131: [POI2008]Sta

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 713  Solved: 196
[Submit][Status]

Description

给出一个N个点的树,找出一个点来,以这个点为根的树时,所有点的深度之和最大

Input

给出一个数字N,代表有N个点.N<=1000000 下面N-1条边.

Output

输出你所找到的点,如果具有多个解,请输出编号最小的那个.

Sample Input

8
1 4
5 6
4 5
6 7
6 8
2 4
3 4

Sample Output

7

HINT

Source

题解:

总算自己yy出一题。。。

刚开始想了一个贪心的算法,每次选取子树大小最小的子树递归下去,这样sum一定会变大?可是怎么证明正确性呢?看了一下样例发现这样做是错的。。。

然后发现sum传到儿子是可以O(1)更新的,于是就两次dfs水过了。。。

从父亲传到儿子,所有儿子的子树内的深度-1,其余深度+1,O(1)

代码:

  1 #include<cstdio>  2   3 #include<cstdlib>  4   5 #include<cmath>  6   7 #include<cstring>  8   9 #include<algorithm> 10  11 #include<iostream> 12  13 #include<vector> 14  15 #include<map> 16  17 #include<set> 18  19 #include<queue> 20  21 #include<string> 22  23 #define inf 1000000000 24  25 #define maxn 1000000+1000 26  27 #define maxm 500+100 28  29 #define eps 1e-10 30  31 #define ll long long 32  33 #define pa pair<int,int> 34  35 #define for0(i,n) for(int i=0;i<=(n);i++) 36  37 #define for1(i,n) for(int i=1;i<=(n);i++) 38  39 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 40  41 using namespace std; 42  43 inline int read() 44  45 { 46  47     int x=0,f=1;char ch=getchar(); 48  49     while(ch<0||ch>9){if(ch==-)f=-1;ch=getchar();} 50  51     while(ch>=0&&ch<=9){x=10*x+ch-0;ch=getchar();} 52  53     return x*f; 54  55 } 56 int n,tot,head[maxn],dep[maxn]; 57 ll s[maxn],sum[maxn]; 58 struct edge{int go,next;}e[2*maxn]; 59 inline void insert(int x,int y){e[++tot].go=y;e[tot].next=head[x];head[x]=tot;} 60 inline void dfs(int x) 61 { 62     s[x]=1; 63     for(int i=head[x],y;i;i=e[i].next) 64      if(!dep[y=e[i].go]) 65      { 66          sum[y]=dep[y]=dep[x]+1; 67          dfs(y); 68          s[x]+=s[y]; 69          sum[x]+=sum[y]; 70      } 71 }      72 inline void dfs2(int x) 73 { 74     for(int i=head[x],y;i;i=e[i].next) 75     if(dep[y=e[i].go]>dep[x]) 76     { 77         sum[y]=sum[x]-s[y]+n-s[y]; 78         dfs2(y); 79     } 80 } 81  82 int main() 83  84 { 85  86     freopen("input.txt","r",stdin); 87  88     freopen("output.txt","w",stdout); 89  90     n=read(); 91     for1(i,n-1){int x=read(),y=read();insert(x,y);insert(y,x);}; 92     dep[1]=sum[1]=1; 93     dfs(1); 94     dfs2(1); 95     int ans=0; 96     for1(i,n)if(sum[i]>sum[ans])ans=i; 97     printf("%d\n",ans);     98  99     return 0;100 101 }
View Code

 

BZOJ1131: [POI2008]Sta